Adaptive non‐local means on local principle neighborhood for noise/artifacts reduction in low‐dose CT images

Purpose Low‐dose CT (LDCT) technique can reduce the x‐ray radiation exposure to patients at the cost of degraded images with severe noise and artifacts. Non‐local means (NLM) filtering has shown its potential in improving LDCT image quality. However, currently most NLM‐based approaches employ a weighted average operation directly on all neighbor pixels with a fixed filtering parameter throughout the NLM filtering process, ignoring the non‐stationary noise nature of LDCT images. In this paper, an adaptive NLM filtering scheme on local principle neighborhoods (PC‐NLM) is proposed for structure‐preserving noise/artifacts reduction in LDCT images. Methods Instead of using neighboring patches directly, in the PC‐NLM scheme, the principle component analysis (PCA) is first applied on local neighboring patches of the target patch to decompose the local patches into uncorrelated principle components (PCs), then a NLM filtering is used to regularize each PC of the target patch and finally the regularized components is transformed to get the target patch in image domain. Especially, in the NLM scheme, the filtering parameter is estimated adaptively from local noise level of the neighborhood as well as the signal‐to‐noise ratio (SNR) of the corresponding PC, which guarantees a “weaker” NLM filtering on PCs with higher SNR and a “stronger” filtering on PCs with lower SNR. The PC‐NLM procedure is iteratively performed several times for better removal of the noise and artifacts, and an adaptive iteration strategy is developed to reduce the computational load by determining whether a patch should be processed or not in next round of the PC‐NLM filtering. Results The effectiveness of the presented PC‐NLM algorithm is validated by experimental phantom studies and clinical studies. The results show that it can achieve promising gain over some state‐of‐the‐art methods in terms of artifact suppression and structure preservation. Conclusions With the use of PCA on local neighborhoods to extract principal structural components, as well as adaptive NLM filtering on PCs of the target patch using filtering parameter estimated based on the local noise level and corresponding SNR, the proposed PC‐NLM method shows its efficacy in preserving fine anatomical structures and suppressing noise/artifacts in LDCT images.

[1]  Thomas W. Parks,et al.  Adaptive principal components and image denoising , 2003, Proceedings 2003 International Conference on Image Processing (Cat. No.03CH37429).

[2]  Jeffrey H Siewerdsen,et al.  PIRPLE: a penalized-likelihood framework for incorporation of prior images in CT reconstruction , 2013, Physics in medicine and biology.

[3]  Tolga Tasdizen,et al.  Principal Neighborhood Dictionaries for Nonlocal Means Image Denoising , 2009, IEEE Transactions on Image Processing.

[4]  Xuanqin Mou,et al.  Statistical Interior Tomography , 2010, IEEE Transactions on Medical Imaging.

[5]  Timothy D. Solberg,et al.  Noise reduction in low-dose cone beam CT by incorporating prior volumetric image information , 2013 .

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

[7]  E. Hoffman,et al.  Ultra-low dose lung CT perfusion regularized by a previous scan. , 2009, Academic radiology.

[8]  Wei Xu,et al.  Efficient low-dose CT artifact mitigation using an artifact-matched prior scan. , 2012, Medical physics.

[9]  David J Brenner,et al.  Cancer risks from CT scans: now we have data, what next? , 2012, Radiology.

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

[11]  D. Brenner,et al.  Computed tomography--an increasing source of radiation exposure. , 2007, The New England journal of medicine.

[12]  David Zhang,et al.  Two-stage image denoising by principal component analysis with local pixel grouping , 2010, Pattern Recognit..

[13]  Huazhong Shu,et al.  Artifact Suppressed Dictionary Learning for Low-Dose CT Image Processing , 2014, IEEE Transactions on Medical Imaging.

[14]  Jing Wang,et al.  Deriving adaptive MRF coefficients from previous normal-dose CT scan for low-dose image reconstruction via penalized weighted least-squares minimization. , 2014, Medical physics.

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

[16]  Gabor Jozsef,et al.  Breast radiotherapy in the prone position primarily reduces the maximum out-of-field measured dose to the ipsilateral lung. , 2012, Medical physics.

[17]  L. Xing,et al.  Iterative image reconstruction for CBCT using edge-preserving prior. , 2008, Medical physics.

[18]  Christine Toumoulin,et al.  Thoracic low-dose CT image processing using an artifact suppressed large-scale nonlocal means , 2012, Physics in medicine and biology.

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

[20]  Klaus Mueller,et al.  Database-assisted low-dose CT image restoration. , 2013, Medical physics.

[21]  Yuanke Zhang,et al.  Statistical Sinogram Smoothing for Low-Dose CT With Segmentation-Based Adaptive Filtering , 2010, IEEE Transactions on Nuclear Science.

[22]  O. Linton,et al.  National conference on dose reduction in CT, with an emphasis on pediatric patients. , 2003, AJR. American journal of roentgenology.

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

[24]  Jing Wang,et al.  Multiscale Penalized Weighted Least-Squares Sinogram Restoration for Low-Dose X-Ray Computed Tomography , 2008, IEEE Transactions on Biomedical Engineering.

[25]  Armando Manduca,et al.  Adaptive nonlocal means filtering based on local noise level for CT denoising. , 2013, Medical physics.

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

[27]  Lei Zhang,et al.  Low-Dose X-ray CT Reconstruction via Dictionary Learning , 2012, IEEE Transactions on Medical Imaging.

[28]  Richard M. Everson,et al.  Inferring the eigenvalues of covariance matrices from limited, noisy data , 2000, IEEE Trans. Signal Process..

[29]  Jing Huang,et al.  Iterative Reconstruction for X-Ray Computed Tomography Using Prior-Image Induced Nonlocal Regularization , 2014, IEEE Transactions on Biomedical Engineering.

[30]  Pierrick Coupé,et al.  MRI noise estimation and denoising using non-local PCA , 2015, Medical Image Anal..

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

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

[33]  Rainer Raupach,et al.  Wavelet Based Noise Reduction in CT-Images Using Correlation Analysis , 2008, IEEE Transactions on Medical Imaging.