finger process.rar_6GZ_指纹

close all name='17.bmp'; img_in=imread(name); imshow(img_in);title('original picture'); img=img_in; [m,n]=size(img); img=double(img); for i=1:m for j=1:n if img(i,j)<240 img(i,j)=255; end end end sobelx=[-1 -2 -1;0 0 0;1 2 1]; sobely=[-1 0 1;-2 0 2;-1 0 1]; xx=filter2(sobelx,img); % figure, imshow(xx, []); yy=filter2(sobely,img); % figure, imshow(yy, []); imgto=abs(xx)+abs(yy); % figure, imshow(imgto);title('sobel'); se=ones(3); for t=1:5 imgto=imdilate(imgto,se); end figure,imshow(imgto);title('dilate'); % figure,imhist(imgto); img_temp=255-img_in; for i=1:m for j=1:n if imgto(i,j)==0 img_temp(i,j)=0; end end end img_temp=255-img_temp; figure,imshow(img_temp);title('temp image'); img1=img_temp; img1=medfilt2(img1); % img2=double(img2); img2=edge(img1,'Sobel',0.25); % figure,imshow(img3); img_edge=edge(img1,'Roberts',0.30); % figure,imshow(img4); img_edge=(img2+img_edge)/2; img_edge=im2bw(img2,0.4); img_edge=im2bw(img_edge,0.4); se=ones(4); for x=1:3 img_edge=imdilate(img_edge,se); end img_edge=medfilt2(img_edge); figure,imshow(img_edge);title('After dilate'); img_out=255-img_in; for i=1:m for j=1:n if img_edge(i,j)==0 img_out(i,j)=0; end end end img_out=255-img_out; figure,imshow(img_out);title('final image');