傅里叶变换.txt

来自「压缩包中共有7个文件」· 文本 代码 · 共 76 行

%二值图象的傅里叶变换
f=zeros(30,30);
f(5:24,13:17)=1;
figure,imshow(f,'notruesize')
F=fft2(f);
figure,mesh(fftshift(abs(F)));%绘制频谱图(量化表现)(此处有问题)fftshift把零频率部分搬移到频谱中间
F2=fftshift(log(abs(F)));%未填充
figure,imshow(F2,[-1 5],'notruesize');%显示频谱图像（颜色表现）
F=fft2(f,256,256);%填充
figure,imshow(fftshift(log(abs(F))),[-1 5]);
colormap(jet);colorbar


%平移特性
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('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;
f=imrotate(f,45,'bilinear','crop');%以图象中心为原点旋转45度,(好像有问题)
imshow(f,'notruesize');
title('图象正向旋转45^{0}');
subplot(224);
F=fftshift(abs(fft2(f)));
imshow(F,[-1 5],'notruesize');
title('图象旋转45^{0}后的傅里叶变换频谱')