Spline based iterative phase retrieval algorithm for X-ray differential phase contrast radiography

Iterative phase contrast CT reconstruction with novel tomographic operator and data-driven prior

Breast cancer remains the most prevalent malignancy in women in many countries around the world, thus calling for better imaging technologies to improve screening and diagnosis. Grating interferometry (GI)-based phase contrast X-ray CT is a promising technique which could make the transition to clinical practice and improve breast cancer diagnosis by combining the high three-dimensional resolution of conventional CT with higher soft-tissue contrast. Unfortunately though, obtaining high-quality images is challenging. Grating fabrication defects and photon starvation lead to high noise amplitudes in the measured data. Moreover, the highly ill-conditioned differential nature of the GI-CT forward operator renders the inversion from corrupted data even more cumbersome. In this paper, we propose a novel regularized iterative reconstruction algorithm with an improved tomographic operator and a powerful data-driven regularizer to tackle this challenging inverse problem. Our algorithm combines the L-BFGS optimization scheme with a data-driven prior parameterized by a deep neural network. Importantly, we propose a novel regularization strategy to ensure that the trained network is non-expansive, which is critical for the convergence and stability analysis we provide. We empirically show that the proposed method achieves high quality images, both on simulated data as well as on real measurements.

Phase retrieval based on deep learning in grating interferometer

Grating interferometry is a promising technique to obtain differential phase contrast images with illumination source of low intrinsic transverse coherence. However, retrieving the phase contrast image from the differential phase contrast image is difficult due to the accumulated noise and artifacts from the differential phase contrast image (DPCI) reconstruction. In this paper, we implemented a deep learning-based phase retrieval method to suppress these artifacts. Conventional deep learning based denoising requires noise/clean image pair, but it is not feasible to obtain sufficient number of clean images for grating interferometry. In this paper, we apply a recently developed neural network called Noise2Noise (N2N) that uses noise/noise image pairs for training. We obtained many DPCIs through combination of phase stepping images, and these were used as input/target pairs for N2N training. The application of the N2N network to simulated and measured DPCI showed that the phase contrast images were retrieved with strongly suppressed phase retrieval artifacts. These results can be used in grating interferometer applications which uses phase stepping method.

Reweighted L1-norm regularized phase retrieval for x-ray differential phase contrast radiograph

Talbot-Lau x-ray grating interferometry greatly decreases the requirements on x-ray sources to realize differential phase contrast imaging and has found many applications in industrial and medical imaging. Phase retrieval from the noisy differential signal is crucial for quantitative analysis, comparison, and fusion with other imaging modalities. In this paper, we introduce a reweighted L1-norm based nonlinear regularization method for the phase retrieval problem. Both simulation and experimental results demonstrated that, comparing with the widely used L1-norm based regularization method and Wiener filter method, the proposed method is more effective both in eliminating the strip noises and in preserving the image detail.

Improved visualization of X-ray phase contrast volumetric data through artifact-free integrated differential images

Artifacts arising when differential phase images are integrated is a common problem to several X-ray phase-based experimental techniques. The combination of noise and insufficient sampling of the high-frequency differential phase signal leads to the formation of streak artifacts in the projections, translating into poor image quality in the tomography slices. In this work, we apply a non-iterative integration algorithm proven to reduce streak artifacts in planar (2D) images to a differential phase tomography scan. We report on how the reduction of streak artifacts in the projections improves the quality of the tomography slices, especially in the directions different from the reconstruction plane. Importantly, the method is compatible with large tomography datasets in terms of computation time.

Fast, non-iterative algorithm for quantitative integration of X-ray differential phase-contrast images

X-ray phase contrast imaging is gaining importance as an imaging tool. However, it is common for X-ray phase detection techniques to be sensitive to the derivatives of the phase. Therefore, the integration of differential phase images is a fundamental step both to access quantitative pixel content and for further analysis such as segmentation. The integration of noisy data leads to artefacts with a severe impact on image quality and on its quantitative content. In this work, an integration method based on the Wiener filter is presented and tested using simulated and real data obtained with the edge illumination differential X-ray phase imaging method. The method is shown to provide high image quality while preserving the quantitative pixel content of the integrated image. In addition, it requires a short computational time making it suitable for large datasets.

Forward model for propagation-based x-ray phase contrast imaging in parallel- and cone-beam geometry

We demonstrate a fast, flexible, and accurate paraxial wave propagation model to serve as a forward model for propagation-based X-ray phase contrast imaging (XPCI) in parallel-beam or cone-beam geometry. This model incorporates geometric cone-beam effects into the multi-slice beam propagation method. It enables rapid prototyping and is well suited to serve as a forward model for propagation-based X-ray phase contrast tomographic reconstructions. Furthermore, it is capable of modeling arbitrary objects, including those that are strongly or multi-scattering. Simulation studies were conducted to compare our model to other forward models in the X-ray regime, such as the Mie and full-wave Rytov solutions.

Subspace-based resolution-enhancing image reconstruction method for few-view differential phase-contrast tomography

It is well known that properly designed image reconstruction methods can facilitate reductions in imaging doses and data-acquisition times in tomographic imaging. The ability to do so is particularly important for emerging modalities, such as differential x-ray phase-contrast tomography (D-XPCT), which are currently limited by these factors. An important application of D-XPCT is high-resolution imaging of biomedical samples. However, reconstructing high-resolution images from few-view tomographic measurements remains a challenging task due to the high-frequency information loss caused by data incompleteness. In this work, a subspace-based reconstruction strategy is proposed and investigated for use in few-view D-XPCT image reconstruction. By adopting a two-step approach, the proposed method can simultaneously recover high-frequency details within a certain region of interest while suppressing noise and/or artifacts globally. The proposed method is investigated by the use of few-view experimental data acquired by an edge-illumination D-XPCT scanner.

Joint absorption and phase retrieval in grating-based x-ray radiography

Given the raw absorption and differential phase-contrast images obtained from a grating-based x-ray radiography, we formulate the joint denoising of the absorption image and retrieval of the non-differential phase image as a regularized inverse problem. The choice of the regularizer is driven by the existing correlation between absorption and differential phase; it leads to the linear combination of a total-variation norm with a total-variation nuclear norm. We then develop the corresponding algorithm to efficiently solve this inverse problem. We evaluate our method using different experiments, including mammography data. We conclude that our method provides useful information in the context of mammography screening and diagnosis.