Accelerated Image Reconstruction Using Ordered Subsets of Projection Data

H. Malcolm Hudson, Richard S. Larkin

Research output: Contribution to journalArticlepeer-review

2789 Citations (Scopus)


We define ordered subset processing For standard algorithms (such as Expectation Maximization, EM) for image restoration from projections. Ordered subsets methods group projection data into an ordered sequence of subsets (or blocks). An iteration of ordered subsets EM is defined as a single pass through all the subsets, in each subset using the current estimate to initialize application of EM with that data subset. This approach is similar in concept to block-Kaczmarz methods introduced by Eggermont et al. [1] for iterative reconstruction. Simultaneous iterative reconstruction (SIRT) and multiplicative algebraic reconstruction (MART) techniques are well known special cases. Ordered subsets EM (OS-EM) provides a restoration imposing a natural positivity condition and with close links to the EM algorithm. OS-EM is applicable in both single photon (SPECT) and positron emission tomography (PET). In simulation studies in SPECT, the OS-EM algorithm provides an order-of-magnitude acceleration over EM, with restoration quality maintained.

Original languageEnglish
Pages (from-to)601-609
Number of pages9
JournalIEEE Transactions on Medical Imaging
Issue number4
Publication statusPublished - 1994


Dive into the research topics of 'Accelerated Image Reconstruction Using Ordered Subsets of Projection Data'. Together they form a unique fingerprint.

Cite this