DocumentCode
3344307
Title
Does OSEM achieve the lowest variance?
Author
Cloquet, Christophe ; Defrise, Michel
Author_Institution
Dept. of Nucl. Med., Vrije Univ. Brussel, Brussels, Belgium
fYear
2011
fDate
23-29 Oct. 2011
Firstpage
2360
Lastpage
2365
Abstract
Maximum Likelihood (ML) reconstruction algorithms are non biased and achieve the lowest variance, called the Cramer-Rao lower bound (CRLB), for an infinite number of counts and iterations. This result is however not true at finite number of counts or iterations. In this study, we concentrate on the two dimensional Ordered Subsets Expectation Maximization (2D OSEM) algorithm with a finite number of counts and iterations, and investigate the question: given its bias, does this algorithm achieve the minimum variance predicted by the Cramer-Rao lower bound? We found a threshold under which the variance significatively exceeds the biased CRLB. We also found that above this threshold, the variance almost equals the biased CRLB, even for a finite number of iterations and in cold regions. A further analysis is needed to investigate the reason of the observed difference, which might indicate that there exists an algorithm with a smaller variance than OSEM for the same bias, or that a higher bound could be found.
Keywords
expectation-maximisation algorithm; image reconstruction; medical image processing; positron emission tomography; 2D ordered subsets expectation maximization algorithm; Cramer-Rao lower bound; biased CRLB; maximum likelihood reconstruction algorithms; positron emission tomography; Irrigation; Cramer-Rao; MLEM; Maximum likelihood; OSEM; Poisson; finite sample size;
fLanguage
English
Publisher
ieee
Conference_Titel
Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), 2011 IEEE
Conference_Location
Valencia
ISSN
1082-3654
Print_ISBN
978-1-4673-0118-3
Type
conf
DOI
10.1109/NSSMIC.2011.6153880
Filename
6153880
Link To Document