方形环向功率谱.m

图像经过二维傅立叶频谱分析后

%   L503频谱分析
f=imread('E:\TUTU\L503.JPG');   %  读取图像；
h=rgb2gray(f);                 %   转化为灰度图；
h1=h(1:300,1:300);             %   截缺取300*360大小的图像；
h11=zeros(500,500);
for i=1:300;
   for j=1:300;
       h11(101:400,101:400)=h1(i,j);   % 将原矩阵扩展成500*500矩阵；
   end
end
e=im2bw(h1);                   %   转化为二值图像                      
F=fft2(e);                      %  快速傅立叶变化（离散函数）
S=abs(F);                      %  获取傅立叶频谱
p = unwrap(S); 
FC=fftshift(S);                  %  将变换的原点转移到频率矩形的中心
figure;
colormap(gray);                % 归一化处理；
FFC=255*FC/max(FC(:));
S2=log10(abs(FFC));              %  使用对数变换进行视觉增强
SS=FFC^2;                      %   
E=sum(sum(SS));               %  
for m=1:150;
    ER(m)=0;
for i=1:size(SS,1);
    for j=1:size(SS,2);
        if  (m-1)<abs(i-0.5-(size(SS,1)/2))&&abs(i-0.5-(size(SS,1)/2))<=m;
            if (m-1)<abs(j-0.5-(size(SS,2/2)))&&abs(j-0.5-(size(SS,2)/2))<=m;
            ER(m)=SS(i,j)+ER(m);    %   正方环向能量谱
            end
        end
     end 
end
end
ER1=sum(ER);




E1=sum(sum(SS(133:147,171:189)));
fh1=E1/E;
E2=sum(sum(SS(126:154,162:198)))-E1;
fh2=E2/E;
E3=sum(sum(SS(119:161,153:207)))-E2-E1;
fh3=E3/E;
E4=sum(sum(SS(112:168,144:216)))-E3-E2-E1;
fh4=E4/E;
E5=sum(sum(SS(105:175,135:225)))-E4-E3-E2-E1;
fh5=E5/E;
E6=sum(sum(SS(98:182,126:234)))-E5-E4-E3-E2-E1;
fh6=E6/E;
E7=sum(sum(SS(91:189,117:243)))-E6-E5-E4-E3-E2-E1;
fh7=E7/E;
E8=sum(sum(SS(84:196,108:252)))-E7-E6-E5-E4-E3-E2-E1;
fh8=E8/E;
E9=sum(sum(SS(77:203,99:261)))-E8-E7-E6-E5-E4-E3-E2-E1;
fh9=E9/E;
E10=sum(sum(SS(70:210,90:270)))-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh10=E10/E;
E11=sum(sum(SS(63:217,81:279)))-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh11=E11/E;
E12=sum(sum(SS(56:224,72:288)))-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh12=E12/E;
E13=sum(sum(SS(49:231,63:297)))-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh13=E13/E;
E14=sum(sum(SS(42:238,54:306)))-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh14=E14/E;
E15=sum(sum(SS(35:245,45:315)))-E14-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh15=E15/E;
E16=sum(sum(SS(28:252,36:324)))-E15-E14-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh16=E16/E;
E17=sum(sum(SS(21:259,27:333)))-E16-E15-E14-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh17=E17/E;
E18=sum(sum(SS(14:266,18:342)))-E17-E16-E15-E14-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh18=E18/E;
E19=sum(sum(SS(7:273,9:351)))-E18-E17-E16-E15-E14-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh19=E19/E;
E20=sum(sum(SS(1:280,1:360)))-E19-E18-E17-E16-E15-E14-E13-E12-E11-E10-E9-E8-E7-E6-E5-E4-E3-E2-E1;
fh20=E20/E;
fh=fh1+fh2+fh3+fh4+fh5+fh6+fh7+fh8+fh9+fh10+fh11+fh12+fh13+fh14+fh15+fh16+fh17+fh18+fh19+fh20;