clear all
close all
clc
%%                            skript BIOMETRIE - snimky LAMP
% vstupem je puvodni sedotovy snimek oka
% jednotlive kroky:
% 1. Detekce vnitrni hranice duhovky
% 2. Normalizece duhovky
% 3. Predzpracovani normalizovaneho obrazu
% 4. Ziskani 1D signalů
% 5. Extrakce priznaku - vystupem 3D matice priznaku
% ---------------------------------------------------------------------------
%% Skript
    I = imread('S2005L01.jpg'); % nacteni vzoroveho snimku
    figure(1); imshow(I);title('original')
        
%% 1. Detekce vnitrni hranice duhovky - fast radial symmetry
    [xs,ys,radius,xp,yp,I] = radial_symmetry (I);
    
    figure(2);
    imshow(I,[]); title('Interval','FontSize', 11);
    hold on
    plot(xs,ys,'gx')
    hold on
    plot(xs+xp,ys+yp,'.g')
    f=figure(2);
    
% pripadne ulozeni je uvedeno v kompletni verzi na DVD

%% 2. Normalizace snimku - Prevod z kartezskych do polarnich
    [normalizovany,R_image,R_step, D_step] = car_pol(xs,ys,radius,I);

    figure(3);
    imshow(normalizovany,[]); title('normalizovany snimek')
    figure(4);
    imshow(R_image,[])
    
% LAMP pripadne ulozeni je uvedeno v kompletni verzi na DVD
    
%% 3. Predzpracovani normalizovaneho obrazu

    [upraveny] = predzpracovani(normalizovany);
    
    figure(5);
    imshow(upraveny,[]); title('upraveny normalizovany')
%LAMP pripadne ulozeni je uvedeno v kompletni verzi na DVD

%% 4. Ziskani 1D signalů

signals1D = extraction1D (upraveny);

%% 5. Extrakce priznaku

m = GaussHermite(signals1D);

%% 6. vyber specifickych uhlu

% LAMP pripadne ulozeni je uvedeno v kompletni verzi na DVD
