[en] An investigation has been conducted to develop and validate techniques for the correction of projection images in SPECT studies of the myocardium subject to misalignment due to voluntary patient motion. The problem is frequently encountered due to the uncomfortable position the patient must assume during the 30 minutes required to obtain a 180 degree set of projection images. The reconstruction of misaligned projections can lead to troublesome artifacts in reconstructed images and degrade the diagnostic potential of the procedure. Significant improvement in the quality of heart reconstructions has been realized with the implementation of an algorithm to provide detection of and correction for patient motion. Normal, involuntary motion is not corrected for, however, since such movement is below the spatial resolution of the thallium imaging system under study. The algorithm is based on a comparison of the positions of an object in a set of projection images to the known, sinusoidal trajectory of an off-axis fixed point in space. Projection alignment, therefore, is achieved by shifting the position of a point or set of points in a projection image to the sinusoid of a fixed position in space

[en] An infinite class of closed-form methods was developed by the authors last year for image reconstruction in 2D SPECT with uniform attenuation. In the work reported here, they extended their approach to develop a class of closed-form methods that compensate for the effects of both uniform attenuation and distance-dependent spatial resolution in 2D SPECT. These methods, which are characterized by an index n that can be assigned any real number, are exact in the absence of noise but propagate noise differently. They implemented this class of methods for SPECT image reconstruction in both computer-simulation and real-data studies. The results demonstrate that this class of methods corrects effectively for the aforementioned effects. Extensive computer simulation studies indicate that the method obtained with n = 2, which they had proved to be the optimal choice of n in 2D SPECT when only attenuation is present, also provides the smallest global image variance among the methods in the class when compensation for both uniform attenuation and distance-dependent spatial resolution is performed

[en] In PET, usually the data are precorrected for accidental coincidence (AC) events by real-time subtraction of the delayed window coincidences. Randoms subtraction compensates in mean for AC events but destroys the Poisson statistics. Furthermore, for transmission tomography the weighted least-squares (WLS) method leads to systematic biases, especially at low count rates. We propose a new open-quotes shiftedclose quotes Poisson (SP) model for precorrected PET data, which properly matches the first and second order moments of the measurement statistics. Using simulations and analytic approximations, we show that estimators based on the open-quotes ordinaryclose quotes Poisson (OP) model for the precorrected data lead to higher standard deviations than the proposed method. Moreover, if one zero-thresholds the data before applying the maximization algorithm, the OP model results in systematic bias. It is shown that the proposed SP model leads to penalized-likelihood estimates free of systematic bias, even for zero-thresholded data. The proposed SP model does not increase the computation requirements compared to OP model and it is robust to errors in the estimates of the AC event rates

[en] One of the most challenging problems in medical imaging is the development of reconstruction algorithms for nonstandard geometries. This work focuses on the application of Fourier analysis to the problem of resampling or rebinning. Conventional resampling methods utilizing some form of interpolation almost always result in a loss of resolution in the tomographic image. Fourier Transform Resampling (FTRS) offers potential improvement because the Modulation Transfer Function (MTF) of the process behaves like an ideal low pass filter. The MTF, however, is nonstationary if the coordinate transformation is nonlinear. FTRS may be viewed as a generalization of the linear coordinate transformations of standard Fourier analysis. Simulated MTF's were obtained by projecting point sources at different transverse positions in the flat fan beam detector geometry. These MTF's were compared to the closed form expression for FIRS. Excellent agreement was obtained for frequencies at or below the estimated cutoff frequency. The resulting FTRS algorithm is applied to simulations with symmetric fan beam geometry, an elliptical orbit and uniform attenuation, with a normalized root mean square error (NRME) of 0.036. Also, a Tc-99m point source study (1 cm dia., placed in air 10 cm from the COR) for a circular fan beam acquisition was reconstructed with a hybrid resampling method. The FWHM of the hybrid resampling method was 11.28 mm and compares favorably with a direct reconstruction (FWHM: 11.03 mm)

[en] Analytical expressions that describe the dependence of slopes and amplitudes of the scatter distribution functions (SDF) on source depth and media density are used to estimate a scatter component in SPECT projection data. Since the ratio of detected scattered to total photons (S/T), SDF amplitude and slope depend strongly on line source length (SL) used to obtain SDFs, we compared estimated scattered components using SDFs, obtained for lengths of 2-21 cm. At 10 cm source depth, S/T changes from 0.19 to 0.36 when SL changes from 2 to 21 cm. Scatter amplitude's dependence on source depth (d) in water was described by 6.38e^-0.186d for a 2 cm and 16.15e^-0.129d for a 21 cm SL. Slope was described by 0.292d^-0.601 for a cm SL and by 0.396d^-0.82 for a 21 cm SL. The estimated scatter components are compared with simulated SPECT projection data obtained with Monte Carlo modeling of six hot spheres placed in a cylindrical water filled phantom. The comparison of estimated with simulated total counts/projection shows very good agreement when approaching SDF for a point source (the % difference varied from 2 to 13% for 2 cm SL). Significant overestimate is seen when source length increases

[en] Smith and Grangeat derived a cone-beam inversion formula that can be applied when a nonplanar orbit satisfying the completeness condition is used. Although Grangeat's inversion formula is mathematically different from Smith's, they have similar overall structures to each other. The contribution of this paper is two-fold. First, based on the derivation of Smith, the authors point out that Grangeat's inversion formula and Smith's can be conveniently described using a single formula (the Smith-Grangeat inversion formula) that is in the form of space-variant filtering followed by cone-beam backprojection. Furthermore, the resulting formula is reformulated for data acquisition systems with a planar detector to obtain a new reconstruction algorithm. Second, the authors make two significant modifications to the new algorithm to reduce artifacts and numerical errors encountered in direct implementation of the new algorithm. As for exactness of the new algorithm, the following fact can be stated. The algorithm based on Grangeat's intermediate function is exact for any complete orbit, whereas that based on Smith's intermediate function should be considered as an approximate inverse excepting the special case where almost every plane in 3-D space meets the orbit. The validity of the new algorithm is demonstrated by simulation studies

[en] Noise equivalent counts are a convenient and effective means to assess PET emission image quality. The authors have extended the method to include the effects of transmission imaging on the statistics of attenuation corrected PET data. The result of the calculations is a noise figure which describes the SNR performance of the elements of the attenuation corrected emission sinogram. The noise figure demonstrates the tradeoff between emission and transmission imaging performance, and can be used to determine optimal partitioning of imaging time between emission and transmission scans. Also, the technique can be used to compare the efficacy of simultaneous transmission/emission imaging techniques and multiple orbiting rod source geometries. Experimental and simulated results from the GE 2,048 PET scanner are used to demonstrate the model. In a sample imaging situation in that system geometry [0.2 μCi/cc activity in a 20-cm flood phantom, 3-mCi orbiting rod source(s)], the dual rod source achieves 80 percent of the noise figure improvement which is available in simultaneous transmission/emission imaging without transmission data filtering, and demonstrates superior performance when a 3-point averaging transmission filter is applied

[en] The interest in fully three-dimensional image reconstruction, especially in positron emission tomography (PET) has significantly increased for the last few years. Taking into account the cross-plane gamma rays in a three-dimensional reconstruction algorithm improves sensitivity. At LETI, the speciality in PET is the time-of-flight (TOF) measurement. The authors present two reconstruction techniques for 3D TOF PET. The first is a backprojection-convolution algorithm. Due to the redundancy in the 3D data set, there exists an infinite number of filters. The authors establish a general condition that characterizes the filters and propose an algorithm with a factorizable filter. However, this first technique requires an acquisition system with revolution symmetry. Thus, they present a second one adapted to a detection geometry with a small number of angular positions. It consists of a multi-image deconvolution algorithm with Wiener filter. (author)

[en] Accidental suppression and effective data rates have been measured for the orbiting transmission source as implemented in the Donner 600-Crystal Positron Tomograph. A mechanical description of the orbiting source and a description of the electronics used to discard scattered and accidental events is included. Since accidental coincidences were the rate-limiting factor in transmission data acquisition, the new method allows us to acquire sufficient transmission data in a shorter time with a more active transmission source

[en] An elementary introduction is given to the basic principles of single photon emission computerized tomography (SPECT) and positron emission tomography (PET) imaging. 15 refs.; 9 figs