二维离散傅里叶变换的旋转性
源码如下:
>> clear
>> %构造原始图像
I=zeros(256,256);
I(28:228,108:148)=1;
subplot(2,2,1),imshow(I);
%求原始图像的傅里叶频谱
J=fft2(I);
F=abs(J);
J1=fftshift(F);
subplot(2,2,2),imshow(J1,[5,50]);
%对原始图像进行旋转
K=imrotate(I,315,'bilinear','crop');
subplot(2,2,3),imshow(K);
%求旋转后图像的傅里叶频谱
K1=fft2(K);
F=abs(K1);
K2=fftshift(F);
subplot(2,2,4),imshow(K2,[5,50]);
%加标题,可以加在subplot之后。
subplot(2,2,1),title('原始图像');
subplot(2,2,2),title('原图像的傅里叶频谱');
subplot(2,2,3),title('旋转后的图像');
subplot(2,2,4),title('旋转后图像的傅里叶频谱');
代码结果图像: