FFT2 图像二维FFT含义解释
一维情况时:
clc; clear;
x = zeros(256,1);
n=100; nw=1;
for i=n:n+nw
x(i) = 1;
end
f = fft(x);
figure(1); subplot(2,1,1); plot(x); axis([0 300 -1 2]);
title(['Signal, width=', num2str(nw)]);
figure(1); subplot(2,1,2); plot(abs(f(1:128)))
title('Spectrum');
矩形脉冲宽度窄时,高频分量多,宽度宽时,高频分量减少
clc; clear;
A = zeros(256, 256);
A = im2uint8(A);
m=100; n=100;
mw=50; nw=1;
for i=m:m+mw
for j=n:n+nw
A(i,j) = 255;
end
end
figure(1); subplot(2,1,1); imshow(A), colorbar
title(['Image, pulse width=', num2str(nw), ', pulse height=', num2str(mw)]);
B = fft2(A);
C = fftshift(B);
figure(1);subplot(2,1,2); imshow((abs(C)), []), colorbar
title('Image spectrum');
只是考察了水平方向的变换,脉冲水平宽度窄时,水平高频分量多,水平宽度宽时,水平高频分量减少
垂直方向可以类推。
目前还不明白斜对角的频谱表示什么含义?
clc; clear;
A = zeros(256, 256);
A = im2uint8(A);
m=100; n=100;
r=1;
for i=1:256
for j=1:256
if( sqrt( (i-m)^2+(j-n)^2 ) <= r )
A(i,j) = 255;
end
end
end
figure(1); subplot(2,1,1); imshow(A), colorbar
title(['Image, pulse center=', num2str(n), ', pulse radius=', num2str(r)]);
B = fft2(A);
C = fftshift(B);
figure(1);subplot(2,1,2); imshow((abs(C)), []), colorbar
title('Image spectrum');
可以看出,FFT2虽然只是在水平和垂直方向的频谱分析,但是也显示了斜方向的频谱情况,如果使用一个椭圆脉冲,应该可以更清楚的看出在斜方向的频谱情况。