Multiscale noise reduction on low-dose CT sinogram by stationary wavelet transform

Low-dose protocol for computed tomography (CT) scans has been gradually used in clinics to lower the radiation exposure for mass screening. However, the high noise during data acquisition (and therefore degraded image quality) impairs diagnostic accuracy. This work explores a multiscale approach to reduce non-stationary Gaussian noise in low-dose CT sinograms by wavelet analysis. To explore the noise property in wavelet domain, statistical analysis on the distribution of wavelet coefficients was performed with different basic functions, using computer simulation. A stationary wavelet transform was chosen to alleviate the Gibbs ringing effect caused by thresholding process with orthogonal basic functions. A Bayesian analysis was applied to estimate the local variance at each decomposed scale so that noise reduction on the wavelet coefficients becomes adaptive to each scale. Both computer simulations and phantom experiments were performed to show the potential of the presented local-adaptive stationary wavelet transform for low-dose CT. Comparing with traditional smoothing filters and wavelet-based thresholding denoising methods, the proposed method demonstrated noticeable improvement on noise reduction and edge preservation of low-dose CT images.

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

[2]  Edward H. Adelson,et al.  Noise removal via Bayesian wavelet coring , 1996, Proceedings of 3rd IEEE International Conference on Image Processing.

[3]  Kannan Ramchandran,et al.  Low-complexity image denoising based on statistical modeling of wavelet coefficients , 1999, IEEE Signal Processing Letters.

[4]  Stéphane Mallat,et al.  A Theory for Multiresolution Signal Decomposition: The Wavelet Representation , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

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

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

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

[8]  B. Silverman,et al.  The Stationary Wavelet Transform and some Statistical Applications , 1995 .

[9]  D. McCauley,et al.  Low-dose CT of the lungs: preliminary observations. , 1990, Radiology.

[10]  Ruola Ning,et al.  Image denoising based on multiscale singularity detection for cone beam CT breast imaging , 2004, IEEE Transactions on Medical Imaging.

[11]  G. A Theory for Multiresolution Signal Decomposition : The Wavelet Representation , 2004 .

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

[13]  Antonin Chambolle,et al.  Nonlinear wavelet image processing: variational problems, compression, and noise removal through wavelet shrinkage , 1998, IEEE Trans. Image Process..