Linear least squares compartmental-model-independent parameter identification in PET

Quality and precision of parametric images created from PET sinogram data by direct reconstruction: proof of concept

We have previously implemented the direct reconstruction of dense kinetic model parameter images ("parametric images") from sinogram data, and compared it to conventional image domain kinetic parameter estimation methods . Although it has been shown that the direct reconstruction algorithm estimates the kinetic model parameters with lower root mean squared error than the conventional image domain techniques, some theoretical obstacles remain. These obstacles include the difficulty of evaluating the accuracy and precision of the estimated parameters. In image domain techniques, the reconstructed time activity curve (TAC) and the model predicted TAC are compared, and the goodness-of-fit is evaluated as a measure of the accuracy and precision of the estimated parameters. This approach cannot be applied to the direct reconstruction technique as there are no reconstructed TACs. In this paper, we propose ways of evaluating the precision and goodness-of-fit of the kinetic model parameters estimated by the direct reconstruction algorithm. Specifically, precision of the estimates requires the calculation of variance images for each parameter, and goodness-of-fit is addressed by reconstructing the difference between the measured and the fitted sinograms. We demonstrate that backprojecting the difference from sinogram space to image space creates error images that can be examined for goodness-of-fit and model selection purposes. The presence of nonrandom structures in the error images may indicate an inadequacy of the kinetic model that has been incorporated into the direct reconstruction algorithm. We introduce three types of goodness-of-fit images. We propose and demonstrate a number-of-runs image as a means of quantifying the adequacy or deficiency of the model. We further propose and demonstrate images of the F statistic and the change in the Akaike Information Criterion as devices for identifying the statistical advantage of one model over another at each voxel. As direct reconstruction to parametric images proliferates, it will be essential for imagers to adopt methods such as those proposed herein to assess the accuracy and precision of their parametric images.

Direct reconstruction of kinetic parameter images from dynamic PET data

Our goal in this paper is the estimation of kinetic model parameters for each voxel corresponding to a dense three-dimensional (3-D) positron emission tomography (PET) image. Typically, the activity images are first reconstructed from PET sinogram frames at each measurement time, and then the kinetic parameters are estimated by fitting a model to the reconstructed time-activity response of each voxel. However, this "indirect" approach to kinetic parameter estimation tends to reduce signal-to-noise ratio (SNR) because of the requirement that the sinogram data be divided into individual time frames. In 1985, Carson and Lange proposed, but did not implement, a method based on the expectation-maximization (EM) algorithm for direct parametric reconstruction. The approach is "direct" because it estimates the optimal kinetic parameters directly from the sinogram data, without an intermediate reconstruction step. However, direct voxel-wise parametric reconstruction remained a challenge due to the unsolved complexities of inversion and spatial regularization. In this paper, we demonstrate and evaluate a new and efficient method for direct voxel-wise reconstruction of kinetic parameter images using all frames of the PET data. The direct parametric image reconstruction is formulated in a Bayesian framework, and uses the parametric iterative coordinate descent (PICD) algorithm to solve the resulting optimization problem. The PICD algorithm is computationally efficient and is implemented with spatial regularization in the domain of the physiologically relevant parameters. Our experimental simulations of a rat head imaged in a working small animal scanner indicate that direct parametric reconstruction can substantially reduce root-mean-squared error (RMSE) in the estimation of kinetic parameters, as compared to indirect methods, without appreciably increasing computation.

Relevance of accurate Monte Carlo modeling in nuclear medical imaging

Monte Carlo techniques have become popular in different areas of medical physics with advantage of powerful computing systems. In particular, they have been extensively applied to simulate processes involving random behavior and to quantify physical parameters that are difficult or even impossible to calculate by experimental measurements. Recent nuclear medical imaging innovations such as single-photon emission computed tomography (SPECT), positron emission tomography (PET), and multiple emission tomography (MET) are ideal for Monte Carlo modeling techniques because of the stochastic nature of radiation emission, transport and detection processes. Factors which have contributed to the wider use include improved models of radiation transport processes, the practicality of application with the development of acceleration schemes and the improved speed of computers. In this paper we present a derivation and methodological basis for this approach and critically review their areas of application in nuclear imaging. An overview of existing simulation programs is provided and illustrated with examples of some useful features of such sophisticated tools in connection with common computing facilities and more powerful multiple-processor parallel processing systems. Current and future trends in the field are also discussed.