Penalized-likelihood sinogram restoration for computed tomography

被引:170
作者
La Riviere, Patrick J. [1 ]
Bian, Junguo [1 ]
Vargas, Phillip A. [1 ]
机构
[1] Univ Chicago, Dept Radiol, Chicago, IL 60637 USA
关键词
computed tomography; penalized likelihood; sinogram restoration;
D O I
10.1109/TMI.2006.875429
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
We formulate computed tomography (CT) sinogram preprocessing as a statistical restoration problem in which the goal is to obtain the best estimate of the line integrals needed for reconstruction from the set of noisy, degraded measurements. CT measurement data are degraded by a number of factors-including beam hardening and off-focal radiation-that produce artifacts in reconstructed images unless properly corrected. Currently, such effects are addressed by a sequence of sinogram-preprocessing steps, including deconvolution corrections for off-focal radiation, that have the potential to amplify noise. Noise itself is generally mitigated through apodization of the reconstruction kernel, which effectively ignores the measurement statistics, Although in high-noise situations adaptive filtering methods that loosely model data statistics are sometimes applied. As an alternative, we present a general imaging model relating the degraded measurements to the sinogram of ideal line integrals and propose to estimate these line integrals by iteratively optimizing a statistically based objective function. We consider three different strategies for estimating the set of ideal line integrals, one based on direct estimation of ideal "monochromatic" line integrals that have been corrected for single-material beam hardening, one based on estimation of ideal "polychromatic" line integrals that can be readily mapped to monochromatic line integrals, and one based on estimation of ideal transmitted intensities, from which ideal, monochromatic line integrals can be readily estimated. The first two approaches involve maximization of a penalized Poisson-likelihood objective function while the third involves minimization of a quadratic penalized weighted least squares (PWLS) objective applied in the transmitted intensity domain. We find that at low exposure levels typical of those being considered for screening CT, the Poisson-likelihood based approaches outperform the PWLS objective as well as a standard approach based on adaptive filtering followed by deconvolution. At higher exposure levels, the approaches all perform similarly.
引用
收藏
页码:1022 / 1036
页数:15
相关论文
共 35 条
[1]  
[Anonymous], 2000, COMPUTED TOMOGRAPHY
[2]   Potential of software-based scatter corrections in cone-beam volume CT [J].
Bertram, M ;
Wiegert, J ;
Rose, G .
MEDICAL IMAGING 2005: PHYSICS OF MEDICAL IMAGING, PTS 1 AND 2, 2005, 5745 :259-270
[3]   BEAM HARDENING IN X-RAY RECONSTRUCTIVE TOMOGRAPHY [J].
BROOKS, RA ;
DICHIRO, G .
PHYSICS IN MEDICINE AND BIOLOGY, 1976, 21 (03) :390-398
[4]  
CARMI R, 2004, NUCL SCI S C REC, V5, P2766
[5]   A MODIFIED EXPECTATION MAXIMIZATION ALGORITHM FOR PENALIZED LIKELIHOOD ESTIMATION IN EMISSION TOMOGRAPHY [J].
DEPIERRO, AR .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1995, 14 (01) :132-137
[6]   Segmentation-free statistical image reconstruction for polyenergetic x-ray computed tomography with experimental validation [J].
Elbakri, IA ;
Fessler, JA .
PHYSICS IN MEDICINE AND BIOLOGY, 2003, 48 (15) :2453-2477
[7]  
ELBAKRI IA, 2003, THESIS U MICHIGAN AN
[8]   Efficient and accurate likelihood for iterative image reconstruction in x-ray computed tomography [J].
Elbakri, LA ;
Fessler, JA .
MEDICAL IMAGING 2003: IMAGE PROCESSING, PTS 1-3, 2003, 5032 :1839-1850
[9]   Ordered subsets algorithms for transmission tomography [J].
Erdogan, H ;
Fessler, JA .
PHYSICS IN MEDICINE AND BIOLOGY, 1999, 44 (11) :2835-2851
[10]   Monotonic algorithms for transmission tomography [J].
Erdogan, H ;
Fessler, JA .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1999, 18 (09) :801-814