Distributed 3-D iterative reconstruction for quantitative SPECT

We describe a distributed three dimensional (3-D) iterative reconstruction library for quantitative single-photon emission computed tomography (SPECT). This library includes 3-D projector-backprojector pairs (PBPs) and distributed 3-D iterative reconstruction algorithms. The 3-D PBPs accurately and efficiently model various combinations of the image degrading factors including attenuation, detector response and scatter response. These PBPs were validated by comparing projection data computed using the projectors with that from direct Monte Carlo (MC) simulations. The distributed 3-D iterative algorithms spread the projection-backprojection operations for all the projection angles over a heterogeneous network of single or multi-processor computers to reduce the reconstruction time. Based on a master/slave paradigm, these distributed algorithms provide dynamic load balancing and fault tolerance. The distributed algorithms were verified by comparing images reconstructed using both the distributed and non-distributed algorithms. Computation times for distributed 3-D reconstructions running on up to 4 identical processors were reduced by a factor approximately 80-90% times the number of the processors participating, compared to those for non-distributed 3-D reconstructions running on a single processor. When combined with faster affordable computers, this library provides an efficient means for implementing accurate reconstruction and compensation methods to improve quality and quantitative accuracy in SPECT images. >

[1]  C E Metz,et al.  The geometric transfer function component for scintillation camera collimators with straight parallel holes. , 1980, Physics in medicine and biology.

[2]  Eric C. Frey,et al.  A fast projector-backprojector pair modeling the asymmetric, spatially varying scatter response function for scatter compensation in SPECT imaging , 1993 .

[3]  N. Pelc,et al.  An attenuated projector-backprojector for iterative SPECT reconstruction. , 1985, Physics in medicine and biology.

[4]  E. Frey,et al.  A practical method for incorporating scatter in a projector-backprojector for accurate scatter compensation in SPECT , 1993 .

[5]  M. Ljungberg,et al.  A Monte Carlo program for the simulation of scintillation camera characteristics. , 1989, Computer methods and programs in biomedicine.

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

[7]  T A Riauka,et al.  Photon propagation and detection in single-photon emission computed tomography--an analytical approach. , 1994, Medical physics.

[8]  Grant T. Gullberg,et al.  Users Manual: Donner Algorithms for Reconstruction Tomography , 1977 .

[9]  E. W. Marchand From Line to Point Spread Function: The General Case , 1965 .

[10]  Max A. Viergever,et al.  Object shape dependent PSF model for SPECT imaging , 1993 .

[11]  Max A. Viergever,et al.  A new phantom for fast determination of the scatter response of a gamma camera , 1993 .

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

[13]  G.T. Gullberg,et al.  Frequency domain implementation of the three-dimensional geometric point response correction in SPECT imaging , 1991, Conference Record of the 1991 IEEE Nuclear Science Symposium and Medical Imaging Conference.

[14]  E. Hoffman,et al.  3-D phantom to simulate cerebral blood flow and metabolic images for PET , 1990 .

[15]  D. Gilland,et al.  Implementation of simultaneous attenuation and detector response correction in SPECT , 1988 .