标签:
1. 先载入一幅灰度图像,如下:
(非灰度图)
2. 利用函数fft2,对其进行快速傅立叶变换, 并利用函数fftshift 将变换后的图像原点移动到频率矩形的中心。
3. 利用abs()函数来得到傅立叶频谱;angle()函数得到相位图;
4. 利用imshow 来可视化图像,观察图像的特点;
一. 结果(图像):
二. 分析说明:
1. 在载入图片的时候要注意图片要存放到该文件的文件夹中。
2. 下载的图片看似是灰度图,实际上并非是灰度图,所以加多了一句代码“A=rgb2gray(A)”
三. 附件程序
A=imread(‘1.jpg‘); %载入图片
A=rgb2gray(A)
B=fftshift(fft2(A)); % 进行傅立叶变换
subplot(231)
imshow(A);
title(‘原始图像‘);
subplot(232)
imshow(abs(B),[ ]);
title(‘原始频谱图‘);
subplot(233)
imshow(log(abs(B)),[ ]);
title(‘取对数后的频谱图‘);
subplot(234)
imshow(angle(B),[ ]);
title(‘相位图‘);
subplot(235)
imshow(real(B),[ ]);
title(‘实部图‘);
subplot(236)
imshow(imag(B),[ ]);
title(‘虚部图‘);
colormap(jet(64)) %给图片上色
转自:http://blog.163.com/luoyun_dreamer/blog/static/215529070201322612858437/
标签:
原文地址:http://www.cnblogs.com/sweeeper/p/4687550.html