%已知一個二值圖像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度后的傅立葉變換頻譜');