数字图像处理实验三

% 将两幅256×256的灰度图像进行二维傅立叶变换，分别得到各自的幅度谱图像和相位谱图像，显示这4幅图像。
% 将两幅图像的相位谱交换，进行反傅立叶变换，显示这两幅重建图像。

% 二维快速傅立叶变换函数fft2( ) 
% 相位提取函数 angle( )
% 幅度提取函数abs ( )

%
clear;clc;
%
im1 = double(imread('1-Debbie4.gif'));
im2 = double(imread('12-Panda.bmp'));
%
im1_f = fft2(im1);
im2_f = fft2(im2);
%
im1_abs = abs(im1_f);
im2_abs = abs(im2_f);
im1_ang = angle(im1_f);
im2_ang = angle(im2_f);
im1_d = fftshift(im1_f);
im2_d = fftshift(im2_f);
%
im1_e = im1_abs.*exp(1i*im2_ang);
im2_e = im2_abs.*exp(1i*im1_ang);
%
im1_r = real(ifft2(im1_e)); 
im2_r = real(ifft2(im2_e)); 


% display
figure
subplot(231)
imshow(mat2gray(im1));
xlabel('Original 1')
subplot(234)
imshow(mat2gray(im2));
xlabel('Original 2')
subplot(232)
imshow(mat2gray(10*log10(abs(im1_d))));
xlabel('Magititute 1')
subplot(235)
imshow(mat2gray(10*log10(abs(im2_d))));
xlabel('Magititute 2')
subplot(233)
imshow(mat2gray(angle(im1_d)));
xlabel('Angle 1')
subplot(236)
imshow(mat2gray(angle(im2_d)));
xlabel('Angle 2')

%
figure
subplot(121)
imshow(mat2gray(im1_r));
xlabel('Rec. Image 1')
subplot(122)
imshow(mat2gray(im2_r));
xlabel('Rec. Image 2')