Fourier correction for spatially variant collimator blurring in SPECT

In single-photon emission computed tomography (SPECT), projection data are acquired by rotating the photon detector around a patient, either in a circular orbit or in a noncircular orbit. The projection data of the desired spatial distribution of emission activity is blurred by the point-response function of the collimator that is used to define the range of directions of gamma-ray photons reaching the detector. The point-response function of the collimator is not spatially stationary, but depends on the distance from the collimator to the point. Conventional methods for deblurring collimator projection data are based on approximating the actual distance-dependent point-response function by a spatially invariant blurring function, so that deconvolution methods can be applied independently to the data at each angle of view. A method is described here for distance-dependent preprocessing of SPECT projection data prior to image reconstruction. Based on the special distance-dependent characteristics of the Fourier coefficients of the sinogram, a spatially variant inverse filter can be developed to process the projection data in all views simultaneously. The algorithm is first derived from Fourier analysis of the projection data from the circular orbit geometry. For circular orbit projection data, experimental results from both simulated data and real phantom data indicate the potential of this method. It is shown that the spatial filtering method can be extended to the projection data from the noncircular orbit geometry. Experiments on simulated projection data from an elliptical orbit demonstrate correction of the spatially variant blurring and distortion in the reconstructed image caused by the noncircular orbit geometry