Improved Compressed Sensing-Based Algorithm for Sparse-View CT Image Reconstruction

In computed tomography (CT), there are many situations where reconstruction has to be performed with sparse-view data. In sparse-view CT imaging, strong streak artifacts may appear in conventionally reconstructed images due to limited sampling rate that compromises image quality. Compressed sensing (CS) algorithm has shown potential to accurately recover images from highly undersampled data. In the past few years, total-variation-(TV-) based compressed sensing algorithms have been proposed to suppress the streak artifact in CT image reconstruction. In this paper, we propose an efficient compressed sensing-based algorithm for CT image reconstruction from few-view data where we simultaneously minimize three parameters: the ℓ 1 norm, total variation, and a least squares measure. The main feature of our algorithm is the use of two sparsity transforms—discrete wavelet transform and discrete gradient transform. Experiments have been conducted using simulated phantoms and clinical data to evaluate the performance of the proposed algorithm. The results using the proposed scheme show much smaller streaking artifacts and reconstruction errors than other conventional methods.

[1]  Michael Unser,et al.  Statistical analysis of functional MRI data in the wavelet domain , 1998, IEEE Transactions on Medical Imaging.

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

[3]  Ishpreet Singh Virk,et al.  Medical Image Denoising In The Wavelet Domain Using Haar And DB3 Filtering , 2012 .

[4]  Xiaoming Huo,et al.  Combined image representation using edgelets and wavelets , 1999, Optics & Photonics.

[5]  Tutut Herawan,et al.  Computational and mathematical methods in medicine. , 2006, Computational and mathematical methods in medicine.

[6]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[7]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[8]  M. Jiang,et al.  Ordered-subset simultaneous algebraic reconstruction techniques (OS-SART) , 2004 .

[9]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[10]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

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

[12]  Harish P. Hiriyannaiah,et al.  X-ray computed tomography for medical imaging , 1997, IEEE Signal Process. Mag..

[13]  Shiqian Ma,et al.  An efficient algorithm for compressed MR imaging using total variation and wavelets , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.

[14]  S. Leng,et al.  High temporal resolution and streak-free four-dimensional cone-beam computed tomography , 2008, Physics in medicine and biology.

[15]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[16]  Vivek Agarwal,et al.  Total Variation Regularization and L-curve method for the selection of regularization parameter , 2003 .

[17]  Per Christian Hansen,et al.  AIR Tools - A MATLAB package of algebraic iterative reconstruction methods , 2012, J. Comput. Appl. Math..

[18]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[19]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[20]  Misha Elena Kilmer,et al.  Choosing Regularization Parameters in Iterative Methods for Ill-Posed Problems , 2000, SIAM J. Matrix Anal. Appl..

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

[22]  E. Miller,et al.  Efficient determination of multiple regularization parameters in a generalized L-curve framework , 2002 .

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

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

[25]  G T Herman,et al.  ART: mathematics and applications. A report on the mathematical foundations and on the applicability to real data of the algebraic reconstruction techniques. , 1973, Journal of theoretical biology.

[26]  Shuqian Luo,et al.  A compressed sensing-based iterative algorithm for CT reconstruction and its possible application to phase contrast imaging , 2011, Biomedical engineering online.

[27]  D. Brenner,et al.  Cancer risks from diagnostic radiology. , 2008, The British journal of radiology.

[28]  Rebecca S Lewis,et al.  Projected cancer risks from computed tomographic scans performed in the United States in 2007. , 2009, Archives of internal medicine.