为什么从MAGNITUDE重建的测试图像之间存在这种强度变化（或交换）的神奇矛盾。
clc;
clear all;
close all;
i1=imread('C:\Users\Admin\Desktop\rough\Capture1.png');
i1=rgb2gray(i1);
f1=fftn(i1);
mag1=abs(f1);
s=log(1+fftshift(f1));
phase1=angle(f1);
r1=ifftshift(ifftn(mag1));
r2=ifftn(exp(1i*phase1));
figure,imshow(i1);
figure,imshow(s,[]);
figure,imshow(uint8(r1));
figure,imshow(r2,[]);
r2=histeq(r2);
r3=histeq(uint8(r2));
figure,imshow(r2);
figure,imshow(r3);
#1 楼
图1.（c）仅显示了从MAGNITUDE光谱重建的测试图像。我们可以说低频像素的强度值比高频像素高。
实际上，这是不正确的。相位值确定图像的正弦分量的偏移。在零相位的情况下，所有正弦波都位于同一位置，您会得到一个对称图像，其结构与原始图像完全没有任何实际相关性。在同一位置居中意味着正弦曲线在该位置最大，这就是为什么在图1.c中部有一个大的白色斑块的原因。
由于相位一致的原理，仅相位重构保留了特征。在边缘和线条的位置，大多数正弦波分量具有相同的相位。请参阅http://homepages.inf.ed.ac.uk/rbf/CVonline/LOCAL_COPIES/OWENS/LECT7/node2.html。http：//www.csse.uwa可单独使用此方法正确地检测线条和边缘。 edu.au/~pk/research/pkpapers/phasecorners.pdf，不考虑幅度。因此，您可以看到相位信息是最重要的。
改变各个分量正弦波的幅度会改变特征的形状。当您执行仅阶段重建时，将所有大小都设置为1，这将更改要素的形状，但不会更改其位置。在许多图像中，低频分量的幅度要比高频分量高，因此仅相位重建的确看起来像是高通滤波器。功能。
您无法添加仅相位和仅幅度的图像来获取原始图像。您可以在傅立叶域中将它们相乘并变换回原始图像。
#2 楼
在您的行mag1=abs(f1); 中，图像的总强度保持不变（通过累加所有像素的强度进行测试）。在傅立叶空间中拒绝相位信息只会导致实际空间中的强度在空间上重新分布，从而使r1的总惯性与i1相同。 />您正在将每个像素的振幅（在傅立叶空间中）标准化为1，因此图像的总强度将发生变化。由于相位携带了图像的大部分空间信息，因此仍然保留了图像的主要特征。
评论
$ \ begingroup $
@geometrical谢谢您的解释。我读过文章，但我对此表示怀疑。先生，您说：“在正弦线的位置，大多数正弦波分量具有相同的相位。”并使用相位一致性方法可以检测到这些。但是，白色大片的低频分量也可以具有相同的相位？因此也应检测这些频率。我也准备了一个代码，就像您在最后一条答案中所说的那样，但是我无法重建原始图像...我要在下一个注释中添加代码。
$ \ endgroup $
–sagar
2014年6月24日6:10
$ \ begingroup $
@geometrical'clc;清除所有关闭所有; i1 = imread（'C：\ Users \ Admin \ Desktop \ rough \ Capture1.png'）; i1 = rgb2gray（i1）;图，imshow（i1）; f1 = fftn（i1）; mag1 = abs（f1）; phase1 = angle（f1）; a1 = fftn（mag1）; a2 = fftn（phase1）; a3 = a1。* a2; a4 = ifftn（a3）;图，imshow（uint8（a4））;'
$ \ endgroup $
–sagar
14年6月24日在6:15
$ \ begingroup $
在大的白色斑块图像中，所有正弦曲线都已偏移为中心具有相同相位（= 0）。相位一致性是关于检测图像中的线条或边缘特征。这是相位对图像结构最重要的另一证明。对于您的代码，我的意思是重建相位和幅度图像。
$ \ endgroup $
–geometrikal
2014年6月24日7:37
$ \ begingroup $
clc;清除所有关闭所有; i1 = imread（'peppers.tif'）; i1 = rgb2gray（i1）;图，imshow（i1）; f1 = fftn（i1）; mag1 = abs（f1）; phase1 = exp（1i * angle（f1））; a1 = ifftn（mag1）; a2 = ifftn（phase1）; a3 = fftn（a1）。* fftn（a2）; a4 = ifftn（a3）;图，imshow（uint8（a4））;
$ \ endgroup $
–geometrikal
14年6月24日在7:38
$ \ begingroup $
抱歉打扰您，但是具有相同相位的低频分量会发生什么。它们也应仅在相位重建中保存。
$ \ endgroup $
–sagar
2014年6月24日13:07