傅立叶变换频谱图怎么看？

I = imread('rice.png'); figure; subplot(1, 2, 1); imshow(I); title('原图');f = fft2(I); f = fftshift(log(1+abs(f))); subplot(1, 2, 2); imshow(f, []); title('频谱图');

f= imread('rice.png'); f=im2double(f); F=fft2(f,256,256); F2=fftshift(F); F3=log(1+abs(F2)); figure(2) subplot(1,2,1); imshow(f); title('源图像');subplot(1,2,2); imshow(F3); title('傅里叶频谱');

上面代码得到的频谱图该怎么看？

图像的log频谱又是指什么？

添加评论

知乎用户，计算成像

Vincent Ku、Scofield Michael、小二 等人赞同

谢邀。

第一张图展示了频谱，注意：
0) fft2() 函数得到的傅里叶变换是低频在四周的。为了获得更好的直观的感受，人们经常用fftshift() 调整结果以使得低频在中央。
1) 频谱是复数，为了用二维图像展示频谱，就取了绝对值；
2) 由于频谱的动态范围太大，而显示器的动态范围有限，所以用log变换降低频谱的动态范围。
3) MATLAB的imshow函数，第二个输入变量是"[]"的时候，imshow会自动拉伸动态范围使得图像对比度提升。

看频谱一般看：
1) 幅度谱。即第一张图的右图，观察图像的频率分布，哪里亮哪里暗。
2) 相位谱。用得不多，一般在需要对比的时候会用到。

第二张图显示的也是频谱，与第一张的唯一区别就是imshow中第二个输入变量不是"[]"，而是缺省。这时候MATLAB会不自动进行动态范围拉伸，而以[0,1]为显示区间。而大部分频谱的绝对值都是大于1的，所以得到的结果大部分都是白色的。

两张图从信号处理的角度看没有任何区别，只是imshow用得不一样，就产生了看上去不一样的结果。