A Fourier-domain algorithm for total-variation regularized phase retrieval in differential X-ray phase contrast imaging

Phase retrieval in differential X-ray phase contrast imaging involves a one dimensional integration step. In the presence of noise, standard integration methods result in image blurring and streak artifacts. This work proposes a regularized integration method which takes the availability of two dimensional data as well as the integration-specific frequency-dependent noise amplification into account. In more detail, a Fourier-domain algorithm is developed comprising a frequency-dependent minimization of the total variation orthogonal to the direction of integration. For both simulated and experimental data, the novel method yielded strong artefact reduction without increased blurring superior to the results obtained by standard integration methods or regularization techniques in the image domain.