Projection space denoising with bilateral filtering and CT noise modeling for dose reduction in CT.

PURPOSE To investigate a novel locally adaptive projection space denoising algorithm for low-dose CT data. METHODS The denoising algorithm is based on bilateral filtering, which smooths values using a weighted average in a local neighborhood, with weights determined according to both spatial proximity and intensity similarity between the center pixel and the neighboring pixels. This filtering is locally adaptive and can preserve important edge information in the sinogram, thus maintaining high spatial resolution. A CT noise model that takes into account the bowtie filter and patient-specific automatic exposure control effects is also incorporated into the denoising process. The authors evaluated the noise-resolution properties of bilateral filtering incorporating such a CT noise model in phantom studies and preliminary patient studies with contrast-enhanced abdominal CT exams. RESULTS On a thin wire phantom, the noise-resolution properties were significantly improved with the denoising algorithm compared to commercial reconstruction kernels. The noise-resolution properties on low-dose (40 mA s) data after denoising approximated those of conventional reconstructions at twice the dose level. A separate contrast plate phantom showed improved depiction of low-contrast plates with the denoising algorithm over conventional reconstructions when noise levels were matched. Similar improvement in noise-resolution properties was found on CT colonography data and on five abdominal low-energy (80 kV) CT exams. In each abdominal case, a board-certified subspecialized radiologist rated the denoised 80 kV images markedly superior in image quality compared to the commercially available reconstructions, and denoising improved the image quality to the point where the 80 kV images alone were considered to be of diagnostic quality. CONCLUSIONS The results demonstrate that bilateral filtering incorporating a CT noise model can achieve a significantly better noise-resolution trade-off than a series of commercial reconstruction kernels. This improvement in noise-resolution properties can be used for improving image quality in CT and can be translated into substantial dose reduction.

[1]  F. J. Anscombe,et al.  THE TRANSFORMATION OF POISSON, BINOMIAL AND NEGATIVE-BINOMIAL DATA , 1948 .

[2]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

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

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

[5]  W A Kalender,et al.  Dose reduction in CT by anatomically adapted tube current modulation. II. Phantom measurements. , 1999, Medical physics.

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

[7]  J M Boone,et al.  Determination of the presampled MTF in computed tomography. , 2001, Medical physics.

[8]  Omer Demirkaya.  Reduction of noise and image artifacts in computed tomography by nonlinear filtration of projection images , 2001, SPIE Medical Imaging.

[9]  Michael Elad,et al.  On the origin of the bilateral filter and ways to improve it , 2002, IEEE Trans. Image Process..

[10]  Rainer Raupach,et al.  Spatial domain filtering for fast modification of the tradeoff between image sharpness and pixel noise in computed tomography , 2003, IEEE Transactions on Medical Imaging.

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

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

[13]  Lucas J. van Vliet,et al.  Separable bilateral filtering for fast video preprocessing , 2005, 2005 IEEE International Conference on Multimedia and Expo.

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

[15]  M. Reiser,et al.  Material differentiation by dual energy CT: initial experience , 2007, European Radiology.

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

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

[18]  Borut Marincek,et al.  Dual-Energy Contrast-Enhanced Computed Tomography for the Detection of Urinary Stone Disease , 2007, Investigative radiology.

[19]  Walter Huda,et al.  Patient radiation doses from adult and pediatric CT. , 2007, AJR. American journal of roentgenology.

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

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

[22]  James C. Williams,et al.  Noninvasive differentiation of uric acid versus non-uric acid kidney stones using dual-energy CT. , 2007, Academic radiology.

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

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