交换两幅图像的幅度谱和相位谱，并重构图像 此博文包含图片

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

图像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相位谱');

 

 

 

实验结果图：

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