Fourier Transforms and Phase Masks

Question

Simon 2014-2-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/119112-fourier-transforms-and-phase-masks

回答： Jesus Eduardo 2014-9-9

I am trying to perform a fourier transform of a phase mask to reconstruct my original image after having performed the fourier transform and extracting its phase. I am able to get a reconstructed image but it seems that the original image is overlayed with a x and y mirror image of the original image.

Does anyone know how I can get rid of this or what I am doing wrong?

My code is as follows:

clear all; close all;
I = imread('batman.gif', 'gif');
I = double(I);
I = I./max(max(I));
avgl=mean(mean(I));
figure; imshow(mat2gray(I));
title('Original Object');
figure;
axis([0 , 101, 0, 1]);
xlabel ('Number of Iterations')
ylabel ('RMSE')
hold on
I1 = I;
for n=1:101; %Iterations to optimize the phase hologram
      H = fftshift((fft2(fftshift(I1))));
      I2 = fftshift(fft2(fftshift(exp(1j.*angle(H)))));
      avg2=mean(mean(abs(I2)));
      I2=(I2./avg2).*avgl;
      rmse=mean(mean((abs(I2)-I).^2))^0.5;
      plot(n,rmse,'o');
      pause(0.00001); %To see the error in each iteration.
      I1=I.*exp(1j*angle(I2));
  end
figure; imshow(H);
colormap gray
ti = get(gca,'TightInset')
set(gca,'Position',[ti(1) ti(2) 1-ti(3)-ti(1) 1-ti(4)-ti(2)]);
set(gca,'units','centimeters')
pos = get(gca,'Position');
ti = get(gca,'TightInset');
set(gcf, 'PaperUnits','centimeters');
set(gcf, 'PaperSize', [pos(3)+ti(1)+ti(3) pos(4)+ti(2)+ti(4)]);
set(gcf, 'PaperPositionMode', 'manual');
set(gcf, 'PaperPosition',[0 0 pos(3)+ti(1)+ti(3) pos(4)+ti(2)+ti(4)]);
imwrite(H, 'SupermanPhaseMask2.gif')
hold off
I2=I2./max(max(abs(I2)));
figure; imshow(abs(I2));
title('Reconstucted Image')