其他分享
首页 > 其他分享> > 【图像识别】基于模板匹配实现指纹识别

【图像识别】基于模板匹配实现指纹识别

作者:互联网

指纹识别技术主要分三个步骤:指纹预处理、特征提取、指纹分类与匹配。

无论是指纹分类还是指纹匹配,都需要提取指纹的有效特征,而特征提取的性能很大程度上要依赖于指纹图像的质量。在实际应用中,由于采集条件和采集设备的因素,采集到的指纹图像质量比较差,容易导致很多问题,影响后续处理的效果。因此,指纹图像的预处理是关系到指纹识别系统性能好坏的一个关键。

预处理主要分以下四个步骤:

    1. 指纹图像灰度归一化和均衡化

      归一化的目的,在于消除指纹采集过程中由于传感器自身的噪声以及因为手指压力不同而造成的灰度差异,将指纹图像的对比度和灰度调整到一个固定的级别上,为后续处理提供一个较为统一的图像规格。均衡化是对图像中像素个数多的灰度级进行展宽,对像素个数少的灰度级进行缩减。

    2. 指纹图像分割

      其目标就是根据特征提取的需要,把指纹图像中质量很差、在后续处理中很难恢复的图像区域与有效区域分开,使后续处理能够集中在有效区域;能提高特征提取的精确度;能大大减少指纹预处理的时间。

    3. 指纹图像二值化

      二值化的目的是把灰度指纹图像变成0、1取值的二值图像

    4. 二值化后处理及细化

      由于灰度滤波的不完全性,而且在二值化过程中有时会引入新的噪声,需要对图像进行滤波处理。采用加权中值滤波的方法,根据前景点的不同方向选用不同的权值模板进行滤波,以便于消除纹线上的孔洞和缺口。二值化后的纹线仍然有一定

      宽度,需要细化为单个像素宽度的骨架。细化算法很多,这里采用骨架提取技术。

    5. clc;
      
      close all;
      
      global immagine n_bands h_bands n_arcs h_radius h_lato n_sectors matrice num_disk
      
      %immagine 双精度类型的灰度图
      
      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);
      
      end
      
      num_disk=8;%8个方向
      
      % 1--> add database
      
      % 0--> recognition
      
      % ok=0;
      
      chos=0;
      
      possibility=6;%有6个菜单,分别为'选择图像并加入数据库','指纹识别','删除数据库','可视化指纹图像','可视化Gabor滤波','退出'
      
       
      
      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 chos~=possibility,
      
          chos=menu('指纹识别系统','选择图像并加入数据库','指纹识别','删除数据库','可视化指纹图像','可视化Gabor滤波','退出');
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          % 计算指纹编码并添加到数据库
      
          if chos==1
      
              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)'},'选在灰度图');%打开一个灰度图
      
                  if namefile~=0 %文件存在
      
                      [img,map]=imread(strcat(pathname,namefile));%读取文件
      
                      selezionato=1;
      
                  else
      
                      disp('选择灰度图');
      
                  end
      
                  if (any(namefile~=0) && (~isgray(img)))%文件不存在或不是灰度图
      
                      disp('选择灰度图');
      
                      selezionato=0;
      
                  end
      
              end
      
              
      
              immagine=double(img);
      
              
      
              if isa(img,'uint8')%如果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);%Gabor滤波
      
                  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
      
              % 增加指纹编号到数据库
      
              if (exist('database.dat')==2)
      
                  load('database.dat','-mat');
      
                  fingerprint_number=fingerprint_number+1;
      
                  data{fingerprint_number,1}=finger_code1;
      
                  data{fingerprint_number,2}=finger_code2;
      
                  save('database.dat','data','fingerprint_number','-append');
      
              else
      
                  fingerprint_number=1;
      
                  data{fingerprint_number,1}=finger_code1;
      
                  data{fingerprint_number,2}=finger_code2;
      
                  save('database.dat','data','fingerprint_number');
      
              end
      
              
      
              message=strcat('指纹增加成功。编号:',num2str(fingerprint_number));
      
              msgbox(message,'指纹编数据库','信息');
      
          end % chos 1
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          % 指纹识别
      
          if chos==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)'},'选择灰度图');
      
                  if namefile~=0
      
                      [img,map]=imread(strcat(pathname,namefile));
      
                      selezionato=1;
      
                  else
      
                      disp('选择灰度图');
      
                  end
      
                  if (any(namefile~=0) && (~isgray(img)))
      
                      disp('选择灰度图');
      
                      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);%归一化
      
              
      
              % 存储每个特征向量d的入口
      
              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('database.dat')==2)
      
                  load('database.dat','-mat');
      
                  %---- 分配内存 -----------------------------------
      
                  %...
      
                  vettore_a=zeros(num_disk*n_sectors,1);
      
                  vettore_b=zeros(num_disk*n_sectors,1);
      
                  best_matching=zeros(fingerprint_number,1);
      
                  valori_rotazione=zeros(n_arcs,1);
      
                  % 开始检查 ---------------------------------------
      
                  for scanning=1:fingerprint_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('与数据库中最相似的指纹是 : ',num2str(posizione_minimo)');% 具有',num2str(distanza_minima),'个实例');
      
                  msgbox(message,'数据库信息','信息');            
      
              else
      
                  message='数据库为空. 不能匹配.';
      
                  msgbox(message,'指纹编号数据库错误','错误');    
      
              end
      
              
      
          end % chos 2
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          % 删除数据库  
      
          if chos==3
      
              clc;
      
              close all;
      
              if (exist('database.dat')==2)
      
                  button = questdlg('你确定要删除数据库吗?');
      
                  if strcmp(button,'Yes')
      
                      delete('database.dat');
      
                      msgbox('数据库删除成功。','数据库删除','信息');
      
                  end
      
              else
      
                  warndlg('数据库为空!','警告')
      
              end
      
          end % chos 3
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          % 可视化指纹图像    
      
          if chos==4
      
              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('选择灰度图');
      
                  end
      
                  if (any(namefile~=0) && (~isgray(img)))
      
                      disp('选择灰度图');
      
                      selezionato=0;
      
                  end
      
              end
      
              figure('name','选择图像');
      
              imshow(img);
      
          end % chos 4
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          %--------------------------------------------------------------------------
      
          % 可视化Gabor滤波
      
          if chos==5
      
              clc;
      
              close all;
      
              figure('name','Gabor滤波');
      
              mesh(gabor2d_sub(0,num_disk));
      
          end % chos 5 
      
      end % end while
      

    6. 完整代码或者代写添加QQ1575304183

标签:end,img,指纹识别,num,灰度,fingerprint,图像识别,disk,模板
来源: https://www.cnblogs.com/ttmatlab/p/14746172.html