clerk=imread('clerk.tif');
lena=imread('lena.tif');

clerk = imresize(clerk,0.5,'nearest')

%
% fo_clerk=fftshift(fft2(clerk));
% fo_lena=fftshift(fft2(lena));
% figure('Name','clerk')
% surf_clerk = surf(abs(log(fo_clerk)));
% figure('Name','lena')
% surf_lena = surf(abs(log(fo_lena)));
%
% figure('Name','phase clerk')
% imagesc( angle(fo_clerk) )
% figure('Name','phase lena')
% imagesc( angle(fo_lena) )

%
%
% fun = @(x) abs(log(fftshift(fft2(x))));
% blk = blkproc(clerk,[5,5],fun)
% figure('Name','FClerk')
% imagesc(blk)
%
% blk = blkproc(lena,[5,5],fun)
% figure('Name','FLena')
% imagesc(blk)


lenafft=fft2(double(lena));
clerkfft=fft2(double(clerk));

amplitude_lena = abs(lenafft);
amplitude_clerk = abs(clerkfft);

%pase swap
freq_lena = amplitude_lena .*exp(i*angle(clerkfft));
freq_clerk = amplitude_clerk .*exp(i*angle(lenafft));


inverselena= ifft2(freq_lena);
inverseclerk= ifft2(freq_clerk);



figure('Name','lena')
imshow(inverselena)
figure('Name','clerk')
imshow(inverseclerk)