MATLAB特征提取代码

1、精品文档for i=1:26f=strcat(D:bishe, nu m2str(i);image=strcat(f,.jpg);PS=imread(image);PS=imresize(PS,300,300,bili nea。;归一化大小PS=rgb2gray(PS);m, n=size(PS);%测量图像尺寸参数GP=zeros(1,256);%预创建存放灰度出现概率的向量for k=0:255GP(k+1)=le ngth(fi nd(PS=k)/(m* n);%计算每级灰度出现的概率，将其存入GP中相应位置end%t方图均衡化S1= zeros(1,256);for i=1:256fo

2、r j=1:iS1(i)=GP(j)+S1(i);endendS2=rou nd(S1*256)+0.5);%图像均衡化%计算Sk%将Sk归到相近级的灰度f=PS;for i=0:255f(fin d(PS=i)=S2(i+1); endfigure,imshow(f);泌缘检测f=edge(f,ca nn y,0.25); imshow(f);%T值法锐化图像%将各个像素归一化后的灰度值赋给这个像素f=double(f);x,y=gradie nt(f); g=sqrt(x.*x+y.*y); i=fin d(g=0.5); g(i)=256;j=fin d(g0)P2=j;if (P2-P

3、1)=1) %判断是否为垂直方向连续的周长像素点Ny=Ny+1;endP仁 P2;endendend% 佥测水平方向连读的周长像素点P仁0;P2=0;Nx=0; % 记录水平方向连续周长像素点的个数for j=1:yfor i=1:xif (BW(i,j)0)P2=i;if (P2-P1)=1) %判断是否为水平方向连续的周长像素点Nx=Nx+1;endP仁 P2;endendendSN=sum(sum(BW); %计算周长像素点的总数Nd=SN-Nx-Ny; %计算奇数码的链码数目H=max(sum(g); % 计算目标的高度W=max(sum(g); % 图象g经矩阵转置后，计算宽度L=s

4、qrt (2)*Nd+Nx+Ny; %计算周长%=形态特征值计算=%A=bwarea(g); %计算目标的面积R=A/(H*W); % 计算矩形度E=mi n( H,W)/max(H,W); %计算伸长度temp1=A,R,E;他取不变矩特征精品文档M,N=size(g);x,y=meshgrid(1:N,1:M);x=x(:);y=y(：);g=g(：);m.mOO=sum(g);if(m.mOO=O)m.mOO=eps;endm.m10=sum(x.*g);m.m01=sum(y.*g);m.m20=sum(x.A2.*g);m.m02=sum(y.A2.*g);xbar=m.m1O/m.

5、mOO;ybar=m.mO1/m.mOO;e.eta11=(m.m11-ybar*m.m10)/m.m00A2;e.eta20=(m.m20-xbar*m.m10)/m.m00A2;e.eta02=(m.m02-ybar*m.m01)/m.mO0A2;e.eta30=(m.m30-3*xbar*m.m20+2*xbae.eta03=(m.m03-3*ybar*m.m02+2*ybae.eta2 仁(m.m21-2*xbar*m.m11-ybar*m.m20+2*xbaphi(2)=(e.eta20-e.eta02)A2+4*e.eta11A2;phi(3)=(e.eta30-3*e.eta12

6、)A2+(3*e.eta21-e.eta03)A2;phi(4)=(e.eta30+e.eta12)A2+(e.eta21+e.eta03)A2;phi(5)=(e.eta30-3*e.eta12)*(e.eta30+e.eta12)*(e.eta30+e.eta12)A2-3*(e.eta21 +e.eta03)A2+(3*e.eta21-e.eta03)*(e.eta21+e.eta03)*(3*(e.eta30+e.eta12)A2-(e.eta21+e.eta03)A2);phi(6)=(e.eta20-e.eta02)*(e.eta30+e.eta12)A2-(e.eta21+e.eta03)A2)+4*e.eta11 *(e.eta30+e.eta12)*(e.eta21+e.eta03);phi(7)=(3*e.eta21-e.eta03)*(e.eta30+e.eta12)*(e.eta30+e.eta12)A2-3*(e.eta21 +e.eta03)A2)+(3*e.eta12-e.eta30)*(e.eta21