10-2.m

matlab图像处理源码

%已知一个二值图像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,800:849) = 1;
imshow(f,'notruesize');
title('X轴方向移动后的图像');
subplot(224)
F=fftshift(abs(fft2(f)));
imshow(F,[-1 5],'notruesize');
title('X轴方向移动后的傅立叶变换频谱');

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(50:349,475:524) = 1;
imshow(f,'notruesize');
title('Y轴方向移动后的图像');
subplot(224)
F=fftshift(abs(fft2(f)));
imshow(F,[-1 5],'notruesize');
title('Y轴方向移动后的傅立叶变换频谱');