Data-Driven Gradient Regularization for Quasi-Newton Optimization in Iterative Grating Interferometry CT Reconstruction.

Grating interferometry CT (GI-CT) is a promising technology that could play an important role in future breast cancer imaging. Thanks to its sensitivity to refraction and small-angle scattering, GI-CT could augment the diagnostic content of conventional absorption-based CT. However, reconstructing GI-CT tomographies is a complex task because of ill problem conditioning and high noise amplitudes. It has previously been shown that combining data-driven regularization with iterative reconstruction is promising for tackling challenging inverse problems in medical imaging. In this work, we present an algorithm that allows seamless combination of data-driven regularization with quasi-Newton solvers, which can better deal with ill-conditioned problems compared to gradient descent-based optimization algorithms. Contrary to most available algorithms, our method applies regularization in the gradient domain rather than in the image domain. This comes with a crucial advantage when applied in conjunction with quasi-Newton solvers: the Hessian is approximated solely based on denoised data. We apply the proposed method, which we call GradReg, to both conventional breast CT and GI-CT and show that both significantly benefit from our approach in terms of dose efficiency. Moreover, our results suggest that thanks to its sharper gradients that carry more high spatial-frequency content, GI-CT can benefit more from GradReg compared to conventional breast CT. Crucially, GradReg can be applied to any image reconstruction task which relies on gradient-based updates.

Towards clinical-dose grating interferometry breast CT with fused intensity-based iterative reconstruction.

X-ray grating interferometry CT (GI-CT) is an emerging imaging modality which provides three complementary contrasts that could increase the diagnostic content of clinical breast CT: absorption, phase, and dark-field. Yet, reconstructing the three image channels under clinically compatible conditions is challenging because of severe ill-conditioning of the tomographic reconstruction problem. In this work we propose to solve this problem with a novel reconstruction algorithm that assumes a fixed relation between the absorption and the phase-contrast channel to reconstruct a single image by automatically fusing the absorption and phase channels. The results on both simulations and real data show that, enabled by the proposed algorithm, GI-CT outperforms conventional CT at a clinical dose.

The choice of an autocorrelation length in dark-field lung imaging.

Respiratory diseases are one of the most common causes of death, and their early detection is crucial for prompt treatment. X-ray dark-field radiography (XDFR) is a promising tool to image objects with unresolved micro-structures such as lungs. Using Talbot-Lau XDFR, we imaged inflated porcine lungs together with Polymethylmethacrylat (PMMA) microspheres (in air) of diameter sizes between 20 and 500 [Formula: see text] over an autocorrelation range of 0.8-5.2 [Formula: see text]. The results indicate that the dark-field extinction coefficient of porcine lungs is similar to that of densely-packed PMMA spheres with diameter of [Formula: see text], which is approximately the mean alveolar structure size. We evaluated that, in our case, the autocorrelation length would have to be limited to [Formula: see text] in order to image [Formula: see text]-thick lung tissue without critical visibility reduction (signal saturation). We identify the autocorrelation length to be the critical parameter of an interferometer that allows to avoid signal saturation in clinical lung dark-field imaging.

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.

Intensity-based iterative reconstruction for helical grating interferometry breast CT with static grating configuration.

Grating interferometry breast computed tomography (GI-BCT) has the potential to provide enhanced soft tissue contrast and to improve visualization of cancerous lesions for breast imaging. However, with a conventional scanning protocol, a GI-BCT scan requires longer scanning time and higher operation complexity compared to conventional attenuation-based CT. This is mainly due to multiple grating movements at every projection angle, so-called phase stepping, which is used to retrieve attenuation, phase, and scattering (dark-field) signals. To reduce the measurement time and complexity and extend the field of view, we have adopted a helical GI-CT setup and present here the corresponding tomographic reconstruction algorithm. This method allows simultaneous reconstruction of attenuation, phase contrast, and scattering images while avoiding grating movements. Experiments on simulated phantom and real initial intensity, visibility and phase maps are provided to validate our method.

INSIDEnet: Interpretable NonexpanSIve Data-Efficient network for denoising in grating interferometry breast CT.

PURPOSE: Breast cancer is the most common malignancy in women. Unfortunately, current breast imaging techniques all suffer from certain limitations: they are either not fully three dimensional, have an insufficient resolution or low soft-tissue contrast. Grating interferometry breast computed tomography (GI-BCT) is a promising X-ray phase contrast modality that could overcome these limitations by offering high soft-tissue contrast and excellent three-dimensional resolution. To enable the transition of this technology to clinical practice, dedicated data-processing algorithms must be developed in order to effectively retrieve the signals of interest from the measured raw data. METHODS: This article proposes a novel denoising algorithm that can cope with the high-noise amplitudes and heteroscedasticity which arise in GI-BCT when operated in a low-dose regime to effectively regularize the ill-conditioned GI-BCT inverse problem. We present a data-driven algorithm called INSIDEnet, which combines different ideas such as multiscale image processing, transform-domain filtering, transform learning, and explicit orthogonality to build an Interpretable NonexpanSIve Data-Efficient network (INSIDEnet). RESULTS: We apply the method to simulated breast phantom datasets and to real data acquired on a GI-BCT prototype and show that the proposed algorithm outperforms traditional state-of-the-art filters and is competitive with deep neural networks. The strong inductive bias given by the proposed model's architecture allows to reliably train the algorithm with very limited data while providing high model interpretability, thus offering a great advantage over classical convolutional neural networks (CNNs). CONCLUSIONS: The proposed INSIDEnet is highly data-efficient, interpretable, and outperforms state-of-the-art CNNs when trained on very limited training data. We expect the proposed method to become an important tool as part of a dedicated plug-and-play GI-BCT reconstruction framework, needed to translate this promising technology to the clinics.