MATLAB特征提取代碼

1、精品文檔for i=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);%測(cè)量圖像尺寸參數(shù)GP=zeros(1,256);%預(yù)創(chuàng)建存放灰度出現(xiàn)概率的向量for k=0:255GP(k+1)=length(find(PS=k)/(m*n);%計(jì)算每級(jí)灰度出現(xiàn)的概率，將其存入GP中相應(yīng)位置end%直方圖均衡化S1=ze

2、ros(1,256);for i=1:256for j=1:iS1(i)=GP(j)+S1(i);%計(jì)算 SkendendS2=round(S1*256)+0.5);%將 Sk 歸到相近級(jí)的灰度%圖像均衡化f=PS;for i=0:255f(find(PS=i)=S2(i+1);%將各個(gè)像素歸一化后的灰度值賦給這個(gè)像素endfigure,imshow(f);%邊緣檢測(cè)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(

3、i)=256;j=find(g<0.5);g(j)=0;imshow(g);title('二值法銳化圖像');%中值濾波g=medfilt2(g);g=dither(g);imshow(g);%提取面積，矩形度，圓形度，拉伸度特征。1歡迎下載精品文檔%g=im2bw(g);x,y=size(g);BW = bwperim(g,8); %檢測(cè)邊緣跟蹤，用于計(jì)算周長(zhǎng)%檢測(cè)垂直方向連讀的周長(zhǎng)像素點(diǎn)%P1=0;P2=0;Ny=0; %記錄垂直方向連續(xù)周長(zhǎng)像素點(diǎn)的個(gè)數(shù)for i=1:xfor j=1:yif (BW(i,j)>0)P2=j;if (P2-P1)=1) %判斷是

4、否為垂直方向連續(xù)的周長(zhǎng)像素點(diǎn)Ny=Ny+1;endP1=P2;endendend%檢測(cè)水平方向連讀的周長(zhǎng)像素點(diǎn)P1=0;P2=0;Nx=0; %記錄水平方向連續(xù)周長(zhǎng)像素點(diǎn)的個(gè)數(shù)for j=1:yfor i=1:xif (BW(i,j)>0)P2=i;if (P2-P1)=1) % 判斷是否為水平方向連續(xù)的周長(zhǎng)像素點(diǎn) Nx=Nx+1;endP1=P2;endendendSN=sum(sum(BW); %計(jì)算周長(zhǎng)像素點(diǎn)的總數(shù)Nd=SN-Nx-Ny; %計(jì)算奇數(shù)碼的鏈碼數(shù)目H=max(sum(g); %計(jì)算目標(biāo)的高度W=max(sum(g'); %圖象 g 經(jīng)矩陣轉(zhuǎn)置后，計(jì)算寬度L=s

5、qrt(2)*Nd+Nx+Ny; %計(jì)算周長(zhǎng)%=形態(tài)特征值計(jì)算=%A=bwarea(g); %計(jì)算目標(biāo)的面積R=A/(H*W); %計(jì)算矩形度E=min(H,W)/max(H,W); %計(jì)算伸長(zhǎng)度temp1=A,R,E;%提取不變矩特征。2歡迎下載精品文檔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

6、.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.m002;e.eta20=(m.m20-xbar*m.m10)/m.m002;e.eta02=(m.m02-ybar*m.m01)/m.m002;e.eta30=(m.m30-3*xbar*m.m20+2*xbar2*m.m10)/m.m002.5;e.eta03=(m.m03-3*ybar*m.m02

7、+2*ybar2*m.m01)/m.m002.5;e.eta21=(m.m21-2*xbar*m.m11-ybar*m.m20+2*xbar2*m.m01)/m.m002.5;e.eta12=(m.m12-2*ybar*m.m11-xbar*m.m02+2*ybar2*m.m10)/m.m002.5;phi(1)=e.eta20+e.eta02;phi(2)=(e.eta20-e.eta02)2+4*e.eta112;phi(3)=(e.eta30-3*e.eta12)2+(3*e.eta21-e.eta03)2;phi(4)=(e.eta30+e.eta12)2+(e.eta21+e.eta

8、03)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