rec.m

Fingerprint source code using matlab

clear;
clc;
close all;

global immagine n_bands h_bands n_arcs h_radius h_lato n_sectors matrice num_disk


n_bands=4;
h_bands=20;
n_arcs=16;
h_radius=12;
h_lato=h_radius+(n_bands*h_bands*2)+16;

if mod(h_lato,2)==0
    h_lato=h_lato-1;
end

n_sectors=n_bands*n_arcs;
matrice=zeros(h_lato);

for ii=1:(h_lato*h_lato)
    matrice(ii)=whichsector(ii,h_radius);
end

num_disk=8;
choix=0;
possibility=6;
messaggio='Insert the number of set: each set determins a class. This set should include a number of images for each person, with some variations in expression and in the lighting.';

while choix~=possibility,
      choix=menu('reconnaissance d"empreinte digitale','selectionner une empreinte','reconnaissance d"empreinte','Exit');
   %=======================================================================
   %=======================================================================
   %=======================================================================
   
        if choix==1,
        clc;
        close all;
        selezionato=0;
        while selezionato==0
            [namefile,pathname]=uigetfile({'*.bmp;*.tif;*.tiff;*.jpg;*.fig;*.jpeg;*.gif','IMAGE Files (*.bmp,*.tif,*.tiff,*.jpg,*.jpeg,*.gif)'},'Chose GrayScale Image');
            if namefile~=0
                [img,map]=imread(strcat(pathname,namefile));
                selezionato=1;
            else
                disp('Select a grayscale image');
            end
            if (any(namefile~=0) && (~isgray(img)))
                disp('Select a grayscale image');
                selezionato=0;
            end
        end
        
        immagine=double(img);
        
        if isa(img,'uint8')
            graylevmax=2^8-1;
        end
        if isa(img,'uint16')
            graylevmax=2^16-1;
        end
        if isa(img,'uint32')
            graylevmax=2^32-1;
        end
        fingerprint = immagine;
     
        N=h_lato;
        
        [BinarizedPrint,XofCenter,YofCenter]=centralizing(fingerprint,0);
        [CroppedPrint]=cropping(XofCenter,YofCenter,fingerprint);
        [NormalizedPrint,vector]=sector_norm(CroppedPrint,0);
        
        for (angle=0:1:num_disk-1)    
            gabor=gabor2d_sub(angle,num_disk);
            ComponentPrint=conv2fft(NormalizedPrint,gabor,'same');
            [disk,vector]=sector_norm(ComponentPrint,1);    
            finger_code1{angle+1}=vector(1:n_sectors);
        end
        
        
        img=imrotate(img,180/(num_disk*2));
        fingerprint=double(img);
        
        [BinarizedPrint,XofCenter,YofCenter]=centralizing(fingerprint,0);
        [CroppedPrint]=cropping(XofCenter,YofCenter,fingerprint);
        [NormalizedPrint,vector]=sector_norm(CroppedPrint,0);
        
        for (angle=0:1:num_disk-1)    
            gabor=gabor2d_sub(angle,num_disk);
            ComponentPrint=conv2fft(NormalizedPrint,gabor,'same');
            [disk,vector]=sector_norm(ComponentPrint,1);    
            finger_code2{angle+1}=vector(1:n_sectors);
        end
        % FingerCode added to database
        if (exist('fp_database.dat')==2)
            load('fp_database.dat','-mat');
            fp_number=fp_number+1;
            data{fp_number,1}=finger_code1;
            data{fp_number,2}=finger_code2;
            save('fp_database.dat','data','fp_number','-append');
        else
            fp_number=1;
            data{fp_number,1}=finger_code1;
            data{fp_number,2}=finger_code2;
            save('fp_database.dat','data','fp_number');
        end
        
        message=strcat('Ajout avec sucess ',num2str(fp_number));
        msgbox(message,'Empreinte Digitale','Aide');
        end
        %==================================================================
        %==================================================================
        %==================================================================
        
        if choix==2,
         clc;
         close all;
         selezionato=0;
         while selezionato==0
              [namefile,pathname]=uigetfile({'*.bmp;*.tif;*.tiff;*.jpg;*.jpeg;*.gif','IMAGE Files (*.bmp,*.tif,*.tiff,*.jpg,*.jpeg,*.gif)'},'Chose GrayScale Image');
              if namefile~=0
                 [img,map]=imread(strcat(pathname,namefile));
                 selezionato=1;
              else
                 disp('Select a grayscale image');
              end
              if (any(namefile~=0) && (~isgray(img)))
                 disp('Select a grayscale image');
                 selezionato=0;
              end
         end
        
          immagine=double(img);
        
         if isa(img,'uint8')
             graylevmax=2^8-1;
         end
         if isa(img,'uint16')
             graylevmax=2^16-1;
         end
         if isa(img,'uint32')
             graylevmax=2^32-1;
         end
         fingerprint = immagine;
         N=h_lato;
         [BinarizedPrint,XofCenter,YofCenter]=centralizing(fingerprint,0);
         [CroppedPrint]=cropping(XofCenter,YofCenter,fingerprint);
         [NormalizedPrint,vector]=sector_norm(CroppedPrint,0);
         vettore_in=zeros(num_disk*n_sectors,1);
         for (angle=0:1:num_disk-1)    
             gabor=gabor2d_sub(angle,num_disk);
             ComponentPrint=conv2fft(NormalizedPrint,gabor,'same');
             [disk,vector]=sector_norm(ComponentPrint,1);    
             finger_code{angle+1}=vector(1:n_sectors);
             vettore_in(angle*n_sectors+1:(angle+1)*n_sectors)=finger_code{angle+1};
         end
         if (exist('fp_database.dat')==2)
            load('fp_database.dat','-mat');
            vettore_a=zeros(num_disk*n_sectors,1);
            vettore_b=zeros(num_disk*n_sectors,1);
            best_matching=zeros(fp_number,1);
            valori_rotazione=zeros(n_arcs,1);
            for scanning=1:fp_number
                fcode1=data{scanning,1};
                fcode2=data{scanning,2};
                for rotazione=0:(n_arcs-1)
                    p1=fcode1;
                    p2=fcode2;
                    % ruoto i valori dentro disco
                    for conta_disco=1:num_disk
                        disco1=p1{conta_disco};
                        disco2=p2{conta_disco};
                        for old_pos=1:n_arcs
                            new_pos=mod(old_pos+rotazione,n_arcs);
                            if new_pos==0
                                new_pos=n_arcs;
                            end
                            for conta_bande=0:1:(n_bands-1)
                                disco1r(new_pos+conta_bande*n_arcs)=disco1(old_pos+conta_bande*n_arcs);
                                disco2r(new_pos+conta_bande*n_arcs)=disco2(old_pos+conta_bande*n_arcs);
                            end
                        end
                        p1{conta_disco}=disco1r;
                        p2{conta_disco}=disco2r;
                    end
                    % ruoto i dischi circolarmente
                    for old_disk=1:num_disk
                        new_disk=mod(old_disk+rotazione,num_disk);
                        if new_disk==0
                            new_disk=num_disk;
                        end
                        pos=old_disk-1;
                        vettore_a(pos*n_sectors+1:(pos+1)*n_sectors)=p1{new_disk};
                        vettore_b(pos*n_sectors+1:(pos+1)*n_sectors)=p2{new_disk};                    
                    end
                    d1=norm(vettore_a-vettore_in);
                    d2=norm(vettore_b-vettore_in);
                    if d1<d2
                        val_minimo=d1;
                    else
                        val_minimo=d2;
                    end
                    valori_rotazione(rotazione+1)=val_minimo;
                end
                [minimo,posizione_minimo]=min(valori_rotazione);
                best_matching(scanning)=minimo;
            end
            [distanza_minima,posizione_minimo]=min(best_matching);
            beep;
            message=strcat('The nearest fingerprint present in DataBase which matchs input fingerprint is  : ',num2str(posizione_minimo),...
                ' with a distance of : ',num2str(distanza_minima));
            msgbox(message,'DataBase Info','help');
            
        else
            message='DataBase is empty. No check is possible.';
            msgbox(message,'FingerCode DataBase Error','warn');    
        end
        
        end 
        
        
        
%         -----------------------------------------------
%         -------------------------
         if choix==7,
        clc;
        close all;
        figure('Name','Gabor Filter');
        mesh(gabor2d_sub(0,num_disk));
        end
end