feature-selection.rar_feature selection_图像孤立点_提取图像曲线_曲线特征提取_特征 m

clear; I=imread('3.tif'); subplot(1,2,1); imshow(I); title('原始图像'); Hist=zeros(1,256); [m,n]=size(I); for i=1:m for j=1:n Hist(I(i,j)+1)=Hist(I(i,j)+1)+1; end end J=zeros(1,9); x=0; for k=30:255 x=x+k*Hist(k+1); end y=0; for k=30:255 y=y+Hist(k+1); end u=round(x/y); for i=1:m for j=1:n if (u-30)<I(i,j)&&I(i,j)<(u+30) I(i,j)=255; else I(i,j)=0; end end end subplot(1,2,2); imshow(I); title('二值化后的图像'); I=double(I); for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==0 J(1)=J(1)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255 J(2)=J(2)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==510 J(3)=J(3)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255*3 J(4)=J(4)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255*4 J(5)=J(5)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255*5 J(6)=J(6)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255*6 J(7)=J(7)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255*7 J(8)=J(8)+1; else end end end for i=2:(m-1) for j=2:(n-1) if (I(i-1,j-1)+I(i-1,j)+I(i-1,j+1)+I(i,j-1)+I(i+1,j-1)+I(i+1,j)+I(i+1,j+1)+I(i,j+1))==255*8 J(9)=J(9)+1; else end end end