Noise properties of low-dose CT projections and noise treatment by scale transformations

Projection data acquired for image reconstruction of low-dose computed tomography (CT) are degraded by many factors. These factors complicate noise analysis on the projection data and render a very challenging task for noise reduction. In this study, we first investigate the noise property of the projection data by analyzing a repeatedly acquired experimental phantom data set, in which the phantom was scanned 900 times at a fixed projection angle. The statistical analysis shows that the noise can be regarded as normally distributed with a nonlinear signal-dependent variance. Based on this observation, we then utilize scale transformations to modulate the projection data so that the data variance can be stabilized to be signal independent. By analyzing the relationship between the data standard deviation and the data mean level, we propose a segmented logarithmic transform for the stabilization of the non-stationary noise. After the scale transformations, the noise variance becomes approximately a constant. A two-dimensional Wiener filter is then designed for an analytical treatment of the noise. Experimental results show that the proposed method has a better noise reduction performance without circular artifacts, by visual judgment, as compared to conventional filters, such as the Harming filter.

[1]  Bartlett Ms.  The use of transformations. , 1947 .

[2]  I Buvat,et al.  Two-dimensional statistical model for regularized backprojection in SPECT. , 1998, Physics in medicine and biology.

[3]  J. Curtiss.  On Transformations Used in the Analysis of Variance , 1943 .

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

[5]  Patrick Dupont,et al.  Maximum-likelihood expectation-maximization reconstruction of sinograms with arbitrary noise distribution using NEC-transformations , 2001, IEEE Transactions on Medical Imaging.

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

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

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

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

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

[11]  M. Bartlett,et al.  The use of transformations. , 1947, Biometrics.

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

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