SNR-weighted sinogram smoothing with improved noise-resolution properties for low-dose x-ray computed tomography

To treat the noise in low-dose x-ray CT projection data more accurately, analysis of the noise properties of the data and development of a corresponding efficient noise treatment method are two major problems to be addressed. In order to obtain an accurate and realistic model to describe the x-ray CT system, we acquired thousands of repeated measurements on different phantoms at several fixed scan angles by a GE high-speed multi-slice spiral CT scanner. The collected data were calibrated and log-transformed by the sophisticated system software, which converts the detected photon energy into sinogram data that satisfies the Radon transform. From the analysis of these experimental data, a nonlinear relation between mean and variance for each datum of the sinogram was obtained. In this paper, we integrated this nonlinear relation into a penalized likelihood statistical framework for a SNR (signal-to-noise ratio) adaptive smoothing of noise in the sinogram. After the proposed preprocessing, the sinograms were reconstructed with unapodized FBP (filtered backprojection) method. The resulted images were evaluated quantitatively, in terms of noise uniformity and noise-resolution tradeoff, with comparison to other noise smoothing methods such as Hanning filter and Butterworth filter at different cutoff frequencies. Significant improvement on noise and resolution tradeoff and noise property was demonstrated.

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

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

[3]  Bruce R. Whiting,et al.  Signal statistics in x-ray computed tomography , 2002, SPIE Medical Imaging.

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

[5]  John G. Proakis,et al.  Probability, random variables and stochastic processes , 1985, IEEE Trans. Acoust. Speech Signal Process..

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

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

[8]  Tianhu Lei,et al.  Statistical approach to X-ray CT imaging and its applications in image analysis. I. Statistical analysis of X-ray CT imaging , 1992, IEEE Trans. Medical Imaging.

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

[10]  Zhengrong Liang,et al.  Detector response restoration in image reconstruction of high resolution positron emission tomography , 1994, IEEE Trans. Medical Imaging.

[11]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

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

[13]  D S Lalush,et al.  A generalized Gibbs prior for maximum a posteriori reconstruction in SPECT. , 1993, Physics in medicine and biology.

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

[15]  P. J. La Riviere Reduction of noise-induced streak artifacts in X-ray computed tomography through penalized-likelihood sinogram smoothing , 2003 .