SPULTRA: Low-Dose CT Image Reconstruction With Joint Statistical and Learned Image Models

Low-dose CT image reconstruction has been a popular research topic in recent years. A typical reconstruction method based on post-log measurements is called penalized weighted-least squares (PWLS). Due to the underlying limitations of the post-log statistical model, the PWLS reconstruction quality is often degraded in low-dose scans. This paper investigates a shifted-Poisson (SP) model based likelihood function that uses the pre-log raw measurements that better represents the measurement statistics, together with a data-driven regularizer exploiting a Union of Learned TRAnsforms (SPULTRA). Both the SP induced data-fidelity term and the regularizer in the proposed framework are nonconvex. The proposed SPULTRA algorithm uses quadratic surrogate functions for the SP induced data-fidelity term. Each iteration involves a quadratic subproblem for updating the image, and a sparse coding and clustering subproblem that has a closed-form solution. The SPULTRA algorithm has a similar computational cost per iteration as its recent counterpart PWLS-ULTRA that uses post-log measurements, and it provides better image reconstruction quality than PWLS-ULTRA, especially in low-dose scans.

[1]  Yoram Bresler,et al.  Structured Overcomplete Sparsifying Transform Learning with Convergence Guarantees and Applications , 2015, International Journal of Computer Vision.

[2]  Yoram Bresler,et al.  Efficient Blind Compressed Sensing Using Sparsifying Transforms with Convergence Guarantees and Application to Magnetic Resonance Imaging , 2015, SIAM J. Imaging Sci..

[3]  Guang-Hong Chen,et al.  Time-Resolved Interventional Cardiac C-arm Cone-Beam CT: An Application of the PICCS Algorithm , 2012, IEEE Transactions on Medical Imaging.

[4]  Wufan Chen,et al.  Variance analysis of x-ray CT sinograms in the presence of electronic noise background. , 2012, Medical physics.

[5]  Jing Wang,et al.  Virtual Colonoscopy Screening with Ultra Low-Dose CT and Less-Stressful Bowel Preparation: A computer simulation study. , 2008, IEEE transactions on nuclear science.

[6]  J. Fessler.  Statistical Image Reconstruction Methods for Transmission Tomography , 2000 .

[7]  Bastian Goldlücke,et al.  Variational Analysis , 2014, Computer Vision, A Reference Guide.

[8]  Jeffrey A. Fessler,et al.  Relaxed Linearized Algorithms for Faster X-Ray CT Image Reconstruction , 2015, IEEE Transactions on Medical Imaging.

[9]  J. Fessler,et al.  Modelling the physics in the iterative reconstruction for transmission computed tomography , 2013, Physics in medicine and biology.

[10]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[11]  M. Shiung,et al.  Development and Validation of a Practical Lower-Dose-Simulation Tool for Optimizing Computed Tomography Scan Protocols , 2012, Journal of computer assisted tomography.

[12]  Jean-Baptiste Thibault,et al.  Modeling and Pre-Treatment of Photon-Starved CT Data for Iterative Reconstruction , 2017, IEEE Trans. Medical Imaging.

[13]  Jeffrey A. Fessler,et al.  3D Forward and Back-Projection for X-Ray CT Using Separable Footprints , 2010, IEEE Transactions on Medical Imaging.

[14]  Yoram Bresler,et al.  $\ell_{0}$ Sparsifying Transform Learning With Efficient Optimal Updates and Convergence Guarantees , 2015, IEEE Transactions on Signal Processing.

[15]  Jeffrey A. Fessler,et al.  PWLS-ULTRA: An Efficient Clustering and Learning-Based Approach for Low-Dose 3D CT Image Reconstruction , 2017, IEEE Transactions on Medical Imaging.

[16]  Jin Liu,et al.  3D Feature Constrained Reconstruction for Low-Dose CT Imaging , 2018, IEEE Transactions on Circuits and Systems for Video Technology.

[17]  Jeffrey A. Fessler,et al.  Modeling Mixed Poisson-Gaussian Noise in Statistical Image Reconstruction for X-Ray CT , 2016 .

[18]  Jean-Baptiste Thibault,et al.  Modeling and Pre-Treatment of Photon-Starved CT Data for Iterative Reconstruction , 2017, IEEE Transactions on Medical Imaging.

[19]  Jaejun Yoo,et al.  Deep Convolutional Framelet Denosing for Low-Dose CT via Wavelet Residual Network , 2018, IEEE Transactions on Medical Imaging.

[20]  Xuanqin Mou,et al.  Tensor-Based Dictionary Learning for Spectral CT Reconstruction , 2017, IEEE Transactions on Medical Imaging.

[21]  Rosemary A. Renaut,et al.  Least squares problems with inequality constraints as quadratic constraints , 2010 .

[22]  Patrick J. La Rivière,et al.  Penalized-likelihood sinogram restoration for computed tomography , 2006, IEEE Transactions on Medical Imaging.

[23]  J. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[24]  Zhaoying Bian,et al.  Optimizing a Parameterized Plug-and-Play ADMM for Iterative Low-Dose CT Reconstruction , 2019, IEEE Transactions on Medical Imaging.

[25]  Jingyan Xu,et al.  Electronic Noise Modeling in Statistical Iterative Reconstruction , 2009, IEEE Transactions on Image Processing.

[26]  Hu Chen,et al.  LEARN: Learned Experts’ Assessment-Based Reconstruction Network for Sparse-Data CT , 2017, IEEE Transactions on Medical Imaging.

[27]  C. Bouman,et al.  A recursive filter for noise reduction in statistical iterative tomographic imaging , 2006, Electronic Imaging.

[28]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 1999, IEEE Transactions on Medical Imaging.

[29]  W P Segars,et al.  Realistic CT simulation using the 4D XCAT phantom. , 2008, Medical physics.

[30]  J. Webster Stayman,et al.  Penalized-Likelihood Reconstruction With High-Fidelity Measurement Models for High-Resolution Cone-Beam Imaging , 2017, IEEE Transactions on Medical Imaging.

[31]  Jeffrey A. Fessler,et al.  Generalizing the Optimized Gradient Method for Smooth Convex Minimization , 2016, SIAM J. Optim..

[32]  Jeffrey A. Fessler,et al.  Efficient and accurate likelihood for iterative image reconstruction in x-ray computed tomography , 2003, SPIE Medical Imaging.

[33]  Paul Kinahan,et al.  Attenuation correction for a combined 3D PET/CT scanner. , 1998, Medical physics.

[34]  Daniel Kolditz,et al.  Iterative reconstruction methods in X-ray CT. , 2012, Physica medica : PM : an international journal devoted to the applications of physics to medicine and biology : official journal of the Italian Association of Biomedical Physics.

[35]  Lin Fu,et al.  Ultra-low dose CT attenuation correction for PET/CT: analysis of sparse view data acquisition and reconstruction algorithms , 2015, Physics in medicine and biology.

[36]  Jong Chul Ye,et al.  Deep Residual Learning for Compressed Sensing CT Reconstruction via Persistent Homology Analysis , 2016, ArXiv.

[37]  Jonas Adler,et al.  Learned Primal-Dual Reconstruction , 2017, IEEE Transactions on Medical Imaging.

[38]  Jeffrey A. Fessler,et al.  Adaptive sparse modeling and shifted-poisson likelihood based approach for low-dosect image reconstruction , 2017, 2017 IEEE 27th International Workshop on Machine Learning for Signal Processing (MLSP).

[39]  Trond Steihaug,et al.  An interior-point trust-region-based method for large-scale non-negative regularization , 2002 .

[40]  Lin Fu,et al.  Comparison Between Pre-Log and Post-Log Statistical Models in Ultra-Low-Dose CT Reconstruction , 2017, IEEE Transactions on Medical Imaging.

[41]  J. Hsieh.  Adaptive streak artifact reduction in computed tomography resulting from excessive x-ray photon noise. , 1998, Medical physics.

[42]  Bruce R. Whiting,et al.  Signal statistics in x-ray computed tomography , 2002, SPIE Medical Imaging.

[43]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[44]  Qianjin Feng,et al.  Low-dose computed tomography image restoration using previous normal-dose scan. , 2011, Medical physics.

[45]  Jeffrey A. Fessler,et al.  Regularization Designs for Uniform Spatial Resolution and Noise Properties in Statistical Image Reconstruction for 3-D X-ray CT , 2015, IEEE Transactions on Medical Imaging.

[46]  Yoram Bresler,et al.  Learning Sparsifying Transforms , 2013, IEEE Transactions on Signal Processing.

[47]  Jean-Baptiste Thibault,et al.  Ultra low dose CT for attenuation correction in PET/CT , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[48]  W. Marsden.  I and J , 2012 .

[49]  Åke Björck,et al.  Numerical methods for least square problems , 1996 .

[50]  Guang-Hong Chen,et al.  Unbiased statistical image reconstruction in low-dose CT , 2019, Medical Imaging.

[51]  Yoram Bresler,et al.  Data-Driven Learning of a Union of Sparsifying Transforms Model for Blind Compressed Sensing , 2015, IEEE Transactions on Computational Imaging.

[52]  Lin Fu,et al.  2.5D dictionary learning based computed tomography reconstruction , 2016, SPIE Defense + Security.

[53]  F. Giannessi.  Variational Analysis and Generalized Differentiation , 2006 .

[54]  P. L. La Riviere.  Penalized-likelihood sinogram smoothing for low-dose CT. , 2005, Medical physics.

[55]  Jeffrey A. Fessler,et al.  Combining Ordered Subsets and Momentum for Accelerated X-Ray CT Image Reconstruction , 2015, IEEE Transactions on Medical Imaging.

[56]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[57]  Zhengrong Liang,et al.  Applications of nonlocal means algorithm in low‐dose X‐ray CT image processing and reconstruction: A review , 2017, Medical physics.

[58]  B. Mordukhovich.  Variational analysis and generalized differentiation , 2006 .

[59]  Jeffrey A. Fessler,et al.  Low dose CT image reconstruction with learned sparsifying transform , 2016, 2016 IEEE 12th Image, Video, and Multidimensional Signal Processing Workshop (IVMSP).

[60]  Jing Wang,et al.  Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography , 2006, IEEE Transactions on Medical Imaging.

[61]  J Trzasko,et al.  Nonconvex prior image constrained compressed sensing (NCPICCS): theory and simulations on perfusion CT. , 2011, Medical physics.

[62]  Bruno De Man,et al.  Comparison Between Pre-Log and Post-Log Statistical Models in Ultra-Low-Dose CT Reconstruction. , 2017, IEEE transactions on medical imaging.

[63]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[64]  Lei Zhang,et al.  Low-Dose X-ray CT Reconstruction via Dictionary Learning , 2012, IEEE Transactions on Medical Imaging.

[65]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Hybrid Poisson/polynomial Objective Functions for Tomographic Image Reconstruction from Transmission Scans , 2022 .

[66]  Quanzheng Li,et al.  Iterative Low-Dose CT Reconstruction With Priors Trained by Artificial Neural Network , 2017, IEEE Transactions on Medical Imaging.