Maximum likelihood SPECT in clinical computation times using mesh-connected parallel computers.

Extending the work of A.W. McCarthy et al. (1988) and M.I. Miller and B. Roysam (1991), the authors demonstrate that a fully parallel implementation of the maximum-likelihood method for single-photon emission computed tomography (SPECT) can be accomplished in clinical time frames on massively parallel systolic array processors. The authors show that for SPECT imaging on 64x64 image grids, with 96 view angles, the single-instruction, multiple data (SIMD) distributed array processor containing 64(2) processors performs the expectation-maximization (EM) algorithm with Good's smoothing at a rate of 1 iteration/1.5 s. This promises for emission tomography fully Bayesian reconstructions including regularization in clinical computation times which are on the order of 1 min/slice. The most important result of the implementations is that the scaling rules for computation times are roughly linear in the number of processors.

[1]  Donald L. Snyder,et al.  Image Reconstruction from List-Mode Data in an Emission Tomography System Having Time-of-Flight Measurements , 1983, IEEE Transactions on Nuclear Science.

[2]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[3]  M.I. Miller,et al.  The role of likelihood and entropy in incomplete-data problems: Applications to estimating point-process intensities and toeplitz constrained covariances , 1987, Proceedings of the IEEE.

[4]  Linda Kaufman,et al.  Implementing and Accelerating the EM Algorithm for Positron Emission Tomography , 1987, IEEE Transactions on Medical Imaging.

[5]  R. A. Gaskins,et al.  Nonparametric roughness penalties for probability densities , 2022 .

[6]  R. Jaszczak,et al.  Inverse Monte Carlo: A Unified Reconstruction Algorithm for SPECT , 1985, IEEE Transactions on Nuclear Science.

[7]  M I Miller,et al.  Bayesian image reconstruction for emission tomography incorporating Good's roughness prior on massively parallel processors. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Robert M. Lewitt,et al.  Accelerated Iterative Reconstruction for Positron Emission Tomography Based on the EM Algorithm for Maximum Likelihood Estimation , 1986, IEEE Transactions on Medical Imaging.

[9]  Edward J. Hoffman,et al.  On The Convergence Of The Maximum Likelihood Estimator Method Of Tomographic Image Reconstruction , 1987, Medical Imaging.

[10]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[11]  Michael I. Miller,et al.  The Use of Sieves to Stabilize Images Produced with the EM Algorithm for Emission Tomography , 1985, IEEE Transactions on Nuclear Science.

[12]  Donald L. Snyder,et al.  A case study in statistical image processing: positron-emission tomography , 1991 .

[13]  K. Lange,et al.  A Theoretical Study of Some Maximum Likelihood Algorithms for Emission and Transmission Tomography , 1987, IEEE Transactions on Medical Imaging.

[14]  F. David Rollo,et al.  Nuclear medicine physics, instrumentation, and agents , 1977 .

[15]  Michael I. Miller,et al.  Bayesian imaging using Good's roughness measure-implementation on a massively parallel processor , 1988, ICASSP-88., International Conference on Acoustics, Speech, and Signal Processing.

[16]  M. Miller,et al.  Maximum-Likelihood Reconstruction for Single-Photon Emission Computed-Tomography , 1985, IEEE Transactions on Nuclear Science.

[17]  M. Miller,et al.  An evaluation of maximum likelihood reconstruction for SPECT. , 1990, IEEE transactions on medical imaging.

[18]  C E Floyd,et al.  Inverse Monte Carlo as a unified reconstruction algorithm for ECT. , 1986, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.