交换两幅图像的幅度谱和相位谱，并重构图像

转载自：http://blog.sina.com.cn/s/blog_909778ea0100y8ju.html
博主：科技狸的博客；
来源：新浪博客；

---

交换两幅图像的幅度谱和相位谱，并重构图像！！！
图像x

图像y

%注意两幅图像的大小要相同，此处均为361*461
%读入图像
x=imread('x.bmp');
y=imread('y.bmp');
% 傅里叶变换
xf=fft2(double(x));
yf=fft2(double(y));
%取幅度和相位
xf1=abs(xf);
xf2=angle(xf);

yf1=abs(yf);
yf2=angle(yf);
%交换相位
xfr=xf1.*cos(yf2)+xf1.*sin(yf2).*i;
yfr=yf1.*cos(xf2)+yf1.*sin(xf2).*i;
%傅里叶反反变换
xr=abs(ifft2(xfr));
yr=abs(ifft2(yfr));
%转换成uint8类型
xr=uint8(xr);
yr=uint8(yr);
%显示
subplot(2,2,1);imshow(x);title('x 原图');
subplot(2,2,2);imshow(y);title('y 原图');
subplot(2,2,3);imshow(xr,[]);title('x 幅度谱 与 y相位谱');
subplot(2,2,4);imshow(yr,[]);title('y 幅度谱 与 x相位谱');

实验结果图：

交换两幅图像的幅度谱和相位谱，并重构图像

实验说明：在双谱重构中，相位谱起了更大作用。