No abstract

A procedure for the determination of the system matrix in single photon emission tomography (SPECT) is described which uses the conjugate gradient reconstruction technique in order to take into account the variable system resolution of a camera equipped with parallel-hole collimators. The procedure involves the acquisition of the system line spread functions (LSF) in the region occupied by the object to be studied. Those data are used to generate a set of weighting factors based on the assumption that the LSFs of the collimated camera are of Gaussian shape with the full width at half maximum (FWHM) linearly dependent on the source depth in the span of image space. The factors are stored on a disc file for subsequent use in the reconstruction process. Afterwards the reconstruction is performed using the conjugate gradient method with the system matrix modified by the incorporation of these precalculated factors in order to take into account the variable geometrical system response. The set of weighting factors is regenerated whenever the acquisition conditions are changed (collimator, radius of rotation). In the case of an ultra high resolution (UHR) collimator 2000 weighting factors need to be calculated. The modification of the system matrix for the geometrical response allows the number of iterations to increase, considerably improving image definition without the appearance of noise artifacts. Moreover, phantom studies show that the number of iterations is less critical because of improved stability in the convergence to the solution. For brain studies of patients 10-15 iterations are usually performed. Studies with a single line source give a value between 7 and 8 mm for the FWHM of the point spread function (PSF) when the conjugate gradient method with modified system matrix is used on data acquired with a UHR collimator, whereas without the modification of the system matrix the result is 9 mm FWHM, if filtered backprojection (FBP) is used with the same filter as in the clinical studies the result is 15 mm FWHM. The results of this work show that proper definition of the system matrix using conjugate gradients influences the quality of the reconstruction remarkably. Nevertheless, further work has to be done in order to assess to what extent the system matrix is ill-conditioned and, eventually, to define a suitable regularization technique.

An accurate model of the nonstationary geometrical response of a camera-collimator system is discussed. The algorithm is compared to three other algorithms that are specialized for region-of-interest evaluation, as well as to the conventional method for summing the reconstructed quantity over the regions of interest. For noise-free data and for regions of accurate shape, least-squares estimates were unbiased within roundoff errors. For noisy data, estimates were still unbiased but precision worsened for regions smaller than resolution: simulating typical statistics of brain perfusion studies performed with a collimated camera, the estimated standard deviation for a 1-cm-square region was 10% with an ultra-high-resolution collimator and 7% with a low-energy all-purpose collimator. Conventional region-of-interest estimates show comparable precision but are heavily biased if filtered backprojection is used for image reconstruction. Using the conjugate-gradient iterative algorithm and the model of nonstationary geometrical response, bias of estimates decreased on increasing the number of iterations, but precision worsened, thus achieving an estimated standard deviation of more than 25% for the same 1-cm region.

Putaminal BP(LS) is an efficient parameter for exploring the correlations between PD severity and rCMRglc cortical changes. The correlation between dopaminergic degeneration and rCMRglc in several prefrontal regions likely represents the cortical functional correlate of the dysfunction in the motor basal ganglia-cortical circuit in PD. This finding suggests focusing on the metabolic course of these areas to follow PD progression and to analyze treatment effects.

A complete theory of camera multihole collimators is presented. The geometrical system response is determined in closed form in frequency space. This closed form accounts for the known efficiency and resolution formulae for parallel beam, fan beam, cone beam and astigmatic collimators as well as for the most frequent hole array patterns and hole shapes. The point spread function in the space domain for a certain collimator and source position can be calculated via a discrete fast Fourier transform. Beside the complete theoretical definition of the response of multihole collimators, this theory allows the definition of accurate models of the geometrical response for SPECT reconstruction and it is suitable for designing new collimators.