PS=imresize(PS,[300,300],'bilinear');%歸一化大小

PS=rgb2gray(PS);

[m,n]=size(PS);?????????????????????? %測量圖像尺寸參數

GP=zeros(1,256);???????????????????? %預創建存放灰度出現概率的向量

for k=0:255

GP(k+1)=length(find(PS==k))/(m*n);? %計算每級灰度出現的概率，將其存入GP中相應位置

end

%直方圖均衡化

S1=zeros(1,256);

for i=1:256

for j=1:i

S1(i)=GP(j)+S1(i);????????????? %計算Sk

end

end

S2=round((S1*256)+0.5);?????????????? %將Sk歸到相近級的灰度

%圖像均衡化

f=PS;

for i=0:255

f(find(PS==i))=S2(i+1);???????? %將各個像素歸一化后的灰度值賦給這個像素

end

figure,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(g<0.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; % 記錄垂直方向連續周長像素點的個數

for i=1:x

for j=1:y

if (BW(i,j)>0)

P2=j;

if ((P2-P1)==1) % 判斷是否為垂直方向連續的周長像素點

Ny=Ny+1;

end

P1=P2;

end

end

end

%檢測水平方向連讀的周長像素點

P1=0;

P2=0;

Nx=0; % 記錄水平方向連續周長像素點的個數

for j=1:y

for i=1:x

if (BW(i,j)>0)

P2=i;

if ((P2-P1)==1) % 判斷是否為水平方向連續的周長像素點

Nx=Nx+1;

end

P1=P2;

end

end

end

SN=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;

end

m.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

《新程序員》：云原生和全面數字化實踐50位技術專家共同創作，文字、視頻、音頻交互閱讀

總結

以上是生活随笔為你收集整理的计算特征数据matlab代码,科学网—MATLAB特征提取代码 - 蒋样明的博文的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。