Search

Presents information on a step being taken towards performing respiration compensated attenuation. Details on the experiment that was done; Information on the hardware gating; Identification of the automated 3-D transmission image segmentation.

Computer simulations, a phantom study and a human study were performed to determine whether a slowly rotating single-photon computed emission tomography (SPECT) system could provide accurate arterial input functions for quantification of myocardial perfusion imaging using kinetic models. The errors induced by data inconsistency associated with imaging with slow camera rotation during tracer injection were evaluated with an approach called SPECT/P (dynamic SPECT from positron emission tomography (PET)) and SPECT/D (dynamic SPECT from database of SPECT phantom projections). SPECT/P simulated SPECT-like dynamic projections using reprojections of reconstructed dynamic (94)Tc-methoxyisobutylisonitrile ((94)Tc-MIBI) PET images acquired in three human subjects (1 min infusion). This approach was used to evaluate the accuracy of estimating myocardial wash-in rate parameters K(1) for rotation speeds providing 180° of projection data every 27 or 54 s. Blood input and myocardium tissue time-activity curves (TACs) were estimated using spatiotemporal splines. These were fit to a one-compartment perfusion model to obtain wash-in rate parameters K(1). For the second method (SPECT/D), an anthropomorphic cardiac torso phantom was used to create real SPECT dynamic projection data of a tracer distribution derived from (94)Tc-MIBI PET scans in the blood pool, myocardium, liver and background. This method introduced attenuation, collimation and scatter into the modeling of dynamic SPECT projections. Both approaches were used to evaluate the accuracy of estimating myocardial wash-in parameters for rotation speeds providing 180° of projection data every 27 and 54 s. Dynamic cardiac SPECT was also performed in a human subject at rest using a hybrid SPECT/CT scanner. Dynamic measurements of (99m)Tc-tetrofosmin in the myocardium were obtained using an infusion time of 2 min. Blood input, myocardium tissue and liver TACs were estimated using the same spatiotemporal splines. The spatiotemporal maximum-likelihood expectation-maximization (4D ML-EM) reconstructions gave more accurate reconstructions than did standard frame-by-frame static 3D ML-EM reconstructions. The SPECT/P results showed that 4D ML-EM reconstruction gave higher and more accurate estimates of K(1) than did 3D ML-EM, yielding anywhere from a 44% underestimation to 24% overestimation for the three patients. The SPECT/D results showed that 4D ML-EM reconstruction gave an overestimation of 28% and 3D ML-EM gave an underestimation of 1% for K(1). For the patient study the 4D ML-EM reconstruction provided continuous images as a function of time of the concentration in both ventricular cavities and myocardium during the 2 min infusion. It is demonstrated that a 2 min infusion with a two-headed SPECT system rotating 180° every 54 s can produce measurements of blood pool and myocardial TACs, though the SPECT simulation studies showed that one must sample at least every 30 s to capture a 1 min infusion input function.

Background: There is a well-recognized need for a new generation of single photon emission computed tomography (SPECT) perfusion tracers with improved myocardial extraction over a wide flow range. Radiotracers that target complex I of the mitochondrial electron transport chain have been proposed as a new class of myocardial perfusion imaging agents. 7-(Z)-[(125)I]iodorotenone ((125)I-ZIROT) has demonstrated superior myocardial extraction and retention characteristics in rats and in isolated perfused rabbit hearts. We sought to fully characterize the biodistribution and myocardial extraction versus flow relationship of (123)I-ZIROT in an intact large-animal model., Methods and Results: The (123)I-ZIROT was administered during adenosine A(2A) agonist-induced hyperemia in 5 anesthetized dogs with critical left anterior descending (LAD) stenoses. When left circumflex (LCx) flow was maximal, (123)I-ZIROT and microspheres were coinjected and the dogs were euthanized 5 minutes later. (123)I-ZIROT biodistribution was evaluated in 2 additional dogs by in vivo planar imaging. At (123)I-ZIROT injection, transmural LAD flow was unchanged from baseline (mean±SEM, 0.90±0.22 versus 0.87±0.11 mL/[min · g]; P=0.92), whereas LCx zone flow increased significantly (mean±SEM, 3.25±0.51 versus 1.00±0.17 mL/[min · g]; P<0.05). Myocardial (123)I-ZIROT extraction tracked regional myocardial flow better than either thallium-201 or (99m)Tc-sestamibi from previous studies using a similar model. Furthermore, the (123)I-ZIROT LAD/LCx activity ratios by ex vivo imaging or well counting (mean±SEM, 0.42±0.08 and 0.45±0.1, respectively) only slightly underestimated the LAD/LCx microsphere flow ratio (0.32±0.09)., Conclusions: The ability of (123)I-ZIROT to more linearly track blood flow over a wide range makes it a promising new SPECT myocardial perfusion imaging agent with potential for improved coronary artery disease detection and better quantitative estimation of the severity of flow impairment.

The goal of this project is to develop radionuclide molecular imaging technologies using a clinical pinhole SPECT/CT scanner to quantify changes in cardiac metabolism using the spontaneously hypertensive rat (SHR) as a model of hypertensive-related pathophysiology. This paper quantitatively compares fatty acid metabolism in hearts of SHR and Wistar-Kyoto normal rats as a function of age and thereby tracks physiological changes associated with the onset and progression of heart failure in the SHR model. The fatty acid analog, (123)I-labeled BMIPP, was used in longitudinal metabolic pinhole SPECT imaging studies performed every seven months for 21 months. The uniqueness of this project is the development of techniques for estimating the blood input function from projection data acquired by a slowly rotating camera that is imaging fast circulation and the quantification of the kinetics of (123)I-BMIPP by fitting compartmental models to the blood and tissue time-activity curves.

The very nature of nuclear medicine, the visual representation of injected radiopharmaceuticals, implies imaging of dynamic processes such as the uptake and wash-out of radiotracers from body organs. For years, nuclear medicine has been touted as the modality of choice for evaluating function in health and disease. This evaluation is greatly enhanced using single photon emission computed tomography (SPECT), which permits three-dimensional (3D) visualization of tracer distributions in the body. However, to fully realize the potential of the technique requires the imaging of in vivo dynamic processes of flow and metabolism. Tissue motion and deformation must also be addressed. Absolute quantification of these dynamic processes in the body has the potential to improve diagnosis. This paper presents a review of advancements toward the realization of the potential of dynamic SPECT imaging and a brief history of the development of the instrumentation. A major portion of the paper is devoted to the review of special data processing methods that have been developed for extracting kinetics from dynamic cardiac SPECT data acquired using rotating detector heads that move as radiopharmaceuticals exchange between biological compartments. Recent developments in multi-resolution spatiotemporal methods enable one to estimate kinetic parameters of compartment models of dynamic processes using data acquired from a single camera head with slow gantry rotation. The estimation of kinetic parameters directly from projection measurements improves bias and variance over the conventional method of first reconstructing 3D dynamic images, generating time-activity curves from selected regions of interest and then estimating the kinetic parameters from the generated time-activity curves. Although the potential applications of SPECT for imaging dynamic processes have not been fully realized in the clinic, it is hoped that this review illuminates the potential of SPECT for dynamic imaging, especially in light of new developments that enable measurement of dynamic processes directly from projection measurements.

The objective of this research was to assess applicability of a technique known as hyperelastic warping for the measurement of local strains in the left ventricle (LV) directly from microPET image data sets. The technique uses differences in image intensities between template (reference) and target (loaded) image data sets to generate a body force that deforms a finite element (FE) representation of the template so that it registers with the target images. For validation, the template image was defined as the end-systolic microPET image data set from a Wistar Kyoto (WKY) rat. The target image was created by mapping the template image using the deformation results obtained from a FE model of diastolic filling. Regression analysis revealed highly significant correlations between the simulated forward FE solution and image derived warping predictions for fiber stretch (R (2) = 0.96), circumferential strain (R (2) = 0.96), radial strain (R (2) = 0.93), and longitudinal strain (R (2) = 0.76) (p < 0.001 for all cases). The technology was applied to microPET image data of two spontaneously hypertensive rats (SHR) and a WKY control. Regional analysis revealed that, the lateral freewall in the SHR subjects showed the greatest deformation compared with the other wall segments. This work indicates that warping can accurately predict the strain distributions during diastole from the analysis of microPET data sets.

Unlabelled: Currently, 2 types of phantoms (physical and computer generated) are used for testing and comparing tomographic reconstruction methods. Data from physical phantoms include all physical effects associated with the detection of radiation. However, with physical phantoms it is difficult to control the number of detected counts, simulate the dynamics of uptake and washout, or create multiple noise realizations of an acquisition. Computer-generated phantoms can overcome some of the disadvantages of physical phantoms, but simulation of all factors affecting the detection of radiation is extremely complex and in some cases impossible. To overcome the problems with both types of phantoms, we developed a physical and computer-generated hybrid phantom that allows the creation of multiple noise realizations of tomographic datasets of the dynamic uptake governed by kinetic models., Methods: The method is phantom and camera specific. We applied it to an anthropomorphic torso phantom with a cardiac insert, using a SPECT system with attenuation correction. First, real data were acquired. For each compartment (heart, blood pool, liver, and background) of the physical phantom, large numbers of short tomographic projections were acquired separately for each angle. Sinograms were built from a database of projections by summing the projections of each compartment of the phantom. The amount of activity in each phantom compartment was regulated by the number of added projections. Sinograms corresponding to various projection times, configurations and numbers of detector heads, numbers of noise realizations, numbers of phantom compartments, and compartment-specific time-activity curves in MBq/cm3 were assembled from the database., Results: The acquisition produced a database of 120 projection angles ranging over 360 degrees . For each angle, 300 projections of 0.5 s each were stored in 128 x 128 matrices for easy access. The acquired database was successful in the generation of static and dynamic sinograms for which the myocardial uptake and washout was governed by a compartment kinetic model., Conclusion: A method has been developed that allows creation of sinograms of physical phantoms with the capacity to control the number of noise realizations, the level of noise, the dynamics of uptake in the phantom compartments, and the acquisition parameters and acquisition modes.

Unlabelled: The goals of this investigation were to assess the accuracy of (18)F-fluorodihydrorotenone ((18)F-FDHR) as a new deposited myocardial flow tracer and to compare the results to those for (201)Tl., Methods: The kinetics of these flow tracers in 22 isolated, erythrocyte- and albumin-perfused rabbit hearts were evaluated over a flow range encountered in patients. The 2 flow tracers plus a vascular reference tracer ((131)I-albumin) were introduced as a bolus through a port just above the aortic cannula. Myocardial extraction, retention, washout, and uptake parameters were computed from the venous outflow curves with the multiple-indicator dilution technique and spectral analysis., Results: The mean +/- SD initial extraction fractions for (18)F-FDHR (0.85 +/- 0.07) and (201)Tl (0.87 +/- 0.05) were not significantly different, although the initial extraction fraction for (18)F-FDHR declined with flow (P < 0.0001), whereas the initial extraction fraction for (201)Tl did not. The washout of (201)Tl was faster (P < 0.001) and more affected by flow (P < 0.05) than was the washout of (18)F-FDHR. Except for the initial extraction fraction, (18)F-FDHR retention was higher (P < 0.001) and less affected by flow (P < 0.05) than was (201)Tl retention. Reflecting its superior retention, the net uptake of (18)F-FDHR was better correlated with flow than was that of (201)Tl at both 1 and 15 min after tracer introduction (P < 0.0001 for both comparisons)., Conclusion: The superior correlation of (18)F-FDHR uptake with flow indicates that it is a better flow tracer than (201)Tl in the isolated rabbit heart. Compared with the other currently available positron-emitting flow tracers ((82)Rb, (13)N-ammonia, and (15)O-water), (18)F-FDHR has the potential of providing excellent image resolution without the need for an on-site cyclotron.

The purpose of this study was to evaluate flow heterogeneity and impaired reflow during reperfusion after 60-min global no-flow ischemia in the isolated rabbit heart. Radiolabeled microspheres were used to measure relative flow in small left ventricular (LV) segments in five ischemia + reperfused hearts and in five nonischemic controls. Relative flow heterogeneity was expressed as relative dispersion (RD) and computed as standard deviation/mean. In postischemic vs. preischemic hearts, RD was increased for the whole LV (0.92 +/- 0.41 vs. 0.37 +/- 0.07, P < 0.05) as well as the subendocardium (Endo) and subepicardium considered separately (1.28 +/- 0.74 vs. 0.30 +/- 0.09 and 0.69 +/- 0.22 vs. 0.38 +/- 0.08; P < 0.05 for both comparisons, respectively) during early reperfusion. During late reperfusion, the increased RD for the whole LV and Endo remained significant (0.70 +/- 0.22 vs. 0.37 +/- 0.07 and 1.06 +/- 0.55 vs. 0.30 +/- 0.09; P < 0.05 for both comparisons, respectively). In addition to the increase in postischemic flow heterogeneity, there were some regions demonstrating severely impaired reflow, indicating that regional ischemia can persist despite restoration of normal global flow. Also, the relationship between regional and global flow was altered by the increased postischemic flow heterogeneity, substantially reducing the significance of measured global LV reflow. These observations emphasize the need to quantify regional flow during reperfusion after sustained no-flow ischemia in the isolated rabbit heart.

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.

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.

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.

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.

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 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.

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.

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.

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.