untitled3.m

一个指纹识别系统

img2=imread('1_1.tif');
img3=double(img2);
[out,xc2,yc2]=centralizing(img3,0);
[CroppedPrint]=cropping(xc2,yc2,img3);
for ( i=1:1:175*175)
    tmp=CroppedPrint(i);
    tmp=double(tmp);
  CroppedPrint(i)=whichsector(i);
     
    if (CroppedPrint(i)==36 | CroppedPrint(i)==37)
        CroppedPrint(i)=tmp/255;
    else 
     
        CroppedPrint(i)=tmp/64;
    end
    
end
num_disk=8
[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);  
finger_code1{angle+1}=vector(1:36);
end