基于MATLAB的指纹识别系统设计

XXV参考文献[1]王玮.自动指纹识别系统关键技术研究[D].重庆:重庆大学,2007.[2]任春晓.自动指纹识别中若干关键算法的研究[D].山东:山东大学,2011.[3]杨高强,李晟,张新鹏.基于方向场的形变指纹检测[J].应用科学学报,2019,37(6):775-782.[4]李俊.低质量指纹图像分割技术研究[D].浙江:浙江理工大学,2013.[5]徐国庆.指纹识别系统中的特征提取与匹配算法的研究和实现[D].电子科技大学,2006.[6]T.S.ONG,T.B.J.ANDREW,N.C.L.DAVID,etal.FingerprintImagesSegmentationUsingTwoStagesCoarsetoFineDiscriminationTechnique[C].2003.[7]郑康亮.基于IP核复用的指纹识别系统的研究与应用[D].广东:广东工业大学,2008.[8]于焱.指纹识别技术及其应用研究[J].长春大学学报,2002,12(3):43-45.[9]王磊.自动指纹识别系统的算法研究及实现[D].黑龙江:哈尔滨工程大学,2008.[10]任丽君.陈蕾.基于嵌入式应用的指纹识别系统设计与实现[J].山西科技.2010,25(4)：32-33.[11]林华.指纹识别系统预处理技术的研究与实现[D].四川:电子科技大学,2006.[12]任毅.指纹识别系统的研究和实现[D].江苏:南京邮电大学,2019.[13]胡涛,林家骐.基于灰度均衡的指纹图像分割算法[J].计算机工程与应用,2007,43(24):205-207.[14]甘树坤,欧宗瑛,魏鸿磊.基于灰度特性的指纹图像分割算法[J].吉林化工学院学报（自然科学版）,2006,23(1):68-71.[15]SPARROW,M.K.,SPARROW,P.J..TopologicalApproachtotheMatchingofSingleFingerprints:DevelopmentofAlgorithmsforUseonLatentFingermarks[R].[16]A.K.Hrechak,J.A.Mchugh.Automatedfingerprintrecognitionusingstructuralmatching.PatternRecognition,2020,23(8):893-904.[17]金国伦.基于Gabor滤波器的指纹图像增强算法的研究[J].数字技术与应用,2016(5):123-124.[18]陈婧,张苏.基于方向场信息和灰度特征的指纹分割算法研究[J].大理大学学报,2016,1(6):30-32.[19]朱树阔.基于DSP的指纹识别系统研究与设计[D].重庆:重庆大学,2010.[20]孙浩程.基于指纹识别的身份认证系统[D].上海:华东理工大学,2015.[21]黄滨.指纹识别二值化方法的研究[D].浙江:浙江大学,2006.[22]闫志锋.指纹识别匹配算法研究[D].黑龙江:哈尔滨工业大学,2008.[23]MaltoniD,MaioD,JainAK,etal.Handbookoffingerprintrecognition[J]．kybemetes,2004（5-6）：1314.[24]李娟,张宇.采用Nios定制指令的嵌入式系统优化设计[J].单片机与嵌入式系统应用用,2004(1):40-42.[25]祝恩.低质量指纹图像的特征提取与识别技术的研究[D].国防科学技术大学,2005.[26]JainA.K,HongL,PankantiS.Biometrics:PromisingFrontiersforEmergingIdentificationMarket.CommunicationsoftheACM,2020:91-98.[27]AnilK.Jain,PatrickFlynn,ArunA.Ross.HandbookofBiometrics[M].Springer,Boston,MA.[28]JulianFierrez,JavierOrtega-Garcia,DoroteoTorreToledano,JoaquinGonzalez-Rodriguez.Biosecbaselinecorpus:Amultimodalbiometricdatabase[J].PatternRecognition,2006,40(4).[29]Nadra,Sugiyanto,GatotMurtiWibowo.DEVELOPMENTOFFUSIONTECHNIQUESBASEDONMATLABCOMPUTER_FUSIONPROGRAMMING:ComparativeStudybetweenMRIfusionandMatlabfusion[J].JournalofMedicalImagingandRadiationSciences,2022,53(4S1).附录附录A扇形区域修饰程序function[sector_num]=whichsector(index)%ModiofiedbyPRTsinghua@%index剪切后图像的index(剪切后的图像是175x175);length=175;x=rem(index,length);y=floor(index/length);x=x-floor(length/2);y=y-floor(length/2);rad=(x*x)+(y*y);ifrad<144%innerestradius=12(144=12*12)sector_num=36;sector_num;returnendifrad>=5184%outtestradius=72(5184=72*72)sector_num=37;sector_num;returnendifx~=0theta=atan(y/x);elseify>0theta=pi/2;elsetheta=-pi/2;endendifx<0theta=theta+pi;elseiftheta<0theta=theta+2*pi;endendiftheta<0theta=theta+2*pi;endr=floor(rad^0.5);ring=floor((r-12)/20);arc=floor(theta/(pi/6));sector_num=ring*12+arc;附录B求Gabor滤波器与指纹图像的卷积程序function[gaborp_2d]=gabor2d_sub(angle,num_disk)%ModifiedbyPRTsinghua@variance=32;k=10;x=cos(angle*pi/num_disk);y=sin(angle*pi/num_disk);w=(2*pi)/k;p=0;m=0;for(i=-16:1:16)p=p+1;sinp(p)=i*y;cosp(p)=i*x;for(j=-16:1:16)m=m+1;x_s(m)=i;y_s(m)=j;endendp=0;for(j=1:1:33)for(i=1:1:33)p=p+1;xx(p)=sinp(i)+cosp(j);yy(p)=cosp(i)-sinp(j);gaborp(p)=1*exp(-((xx(p)*xx(p))+(yy(p)*yy(p)))/variance)*cos(w*xx(p));if(i-16)*(i-16)+(j-16)*(j-16)<=16*16gaborp_2d(i,j)=gaborp(p);elsegaborp_2d(i,j)=0;endendEnd附录C指纹中心值的计算标记程序function[Outputprint,XofCenter,YofCenter]=centralizing(fingerprint,ctrl)% modifiedbyPRTsinghua@imgN=size(fingerprint,1);imgM=size(fingerprint,2);image=wiener2(fingerprint,[33]);[Gx,Gy]=gradient(image);orientnum=wiener2(2.*Gx.*Gy,[33]);orientden=wiener2((Gx.^2)-(Gy.^2),[33]);W=8;ll=9;orient=zeros(imgN/W,imgM/W);%points=(imgN/W)*(imgM/W);fori=1:1:pointsx=floor((i-1)/(imgM/W))*W+1;y=mod(i-1,(imgN/W))*W+1;numblock=orientnum(y:y+W-1,x:x+W-1);denblock=orientden(y:y+W-1,x:x+W-1);somma_num=sum(sum(numblock));somma_denom=sum(sum(denblock));ifsomma_denom~=0inside=somma_num/somma_denom;angle=0.5*atan(inside);elseangle=pi/2;end%eachblockifangle<0ifsomma_num<0angle=angle+pi/2;elseangle=angle+pi;endelseifsomma_num>0angle=angle+pi/2;endendorient(1+(y-1)/W,1+(x-1)/W)=angle;endbinarize=(orient<pi/2);[bi,bj]=find(binarize);xdir=zeros(W,W);ydir=zeros(W,W);fork=1:1:size(bj,1)i=bj(k);j=bi(k);iforient(j,i)<pi/2x=fix(ll*cos(orient(j,i)-pi/2)/(W/2));y=fix(ll*sin(orient(j,i)-pi/2)/(W/2));xdir(j,i)=i-x;ydir(j,i)=j-y;endendbinarize2=zeros(imgN/W,imgM/W);fori=1:1:size(bj,1)x=bj(i);y=bi(i);if~(xdir(y,x)<1|ydir(y,x)<1|xdir(y,x)>imgM/W|ydir(y,x)>imgN/W)whilebinarize(ydir(y,x),xdir(y,x))>0xtemp=xdir(y,x);ytemp=ydir(y,x);ifxtemp<1|ytemp<1|xtemp>imgM/W|ytemp>imgN/Wbreak;endx=xtemp;y=ytemp;ifxdir(y,x)<1|ydir(y,x)<1|xdir(y,x)>imgM/W|ydir(y,x)>imgN/Wifx-1>0whilebinarize(y,x-1)>0x=x-1;ifx-1<1break;endendendbreak;endendendbinarize2(y,x)=binarize2(y,x)+1;end[temp,y]=max(binarize2(1:end-7,:));[temp2,x]=max(temp);angle=orient(y(x),x)-pi/2;XofCenter=round(x*W-(W/2)-(ll/2)*cos(angle));YofCenter=round(y(x)*W-(W/2)-(ll/2)*sin(angle));Outputprint=binarize2;附录D裁剪图像程序function[CroppedPrint]=cropping(XofCenter,YofCenter,CentralizedPrint)%ModifiedbyPRTsinghua@N=175;M=size(CentralizedPrint,1);imgN=size(CentralizedPrint,1);imgM=size(CentralizedPrint,2);if(YofCenter+30)<=MYofCenter=YofCenter+20;elseYofCenter=M;endX=XofCenter-floor(N/2);Y=YofCenter-floor(N/2);%命令行中能显示中心点%%||%||%|.|y=row%||%||%%x=column%XofCenter%矩阵列%YofCenter%矩阵行%%ifA=123%456%789%B=A(1:2,2:3)%=23%56%createsBbyextractingthefirsttwosandlasttwo%columnsofA%if(YofCenter-floor(N/2)<1)||(YofCenter+floor(N/2)>imgN)||(XofCenter-floor(N/2)<1)||(XofCenter+floor(N/2)>imgM)message='Croppingerror:whentheinputimageiscroppedanerroroccurs:apossibleerrorduringcenterpointdetermination.';msgbox(message,'CroppingError','warn');CroppedPrint=zeros(175);return;elseCroppedPrint=CentralizedPrint(YofCenter-floor(N/2):YofCenter+floor(N/2),XofCenter-floor(N/2):XofCenter+floor(N/2));End附录E扇形归一化输入图像程序function[disk,vector]=sector_norm(image,mode,mix)%ModifiedbyPRTsinghua@%N=175sizeofcroppedimage(175x175)N=175;%NumberofsectorsM=38;size_m=N*N;mean_s=zeros(M,1);varn_s=zeros(M,1);num_s=zeros(M,1);image1=zeros(175,175);Mo=50;Vo=50;for(i=1:1:size_m)tmp=whichsector(i);tmp=tmp+1;if(tmp>=1)mean_s(tmp)=mean_s(tmp)+image(i);num_s(tmp)=num_s(tmp)+1;endendfor(i=1:1:M)mean_s(i)=mean_s(i)/num_s(i);endfor(i=1:1:size_m)tmp=whichsector(i);tmp=tmp+1;if(tmp>=1)varn_s(tmp)=varn_s(tmp)+(image(i)-mean_s(tmp))^2;endendfor(i=1:1:M)varn_s(i)=varn_s(i)/num_s