Stability estimates for linearized near-field phase retrieval in X-ray phase contrast imaging
Provides theoretical stability guarantees for a widely used phase retrieval model in X-ray imaging, addressing a fundamental mathematical challenge for practitioners.
The paper shows that the linearized near-field phase retrieval problem in X-ray phase contrast imaging is well-posed under compact support, with Lipschitz stability constants that grow exponentially with Fresnel number generally, but algebraically for homogeneous objects or two-distance measurements.
Propagation-based X-ray phase contrast enables nanoscale imaging of biological tissue by probing not only the attenuation, but also the real part of the refractive index of the sample. Since only intensities of diffracted waves can be measured, the main mathematical challenge consists in a phase-retrieval problem in the near-field regime. We treat an often used linearized version of this problem known as contract transfer function model. Surprisingly, this inverse problem turns out to be well-posed assuming only a compact support of the imaged object. Moreover, we establish bounds on the Lipschitz stability constant. In general this constant grows exponentially with the Fresnel number of the imaging setup. However, both for homogeneous objects, characterized by a fixed ratio of the induced refractive phase shifts and attenuation, and in the case of measurements at two distances, a much more favorable algebraic dependence on the Fresnel number can be shown. In some cases we establish order optimality of our estimates.