Stability estimates for linearized near-field phase retrieval in X-ray phase contrast imaging

TL;DR

This paper analyzes the stability of linearized phase retrieval in X-ray phase contrast imaging, showing it is well-posed under certain conditions and providing bounds on stability constants that depend on the imaging setup.

Contribution

It establishes the well-posedness of the linearized phase retrieval problem and derives bounds on the Lipschitz stability constant, including cases with improved algebraic dependence.

Findings

The inverse problem is well-posed with compact support assumptions.

Stability constants grow exponentially with the Fresnel number generally.

For homogeneous objects and measurements at two distances, stability bounds are more favorable.

Abstract

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.