在做数字图像处理实验的时候,需要将灰度图像做fft变换,显示其频谱图。代码如下:
clc, clear, close all
figure(1);
grayimg =imread('grayimg.jpg');
imshow(grayimg);
figure(2);
B=fftshift(fft2(grayimg));
imshow(log(abs(B)),[]);
colormap(jet(64));
colorbar;
按道理应该能得到一个彩色的频谱图,然而我得到的却是:
我的版本是R2016b,相同的代码在同学的R2018b上据说没遇到这个问题。
解决方法
Google到的一个相关讨论帖子
方法1
R2016b上可以,R2017a据说不可
就是把colormap(jet(64))
改成了colormap(gca, jet(64))
clc, clear, close all
figure(1);
grayimg =imread('grayimg.jpg');
imshow(grayimg);
figure(2);
B=fftshift(fft2(grayimg));
imshow(log(abs(B)),[]);
colormap(gca, jet(64));
colorbar;
方法2
R2016b和R2017a都可
clc, clear, close all
figure(1);
grayimg =imread('grayimg.jpg');
imshow(grayimg);
figure(2);
B=fftshift(fft2(grayimg));
AxesH = axes;
imshow(log(abs(B)),[], 'Parent', AxesH);
colormap(AxesH, jet(64));
colorbar;
最后结果: