Sparse-view x-ray CT reconstruction via total generalized variation regularization

Sparse-view CT reconstruction algorithms via total variation (TV) optimize the data iteratively on the basis of a noise- and artifact-reducing model, resulting in significant radiation dose reduction while maintaining image quality. However, the piecewise constant assumption of TV minimization often leads to the appearance of noticeable patchy artifacts in reconstructed images. To obviate this drawback, we present a penalized weighted least-squares (PWLS) scheme to retain the image quality by incorporating the new concept of total generalized variation (TGV) regularization. We refer to the proposed scheme as 'PWLS-TGV' for simplicity. Specifically, TGV regularization utilizes higher order derivatives of the objective image, and the weighted least-squares term considers data-dependent variance estimation, which fully contribute to improving the image quality with sparse-view projection measurement. Subsequently, an alternating optimization algorithm was adopted to minimize the associative objective function. To evaluate the PWLS-TGV method, both qualitative and quantitative studies were conducted by using digital and physical phantoms. Experimental results show that the present PWLS-TGV method can achieve images with several noticeable gains over the original TV-based method in terms of accuracy and resolution properties.

[1]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[2]  Grant T. Gullberg,et al.  Total variation regulated EM algorithm , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

[3]  Gengsheng L. Zeng,et al.  Total variation regulated EM algorithm [SPECT reconstruction] , 1999 .

[4]  M. Persson,et al.  Total variation norm for three-dimensional iterative reconstruction in limited view angle tomography , 2001, Physics in medicine and biology.

[5]  A. Bovik,et al.  A universal image quality index , 2002, IEEE Signal Processing Letters.

[6]  Jeffrey A. Fessler,et al.  Edge-preserving tomographic reconstruction with nonlocal regularization , 2002 .

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

[8]  R. Lewis,et al.  Medical phase contrast x-ray imaging: current status and future prospects. , 2004, Physics in medicine and biology.

[9]  Dan Xu,et al.  Detection of gamma ray polarization using a 3D position sensitive CdZnTe detector , 2004, IEEE Symposium Conference Record Nuclear Science 2004..

[10]  M. Kalra,et al.  Techniques and applications of automatic tube current modulation for CT. , 2004, Radiology.

[11]  Hongbing Lu,et al.  Nonlinear sinogram smoothing for low-dose X-ray CT , 2004 .

[12]  Fang Xu,et al.  Accelerating popular tomographic reconstruction algorithms on commodity PC graphics hardware , 2005, IEEE Transactions on Nuclear Science.

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

[14]  M. Knaup,et al.  Statistical Cone-Beam CT Image Reconstruction using the Cell Broadband Engine , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[15]  C. McCollough,et al.  CT dose reduction and dose management tools: overview of available options. , 2006, Radiographics : a review publication of the Radiological Society of North America, Inc.

[16]  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.

[17]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[18]  A. Einstein,et al.  Estimating risk of cancer associated with radiation exposure from 64-slice computed tomography coronary angiography. , 2007, JAMA.

[19]  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.

[20]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[21]  Zhengrong Liang,et al.  An experimental study on the noise properties of x-ray CT sinogram data in Radon space , 2008, Physics in medicine and biology.

[22]  L. Xing,et al.  Iterative image reconstruction for CBCT using edge-preserving prior. , 2008, Medical physics.

[23]  Hengyong Yu,et al.  Compressed sensing based interior tomography , 2009, Physics in medicine and biology.

[24]  Jie Tang,et al.  Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms , 2009, Physics in medicine and biology.

[25]  Zhengrong Liang,et al.  Recent development of Low-dose X-ray Cone-beam computed tomography , 2010 .

[26]  Peter Cloetens,et al.  Low-dose x-ray phase-contrast and absorption CT using equally sloped tomography , 2010, Physics in medicine and biology.

[27]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[28]  Qianjin Feng,et al.  Generalized Gibbs priors based positron emission tomography reconstruction. , 2010, Computers in biology and medicine.

[29]  Karl Kunisch,et al.  Total Generalized Variation , 2010, SIAM J. Imaging Sci..

[30]  Xiaochuan Pan,et al.  Evaluation of sparse-view reconstruction from flat-panel-detector cone-beam CT , 2010, Physics in medicine and biology.

[31]  W. Segars,et al.  4D XCAT phantom for multimodality imaging research. , 2010, Medical physics.

[32]  Ge Wang The meaning of interior tomography , 2013, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

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

[34]  T. Pock,et al.  Second order total generalized variation (TGV) for MRI , 2011, Magnetic resonance in medicine.

[35]  Steve B. Jiang,et al.  Low-dose CT reconstruction via edge-preserving total variation regularization. , 2010, Physics in medicine and biology.

[36]  C. McCollough,et al.  Optimal tube potential for radiation dose reduction in pediatric CT: principles, clinical implementations, and pitfalls. , 2011, Radiographics : a review publication of the Radiological Society of North America, Inc.

[37]  Jing Wang,et al.  Inverse determination of the penalty parameter in penalized weighted least-squares algorithm for noise reduction of low-dose CBCT. , 2011, Medical physics.

[38]  Steve B. Jiang,et al.  Low-dose CT reconstruction via edge-preserving total variation regularization , 2010, Physics in medicine and biology.

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

[40]  Xiao Han,et al.  Optimization-based reconstruction of sparse images from few-view projections , 2012, Physics in medicine and biology.

[41]  Jin Sung Kim,et al.  Fast compressed sensing-based CBCT reconstruction using Barzilai-Borwein formulation for application to on-line IGRT. , 2012, Medical physics.

[42]  Zhengrong Liang,et al.  Adaptive-weighted total variation minimization for sparse data toward low-dose x-ray computed tomography image reconstruction , 2012, Physics in medicine and biology.

[43]  Jie Tang,et al.  Prior image constrained compressed sensing: implementation and performance evaluation. , 2011, Medical physics.

[44]  Lei Zhu,et al.  Accelerated barrier optimization compressed sensing (ABOCS) reconstruction for cone-beam CT: Phantom studies. , 2012, Medical physics.

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

[46]  Zhengrong Liang,et al.  Iterative image reconstruction for cerebral perfusion CT using a pre-contrast scan induced edge-preserving prior , 2012, Physics in medicine and biology.

[47]  Yiqiu Dong,et al.  Spatially dependent regularization parameter selection in total generalized variation models for image restoration , 2013, Int. J. Comput. Math..

[48]  Xin Jin,et al.  A limited-angle CT reconstruction method based on anisotropic TV minimization , 2013, Physics in medicine and biology.

[49]  Emil Y. Sidky,et al.  Quantifying Admissible Undersampling for Sparsity-Exploiting Iterative Image Reconstruction in X-Ray CT , 2011, IEEE Transactions on Medical Imaging.

[50]  Zhaoying Bian,et al.  Iterative Image Reconstruction for Sparse-View CT Using Normal-Dose Image Induced Total Variation Prior , 2013, PloS one.

[51]  Jakob S Jørgensen,et al.  Few-view single photon emission computed tomography (SPECT) reconstruction based on a blurred piecewise constant object model. , 2013, Physics in medicine and biology.

[52]  Guang-Hong Chen,et al.  Characterization of statistical prior image constrained compressed sensing (PICCS): II. Application to dose reduction. , 2013, Medical physics.

[53]  Jing Huang,et al.  SR-NLM: A sinogram restoration induced non-local means image filtering for low-dose computed tomography , 2013, Comput. Medical Imaging Graph..

[54]  Jianhua Ma,et al.  Total Variation-Stokes Strategy for Sparse-View X-ray CT Image Reconstruction , 2014, IEEE Transactions on Medical Imaging.