Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT

David S. Lalush, Seemeen S. Karimi, Benjamin Tsui

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

We evaluate fast reconstruction algorithms including ordered subsets- EM (OS-EM) and Rescaled Block Iterative EM (RBI-EM) in fully 3D SPECT applications on the basis of their convergence and resolution recovery properties as iterations proceed. Using a 3D computer-simulated phantom consisting of 3D Gaussian objects, we simulated projection data that includes only the effects of sampling and detector response of a parallel-hole collimator. Reconstructions were performed using each of the three algorithms (ML-EM, OS-EM, and RBI-EM) modeling the 3D detector response in the projection function. Resolution recovery was evaluated by fitting Gaussians to each of the four objects in the iterated image estimates at selected intervals. Results show that OS-EM and RBI-EM behave identically in this case; their resolution recovery results are virtually indistinguishable. Their resolution behavior appears to be very similar to that of ML-EM, but accelerated by a factor of twenty. For all three algorithms, smaller objects take more iterations to converge. Next, we consider the effect noise has on convergence. For both noise-free and noisy data, we evaluate the log likelihood function at each subiteration of OS-EM and RBI-EM, and at each iteration of ML-EM. With noisy data, both OS-EM and RBI-EM give results for which the log-likelihood function oscillates. Especially for 180-degree acquisitions, RBI-EM oscillates less than OS-EM. Both OS-EM and RBI-EM appear to converge to solutions, but not to the ML solution. We conclude that both OS-EM and RBI-EM can be effective algorithms for fully 3D SPECT reconstruction. Both recover resolution similarly to ML-EM, only more quickly.

Original languageEnglish (US)
Title of host publicationIEEE Nuclear Science Symposium & Medical Imaging Conference
PublisherIEEE
Pages1618-1622
Number of pages5
Volume3
StatePublished - 1996
Externally publishedYes
EventProceedings of the 1996 IEEE Nuclear Science Symposium. Part 1 (of 3) - Anaheim, CA, USA
Duration: Nov 2 1996Nov 9 1996

Other

OtherProceedings of the 1996 IEEE Nuclear Science Symposium. Part 1 (of 3)
CityAnaheim, CA, USA
Period11/2/9611/9/96

Fingerprint

Set theory
Detectors
Recovery
Sampling

ASJC Scopus subject areas

  • Engineering(all)

Cite this

Lalush, D. S., Karimi, S. S., & Tsui, B. (1996). Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT. In IEEE Nuclear Science Symposium & Medical Imaging Conference (Vol. 3, pp. 1618-1622). IEEE.

Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT. / Lalush, David S.; Karimi, Seemeen S.; Tsui, Benjamin.

IEEE Nuclear Science Symposium & Medical Imaging Conference. Vol. 3 IEEE, 1996. p. 1618-1622.

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Lalush, DS, Karimi, SS & Tsui, B 1996, Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT. in IEEE Nuclear Science Symposium & Medical Imaging Conference. vol. 3, IEEE, pp. 1618-1622, Proceedings of the 1996 IEEE Nuclear Science Symposium. Part 1 (of 3), Anaheim, CA, USA, 11/2/96.
Lalush DS, Karimi SS, Tsui B. Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT. In IEEE Nuclear Science Symposium & Medical Imaging Conference. Vol. 3. IEEE. 1996. p. 1618-1622
Lalush, David S. ; Karimi, Seemeen S. ; Tsui, Benjamin. / Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT. IEEE Nuclear Science Symposium & Medical Imaging Conference. Vol. 3 IEEE, 1996. pp. 1618-1622
@inproceedings{97fac5d05aeb4e71a57c029fadfa9e92,
title = "Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT",
abstract = "We evaluate fast reconstruction algorithms including ordered subsets- EM (OS-EM) and Rescaled Block Iterative EM (RBI-EM) in fully 3D SPECT applications on the basis of their convergence and resolution recovery properties as iterations proceed. Using a 3D computer-simulated phantom consisting of 3D Gaussian objects, we simulated projection data that includes only the effects of sampling and detector response of a parallel-hole collimator. Reconstructions were performed using each of the three algorithms (ML-EM, OS-EM, and RBI-EM) modeling the 3D detector response in the projection function. Resolution recovery was evaluated by fitting Gaussians to each of the four objects in the iterated image estimates at selected intervals. Results show that OS-EM and RBI-EM behave identically in this case; their resolution recovery results are virtually indistinguishable. Their resolution behavior appears to be very similar to that of ML-EM, but accelerated by a factor of twenty. For all three algorithms, smaller objects take more iterations to converge. Next, we consider the effect noise has on convergence. For both noise-free and noisy data, we evaluate the log likelihood function at each subiteration of OS-EM and RBI-EM, and at each iteration of ML-EM. With noisy data, both OS-EM and RBI-EM give results for which the log-likelihood function oscillates. Especially for 180-degree acquisitions, RBI-EM oscillates less than OS-EM. Both OS-EM and RBI-EM appear to converge to solutions, but not to the ML solution. We conclude that both OS-EM and RBI-EM can be effective algorithms for fully 3D SPECT reconstruction. Both recover resolution similarly to ML-EM, only more quickly.",
author = "Lalush, {David S.} and Karimi, {Seemeen S.} and Benjamin Tsui",
year = "1996",
language = "English (US)",
volume = "3",
pages = "1618--1622",
booktitle = "IEEE Nuclear Science Symposium & Medical Imaging Conference",
publisher = "IEEE",

}

TY - GEN

T1 - Convergence and resolution recovery of block-iterative EM algorithms modeling 3D detector response in SPECT

AU - Lalush, David S.

AU - Karimi, Seemeen S.

AU - Tsui, Benjamin

PY - 1996

Y1 - 1996

N2 - We evaluate fast reconstruction algorithms including ordered subsets- EM (OS-EM) and Rescaled Block Iterative EM (RBI-EM) in fully 3D SPECT applications on the basis of their convergence and resolution recovery properties as iterations proceed. Using a 3D computer-simulated phantom consisting of 3D Gaussian objects, we simulated projection data that includes only the effects of sampling and detector response of a parallel-hole collimator. Reconstructions were performed using each of the three algorithms (ML-EM, OS-EM, and RBI-EM) modeling the 3D detector response in the projection function. Resolution recovery was evaluated by fitting Gaussians to each of the four objects in the iterated image estimates at selected intervals. Results show that OS-EM and RBI-EM behave identically in this case; their resolution recovery results are virtually indistinguishable. Their resolution behavior appears to be very similar to that of ML-EM, but accelerated by a factor of twenty. For all three algorithms, smaller objects take more iterations to converge. Next, we consider the effect noise has on convergence. For both noise-free and noisy data, we evaluate the log likelihood function at each subiteration of OS-EM and RBI-EM, and at each iteration of ML-EM. With noisy data, both OS-EM and RBI-EM give results for which the log-likelihood function oscillates. Especially for 180-degree acquisitions, RBI-EM oscillates less than OS-EM. Both OS-EM and RBI-EM appear to converge to solutions, but not to the ML solution. We conclude that both OS-EM and RBI-EM can be effective algorithms for fully 3D SPECT reconstruction. Both recover resolution similarly to ML-EM, only more quickly.

AB - We evaluate fast reconstruction algorithms including ordered subsets- EM (OS-EM) and Rescaled Block Iterative EM (RBI-EM) in fully 3D SPECT applications on the basis of their convergence and resolution recovery properties as iterations proceed. Using a 3D computer-simulated phantom consisting of 3D Gaussian objects, we simulated projection data that includes only the effects of sampling and detector response of a parallel-hole collimator. Reconstructions were performed using each of the three algorithms (ML-EM, OS-EM, and RBI-EM) modeling the 3D detector response in the projection function. Resolution recovery was evaluated by fitting Gaussians to each of the four objects in the iterated image estimates at selected intervals. Results show that OS-EM and RBI-EM behave identically in this case; their resolution recovery results are virtually indistinguishable. Their resolution behavior appears to be very similar to that of ML-EM, but accelerated by a factor of twenty. For all three algorithms, smaller objects take more iterations to converge. Next, we consider the effect noise has on convergence. For both noise-free and noisy data, we evaluate the log likelihood function at each subiteration of OS-EM and RBI-EM, and at each iteration of ML-EM. With noisy data, both OS-EM and RBI-EM give results for which the log-likelihood function oscillates. Especially for 180-degree acquisitions, RBI-EM oscillates less than OS-EM. Both OS-EM and RBI-EM appear to converge to solutions, but not to the ML solution. We conclude that both OS-EM and RBI-EM can be effective algorithms for fully 3D SPECT reconstruction. Both recover resolution similarly to ML-EM, only more quickly.

UR - http://www.scopus.com/inward/record.url?scp=0030348705&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0030348705&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:0030348705

VL - 3

SP - 1618

EP - 1622

BT - IEEE Nuclear Science Symposium & Medical Imaging Conference

PB - IEEE

ER -