System matrix analysis for sparse-view iterative image reconstruction in X-ray CT.

Iterative image reconstruction (IIR) with sparsity-exploiting methods, such as total variation (TV) minimization, used for investigations in compressive sensing (CS) claim potentially large reductions in sampling requirements. Quantifying this claim for computed tomography (CT) is non-trivial, as both the singularity of undersampled reconstruction and the sufficient view number for sparse-view reconstruction are ill-defined. In this paper, the singular value decomposition method is used to study the condition number and singularity of the system matrix and the regularized matrix. An estimation method of the empirical lower bound is proposed, which is helpful for estimating the number of projection views required for exact reconstruction. Simulation studies show that the singularity of the system matrices for different projection views is effectively reduced by regularization. Computing the condition number of a regularized matrix is necessary to provide a reference for evaluating the singularity and recovery potential of reconstruction algorithms using regularization. The empirical lower bound is helpful for estimating the projections view number with a sparse reconstruction algorithm.

[1]  A. Ramm,et al.  The RADON TRANSFORM and LOCAL TOMOGRAPHY , 1996 .

[2]  Aleksandra Pizurica,et al.  Split-Bregman-based sparse-view CT reconstruction , 2011 .

[3]  Xin Jin,et al.  Experimental studies on few-view reconstruction for high-resolution micro-CT. , 2013, Journal of X-ray science and technology.

[4]  E.J. Candes,et al.  An Introduction To Compressive Sampling , 2008, IEEE Signal Processing Magazine.

[5]  H. Tuy AN INVERSION FORMULA FOR CONE-BEAM RECONSTRUCTION* , 1983 .

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

[7]  Melanie Grunwald,et al.  Foundations Of Image Science , 2016 .

[8]  Emil Y. Sidky,et al.  Algorithm-Enabled Low-Dose Micro-CT Imaging , 2011, IEEE Transactions on Medical Imaging.

[9]  Emil Y. Sidky,et al.  Analysis of discrete-to-discrete imaging models for iterative tomographic image reconstruction and compressive sensing , 2011 .

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

[11]  Emil Y. Sidky,et al.  Quantitative study of undersampled recoverability for sparse images in computed tomography , 2012 .

[12]  Vijay K. Madisetti,et al.  The fast discrete Radon transform , 1992, [Proceedings] ICASSP-92: 1992 IEEE International Conference on Acoustics, Speech, and Signal Processing.

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

[14]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

[15]  William H Press,et al.  Discrete Radon transform has an exact, fast inverse and generalizes to operations other than sums along lines , 2006, Proceedings of the National Academy of Sciences.

[16]  Gangrong Qu Singularities of the radon transform of a class of piecewise smooth functions in R2 , 2011 .

[17]  M. Defrise,et al.  An algorithm for total variation regularization in high-dimensional linear problems , 2011 .

[18]  Gregory Beylkin,et al.  Discrete radon transform , 1987, IEEE Trans. Acoust. Speech Signal Process..

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

[20]  P. Cochat,et al.  Et al , 2008, Archives de pediatrie : organe officiel de la Societe francaise de pediatrie.

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

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

[23]  G. Beylkin The inversion problem and applications of the generalized radon transform , 1984 .

[24]  Vijay K. Madisetti,et al.  The fast discrete Radon transform. I. Theory , 1993, IEEE Trans. Image Process..

[25]  D. Sorensen Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[26]  M. Vannier,et al.  Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? , 2009, Inverse problems.

[27]  李磊,et al.  Image reconstruction based on total-variation minimization and alternating direction method in linear scan computed tomography , 2013 .

[28]  Liang Li,et al.  A few-view reweighted sparsity hunting (FRESH) method for CT image reconstruction. , 2013, Journal of X-ray science and technology.

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

[30]  F. Natterer The Mathematics of Computerized Tomography , 1986 .

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