Maximum likelihood reconstruction for spect with monte carlo modeling: Asymptotic behavior


Journal Article

Inverse Monte Carlo (IMOC) reconstruction for SPECT uses a Maximum Likelihood (EM) estimator with detection probabilities generated as Monte Carlo solutions to the photon transport equation (PTE). To evaluate the behavior of the iterative EM algorithm, reconstructions from experimental projection data for up to 1000 iterations were examined. Compensation for scatter and attenuation was achieved by including those effects in the PTE. For uniform activity distributions, noise increased monotonically with iteration. With line sources included, both noise and resolution improved between 1 and 30 iterations after which resolution was slightly improved at the expense of noise. Contrast (for a nonactive region surrounded by activity) improved from 0.32 at 5 iterations to 0.90 at 200 iterations, and 0.97 at 1000. Uncertainty in the measurement increased due to increased noise in the active region. Copyright © 1987 by The Institute of Electrical and Electronics Engineers, Inc.

Full Text

Duke Authors

Cited Authors

  • Floyd, CE; Jaszczak, RJ; Coleman, RE

Published Date

  • January 1, 1987

Published In

Volume / Issue

  • 34 / 1

Start / End Page

  • 285 - 287

Electronic International Standard Serial Number (EISSN)

  • 1558-1578

International Standard Serial Number (ISSN)

  • 0018-9499

Digital Object Identifier (DOI)

  • 10.1109/TNS.1987.4337349

Citation Source

  • Scopus