图像的傅里叶变换包含了幅度谱和相位谱,分别对幅度谱和相位谱做二维傅里叶逆变换可得到,相位谱在重构过程中的重要性。
现在需要用幅度谱和四分之一的相位谱来恢复原图。
参考文献见“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),[])
上面的程序只能做到恢复原图的四分之一,而不能达到文献中的情况
有没有哪位能帮忙看看问题出现在哪里了呢