图像的四分之一相位重构

miki321

图像的傅里叶变换包含了幅度谱和相位谱，分别对幅度谱和相位谱做二维傅里叶逆变换可得到，相位谱在重构过程中的重要性。
现在需要用幅度谱和四分之一的相位谱来恢复原图。
参考文献见“Image  Representation and  Reconstruction from Spectral  Amplitude or Phase”。
具体流程如下

文献中实现如下：a为原图，b为采用上述流程恢复的图像

我编写的程序如下：
clc
clear
close all
L=imread('lena1.bmp');
[N,M]=size(L);
x=double(L);
X=fft2(x);
yk=zeros(N,M);
for k=0:50
    yk(1:N/2,1:M/2)=x(1:N/2,1:M/2);
    Yk=fft2(yk);
    Yk1=abs(X).*exp(j*angle(Yk));
    yk1=ifft2(Yk1);
    yk=yk1;
end
   imshow(abs(yk),[])
上面的程序只能做到恢复原图的四分之一，而不能达到文献中的情况
有没有哪位能帮忙看看问题出现在哪里了呢