Massively parallel computers for 3D single-photon-emission computed tomography.
Since the introduction of the expectation-maximization (EM) algorithm for generating maximum-likelihood (ML) and maximum a posteriori (MAP) estimates in emission tomography, there have been many investigators applying the ML method. However, almost all of the previous work has been restricted to two-dimensional (2D) reconstructions. The major focus and contribution of this paper is to demonstrate a fully three-dimensional (3D) implementation of the MAP method for single-photon-emission computed tomography (SPECT). The 3D reconstruction exhibits an improvement in resolution when compared to the generation of the series of separate 2D slice reconstructions. As has been noted, the iterative EM algorithm for 2D reconstruction is highly computational; the 3D algorithm is far worse. To accommodate the computational complexity, we have extended our previous work in the 2D arena and demonstrate an implementation on the class of massively parallel processors of the 3D algorithm. Using a 16000 processor MasPar machine, the algorithm is demonstrated to execute at 1.24 s/EM iteration for the entire 64 x 64 x 64 cube of 64 planar measurements obtained from the Siemens Orbiter rotating camera operating in the high-resolution mode.
Butler
C S
CS
Electronic Systems and Signals Research Laboratory, Department of Electrical Engineering, Washington University, Campus Box 1161, St Louis, MO 63110, USA.
Miller
M I
MI
Miller
T R
TR
Wallis
J W
JW
D000465:Algorithms; D001921:Brain / Q000473:pathology; D003201:Computers; D006801:Humans; D007091:Image Processing, Computer-Assisted; D021621:Imaging, Three-Dimensional / Q000295:instrumentation / Q000379:methods*; D016013:Likelihood Functions; D015233:Models, Statistical; D009010:Monte Carlo Method; D012984:Software; D015899:Tomography, Emission-Computed, Single-Photon / Q000295:instrumentation* / Q000379:methods*