ffingercode.m

介绍指纹编码算法

function Ffingercode

global img img2  img3 xc yc ComponentPrint graylevmax xc2 yc2  gabor1  gabor2  gabor3  gabor4  gabor5  gabor6  gabor7  gabor8 z1 z1x z1y ComponentPrint1 ComponentPrint2 ComponentPrint3 ComponentPrint4 ComponentPrint5 ComponentPrint6 ComponentPrint7 ComponentPrint8
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;


[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);    
    img20 = double(ComponentPrint)/256;
    
    finger_code1{angle+1}=vector(1:36);
end


img21=imrotate(img20,22.5/2);
imgN=size(img21,1);
imgM=size(img21,2);
modN=mod(imgN,8);
modM=mod(imgM,8);
img22=double(img21(modN+1:imgN,modM+1:imgM));

[out,xc3,yc3]=centralizing(img22,0);
[CroppedPrint]=cropping(xc3,yc3,img22);
[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);    
    ComponentPrint9=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;
    ComponentPrint10=ComponentPrint9(px+1:px+z1x,py+1:py+z1y);        
    [disk,vector]=sector_norm(ComponentPrint10,1,0);    
    %img = double(ComponentPrint)/graylevmax;
    %img1 = double(disk)/51200;
    finger_code2{angle+1}=vector(1:36);
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('已成功加入 ',num2str(fp_number),'号模板');        
        set(text1,'visible','on','string',message);
        %模板库信息显示
        message2=strcat('系统现有 ',num2str(fp_number),'个模板');
        set(text2,'visible','on','string',message2);