Tomographic image reconstruction via estimation of sparse unidirectional gradients

Since computed tomography (CT) was developed over 35 years ago, new mathematical ideas and computational algorithms have been continuingly elaborated to improve the quality of reconstructed images. In recent years, a considerable effort can be noticed to apply the sparse solution of underdetermined system theory to the reconstruction of CT images from undersampled data. Its significance stems from the possibility of obtaining good quality CT images from low dose projections. Among diverse approaches, total variation (TV) minimizing 2D gradients of an image, seems to be the most popular method. In this paper, a new method for CT image reconstruction via sparse gradients estimation (SGE), is proposed. It consists in estimating 1D gradients specified in four directions using the iterative reweighting algorithm. To investigate its properties and to compare it with TV and other related methods, numerical simulations were performed according to the Monte Carlo scheme, using the Shepp-Logan and more realistic brain phantoms scanned at 9-60 directions in the range from 0 to 179°, with measurement data disturbed by additive Gaussians noise characterized by the relative level of 0.1%, 0.2%, 0.5%, 1%, 2% and 5%. The accuracy of image reconstruction was assessed in terms of the relative root-mean-square (RMS) error. The results show that the proposed SGE algorithm has returned more accurate images than TV for the cases fulfilling the sparsity conditions. Particularly, it preserves sharp edges of regions representing different tissues or organs and yields images of much better quality reconstructed from a small number of projections disturbed by relatively low measurement noise.

[1]  Dong Liang,et al.  Highly Undersampled Magnetic Resonance Image Reconstruction Using Two-Level Bregman Method With Dictionary Updating , 2013, IEEE Transactions on Medical Imaging.

[2]  Li Liu,et al.  Nonparametric optimization of constrained total variation for tomography reconstruction , 2013, Comput. Biol. Medicine.

[3]  Daniele Funaro,et al.  Domain decomposition methods for pseudo spectral approximations , 1987 .

[4]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[5]  A Tikhonov,et al.  Solution of Incorrectly Formulated Problems and the Regularization Method , 1963 .

[6]  Michael P. Friedlander,et al.  Probing the Pareto Frontier for Basis Pursuit Solutions , 2008, SIAM J. Sci. Comput..

[7]  A. Gajos,et al.  Compressive sensing of signals generated in plastic scintillators in a novel J-PET instrument , 2015, 1503.05188.

[8]  Wotao Yin,et al.  Bregman Iterative Algorithms for (cid:2) 1 -Minimization with Applications to Compressed Sensing ∗ , 2008 .

[9]  Robert Cierniak,et al.  Regularized Image Reconstruction from Projections Method , 2014, 2014 International Conference on Engineering and Telecommunication.

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

[11]  Yves Goussard,et al.  GPU-accelerated regularized iterative reconstruction for few-view cone beam CT. , 2015, Medical physics.

[12]  Ge Wang,et al.  Simultaneous CT-MRI Reconstruction for Constrained Imaging Geometries Using Structural Coupling and Compressive Sensing , 2016, IEEE Transactions on Biomedical Engineering.

[13]  David A. Clausi,et al.  Sparse Reconstruction of Breast MRI Using Homotopic $L_0$ Minimization in a Regional Sparsified Domain , 2013, IEEE Transactions on Biomedical Engineering.

[14]  Predrag Sukovic,et al.  Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography , 2000, IEEE Transactions on Medical Imaging.

[15]  Armando Manduca,et al.  Highly Undersampled Magnetic Resonance Image Reconstruction via Homotopic $\ell_{0}$ -Minimization , 2009, IEEE Transactions on Medical Imaging.

[16]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

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

[18]  Alireza Entezari,et al.  Gradient-based sparse approximation for computed tomography , 2015, 2015 IEEE 12th International Symposium on Biomedical Imaging (ISBI).

[19]  Huazhong Shu,et al.  Median prior constrained TV algorithm for sparse view low-dose CT reconstruction , 2015, Comput. Biol. Medicine.

[20]  Stephen J. Wright,et al.  Computational Methods for Sparse Solution of Linear Inverse Problems , 2010, Proceedings of the IEEE.

[21]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation. , 2010, Medical physics.

[22]  Adam G. Polak,et al.  An Error-Minimizing Approach to Regularization in Indirect Measurements , 2010, IEEE Transactions on Instrumentation and Measurement.

[23]  Bhaskar D. Rao,et al.  Sparse solutions to linear inverse problems with multiple measurement vectors , 2005, IEEE Transactions on Signal Processing.

[24]  Yi Sun,et al.  Optimization for limited angle tomography in medical image processing , 2011, Pattern Recognit..

[25]  Klaas Paul Pruessmann,et al.  Realistic Analytical Phantoms for Parallel Magnetic Resonance Imaging , 2012, IEEE Transactions on Medical Imaging.

[26]  J. C. Ye,et al.  Projection reconstruction MR imaging using FOCUSS , 2007, Magnetic resonance in medicine.

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

[28]  Junyan Rong,et al.  CT reconstruction from few-views with anisotropic edge-guided total variance , 2016 .

[29]  Jean-Luc Starck,et al.  Sparse Solution of Underdetermined Systems of Linear Equations by Stagewise Orthogonal Matching Pursuit , 2012, IEEE Transactions on Information Theory.

[30]  Janusz Mroczka,et al.  Performance analysis of regularization algorithms used for image reconstruction in computed tomography , 2013 .

[31]  Søren Holdt Jensen,et al.  Implementation of an optimal first-order method for strongly convex total variation regularization , 2011, ArXiv.

[32]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[33]  Michael Unser,et al.  Joint image reconstruction and segmentation using the Potts model , 2014, 1405.5850.

[34]  Jie Tian,et al.  Sparse reconstruction for quantitative bioluminescence tomography based on the incomplete variables truncated conjugate gradient method. , 2010, Optics express.

[35]  Bhaskar D. Rao,et al.  Extension of SBL Algorithms for the Recovery of Block Sparse Signals With Intra-Block Correlation , 2012, IEEE Transactions on Signal Processing.

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

[37]  Rick Chartrand,et al.  Exact Reconstruction of Sparse Signals via Nonconvex Minimization , 2007, IEEE Signal Processing Letters.

[38]  Marc Moonen,et al.  Enhancing sparsity in linear prediction of speech by iteratively reweighted 1-norm minimization , 2010, 2010 IEEE International Conference on Acoustics, Speech and Signal Processing.

[39]  Bhaskar D. Rao,et al.  Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm , 1997, IEEE Trans. Signal Process..

[40]  Hiêp Quang Luong,et al.  Augmented Lagrangian based reconstruction of non-uniformly sub-Nyquist sampled MRI data , 2011, Signal Process..

[41]  Stephen P. Boyd,et al.  Enhancing Sparsity by Reweighted ℓ1 Minimization , 2007, 0711.1612.

[42]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

[43]  Hiroyuki Kudo,et al.  Sparsity-constrained three-dimensional image reconstruction for C-arm angiography , 2015, Comput. Biol. Medicine.

[44]  Tie Zhang,et al.  Three penalized EM-type algorithms for PET image reconstruction , 2012, Comput. Biol. Medicine.

[45]  Arthur E. Hoerl,et al.  Application of ridge analysis to regression problems , 1962 .

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

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

[48]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[49]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

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

[51]  M. Lai,et al.  An Unconstrained $\ell_q$ Minimization with $0q\leq1$ for Sparse Solution of Underdetermined Linear Systems , 2011 .

[52]  K. Tanabe.  Projection method for solving a singular system of linear equations and its applications , 1971 .

[53]  Luminita A. Vese,et al.  Expectation maximization and total variation-based model for computed tomography reconstruction from undersampled data , 2011, Medical Imaging.

[54]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

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

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

[57]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

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

[59]  Wotao Yin,et al.  Iteratively reweighted algorithms for compressive sensing , 2008, 2008 IEEE International Conference on Acoustics, Speech and Signal Processing.

[60]  Antonio Soriano,et al.  Expectation maximization (EM) algorithms using polar symmetries for computed tomography (CT) image reconstruction , 2013, Comput. Biol. Medicine.

[61]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .