[论文解读] Regularized Newton methods for simultaneous Radon inversion and phase retrieval in phase contrast tomography

Simon Maretzke|arXiv (Cornell University)|Feb 17, 2015

Advanced X-ray Imaging Techniques被引用 5

一句话总结

本文提出了一种正则化牛顿方法，用于在相位对比断层扫描中同时实现断层重建与相位恢复，利用近场和远场传播模型。理论证明仅使用单一探测器距离即可唯一恢复折射和吸收特性，显著提升了弱吸收样品（如生物细胞）在三维纳米尺度成像中的稳定性与准确性。

ABSTRACT

Promoted by the advent of coherent synchrotron light sources, phase contrast tomography allows to resolve three-dimensional variations of an unknown sample's complex refractive index from scattering intensities recorded at different incident angles of an X-ray beam. By diffractive free-space propagation of the transmitted wave field, this method is sensitive not only to absorption but also to refractive phase shifts induced by the specimen, permitting three-dimensional nanoscale imaging of quasi-transparent samples such as biological cells. However, the reconstruction of the specimen structure from the observed data constitutes an algorithmically challenging nonlinear ill-posed inverse problem, mainly due to the characteristic loss of phase information in the detection of the wave field. In this work, regularized Newton methods are developed for the solution of this tomographic phase retrieval problem, based on a detailed analysis of its mathematical structure. We consider both the near-field- or Fresnel regime characterized by a moderate propagation length between sample and detector and the far-field limit of large detector distances, where propagation is governed by the Fourier transform. In the former setting, excellent numerical reconstructions are obtained via the chosen Newton-type approach, supplemented by novel theoretical results stating that measurements from a single detector distance are sufficient to uniquely recover both refraction and absorption of a sample. The proposed algorithm simultaneously performs tomographic- and phase reconstruction, which is found to stabilize the latter by exploiting correlations between the diffraction patterns recorded under different incident angles.

研究动机与目标

• 解决由于强度测量中缺失相位信息而引起的相位对比断层扫描中的非线性、不适定逆问题。
• 开发一种稳定且高效的重建方法，同时恢复准透明样品的复折射率与断层结构。
• 证明单一探测器距离足以唯一恢复吸收与折射分量，挑战了以往认为需要多距离的假设。
• 利用不同入射角度之间的角相关性以稳定相位恢复并提升重建质量。
• 为Fresnel与远场传播模型背景下正则化牛顿方法提供数学上严谨的框架。

提出的方法

• 基于傍轴波方程和近场条件下的Fresnel传播，建立前向模型，其中强度在自由空间传播后由探测器测量。
• 推导前向算子的Fréchet导数，以支持牛顿型迭代反演，确保局部二次收敛。
• 应用Tikhonov型正则化以稳定不适定逆问题的求解，惩罚项促进折射率的平滑性。
• 利用圆柱坐标系下的傅里叶变换与Radon变换，分析近场与远场区域中前向算子的数学结构。
• 实现牛顿方法的离散化版本，结合伴随算子与高效的数值求解器，用于大规模三维重建。
• 将多角度入射数据整合到联合反演框架中，以利用角度相关性并提升相位恢复的稳定性。

实验结果

研究问题

• RQ1 在相位对比断层扫描中，单一探测器距离是否足以唯一确定样品复折射率的折射与吸收分量？
• RQ2 与顺序方法相比，同时进行断层与相位重建在相位恢复的稳定性与准确性方面有何提升？
• RQ3 近场与远场区域中前向算子的数学结构是什么？其如何影响牛顿型方法的收敛性？
• RQ4 不同入射角度之间的角相关性在多大程度上增强了相位分量的重建？
• RQ5 如何最优地将正则化与牛顿迭代结合，以稳定非线性逆问题的求解？

主要发现

• 通过前向算子的单射性理论分析，证明单一探测器距离足以唯一恢复样品复折射率的折射与吸收分量。
• 所提出的正则化牛顿方法在近场区域实现了出色的数值重建，展现出在弱吸收样品三维成像中高保真度与高稳定性。
• 通过利用不同入射角度下衍射图案的相关性，同时进行断层与相位重建显著稳定了相位恢复。
• 该方法在实际噪声水平下实现收敛，正则化有效抑制了伪影并提升了重建图像的分辨率。
• 理论分析证实，在温和条件下，近场区域中前向算子是单射的，从而证明了解的唯一性。
• 远场极限由傅里叶变换主导，从而可建立一种简化但精确的模型，适用于大探测器距离。