L0 constrained sparse reconstruction for multi-slice helical CT reconstruction

In this paper, we present a Bayesian maximum a posteriori method for multi-slice helical CT reconstruction based on an L0-norm prior. It makes use of a very low number of projections. A set of surrogate potential functions is used to successively approximate the L0-norm function while generating the prior and to accelerate the convergence speed. Simulation results show that the proposed method provides high quality reconstructions with highly sparse sampled noise-free projections. In the presence of noise, the reconstruction quality is still significantly better than the reconstructions obtained with L1-norm or L2-norm priors.

[1]  Alfred O. Hero,et al.  Ieee Transactions on Image Processing: to Appear Penalized Maximum-likelihood Image Reconstruction Using Space-alternating Generalized Em Algorithms , 2022 .

[2]  Hakan Erdogan,et al.  Monotonic algorithms for transmission tomography , 2002, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[3]  Justin Romberg,et al.  Practical Signal Recovery from Random Projections , 2005 .

[4]  H Rusinek,et al.  Pulmonary nodule detection: low-dose versus conventional CT. , 1998, Radiology.

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

[6]  Ken D. Sauer,et al.  Parallel computation of sequential pixel updates in statistical tomographic reconstruction , 1995, Proceedings., International Conference on Image Processing.

[7]  Jean-Baptiste Thibault,et al.  High Quality Iterative Image Reconstruction For Multi-Slice Helical CT , 2013 .

[8]  Ariela Sofer,et al.  A data-parallel algorithm for iterative tomographic image reconstruction , 1999, Proceedings. Frontiers '99. Seventh Symposium on the Frontiers of Massively Parallel Computation.

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

[10]  P. L. La Riviere.  Penalized-likelihood sinogram smoothing for low-dose CT. , 2005, Medical physics.

[11]  Alfred O. Hero,et al.  Space-alternating generalized expectation-maximization algorithm , 1994, IEEE Trans. Signal Process..

[12]  Matthijs Oudkerk,et al.  Coronary angiography with multi-slice computed tomography , 2001, The Lancet.

[13]  Jeffrey A. Fessler,et al.  A paraboloidal surrogates algorithm for convergent penalized-likelihood emission image reconstruction , 1998, 1998 IEEE Nuclear Science Symposium Conference Record. 1998 IEEE Nuclear Science Symposium and Medical Imaging Conference (Cat. No.98CH36255).

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

[15]  Lee F Rogers,et al.  Dose reduction in CT: how low can we go? , 2002, AJR. American journal of roentgenology.

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

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

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

[19]  M. L. Fernanda de la Cruz Rodríguez.  [Coronary angiography]. , 1983, Revista de enfermeria.

[20]  T. Slovis,et al.  CT and computed radiography: the pictures are great, but is the radiation dose greater than required? , 2002, AJR. American journal of roentgenology.

[21]  Emmanuel J. Candès,et al.  SPARSE SIGNAL AND IMAGE RECOVERY FROM COMPRESSIVE SAMPLES , 2007, 2007 4th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[22]  J. Hsieh,et al.  Nonstationary noise characteristics of the helical scan and its impact on image quality and artifacts. , 1997, Medical physics.

[23]  K S Lee,et al.  Low-Dose, Volumetric Helical CT: Image Quality, Radiation Dose, and Usefulness for Evaluation of Bronchiectasis , 2000, Investigative radiology.

[24]  T Nielsen,et al.  Cardiac cone-beam CT volume reconstruction using ART. , 2005, Medical physics.

[25]  Zhengrong Liang,et al.  Noise properties of low-dose CT projections and noise treatment by scale transformations , 2001, 2001 IEEE Nuclear Science Symposium Conference Record (Cat. No.01CH37310).

[26]  W. Bautz,et al.  Detection of Coronary Artery Stenoses With Thin-Slice Multi-Detector Row Spiral Computed Tomography and Multiplanar Reconstruction , 2003, Circulation.

[27]  Hakan Erdogan,et al.  Ordered subsets algorithms for transmission tomography. , 1999, Physics in medicine and biology.

[28]  D. Donoho.  For most large underdetermined systems of linear equations the minimal 𝓁1‐norm solution is also the sparsest solution , 2006 .

[29]  M. Reiser,et al.  Detection of Coronary Artery Stenoses With Multislice Helical CT Angiography , 2002, Journal of computer assisted tomography.

[30]  Ping Xia,et al.  Low-dose megavoltage cone-beam CT for radiation therapy. , 2005, International journal of radiation oncology, biology, physics.

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

[32]  Z. Liang,et al.  Noise reduction for low-dose single-slice Helical CT sinogram , 2006, IEEE Symposium Conference Record Nuclear Science 2004..

[33]  Patrick J. La Riviere.  Penalized‐likelihood sinogram smoothing for low‐dose CT , 2005 .

[34]  S. Manglos,et al.  Transmission maximum-likelihood reconstruction with ordered subsets for cone beam CT. , 1995, Physics in medicine and biology.

[35]  Xiaolin Wu.  Adaptive binary vector quantization using Hamming codes , 1995, Proceedings., International Conference on Image Processing.

[36]  Xiaochuan Pan,et al.  Image reconstruction on PI-lines by use of filtered backprojection in helical cone-beam CT. , 2004, Physics in medicine and biology.

[37]  D. McCauley,et al.  Low-dose CT of the lungs: preliminary observations. , 1990, Radiology.

[38]  B. De Man,et al.  A comparison between Filtered Backprojection, Post-Smoothed Weighted Least Squares, and Penalized Weighted Least Squares for CT reconstruction , 2006, 2006 IEEE Nuclear Science Symposium Conference Record.

[39]  Ken D. Sauer,et al.  Parallelizable Bayesian tomography algorithms with rapid, guaranteed convergence , 2000, IEEE Trans. Image Process..