Computational burden is a major concern when an iterative algorithm is used to reconstruct a three-dimensional (3-D) image with attenuation, detector response, and scatter corrections. Most of the computation time is spent executing the projector and backprojector of an iterative algorithm. Usually, the projector and the backprojector are transposed operators of each other. The projector should model the imaging geometry and physics as accurately as possible. Some researchers have used backprojectors that are computationally less expensive than the projectors to reduce computation time. This paper points out that valid backprojectors should satisfy a condition that the projector/backprojector matrix must not contain negative eigenvalues. This paper also investigates the effects when unmatched projector/backprojector pairs are used.

A visualization algorithm was developed to demonstrate the complex 3-dimensional orientation of the fibers and sheets in human myocardium.

The point response resolution of a gamma camera deteriorates with increased distance from the face of the collimator. This results in reconstruction artifacts that are seen as shape distortion and density non-uniformity. For parallel, fan, and cone beam geometries, iterative reconstruction algorithms have been developed which eliminate these artifacts by incorporating the three-dimensional spatially varying geometric point response into models for the projection and backprojection operations which also model photon attenuation. The algorithms have been tested on an IBM 3090-600s supercomputer. The iterative EM reconstruction algorithm is 50 times longer with geometric response and photon attenuation models than without modeling these physical effects. We have demonstrated an improvement in image quality in the reconstruction of projection data collected from a single photon emission computed tomography (SPECT) imaging system. Using phantom experiments, we observe that the modeling of the spatial system response imposes a smoothing without loss of resolution.

Echo-planar imaging (EPI) is very sensitive to patient-induced field inhomogeneity caused by susceptibility changes between different anatomical regions. This results in geometric and intensity distortions in the image, especially near tissue/air and tissue/bone interfaces. A new approach is presented to reduce geometric and intensity distortions in EPI. A phase-encoded multireference scan is used to estimate the amplitude and phase errors in the measured signals due to the field inhomogeneity. The EPI data is corrected using both the amplitude and phase of the measured errors. This technique has been evaluated using EPI pulse sequences implemented with conventional gradients and implemented with imaging systems that have special resonating gradients and fast analog to digital converters. The results in both phantom and human studies show that in the absence of object motion the new correction technique can effectively reduce the geometric and intensity distortions.

The Compton camera can collect SPECT data with high efficiency due to electronic collimation. The data acquired from a Compton camera are projections of source activity along cones and are approximated in this paper by cone-surface integrals. This paper proposes the use of an orthogonal spherical expansion to convert the cone-surface integrals into plane integrals. The conversion technique is efficient. Once the plane integrals are obtained, a 3D image can be reconstructed by the 3D Radon inversion formula. The algorithm is implemented and computer simulations are used to demonstrate the efficiency and accuracy of the proposed reconstruction algorithm.

No abstract

The very nature of nuclear medicine, the visual representation of injected radiopharmaceuticals, implies imaging of dynamic processes such as the uptake and wash-out of radiotracers from body organs. For years, nuclear medicine has been touted as the modality of choice for evaluating function in health and disease. This evaluation is greatly enhanced using single photon emission computed tomography (SPECT), which permits three-dimensional (3D) visualization of tracer distributions in the body. However, to fully realize the potential of the technique requires the imaging of in vivo dynamic processes of flow and metabolism. Tissue motion and deformation must also be addressed. Absolute quantification of these dynamic processes in the body has the potential to improve diagnosis. This paper presents a review of advancements toward the realization of the potential of dynamic SPECT imaging and a brief history of the development of the instrumentation. A major portion of the paper is devoted to the review of special data processing methods that have been developed for extracting kinetics from dynamic cardiac SPECT data acquired using rotating detector heads that move as radiopharmaceuticals exchange between biological compartments. Recent developments in multi-resolution spatiotemporal methods enable one to estimate kinetic parameters of compartment models of dynamic processes using data acquired from a single camera head with slow gantry rotation. The estimation of kinetic parameters directly from projection measurements improves bias and variance over the conventional method of first reconstructing 3D dynamic images, generating time–activity curves from selected regions of interest and then estimating the kinetic parameters from the generated time–activity curves. Although the potential applications of SPECT for imaging dynamic processes have not been fully realized in the clinic, it is hoped that this review illuminates the potential of SPECT for dynamic imaging, especially in light of new developments that enable measurement of dynamic processes directly from projection measurements.

A method is presented for estimating the geometrical parameters for a cone beam detector geometry from the coordinates of the centroid of a projected point source sampled over 360 degrees. Nonlinear expressions are derived for the coordinates of the centroids in terms of the geometrical parameters which include: the two-dimensional coordinates of the projection of the center of rotation onto the detector image plane; the focal length; the distance from the focal point to the center of rotation; and the spatial coordinates of the point source itself. Experimental data were obtained using a rotating gamma camera with a symmetrically converging collimator. The Marquardt algorithm was used to estimate the parameters for this particular cone beam geometry. The method was able to estimate the geometrical parameters and evaluate the accuracy of the collimator construction.