掌纹识别 定位分割 形态学腐蚀膨胀算法 matlab源码

《掌纹识别 定位分割 形态学腐蚀膨胀算法 matlab源码》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、掌纹识别定位分割形态学腐蚀膨胀算法matlab源码掌纹识别定位分割形态学腐蚀膨胀算法matlab源码西南交通大学温长芝2008年10月02日星期四14:09%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Function:掌纹图像预处理程序之形态学腐蚀膨胀算法%Author:温长芝%Date:2006/12/7%Copyright:西南交通大学信号与信息处理四川省重点实验室%%%%%%%%%%%%%%%%%%%%

2、%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%第二种算法采用的是形态学中的腐蚀和膨胀运算。首先，对图像进行低通滤波，消除细小噪声的干扰。%在滤波时，没有采用传统的中值滤波器，而是采用顺序统计滤波器ordfilt2。%imm=ordfilt2(im,n,ones(N,M));%上式表示，将ones(N,M)范围内的像素值依次排序，用第n个像素点的值代替该点的像素值，%此种滤波器与中值滤波器一样不会使图像的边缘变的模糊，在滤除噪声的同时很好的保留了图

3、像的边缘。%之后，设定门限值，将掌纹图像转化为二值图像。由经验得到门限值应为15（或20）。%对二值化后的图像先进行膨胀运算、再进行腐蚀运算，将得到的结果与原图像相减，图像上所剩的点即为角点。%先膨胀运算，再腐蚀运算，SE表示在进行膨胀与腐蚀时所用的结构元素。%SE的确定是此次角点检测的关键所在。SE的形状应与要检测的目标的形状相似。%知要检测的目标的形状为弧形，所以SE的形状也应该为弧形。%找到两角点的坐标，并将其旋转到一条垂直线上，旋转后的直线作为坐标系的纵坐标。%在图像的中心位置切割大小为128*128的图

4、像。%这种方法切割出来的掌纹图像的旋转平移较小，而且算法简单，计算量小，比较实用。%clear,closeall;clc;fori=6forj=1%im_1=imread(strcat('PolyU_',num2str(i),'_F_',num2str(j),'.bmp'));im_2=imread(strcat('PolyU_',num2str(i),'_S_',num2str(j),'.bmp'));%输入图像%out_1=Cropedimage(im_1);%out_2=Cropedimage(im_2);

5、SE=[111111110011111100001111000000110000000100000011000011110011111111111111];%定义的结构元素im=im_2;subplot(3,3,1);imshow(im);title('(a)原始掌纹图像');imm=ordfilt2(im,200,ones(10,40));%中值滤波imwrite(imm,'ordfilt2.bmp');subplot(3,3,3);imshow(imm);title('(b)顺序统计滤波后的图像');imag

6、e(find(imm<15))=0;%二值化image(find(imm>15))=255;image=reshape(image,284,384);imwrite(image,'binaryimage.bmp');subplot(3,3,2);imshow(image);title('(c)二值化后的图像');image1=imdilate(image,SE);%进行膨胀腐蚀运算subplot(3,3,5);imshow(image1);title('(d)膨胀后的图像');image1=imerode(ima

7、ge1,SE);subplot(3,3,8);imshow(image1);title('(e)腐蚀后的图像');imwrite(image1,'dilate_erode.bmp')image1=image1-image;%与原图像相减得到角点image1(find(image1<1))=0;image1(find(image1>1))=255;subplot(3,3,6);imshow(image1);title('(f)相减得到角点的图像');imwrite(image1,'keypoints.bmp');[

8、Y1,X1]=find(image1(1:100,1:150));[Y2,X2]=find(image1(150:250,1:150));%x1=max(X1);y1=max(Y1);x2=max(X2);y2=max(Y2)+150;angle=atan((x1-x2)/(y2-y1))*360/2/pi;%计算图像应旋转的角度im=imrotate(im,angle,