Feature constrained compressed sensing CT image reconstruction from incomplete data via robust principal component analysis of the database

In computed tomography (CT), incomplete data problems such as limited angle projections often cause artifacts in the reconstruction results. Additional prior knowledge of the image has shown the potential for better results, such as a prior image constrained compressed sensing algorithm. While a pre-full-scan of the same patient is not always available, massive well-reconstructed images of different patients can be easily obtained from clinical multi-slice helical CTs. In this paper, a feature constrained compressed sensing (FCCS) image reconstruction algorithm was proposed to improve the image quality by using the prior knowledge extracted from the clinical database. The database consists of instances which are similar to the target image but not necessarily the same. Robust principal component analysis is employed to retrieve features of the training images to sparsify the target image. The features form a low-dimensional linear space and a constraint on the distance between the image and the space is used. A bi-criterion convex program which combines the feature constraint and total variation constraint is proposed for the reconstruction procedure and a flexible method is adopted for a good solution. Numerical simulations on both the phantom and real clinical patient images were taken to validate our algorithm. Promising results are shown for limited angle problems.

[1]  D. G. Grant Tomosynthesis: a three-dimensional radiographic imaging technique. , 1972, IEEE transactions on bio-medical engineering.

[2]  R. Gerchberg Super-resolution through Error Energy Reduction , 1974 .

[3]  Albert Macovski,et al.  Iterative Reconstruction-Reprojection: An Algorithm for Limited Data Cardiac-Computed Tomography , 1982, IEEE Transactions on Biomedical Engineering.

[4]  M. Turk,et al.  Eigenfaces for Recognition , 1991, Journal of Cognitive Neuroscience.

[5]  Gyula Faigel,et al.  X-Ray Holography , 1999 .

[6]  James T Dobbins,et al.  Digital x-ray tomosynthesis: current state of the art and clinical potential. , 2003, Physics in medicine and biology.

[7]  Alejandro F. Frangi,et al.  Two-dimensional PCA: a new approach to appearance-based face representation and recognition , 2004 .

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

[9]  B. Wilson,et al.  Volume CT with a flat-panel detector on a mobile, isocentric C-arm: pre-clinical investigation in guidance of minimally invasive surgery. , 2005, Medical physics.

[10]  Daoqiang Zhang,et al.  (2D)2PCA: Two-directional two-dimensional PCA for efficient face representation and recognition , 2005, Neurocomputing.

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

[12]  Takeo Kanade,et al.  Robust L/sub 1/ norm factorization in the presence of outliers and missing data by alternative convex programming , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[13]  Emmanuel J. Candès,et al.  Signal recovery from random projections , 2005, IS&T/SPIE Electronic Imaging.

[14]  Chris H. Q. Ding,et al.  R1-PCA: rotational invariant L1-norm principal component analysis for robust subspace factorization , 2006, ICML.

[15]  Li Zhang,et al.  An exact reconstruction algorithm in variable pitch helical cone-beam CT when PI-line exists , 2006 .

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

[17]  Li Zhang,et al.  Volumetric imaging from a multisegment straight-line trajectory and a practical reconstruction algorithm , 2007 .

[18]  E. Roessl,et al.  K-edge imaging in x-ray computed tomography using multi-bin photon counting detectors , 2007, Physics in medicine and biology.

[19]  Rama Chellappa,et al.  Super-Resolution of Face Images Using Kernel PCA-Based Prior , 2007, IEEE Transactions on Multimedia.

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

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

[22]  Nojun Kwak,et al.  Principal Component Analysis Based on L1-Norm Maximization , 2008, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[23]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[24]  Jie Tang,et al.  Low radiation dose C-arm cone-beam CT based on prior image constrained compressed sensing (PICCS): including compensation for image volume mismatch between multiple data acquisitions , 2009, Medical Imaging.

[25]  Xuelong Li,et al.  L1-Norm-Based 2DPCA , 2010, IEEE Transactions on Systems, Man, and Cybernetics, Part B (Cybernetics).

[26]  Hengyong Yu,et al.  Can interior tomography outperform lambda tomography? , 2010, Proceedings of the National Academy of Sciences.

[27]  Hongkai Zhao,et al.  Robust principal component analysis-based four-dimensional computed tomography , 2011, Physics in medicine and biology.

[28]  Yi Ma,et al.  Robust principal component analysis? , 2009, JACM.

[29]  S. Osher,et al.  Multi-energy CT based on a prior rank, intensity and sparsity model (PRISM) , 2011, Inverse problems.

[30]  Christoph Hoeschen,et al.  Investigating tomographic reconstruction with a priori geometrical information. , 2012, Journal of X-ray science and technology.

[31]  Liang Li,et al.  A reweighted total variation minimization method for few view CT reconstruction in the instant CT , 2012, 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference Record (NSS/MIC).

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

[33]  Liang Li,et al.  Anisotropic total variation minimization method for limited-angle CT reconstruction , 2012, Optics & Photonics - Optical Engineering + Applications.

[34]  G. Sapiro,et al.  A collaborative framework for 3D alignment and classification of heterogeneous subvolumes in cryo-electron tomography. , 2013, Journal of structural biology.