1. Iterative phase contrast CT reconstruction with novel tomographic operator and data-driven prior
- Author
-
van Gogh, Stefano, Mukherjee, Subhadip, Xu, Jinqiu, Wang, Zhentian, Rawlik, Michał, Varga, Zsuzsanna, Alaifari, Rima, Schönlieb, Carola-Bibiane, Stampanoni, Marco, van Gogh, Stefano [0000-0001-5856-1124], Mukherjee, Subhadip [0000-0002-7957-8758], Wang, Zhentian [0000-0002-1646-1405], Rawlik, Michał [0000-0002-1232-4497], Apollo - University of Cambridge Repository, and University of Zurich
- Subjects
Medicine and health sciences ,FOS: Computer and information sciences ,1000 Multidisciplinary ,Computer and information sciences ,Biology and life sciences ,Phantoms, Imaging ,FOS: Physical sciences ,610 Medicine & health ,Breast Neoplasms ,Physical sciences ,Research and analysis methods ,10049 Institute of Pathology and Molecular Pathology ,Image Processing, Computer-Assisted ,Humans ,Female ,Tomography, X-Ray Computed ,Tomography ,Algorithms ,Research Article - Abstract
Funder: Swisslos Lottery Fund of canton Aargau, Funder: ETH Doc.Mobility Fellowship, Funder: Promedica Stiftung, 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.
- Published
- 2022
- Full Text
- View/download PDF