Effects of temporal modelling on the statistical uncertainty of spatiotemporal distributions estimated directly from dynamic SPECT projections.

Artefacts can result when reconstructing a dynamic image sequence from inconsistent single photon emission computed tomography (SPECT) projection data acquired by a slowly rotating gantry. The artefacts can lead to biases in kinetic parameters estimated from time-activity curves generated by overlaying volumes of interest on the images. Insufficient sampling and truncation of projections by cone-beam collimators can cause additional artefacts. To overcome these sources of bias in conventional image based dynamic data analysis, we have been investigating the estimation of time-activity curves and kinetic model parameters directly from dynamic SPECT projection data by modelling the spatial and temporal distribution of the radiopharmaceutical throughout the projected field of view. In the present work, we perform Monte Carlo simulations to study the effects of the temporal modelling on the statistical variability of the reconstructed spatiotemporal distributions. The simulations utilize fast methods for fully four-dimensional (4D) direct estimation of spatiotemporal distributions and their statistical uncertainties, using a spatial segmentation and temporal B-splines. The simulation results suggest that there is benefit in modelling higher orders of temporal spline continuity. In addition, the accuracy of the time modelling can be increased substantially without unduly increasing the statistical uncertainty, by using relatively fine initial time sampling to capture rapidly changing activity distributions.

Kinetic analysis of 125I-iodorotenone as a deposited myocardial flow tracer: comparison with 99mTc-sestamibi.

UNLABELLED: The goal of this investigation was to assess the accuracy of 7'-Z-[125I]iodorotenone (125I-iodorotenone) as a new deposited myocardial flow tracer and compare the results with those for 99mTc-sestamibi. METHODS: The kinetics of these two flow tracers were evaluated in 25 isolated, erythrocyte- and albumin-perfused rabbit hearts over a flow range relevant to patients. The two flow tracers and a vascular reference tracer (131I-albumin) were introduced simultaneously as a compact bolus through a port just above the aortic cannula in the absence of tracer recirculation. Myocardial extraction, retention, washout, and uptake parameters were computed from the venous outflow curves using the multiple-indicator dilution technique and spectral analysis. RESULTS: The extraction of 125I-iodorotenone was much higher than the extraction of 99mTc-sestamibi (0.84 +/- 0.05 vs. 0.48 +/- 0.10, respectively, P < 0.001). 125I-iodorotenone extraction was also less affected by flow than was 99mTc-sestamibi (P < 0.001). Net retention of 125I-iodorotenone was significantly greater than 99mTc-sestamibi net retention at 1 min (0.77 +/- 0.08 vs. 0.41 +/- 0.11, respectively, P < 0.001) and 26 min (0.46 +/- 0.13 vs. 0.27 +/- 0.11, respectively, P < 0.001) after tracer injection. Flow had less effect on 125I-iodorotenone net retention than on 99mTc-sestamibi net retention 1 min after tracer injection (P < 0.04). However, at 26 min, flow had an equivalent effect on the retention of both flow tracers (P < 0.4). The relationship between 125I-iodorotenone and 99mTc-sestamibi washout was complex and depended on elapsed time after isotope introduction and perfusion rate. Reflecting the favorable extraction and retention characteristics of 125I-iodorotenone, both its maximum myocardial uptake and its 26-min uptake were more closely related to flow than were those of 99mTc-sestamibi (P < 0.001 for both comparisons). CONCLUSION: The extraction and retention of 125I-iodorotenone were greater than those of 99mTc-sestamibi, making 125I-iodorotenone the superior flow tracer in the isolated rabbit heart.

Direct least-squares estimation of spatiotemporal distributions from dynamic SPECT projections using a spatial segmentation and temporal B-splines.

Artifacts can result when reconstructing a dynamic image sequence from inconsistent, as well as insufficient and truncated, cone beam single photon emission computed tomography (SPECT) projection data acquired by a slowly rotating gantry. The artifacts can lead to biases in kinetic model parameters estimated from time-activity curves generated by overlaying volumes of interest on the images. However, the biases in time-activity curve estimates and subsequent kinetic parameter estimates can be reduced significantly by first modeling the spatial and temporal distribution of the radiopharmaceutical throughout the projected field of view, and then estimating the time-activity curves directly from the projections. This approach is potentially useful for clinical SPECT studies involving slowly rotating gantries, particularly those using a single-detector system or body contouring orbits with a multidetector system. We have implemented computationally efficient methods for fully four-dimensional (4-D) direct estimation of spatiotemporal distributions from dynamic SPECT projection data. Temporal B-splines providing various orders of temporal continuity, as well as various time samplings, were used to model the time-activity curves for segmented blood pool and tissue volumes in simulated cone beam and parallel beam cardiac data acquisitions. Least-squares estimates of time-activity curves were obtained quickly using a workstation. Given faithful spatial modeling, accurate curve estimates were obtained using cubic, quadratic, or linear B-splines and a relatively rapid time sampling during initial tracer uptake. From these curves, kinetic parameters were estimated accurately for noiseless data and with some bias for noisy data. A preliminary study of spatial segmentation errors showed that spatial model mismatch adversely affected quantitative accuracy, but also resulted in structured errors (projected model versus raw data) that were easily detected in our simulations. This suggests iterative refinement of the spatial model to reduce structured errors as an area of future research.

List-mode maximum-likelihood reconstruction applied to positron emission mammography (PEM) with irregular sampling.

We present a preliminary study of list-mode likelihood reconstruction of images for a rectangular positron emission tomograph (PET) specifically designed to image the human breast. The prospective device consists of small arrays of scintillation crystals for which depth of interaction is estimated. Except in very rare instances, the number of annihilation events detected is expected to be far less than the number of distinguishable events. If one were to histogram the acquired data, most histogram bins would remain vacant. Therefore, it seems natural to investigate the efficacy of processing events one at a time rather than processing the data in histogram format. From a reconstruction perspective, the new tomograph presents a challenge in that the rectangular geometry leads to irregular radial and angular sampling, and the field of view extends completely to the detector faces. Simulations are presented that indicate that the proposed tomograph can detect 8-mm-diameter spherical tumors with a tumor-to-background tracer density ratio of 3:1 using realistic image acquisition parameters. Spherical tumors of 4-mm diameter are near the limit of detectability with the image acquisition parameters used. Expressions are presented to estimate the loss of image contrast due to Compton scattering.

Estimating glucose metabolism using glucose analogs and two tracer kinetic models in isolated rabbit heart.

The purpose of this investigation was to 1) evaluate the relative accuracy of the Sokoloff and Patlak tracer kinetic models in estimating glucose metabolic rate (GMR) in the presence and absence of insulin; 2) evaluate the effect of nutritional state on the lumped constant (LC); and 3) compare the kinetics of 2-fluoro-2-deoxy-D-[14C]glucose (FDG) and 2-deoxy-D-[3H]glucose (DG) membrane transport and phosphorylation. The experimental preparation was the isolated, red blood cell-albumin-perfused rabbit heart. Our results showed that both tracer kinetic models provided GMR estimates that correlated well with the Fick method (for FDG, R = 0. 84 and 0.91 for the Sokoloff and Patlak models, respectively); nutritional state did not affect the LC; and FDG and DG have different transport and/or phosphorylation parameters. We also observed that 1) the addition of a fourth compartment to the Sokoloff model reduced the mean squared error between measured and modeled data by a factor of 7.4; 2) a longer time (21.8 min) was required to obtain a linear phase of the Patlak plot than is allowed in clinical studies; and 3) accurate GMR estimates were obtained only by using different LCs reflecting insulin's presence or absence. Our results indicate potential sources of error in the use of FDG and positron emission tomography to quantify GMR in patients.

Kinetic parameter estimation from SPECT cone-beam projection measurements.

Kinetic parameters are commonly estimated from dynamically acquired nuclear medicine data by first reconstructing a dynamic sequence of images and subsequently fitting the parameters to time-activity curves generated from regions of interest overlaid upon the image sequence. Biased estimates can result from images reconstructed using inconsistent projections of a time-varying distribution of radiopharmaceutical acquired by a rotating SPECT system. If the SPECT data are acquired using cone-beam collimators wherein the gantry rotates so that the focal point of the collimators always remains in a plane, additional biases can arise from images reconstructed using insufficient, as well as truncated, projection samples. To overcome these problems we have investigated the estimation of kinetic parameters directly from SPECT cone-beam projection data by modelling the data acquisition process. To accomplish this it was necessary to parametrize the spatial and temporal distribution of the radiopharmaceutical within the SPECT field of view. In a simulated chest image volume, kinetic parameters were estimated for simple one-compartment models for four myocardial regions of interest. Myocardial uptake and washout parameters estimated by conventional analysis of noiseless simulated cone-beam data had biases ranging between 3-26% and 0-28%, respectively. Parameters estimated directly from the noiseless projection data were unbiased as expected, since the model used for fitting was faithful to the simulation. Statistical uncertainties of parameter estimates for 10,000,000 events ranged between 0.2-9% for the uptake parameters and between 0.3-6% for the washout parameters.

A methodology for specifying PET VOI's using multimodality techniques.

Volume-of-interest (VOI) extraction for radionuclide and anatomical measurements requires correct identification and delineation of the anatomical feature being studied. We have developed a toolset for specifying three-dimensional (3-D) VOI's on a multislice positron emission tomography (PET) dataset. The software is particularly suited for specifying cerebral cortex VOI's which represent a particular gyrus or deep brain structure. A registered 3-D magnetic resonance image (MRI) dataset is used to provide high-resolution anatomical information, both as oblique two-dimensional (2-D) sections and as volume renderings of a segmented cortical surface. VOI's are specified indirectly in two dimensions by drawing a stack of 2-D regions on the MRI data. The regions are tiled together to form closed triangular mesh surface models, which are subsequently transformed into the observation space of the PET scanner. Quantification by this method allows calculation of radionuclide activity in the VOI's, as well as their statistical uncertainties and correlations. The methodology for this type of analysis and validation results are presented.

Kinetic analysis of rubidium and thallium as deposited myocardial blood flow tracers in isolated rabbit heart.

Evaluation of myocardial perfusion with tracers such as thallium and rubidium is based on the assumption that tissue tracer content is proportional to flow. The purpose of this study was to evaluate the relationship between flow and tissue tracer content of 201Tl and 83Rb in the isolated perfused rabbit heart. 83Rb (86-day half-life), an isotope that is not used clinically, was used as a subsitute for 82Rb (76-s half-life) to improve the accuracy and precision of data acquisition. The multiple indicator-dilution technique was employed with two independent computational approaches. The first approach explicitly deconvolved 201Tl and 83Rb venous concentration curves by the intravascular reference tracer curve. The second approach used a conventional analysis. Both approaches showed that there was more early washout of 83Rb than 201Tl and that the heart retained 201Tl better than 83Rb within 2 min after isotope introduction. These data indicate that 201Tl is a better perfusion tracer than 83Rb in the isolated rabbit heart.

Strategies for extraction of quantitative data from volumetric dynamic cardiac positron emission tomography data.

The ability of positron emission tomography (PET) to serve as a useful myocardial perfusion indicator is well established. We describe a methodology for obtaining reliable quantitative kinetic parameters from dynamic cardiac PET data. Reconstructed images of the myocardium are subdivided into three-dimensional volumes of interest which are used to obtain quantitative measures of myocardial perfusion over physiologically meaningful anatomical regions. The quantitation technique rigorously models the uncertainty of estimated parameters while compensating for effects such as patient motion and partial volumes to arrive at model parameters with well-established confidence intervals.