求助 如何正确查看图像频谱

花如月

[X,Y]=meshgrid(1:79,1:75);
M=180*cos(sqrt(2)*pi*10.03*(X+Y));
M是一个规则的线条图，很多资料都说“规则线条的幅值谱，只有几个亮点”。问题是如何得出这些亮点图呢？


附图是一个文献里的结果，貌似也不是多复杂的问题。不过一直没搞定，我绘制的图不对，大家帮忙看看

zhlong

呵呵，二维的频谱没搞过。

zhai1

FFT2对图像进行二维快速傅立叶变换，得到F(U,V);其绝对值就是幅度图，由于频谱中值得动态范围较大，一般要经过对数变换才能显示图像的频谱图（幅度图）

yangzj

应该是显示的问题。

[X,Y]=meshgrid(1:79,1:75);
M=180*cos(sqrt(2)*pi*10.03*(X+Y));
Mf=fft2(M);
Mf=abs(Mf);
figure();
subplot(121);imshow(M,[]);
subplot(122);imshow(Mf,[]);

nim001

很久没有弄过图像的二维频谱了,但还记得一点儿.
图像经过FFT2运算之后,其中间的部分为图像高频部分,就是图像中比较亮，或者棱角部分；外侧部分为低频部分，就是图像的细节，
要正确查看图像频谱，一般把数字图像用FFT2之后，都会做一个频谱移动，把低频部分移到中间，而高频部分移到外侧，采用FFTSHIFT命令就可以实现。然后就可以观察其频谱图了，或者对图像进行滤波等操作了

花如月

好久没来这里了，多谢4楼和5楼的帮忙

jhqing

回复 nim001 的帖子

请问，得到的二维频谱图跟原图像有什么对应关系，x，y坐标代表什么，