Regularizer Performance for SparseCT Image Reconstruction With Practical Subsampling

Matthew Muckley, Baiyu Chen, Thomas Vahle, Aaron Sodickson, Florian Knoll, Daniel K. Sodickson and Ricardo Otazo

DOI:10.12059/Fully3D.2017-11-3108005

Published in:Fully3D 2017 Proceedings

Pages:572-575

Keywords:
Low-dose X-ray computed tomography (CT) is a major area of research due to the diagnostic capability of CT counterbalanced by the risk of radiation exposure. The standard method for reducing the dose is to decrease the tube current, but this negatively impacts image quality at high dose reduction factors due to photon starvation effects. We investigate an alternative paradigm, called SparseCT, in which a multi-slit collimator (MSC) is placed between the source and the patient. Interrupting the X-ray beam in this way reduces radiation dose to the patient and produces undersampled data from which the image can be estimated using sparse image reconstruction techniques. However, the MSC introduces a number of new considerations, including penumbra effects that require wider slits and more spacing between slits in order to optimize dose efficiency and beam separation. These design choices reduce incoherence in the undersampling scheme and affect the performance of standard edge-preserving/sparsity-promoting regularizers. Here, we simulate these effects in the ideal setting where penumbra effects do not increase the noise in the sinogram. For modest 4-fold retrospective undersampling factors, we observe a slight degradation in an abdominal scan between incoherent-optimal and practical MSC designs for SparseCT. In the future, we plan to use the simulation to inform MSC design prior to fabrication.
Matthew Muckley
New York University School of Medicine, USA
Baiyu Chen
New York University School of Medicine, USA
Thomas Vahle
Siemens Healthcare GmbH, Germany
Aaron Sodickson
Harvard Medical School, USA
Florian Knoll
New York University School of Medicine, USA
Daniel K. Sodickson
New York University School of Medicine, USA
Ricardo Otazo
New York University School of Medicine, USA
  1. J.-B. Thibault, K. Sauer, C. Bouman, and J. Hsieh, “A three-dimensional statistical approach to improved image quality for multi-slice helical CT,” Med. Phys., vol. 34, no. 11, pp. 4526–44, Nov. 2007.
  2. J. Luo, H. Eri, A. Can, S. Ramani, L. Fu, and B. De Man, “2.5D dictionary learning based computed tomography reconstruction,” in spie9847, 2016, p. 98470L.
  3. T. Koesters, F. Knoll, A. Sodickson, D. Sodickson, and R. Otazo, “SparseCT: Interrupted-beam acquisition and sparse reconstruction for radiation dose reduction,” in SPIE Med. Im. 2017, 2017, to appear.
  4. E. J. Candes and M. B. Wakin, “An introduction to compressive sampling,” IEEE Sig. Proc. Mag., vol. 25, no. 2, pp. 21–30, Mar. 2008.
  5. B. Chen, M. J. Muckley, T. O’Donnell, A. Sodickson, T. Flohr, K. Stierstorfer, B. Schmidt, F. Knoll, A. Primak, D. Faul, R. Otazo, and D. Sodickson, “Realistic undersampling model for compressed sensing using a multi-slit collimator,” in Proc. Intl. Mtg. on Fully 3D Image Recon. in Rad. and Nuc. Med, 2017, submitted.
  6. H. Chandarana, A. M. Doshi, A. Shanbhogue, J. S. Babb, M. T. Bruno, T. Zhao, E. Raithel, M. O. Zenge, G. Li, and R. Otazo, “Three-dimensional MR cholangiopancreatography in a breath hold with sparsity-based reconstruction of highly undersampled data,” Radiology, pp. 585–94, 2016.
  7. A. R. De Pierro, “A modified expectation maximization algorithm for penalized likelihood estimation in emission tomography,” IEEE Trans. Med. Imag., vol. 14, no. 1, pp. 132–7, Mar. 1995.
  8. M. J. Muckley, D. C. Noll, and J. A. Fessler, “Fast parallel MR image reconstruction via B1-based, adaptive restart, iterative soft thresholding algorithms (BARISTA),” IEEE Trans. Med. Imag., vol. 34, no. 2, pp. 578–88, Feb. 2015.
  9. A. Beck and M. Teboulle, “A fast iterative shrinkage-thresholding algorithm for linear inverse problems,” SIAM J. Imaging Sci., vol. 2, no. 1, pp. 183–202, 2009.
  10. D. Kim, D. Pal, J.-B. Thibault, and J. A. Fessler, “Accelerating ordered subsets image reconstruction for X-ray CT using spatially non-uniform optimization transfer,” IEEE Trans. Med. Imag., vol. 32, no. 11, pp. 1965–78, Nov. 2013.