实验截图:
- 傅里叶变换的谱和相角:
实验代码:
img = imread('erciyuan.jpg');
img = rgb2gray(img);
img1 = fft2(img);
img1 = fftshift(img1);%快速傅里叶变换
img2 = log(abs(img1)+1);
img3 = angle(img1)*180/pi;%根据快速傅里叶变换的结果得到相角
subplot(1,3,1);imshow(img,[]);title('原图');
subplot(1,3,2);imshow(img2,[]);title('原图进行傅里叶变换的谱');
subplot(1,3,3);imshow(img3,[]);title('原图进行傅里叶变换的相角');
结果分析:
- 我们根据上一个实验,将快速傅里叶变换后的零频分量移动到频谱中心,然后将其取绝对值后,就能得到快速傅里叶变换的相位谱,接着再对其取角度值就能得到快速傅里叶变换的相角谱。