Adaptive non-local means filtering based on local noise level for CT denoising

Radiation dose from CT scans is an increasing health concern in the practice of radiology. Higher dose scans can produce clearer images with high diagnostic quality, but may increase the potential risk of radiation-induced cancer or other side effects. Lowering radiation dose alone generally produces a noisier image and may degrade diagnostic performance. Recently, CT dose reduction based on non-local means (NLM) filtering for noise reduction has yielded promising results. However, traditional NLM denoising operates under the assumption that image noise is spatially uniform noise, while in CT images the noise level varies significantly within and across slices. Therefore, applying NLM filtering to CT data using a global filtering strength cannot achieve optimal denoising performance. In this work, we have developed a technique for efficiently estimating the local noise level for CT images, and have modified the NLM algorithm to adapt to local variations in noise level. The local noise level estimation technique matches the true noise distribution determined from multiple repetitive scans of a phantom object very well. The modified NLM algorithm provides more effective denoising of CT data throughout a volume, and may allow significant lowering of radiation dose. Both the noise map calculation and the adaptive NLM filtering can be performed in times that allow integration with the clinical workflow.

[1]  Xiaochuan Pan,et al.  Image reconstruction with shift-variant filtration and its implication for noise and resolution properties in fan-beam computed tomography. , 2003, Medical physics.

[2]  Fang Xu,et al.  Accelerating popular tomographic reconstruction algorithms on commodity PC graphics hardware , 2005, IEEE Transactions on Nuclear Science.

[3]  F. Beekma,et al.  Ordered subset reconstruction for x-ray CT. , 2001, Physics in medicine and biology.

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

[5]  Cynthia M. McCollough,et al.  Projection space denoising with bilateral filtering and CT noise modeling for dose reduction in CT. , 2009, Medical physics.

[6]  Adam Wunderlich,et al.  Image covariance and lesion detectability in direct fan-beam x-ray computed tomography , 2008, Physics in medicine and biology.

[7]  Klaus Mueller,et al.  IOP PUBLISHING PHYSICS IN MEDICINE AND BIOLOGY , 2007 .

[8]  L. Tanoue.  Computed Tomography — An Increasing Source of Radiation Exposure , 2009 .

[9]  Armando Manduca,et al.  Feasibility of dose reduction using novel denoising techniques for low kV (80 kV) CT enterography: optimization and validation. , 2010, Academic radiology.

[10]  C. McCollough,et al.  Optimal tube potential for radiation dose reduction in pediatric CT: principles, clinical implementations, and pitfalls. , 2011, Radiographics : a review publication of the Radiological Society of North America, Inc.

[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]  Rainer Raupach,et al.  Effect of nonlinear three-dimensional optimized reconstruction algorithm filter on image quality and radiation dose: validation on phantoms. , 2008, Medical physics.

[13]  W. Kalender,et al.  Generalized multi-dimensional adaptive filtering for conventional and spiral single-slice, multi-slice, and cone-beam CT. , 2001, Medical physics.

[14]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[15]  J. S. Kole,et al.  Evaluation of accelerated iterative X-ray CT image reconstruction using floating point graphics hardware , 2004 .

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

[17]  Jeffrey A. Fessler,et al.  Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction , 1997, IEEE Transactions on Medical Imaging.

[18]  Xi Chen,et al.  A Robust and Fast Non-Local Means Algorithm for Image Denoising , 2008, Journal of Computer Science and Technology.

[19]  Roberto Manduchi,et al.  Bilateral filtering for gray and color images , 1998, Sixth International Conference on Computer Vision (IEEE Cat. No.98CH36271).

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

[21]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[22]  A. Bruckstein,et al.  K-SVD : An Algorithm for Designing of Overcomplete Dictionaries for Sparse Representation , 2005 .

[23]  K. Stierstorfer,et al.  Weighted FBP--a simple approximate 3D FBP algorithm for multislice spiral CT with good dose usage for arbitrary pitch. , 2004, Physics in medicine and biology.

[24]  Michael Knaup,et al.  Hyperfast parallel-beam and cone-beam backprojection using the cell general purpose hardware. , 2007, Medical physics.

[25]  Jérôme Darbon,et al.  Fast nonlocal filtering applied to electron cryomicroscopy , 2008, 2008 5th IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[26]  Jing Wang,et al.  Nonlinear sinogram smoothing for low-dose X-ray CT , 2004, IEEE Transactions on Nuclear Science.

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

[28]  Yan Zhang,et al.  Investigation of image noise in cone-beam CT imaging due to photon counting statistics with the Feldkamp algorithm by computer simulations , 2008 .

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

[30]  Jeffrey A. Fessler,et al.  Ieee Transactions on Image Processing: to Appear Globally Convergent Algorithms for Maximum a Posteriori Transmission Tomography , 2022 .

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

[32]  M. Defrise,et al.  Iterative reconstruction for helical CT: a simulation study. , 1998, Physics in medicine and biology.

[33]  Lifeng Yu,et al.  Pilot multi-reader study demonstrating potential for dose reduction in dual energy hepatic CT using non-linear blending of mixed kV image datasets , 2011, European Radiology.

[34]  W A Kalender,et al.  Dose reduction in CT by anatomically adapted tube current modulation. I. Simulation studies. , 1999, Medical physics.