HIO-algorithm 使用HIO算法实图像现相位恢复

使用HIO算法实图像现相位恢复，包含源代码，简单实用，适于初学者学习。

x=imread（‘lena.bmp‘）;
x=rgb2gray（x）;
x=double（x）;
mask=zeros（size（x））;
mask（200:300200:300）=1;
x1=x.*mask;
xm=uint8（x1）;
xf=fft2（x1）;
absfx=abs（xf）;
y=randn（size（x））;
iteration=400;
beta=0.97;
i=1:iteration;
for iter=1:iteration
    y=myHIO（y absfx mask beta）;
%          y=myER（y absfx mask）;
          x1=uint8（x1）;
          y=uint8（real（y））;
          figure（1）
          subplot（121） imshow（x1） axis image title（‘ER‘）;
          subplot（122） imshow（y） axis image title（iter）;
          y=double（y）;           
end
PSNR = Psnr（x1y）;
y=uint8（real（y））;
imwrite（y‘poisson HIO 400.bmp‘‘bmp‘）;
figure（2）
subplot（121）imshow（xm[]）title（‘原始图像‘）;
subplot（122）imshow（y[]）title（‘相位恢复后图象‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-03-22 22:53  HIO algorithm\
     文件      786486  2012-04-09 19:55  HIO algorithm\lena.bmp
     文件         750  2015-03-22 22:56  HIO algorithm\main.m
     文件         106  2014-04-21 17:34  HIO algorithm\myHIO.m
     文件          88  2014-04-21 17:21  HIO algorithm\Pm1.m
     文件          50  2014-06-10 16:16  HIO algorithm\Ps1.m
     文件         398  2014-05-05 09:50  HIO algorithm\Psnr.m