10-3.m

来自「MATLAB数字图象处理的主要源码.包括图象去噪变换增强恢复识别压缩的各种算法」· M 代码 · 共 24 行

%已知一个二值图像1000×1000，中间为一个50×300的白色区域，其它为黑色。
%下面对这个图像来研究傅立叶变换的旋转性质
f = zeros(1000,1000);
f(350:649,475:524) = 1;
%图像定义
subplot(221)
imshow(f,'notruesize');
title('原始图像');
subplot(222)
F=fftshift(abs(fft2(f)));
imshow(F,[-1 5],'notruesize');
title('原始图像的傅立叶变换频谱');
subplot(223)
f = zeros(1000,1000);
f(350:649,475:524) = 1;
f = imrotate(f,45,'bilinear','crop');
%以图像中心为原点旋转45度
imshow(f,'notruesize');
title('图像正向旋转45度');
subplot(224)
F=fftshift(abs(fft2(f)));
imshow(F,[-1 5],'notruesize');
title('图像旋转45度后的傅立叶变换频谱');