我的问题可能很简单，但我确实需要确保我正在进行的流程是有效的。
MRI扫描仪获取k空间数据S，然后采用3D逆傅立叶变换重建空间域中的实际图像，其形式为M.exp(j.P)，其中{ {1}}是幅度图像，M是相位图像。
现在想象我所拥有的是幅度和相位图像，我想从它们生成P。
现在我检查了相位图像，其值介于0到4096之间;所以，首先我必须将其标准化为S范围。然后，我可以将[-pi,pi]计算为S的3D傅里叶变换。
在这个过程中还有什么我应该考虑的吗？
现在，如何从复杂/原始图像I = M.*exp(j.*P)计算相位图像P？我尝试了MATLAB的内置函数I;所以angle，但我注意到P1 = angle(I)和P1不相等，即使P不等于exp(j*P1)。但是，当我忽略幅度图像exp(j*P)并从中计算I1 = exp(j*P)（P2）时，P2 = angle(I1)等于exp(j*P2) !!我认为原因是exp(j*P)的零值，但是当我有M时，我不知道重新生成P的正确方法是什么。