Adaptive noise reduction toward low-dose computed tomography

An efficient noise treatment scheme has been developed to achieve low-dose CT diagnosis based on currently available CT hardware and image reconstruction technologies. The scheme proposed includes two main parts: filtering in sinogram domain and smoothing in image domain. The acquired projection sinograms were first treated by our previously proposed Karhunen-Loeve (K-L) domain penalized weighted least-square (PWLS) filtering, which fully utilizes the prior statistical noise property and three-dimensional (3D) spatial information for an accurate restoration of the low-dose projections. To treat the streak artifacts due to photon starvation, we also incorporated an adaptive filtering into our PWLS framework, which selectively smoothes those channels contributing most to the streak artifacts. After the sinogram filtering, the image was reconstructed by the conventional filtered backprojection (FBP) method. The image is assumed as piecewise regions each has a unique texture. Therefore, an edge-preserving smoothing (EPS) with locally-adaptive parameters to the noise variation was applied for further noise reduction in image domain. Experimental phantom projections acquired by a GE spiral computed tomography (CT) scanner under 10 mAs tube current were used to evaluate the proposed smoothing scheme. The reconstructed imaged demonstrated that the smoothing scheme with appropriate control parameters provides a significant improvement on noise suppression without sacrificing the spatial resolution.

[1]  S. Riederer,et al.  The noise power spectrum in computed X-ray tomography. , 1978, Physics in medicine and biology.

[2]  B. R. Hunt,et al.  Karhunen-Loeve multispectral image restoration, part I: Theory , 1984 .

[3]  Jeffrey A. Fessler Tomographic Reconstruction Using Information-Weighted Spline Smoothing , 1993, IPMI.

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

[5]  J. Polzehl,et al.  Adaptive weights smoothing with applications to image restoration , 1998 .

[6]  Miles N. Wernick,et al.  Fast optimal pre-reconstruction filters for dynamic PET , 1997, 1997 IEEE Nuclear Science Symposium Conference Record.

[7]  J. Marron,et al.  Edge-Preserving Smoothers for Image Processing , 1998 .

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

[9]  M. Wernick,et al.  Fast spatio-temporal image reconstruction for dynamic PET , 1999, IEEE Transactions on Medical Imaging.

[10]  J. B. Lee,et al.  Enhancement of high spectral resolution remote-sensing data by a noise-adjusted principal components transform , 1990 .

[11]  Keishi Kitamura,et al.  Noise reduction in PET attenuation correction using non-linear Gaussian filters , 1999, 1999 IEEE Nuclear Science Symposium. Conference Record. 1999 Nuclear Science Symposium and Medical Imaging Conference (Cat. No.99CH37019).

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

[13]  Wenyuan Xu,et al.  Behavioral analysis of anisotropic diffusion in image processing , 1996, IEEE Trans. Image Process..

[14]  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).

[15]  Maximilian F. Reiser,et al.  Noise/dose reduction and image improvements in screening virtual colonoscopy with tube currents of 20 mAs with nonlinear Gaussian filter chains , 2002, SPIE Medical Imaging.

[16]  Jörg Weule,et al.  Non-Linear Gaussian Filters Performing Edge Preserving Diffusion , 1995, DAGM-Symposium.

[17]  Fred Godtliebsen,et al.  M-Smoother with local Linear Fit , 2002 .

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

[19]  M F Kijewski,et al.  The noise power spectrum of CT images. , 1987, Physics in medicine and biology.

[20]  K Sauer,et al.  Nonstationary filtering of transmission tomograms in high photon counting noise. , 1991, IEEE transactions on medical imaging.

[21]  J. Fessler,et al.  Spatial resolution properties of penalized-likelihood image reconstruction: space-invariant tomographs , 1996, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[22]  Jeffrey A. Fessler,et al.  Nonparametric fixed-interval smoothing with vector splines , 1991, IEEE Trans. Signal Process..

[23]  G. Winkler,et al.  Noise Reduction in Images: Some Recent Edge-Preserving Methods , 1998 .

[24]  Xiaochuan Pan,et al.  Nonparametric regression sinogram smoothing using a roughness-penalized Poisson likelihood objective function , 2000, IEEE Transactions on Medical Imaging.

[25]  Jeffrey A. Fessler,et al.  Three-dimensional non-local edge-preserving regularization for PET transmission reconstruction , 2000, 2000 IEEE Nuclear Science Symposium. Conference Record (Cat. No.00CH37149).

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

[27]  Zhengrong Liang,et al.  Combined transformation of ordering SPECT sinograms for signal extraction from measurements with Poisson noise , 2001, SPIE Medical Imaging.