Cited by 4,429Medical Imaging Techniques and ApplicationsNumerical methods in inverse problemsImage and Signal Denoising Methods

Maximum Likelihood Reconstruction for Emission Tomography

Shepp, L. A., Vardi, Y.·1982·IEEE Transactions on Medical Imaging

Abstract

Previous models for emission tomography (ET) do not distinguish the physics of ET from that of transmission tomography. We give a more accurate general mathematical model for ET where an unknown emission density lambda = lambda(x, y, z) generates, and is to be reconstructed from, the number of counts n(*)(d) in each of D detector units d. Within the model, we give an algorithm for determining an estimate lambdainsertion mark of lambda which maximizes the probability p(n(*)|lambda) of observing the actual detector count data n(*) over all possible densities lambda. Let independent Poisson variables n(b) with unknown means lambda(b), b = 1, ..., B represent the number of unobserved emissions in each of B boxes (pixels) partitioning an object containing an emitter. Suppose each emission in box b is detected in detector unit d with probability p(b, d), d = 1, ..., D with p(b,d) a one-step transition matrix, assumed known. We observe the total number n(*) = n(*)(d) of emissions in each detector unit d and want to estimate the unknown lambda = lambda(b), b = 1, ..., B. For each lambda, the observed data n(*) has probability or likelihood p(n(*)|lambda). The EM algorithm of mathematical statistics starts with an initial estimate lambda(0) and gives the following simple iterative procedure for obtaining a new estimate lambdainsertion mark(new), from an old estimate lambdainsertion mark(old), to obtain lambdainsertion mark(k), k = 1, 2, ..., lambdainsertion mark(new)(b)= lambdainsertion mark(old)(b)Sum of (n(*)p(b,d) from d=1 to D/Sum of lambdainsertion mark()old(b('))p(b('),d) from b(')=1 to B), b=1,...B.

Integrity signals

Registry

Crossref · 10.1109/tmi.1982.4307558

Cites

9 of 9 references verified0 citation errors · 0 not checkable

Statistics

Recomputing p-values needs the paper’s full text.

Claims the paper makes

This paper’s full text isn’t openly available, so its claims and self-check weren’t run. The registry and citation checks still apply.

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Cites

Flagged first, then by citations

Cited by

Sorted by citations

Authors
Shepp, L. A., & Vardi, Y.
Affiliations
AT&T (United States) (US)
Venue
IEEE Transactions on Medical Imaging
Published
1982-10
Volume
1
Issue
2
Pages
113-122
ISSN
  • 0278-0062 (ISSN-L)
  • 1558-254X
OpenAlex
W2069629287
Cited by
4,429
References
9
Citations by year
2026: 74 · 2025: 101 · 2024: 122 · 2023: 101 · 2022: 117 · 2021: 140 · 2020: 154 · 2019: 153 · 2018: 137 · 2017: 142 · 2016: 152 · 2015: 164 · 2014: 167 · 2013: 182 · 2012: 176
FWCI
10.00
Type
journal article
Language
en
Open access
closed
Locations
  • IEEE Transactions on Medical Imaging (publishedVersion)
  • PubMed (publishedVersion)
Topics
  • Health Sciences › Medicine › Radiology, Nuclear Medicine and Imaging › Medical Imaging Techniques and Applications
  • Physical Sciences › Mathematics › Mathematical Physics › Numerical methods in inverse problems
  • Physical Sciences › Computer Science › Computer Vision and Pattern Recognition › Image and Signal Denoising Methods
Keywords
emission tomography, maximum likelihood reconstruction, expectation-maximization algorithm, Poisson likelihood, tomographic reconstruction

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review