An Expectation Maximization Method for Joint Estimation of Emission Activity Distribution and Photon Attenuation Map in PET

被引:8
作者
Mihlin, Alexander [1 ,2 ]
Levin, Craig S. [1 ,2 ,3 ,4 ]
机构
[1] Stanford Univ, Dept Elect Engn, Stanford, CA 94305 USA
[2] Stanford Univ, Dept Radiol, Stanford, CA 94305 USA
[3] Stanford Univ, Dept Phys, Stanford, CA 94305 USA
[4] Stanford Univ, Dept Bioengn, Stanford, CA 94305 USA
关键词
OF-FLIGHT PET; IMAGE-RECONSTRUCTION; COMPUTED-TOMOGRAPHY; EM ALGORITHM; TRANSMISSION;
D O I
10.1109/TMI.2016.2602339
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A maximum likelihood expectation maximization (MLEM) method is proposed for joint estimation of emission activity distribution and photon attenuation map from positron emission tomography (PET) emission data alone. The method is appealing since: (i) it guarantees monotonic likelihood increase to a local extremum, (ii) does not require arbitrary parameters, and (iii) guarantees the positivity of the estimated distributions. Moreover, we propose a discrete Poisson data acquisition model and numerical algorithm for: (i) efficient graphics processing unit (GPU) based formulation, and (ii) a closed form exact solution for the MLEM update equations, which is essential for accurate and robust estimation. Numerical experiments indicate that in the presence of noise, joint EMAA estimation converges to the true emission activity distribution with root mean square errors of 4% and 0.5% respectively in estimation of lung-and myocardial emission activity distributions for a computational XCAT thorax phantom.
引用
收藏
页码:214 / 224
页数:11
相关论文
共 21 条
[1]  
De Pierro AR, 2006, COMPUT APPL MATH, V25, P205, DOI 10.1590/S1807-03022006000200006
[2]   Time-of-flight PET data determine the attenuation sinogram up to a constant [J].
Defrise, Michel ;
Rezaei, Ahmadreza ;
Nuyts, Johan .
PHYSICS IN MEDICINE AND BIOLOGY, 2012, 57 (04) :885-899
[3]   MAXIMUM LIKELIHOOD FROM INCOMPLETE DATA VIA EM ALGORITHM [J].
DEMPSTER, AP ;
LAIRD, NM ;
RUBIN, DB .
JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-METHODOLOGICAL, 1977, 39 (01) :1-38
[4]   Grouped-coordinate ascent algorithms for penalized-likelihood transmission image reconstruction [J].
Fessler, JA ;
Ficaro, EP ;
Clinthorne, NH ;
Lange, K .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1997, 16 (02) :166-175
[5]  
Frodesen A.G., 1979, Probability and Statistics in Particle Physics
[6]   ACCELERATED IMAGE-RECONSTRUCTION USING ORDERED SUBSETS OF PROJECTION DATA [J].
HUDSON, HM ;
LARKIN, RS .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1994, 13 (04) :601-609
[7]   An EM algorithm for estimating SPECT emission and transmission parameters from emission data only [J].
Krol, A ;
Bowsher, JE ;
Manglos, SH ;
Feiglin, DH ;
Tornai, MP ;
Thomas, ED .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2001, 20 (03) :218-232
[8]  
LANGE K, 1984, J COMPUT ASSIST TOMO, V8, P306
[9]  
Mihlin A., 2013, P IEEE NUCL SCI S ME, P1
[10]   Iterative reconstruction for helical CT: a simulation study [J].
Nuyts, J ;
De Man, B ;
Dupont, P ;
Defrise, M ;
Suetens, P ;
Mortelmans, L .
PHYSICS IN MEDICINE AND BIOLOGY, 1998, 43 (04) :729-737