fcheck.m

一个指纹识别系统

function Fcheck
global   img img3
global text1 text2 text3 push1 push2 push3 push4 push5 push6 push7 push8 push9 push10 push11 push12 push13
global push14 push15 push16 push17 push18 push19 push20 push21 push22 push23 push24 push25 push26 push27 push28 push29 push30 push31 push32



N=175;
num_disk=8;


[out,xc2,yc2]=centralizing(img3,0);
[CroppedPrint]=cropping(xc2,yc2,img3);

[NormalizedPrint,vector]=sector_norm(CroppedPrint,0,1);

for (angle=0:1:num_disk-1)    
    gabor=gabor2d_sub(angle,num_disk);
    z2=gabor;
    z1=NormalizedPrint;
    z1x=size(z1,1);
    z1y=size(z1,2);
    z2x=size(z2,1);
    z2y=size(z2,2);    
    ComponentPrint=real(ifft2(fft2(z1,z1x+z2x-1,z1y+z2y-1).*fft2(z2,z1x+z2x-1,z1y+z2y-1)));    
    px=((z2x-1)+mod((z2x-1),2))/2;
    py=((z2y-1)+mod((z2y-1),2))/2;
    ComponentPrint=ComponentPrint(px+1:px+z1x,py+1:py+z1y);        
    [disk,vector]=sector_norm(ComponentPrint,1,0);    
    %img = double(ComponentPrint)/graylevmax;
    %img1 = double(disk)/51200;
    finger_code{angle+1}=vector(1:36);
end


% 输入指纹的编码已经计算，检查数据库
if (exist('fp_database.dat')==2)
    load('fp_database.dat','-mat');
    %---- 申请内存 -----------------------------------
    ruoto1=zeros(36,1);
    ruoto2=zeros(36,1);
    vettore_d1=zeros(12,1);
    vettore_d2=zeros(12,1);
    best_matching=zeros(fp_number,1);
    % 开始检查 ---------------------------------------
    for scanning=1:fp_number
        fcode1=data{scanning,1};
        fcode2=data{scanning,2};
        for rotazione=0:1:11
            d1=0;
            d2=0;
            for disco=1:8
                f1=fcode1{disco};
                f2=fcode2{disco};
                % ora ruoto f1 ed f2  della rotazione ciclica ----------
                for old_pos=1:12
                    new_pos=mod(old_pos+rotazione,12);
                    if (new_pos==0)
                        new_pos=12;
                    end
                    ruoto1(new_pos)=f1(old_pos);
                    ruoto1(new_pos+12)=f1(old_pos+12);
                    ruoto1(new_pos+24)=f1(old_pos+24);
                    ruoto2(new_pos)=f2(old_pos);
                    ruoto2(new_pos+12)=f2(old_pos+12);
                    ruoto2(new_pos+24)=f2(old_pos+24);
                end
                %-------------------------------------------------------
                d1=d1+norm(finger_code{disco}-ruoto1);
                d2=d2+norm(finger_code{disco}-ruoto2);                
            end
            vettore_d1(rotazione+1)=d1;
            vettore_d2(rotazione+1)=d2;
        end
        [min_d1,pos_min_d1]=min(vettore_d1);
        [min_d2,pos_min_d2]=min(vettore_d2);
        if min_d1<min_d2
            minimo=min_d1;
        else
            minimo=min_d2;
        end
        best_matching(scanning)=minimo;
    end
    [distanza_minima,posizione_minimo]=min(best_matching);
    beep;
    message3=strcat('数据库中与输入指纹最匹配的是：',num2str(posizione_minimo),...
                        ' 其距离为： ',num2str(distanza_minima));
   
    %-------------------------------------------------------
    
else
    message3='数据库为空，无法检查';
     
end

            set(text3,'visible','on','string',message3);