Model-Based Iterative Reconstruction for Propagation-Based Phase-Contrast X-Ray CT including Models for the Source and the Detector

High resolution propagation-based lung imaging at clinically relevant X-ray dose levels

Absorption-based clinical computed tomography (CT) is the current imaging method of choice in the diagnosis of lung diseases. Many pulmonary diseases are affecting microscopic structures of the lung, such as terminal bronchi, alveolar spaces, sublobular blood vessels or the pulmonary interstitial tissue. As spatial resolution in CT is limited by the clinically acceptable applied X-ray dose, a comprehensive diagnosis of conditions such as interstitial lung disease, idiopathic pulmonary fibrosis or the characterization of small pulmonary nodules is limited and may require additional validation by invasive lung biopsies. Propagation-based imaging (PBI) is a phase sensitive X-ray imaging technique capable of reaching high spatial resolutions at relatively low applied radiation dose levels. In this publication, we present technical refinements of PBI for the characterization of different artificial lung pathologies, mimicking clinically relevant patterns in ventilated fresh porcine lungs in a human-scale chest phantom. The combination of a very large propagation distance of 10.7 m and a photon counting detector with [Formula: see text] pixel size enabled high resolution PBI CT with significantly improved dose efficiency, measured by thermoluminescence detectors. Image quality was directly compared with state-of-the-art clinical CT. PBI with increased propagation distance was found to provide improved image quality at the same or even lower X-ray dose levels than clinical CT. By combining PBI with iodine k-edge subtraction imaging we further demonstrate that, the high quality of the calculated iodine concentration maps might be a potential tool for the analysis of lung perfusion in great detail. Our results indicate PBI to be of great value for accurate diagnosis of lung disease in patients as it allows to depict pathological lesions non-invasively at high resolution in 3D. This will especially benefit patients at high risk of complications from invasive lung biopsies such as in the setting of suspected idiopathic pulmonary fibrosis (IPF).

Suppressing multi-material and streak artifacts with an accelerated 3D iterative image reconstruction algorithm for in-line X-ray phase-contrast computed tomography

In-line X-ray phase-contrast computed tomography typically contains two independent procedures: phase retrieval and computed tomography reconstruction, in which multi-material and streak artifacts are two important problems. To address these problems simultaneously, an accelerated 3D iterative image reconstruction algorithm is proposed. It merges the above-mentioned two procedures into one step, and establishes the data fidelity term in raw projection domain while introducing 3D total variation regularization term in image domain. Specifically, a transport-of-intensity equation (TIE)-based phase retrieval method is updated alternately for different areas of the multi-material sample. Simulation and experimental results validate the effectiveness and efficiency of the proposed algorithm.

Improving spatial resolution with an edge-enhancement model for low-dose propagation-based X-ray phase-contrast computed tomography

Propagation-based X-ray phase-contrast computed tomography (PB-PCCT) has been increasingly popular for distinguishing low contrast tissues. Phase retrieval is an important step to quantitatively obtain the phase information before the tomographic reconstructions, while typical phase retrieval methods in PB-PCCT, such as homogenous transport of intensity equation (TIE-Hom), are essentially low-pass filters and thus improve the signal to noise ratio at the expense of the reduced spatial resolution of the reconstructed image. To improve the reconstructed spatial resolution, measured phase contrast projections with high edge enhancement and the phase projections retrieved by TIE-Hom were weighted summed and fed into an iterative tomographic algorithm within the framework of the adaptive steepest descent projections onto convex sets (ASD-POCS), which was employed for suppressing the image noise in low dose reconstructions because of the sparse-view scanning strategy or low exposure time for single phase contrast projection. The merging strategy decreases the accuracy of the linear model of PB-PCCT and would finally lead to the reconstruction failure in iterative reconstructions. Therefore, the additive median root prior is also introduced in the algorithm to partly increase the model accuracy. The reconstructed spatial resolution and noise performance can be flexibly balanced by a pair of antagonistic hyper-parameters. Validations were performed by the established phase-contrast Feldkamp-Davis-Kress, phase-retrieved Feldkamp-Davis-Kress, conventional ASD-POCS and the proposed enhanced ASD-POCS with a numerical phantom dataset and experimental biomaterial dataset. Simulation results show that the proposed algorithm outperforms the conventional ASD-POCS in spatial evaluation assessments such as root mean square error (a ratio of 9.78%), contrast to noise ratio (CNR) (a ratio of 7.46%), and also frequency evaluation assessments such as modulation transfer function (a ratio of 66.48% of MTF_50% (50% MTF value)), noise power spectrum (a ratio of 35.25% of f_50% (50% value of the Nyquist frequency)) and noise equivalent quanta (1-2 orders of magnitude at high frequencies). Experimental results again confirm the superiority of proposed strategy relative to the conventional one in terms of edge sharpness and CNR (an average increase of 67.35%).