An adaptive regularization method for low-dose CT reconstruction from CT transmission data in Poisson–Gaussian noise

Abstract Much research has been devoted to the problem of low-dose computerized tomography (CT) image reconstruction with a scan protocol by both lowering the X-ray tube current (low-mAs) and reducing the total number of projection views(sparse view). However, the CT transmission data may be severely corrupted by X-ray quanta noise and system electronic noise. Recently, a nonlocal means (NLM) regularization method for low-dose CT reconstruction was proposed. Although this method is effective in suppressing both Poisson–Gaussian noise and artifacts, it has two disadvantages: the heavy computational burden and blurred edge information in CT reconstruction. This paper proposes an adaptive patch-wise regularization method for low-dose CT reconstruction from available CT transmission data in Poisson–Gaussian noise. The proposed cost function includes a penalized weighted least-square term and two patch-wise regularization terms which are combined with a novel adaptive regularization parameter. The two regularization terms take advantage of image redundant information across different scales. By exploiting both global and local structure information, the reconstruction accuracy is enhanced. In addition, we select the adaptive regularization parameter based on the texture of the image patch so that the edge information is further maintained. Furthermore, our algorithm updates a patch rather than a pixel, the computation burden is greatly reduced. Finally, experiment results show that the proposed adaptive patch-wise regularization method for low-dose CT reconstruction is superior to several conventional regularization-based approaches in terms of computation efficiency and resolution quality.

[1]  Michael Elad,et al.  Multi-Scale Patch-Based Image Restoration , 2016, IEEE Transactions on Image Processing.

[2]  Zhengrong Liang,et al.  Analytical noise treatment for low-dose CT projection data by penalized weighted least-square smoothing in the K-L domain , 2002, SPIE Medical Imaging.

[3]  Patrick J. La Rivière,et al.  Reduction of noise-induced streak artifacts in X-ray computed tomography through spline-based penalized-likelihood sinogram smoothing , 2005, IEEE Transactions on Medical Imaging.

[4]  Yi Liu,et al.  The statistical sinogram smoothing via adaptive-weighted total variation regularization for low-dose X-ray CT , 2014 .

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

[6]  Linghong Zhou,et al.  Iterative image reconstruction using modified non-local means filtering for limited-angle computed tomography. , 2016, Physica medica : PM : an international journal devoted to the applications of physics to medicine and biology : official journal of the Italian Association of Biomedical Physics.

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

[8]  Hong Shangguan,et al.  The adaptive sinogram restoration algorithm based on anisotropic diffusion by energy minimization for low-dose X-ray CT , 2014 .

[9]  Jiang Hsieh,et al.  Computed Tomography: Principles, Design, Artifacts, and Recent Advances, Fourth Edition , 2022 .

[10]  Dai-Qiang Chen,et al.  Spatially adapted regularization parameter selection based on the local discrepancy function for Poissonian image deblurring , 2011 .

[11]  Jing Wang,et al.  Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography , 2006, IEEE Transactions on Medical Imaging.

[12]  Jong Chul Ye,et al.  Sparse-View Spectral CT Reconstruction Using Spectral Patch-Based Low-Rank Penalty , 2015, IEEE Transactions on Medical Imaging.

[13]  Steve B. Jiang,et al.  Low-dose 4DCT reconstruction via temporal nonlocal means. , 2010, Medical physics.

[14]  Jing Huang,et al.  Penalized weighted least-squares approach for multienergy computed tomography image reconstruction via structure tensor total variation regularization , 2016, Comput. Medical Imaging Graph..

[15]  Aleksandra Pizurica,et al.  Iterative CT Reconstruction Using Shearlet-Based Regularization , 2013, IEEE Transactions on Nuclear Science.

[16]  Zhengrong Liang,et al.  Assessment of prior image induced nonlocal means regularization for low‐dose CT reconstruction: Change in anatomy , 2017, Medical physics.

[17]  R. Siddon Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[18]  Zhengrong Liang,et al.  Iterative image reconstruction for cerebral perfusion CT using a pre-contrast scan induced edge-preserving prior , 2012, Physics in medicine and biology.

[19]  Jing Huang,et al.  Spectral CT Image Restoration via an Average Image-Induced Nonlocal Means Filter , 2016, IEEE Transactions on Biomedical Engineering.

[20]  Jing Huang,et al.  Sparse angular CT reconstruction using non-local means based iterative-correction POCS , 2011, Comput. Biol. Medicine.

[21]  Jianhua Ma,et al.  Statistical image reconstruction for low-dose CT using nonlocal means-based regularization. Part II: An adaptive approach , 2015, Comput. Medical Imaging Graph..

[22]  Andrés Almansa,et al.  A TV Based Restoration Model with Local Constraints , 2008, J. Sci. Comput..

[23]  Jeffrey A. Fessler,et al.  Efficient and accurate likelihood for iterative image reconstruction in x-ray computed tomography , 2003, SPIE Medical Imaging.

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

[25]  Bruce R Whiting,et al.  Statistical reconstruction for x-ray computed tomography using energy-integrating detectors , 2007, Physics in medicine and biology.

[26]  Raanan Fattal,et al.  Image and video upscaling from local self-examples , 2011, TOGS.

[27]  Junyan Rong,et al.  Sparse-view reconstruction from restored low-dose CT projections , 2013, 2013 IEEE Nuclear Science Symposium and Medical Imaging Conference (2013 NSS/MIC).

[28]  Zhaoying Bian,et al.  A Simple Low-Dose X-Ray CT Simulation From High-Dose Scan , 2015, IEEE Transactions on Nuclear Science.

[29]  Gang Liu,et al.  High-order TVL1-based images restoration and spatially adapted regularization parameter selection , 2014, Comput. Math. Appl..

[30]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[31]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[32]  Michael Hintermüller,et al.  Expected absolute value estimators for a spatially adapted regularization parameter choice rule in L1-TV-based image restoration , 2010 .

[33]  Wufan Chen,et al.  Variance analysis of x-ray CT sinograms in the presence of electronic noise background. , 2012, Medical physics.

[34]  Jing Wang,et al.  Statistical image reconstruction for low-dose CT using nonlocal means-based regularization , 2014, Comput. Medical Imaging Graph..

[35]  Jian Yu,et al.  A Dictionary Learning Approach for Poisson Image Deblurring , 2013, IEEE Transactions on Medical Imaging.

[36]  Jianhua Ma,et al.  Iterative reconstruction for dual energy CT with an average image-induced nonlocal means regularization , 2017, Physics in medicine and biology.

[37]  J. O’Sullivan,et al.  Properties of preprocessed sinogram data in x-ray computed tomography. , 2006, Medical physics.

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

[39]  Jingyan Xu,et al.  Electronic Noise Modeling in Statistical Iterative Reconstruction , 2009, IEEE Transactions on Image Processing.

[40]  Eric L. Miller,et al.  Tensor-Based Formulation and Nuclear Norm Regularization for Multienergy Computed Tomography , 2013, IEEE Transactions on Image Processing.

[41]  Jean-Michel Morel,et al.  A non-local algorithm for image denoising , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[42]  Lee-Tzuu Chang,et al.  A Method for Attenuation Correction in Radionuclide Computed Tomography , 1978, IEEE Transactions on Nuclear Science.

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

[44]  Qianjin Feng,et al.  Low-dose computed tomography image restoration using previous normal-dose scan. , 2011, Medical physics.

[45]  Jean-Michel Morel,et al.  A Review of Image Denoising Algorithms, with a New One , 2005, Multiscale Model. Simul..

[46]  Patrick J. La Rivière,et al.  Penalized-likelihood sinogram restoration for computed tomography , 2006, IEEE Transactions on Medical Imaging.

[47]  Emma E. Regentova,et al.  Sparse-View CT Reconstruction Using Curvelet and TV-Based Regularization , 2016, ICIAR.

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

[49]  J. Hsieh Adaptive streak artifact reduction in computed tomography resulting from excessive x-ray photon noise. , 1998, Medical physics.

[50]  Dae-Hong Kim,et al.  Feasibility of sinogram reconstruction based on inpainting method with decomposed sinusoid-like curve (S-curve) using total variation (TV) noise reduction algorithm in computed tomography (CT) imaging system: A simulation study , 2018 .

[51]  Robert M. Haralick,et al.  Textural Features for Image Classification , 1973, IEEE Trans. Syst. Man Cybern..

[52]  C. E. Honeycutt,et al.  Image analysis techniques and gray-level co-occurrence matrices (GLCM) for calculating bioturbation indices and characterizing biogenic sedimentary structures , 2008, Comput. Geosci..

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

[54]  Sudeep D. Thepade,et al.  Image Retrieval using Texture Features extracted from GLCM, LBG and KPE , 2010 .

[55]  Michal Irani,et al.  Super-resolution from a single image , 2009, 2009 IEEE 12th International Conference on Computer Vision.

[56]  Jeffrey A. Fessler Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[57]  Yiqiu Dong,et al.  A Multi-Scale Vectorial Lτ-TV Framework for Color Image Restoration , 2011, International Journal of Computer Vision.