MATLAB特征提取代码

MATLAB特征提取代码fori=1:26f=strcat('D:\bishe\',num2str(i));image=strcat(f,'.jpg');PS=imread(image);PS=imresize(PS,[300,300],'bilinear');%归一化大小PS=rgb2gray(PS);[m,n]=size(PS);%测量图像尺寸参数GP=zeros(1,256);%预创建存放灰度出现概率的向量fork=0:255GP(k+1)=length(find(PS==k))/(m*n);%计算每级灰度出现的概率，将其存入GP中相应位置end%直方图均衡化S1=zeros(1,256);fori=1:256forj=1:iS1(i)=GP(j)+S1(i);%计算SkendendS2=round((S1*256)+0.5);%将Sk归到相近级的灰度%图像均衡化f=PS;fori=0:255f(find(PS==i))=S2(i+1);%将各个像素归一化后的灰度值赋给这个像素endfigure,imshow(f);%边缘检测f=edge(f,'canny',0.25);imshow(f);%二值法锐化图像f=double(f);[x,y]=gradient(f);g=sqrt(x.*x+y.*y);i=find(g=0.5);g(i)=256;j=find(g0.5);g(j)=0;imshow(g);title('二值法锐化图像');%中值滤波g=medfilt2(g);g=dither(g);imshow(g);%提取面积，矩形度，圆形度，拉伸度特征%g=im2bw(g);[x,y]=size(g);BW=bwperim(g,8);%检测边缘跟踪，用于计算周长%检测垂直方向连读的周长像素点%P1=0;P2=0;Ny=0;%记录垂直方向连续周长像素点的个数fori=1:xforj=1:yif(BW(i,j)0)P2=j;if((P2-P1)==1)%判断是否为垂直方向连续的周长像素点Ny=Ny+1;endP1=P2;endendend%检测水平方向连读的周长像素点P1=0;P2=0;Nx=0;%记录水平方向连续周长像素点的个数forj=1:yfori=1:xif(BW(i,j)0)P2=i;if((P2-P1)==1)%判断是否为水平方向连续的周长像素点Nx=Nx+1;endP1=P2;endendendSN=sum(sum(BW));%计算周长像素点的总数Nd=SN-Nx-Ny;%计算奇数码的链码数目H=max(sum(g));%计算目标的高度W=max(sum(g'));%图象g经矩阵转置后，计算宽度L=sqrt(2)*Nd+Nx+Ny;%计算周长%====形态特征值计算===%A=bwarea(g);%计算目标的面积R=A/(H*W);%计算矩形度E=min(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.m00=sum(g);if(m.m00==0)m.m00=eps;endm.m10=sum(x.*g);m.m01=sum(y.*g);m.m11=sum(x.*y.*g);m.m20=sum(x.^2.*g);m.m02=sum(y.^2.*g);m.m30=sum(x.^3.*g);m.m03=sum(y.^3.*g);m.m12=sum(x.*y.^2.*g);m.m21=sum(x.^2.*y.*g);xbar=m.m10/m.m00;ybar=m.m01/m.m00;e.eta11=(m.m11-ybar*m.m10)/m.m00^2;e.eta20=(m.m20-xbar*m.m10)/m.m00^2;e.eta02=(m.m02-ybar*m.m01)/m.m00^2;e.eta30=(m.m30-3*xbar*m.m20+2*xbar^2*m.m10)/m.m00^2.5;e.eta03=(m.m03-3*ybar*m.m02+2*ybar^2*m.m01)/m.m00^2.5;e.eta21=(m.m21-2*xbar*m.m11-ybar*m.m20+2*xbar^2*m.m01)/m.m00^2.5;e.eta12=(m.m12-2*ybar*m.m11-xbar*m.m02+2*ybar^2*m.m10)/m.m00^2.5;phi(1)=e.eta20+e.eta02;phi(2)=(e.eta20-e.eta02)^2+4*e.eta11^2;phi(3)=(e.eta30-3*e.eta12)^2+(3*e.eta21-e.eta03)^2;phi(4)=(e.eta30+e.eta12)^2+(e.eta21+e.eta03)^2;phi(5)=(e.eta30-3*e.eta12)*(e.eta30+e.eta12)*((e.eta30+e.eta12)^2-3*(e.eta21+e.eta03)^2+(3*e.eta21-e.eta03)*(e.eta21+e.eta03)*(3*(e.eta30+e.eta12)^2-(e.eta21+e.eta03)^2));phi(6)=(e.eta20-e.eta02)*((e.eta30+e.eta12)^2-(e.eta21+e.eta03)^2)+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)^2-3*(e.eta21+e.eta03)^2)+(3*e.eta12-e.eta30)*(e.eta21+e.eta03)*(3*(e.eta30+e.eta12)^2-(e.eta21+e.eta03)^2);temp2=abs(log(phi));%包含七个特征值temp=[temp1,temp2]end