test_img111213.m

对于普通图象的一些频域的分析小程序,可以对图象处理参考作用.

clear;
I=imread('img21.bmp');
I1=rgb2gray(I);
figure;imshow(I1);
title('img21')
J1=fft2(I1);
K1=fftshift(J1);

figure;imshow(log(abs(K1)),[]);
title('频谱img21');
G2=imread('img22.bmp');
I2=rgb2gray(G2);
figure;imshow(I2);
title('img22')
J2=fft2(I2);
K2=fftshift(J2);
figure;imshow(log(abs(K2)),[]);
title('频谱img22');
G3=imread('img23.bmp');
I3=rgb2gray(G3);
figure;imshow(I3);
title('img23')
J3=fft2(I3);
K3=fftshift(J3);
figure;imshow(log(abs(K3)),[]);
title('频谱img23');
%以下实现功率谱的比较。实验时更换K和I下标达到对不同图像处理的目的

a=0;
b=0;
T=abs(K3)^2;
[n,m]=size(I3);

for i=1:n
    for j=1:m
        if (n/2-i)^2+(m/2-j)^2<80^2  %确定圆的范围
            a=a+T(i,j);
        end
    b=b+T(i,j);    
    end
end
e=a/b;