Abstract
Three-dimensional (3D) mode imaging is the current standard for PET/CT systems. Dynamic imaging for quantification of myocardial blood flow with short-lived tracers, such as 82Rb-chloride, requires accuracy to be maintained over a wide range of isotope activities and scanner counting rates. We proposed new performance standard measurements to characterize the dynamic range of PET systems for accurate quantitative imaging. Methods: 82Rb or 13N-ammonia (1,100–3,000 MBq) was injected into the heart wall insert of an anthropomorphic torso phantom. A decaying isotope scan was obtained over 5 half-lives on 9 different 3D PET/CT systems and 1 3D/2-dimensional PET-only system. Dynamic images (28 × 15 s) were reconstructed using iterative algorithms with all corrections enabled. Dynamic range was defined as the maximum activity in the myocardial wall with less than 10% bias, from which corresponding dead-time, counting rates, and/or injected activity limits were established for each scanner. Scatter correction residual bias was estimated as the maximum cavity blood–to–myocardium activity ratio. Image quality was assessed via the coefficient of variation measuring nonuniformity of the left ventricular myocardium activity distribution. Results: Maximum recommended injected activity/body weight, peak dead-time correction factor, counting rates, and residual scatter bias for accurate cardiac myocardial blood flow imaging were 3–14 MBq/kg, 1.5–4.0, 22–64 Mcps singles and 4–14 Mcps prompt coincidence counting rates, and 2%–10% on the investigated scanners. Nonuniformity of the myocardial activity distribution varied from 3% to 16%. Conclusion: Accurate dynamic imaging is possible on the 10 3D PET systems if the maximum injected MBq/kg values are respected to limit peak dead-time losses during the bolus first-pass transit.
PET imaging in 3-dimensional (3D) mode has become the standard for new whole-body scanners. The increased sensitivity allows for reduction of injected activity to the patient while maintaining excellent image quality; however, random and scattered photon counts are increased, requiring systems with high counting rate capability and accurate corrections for these physical effects. Current PET instrumentation and National Electrical Manufacturers Association (NEMA) performance evaluation methods (1) have been developed primarily to optimize whole-body oncology imaging with 18F-FDG. However, dynamic PET imaging for myocardial blood flow (MBF) quantification with short-lived tracers, such as 82Rb, 15O-water, or 13N-ammonia, requires high counting rates and correction accuracy to be maintained over a wide range of measured activities (2). An ideal PET system should allow for conventional relative myocardial perfusion imaging (MPI) of tracer retention without compromising accuracy of first-pass dynamic data (3). Routine MBF imaging is clinically feasible with the 76-s half-life generator-produced tracer 82Rb, resulting in accurate (4,5) and reproducible measurements (3,6–8), as validated against 13N-ammonia and 15O-water standards (9–12).
We proposed methods to evaluate the dynamic operating range of 3D PET systems for quantitative imaging of MBF. Patient imaging protocols were implemented and used to confirm the predicted operating range.
MATERIALS AND METHODS
Phantom Scans
Image Acquisition
82Rb decaying isotope scans were obtained over 5 half-lives using an anthropomorphic torso phantom (model ECT/TOR/P; Data Spectrum Corp.), approximating a small male upper torso (38 × 26 cm) (13) on 8 different 3D PET/CT systems and 1 3D/2-dimensional (2D) PET-only system. An 13N-ammonia scan was obtained on 1 other 3D PET/CT camera (Supplemental Table 1; supplemental materials are available at http://jnm.snmjournals.org). The phantom contained a myocardial heart cavity and wall insert (model ECT/CAR/I; Data Spectrum Corp.), lungs, spine, and liver chamber (Figs. 1A and 1B). Liver and body cavities were filled with water to mimic soft-tissue attenuation. The phantom was placed in the prone position in the PET field of view to facilitate infusion directly into the myocardial wall, and 1,500–3,000 MBq (40–80 mCi) of 82Rb or 1,100 MBq of 13N-ammonia were infused. 82Rb was infused either as a 30-s square-wave with saline push or as a 50 mL/min bolus (Supplemental Table 1). A list-mode PET acquisition was started immediately after completion of tracer infusion, simulating the localized activity and high counting rate observed during tracer first-pass transit through the heart. For attenuation correction, the PET scan was followed by a low-dose CT scan on the PET/CT systems or a 4-min transmission scan on the PET-only camera.
Image Reconstruction
Dynamic images (28 frames × 15 s) were reconstructed using vendor-supplied Fourier rebinning-filtered backprojection or iterative expectation-maximization algorithms (14), with an 8- or 12-mm Hann or Gaussian postprocessing filter and all corrections enabled for isotope decay, attenuation, scatter, randoms, prompt-gammas, detector efficiency, and dead-time, according to routine clinical practice at each institution. Most systems had explicit prompt γ-correction enabled during reconstruction (Supplemental Table 1); others used a 50-cm CT attenuation correction field of view to minimize the contribution of prompt γ-photons to the 3D coincidence background (15).
Quantitative Analysis
Reconstructed image time–activity curves were analyzed to determine the dynamic operating range for which quantitative accuracy was maintained. Total injected activity time–activity curves were measured using Inveon Research Workplace software (Siemens) (Figs. 1C and 1D). A spheric volume of interest (10-cm diameter) encapsulating the activity in the heart insert was drawn (Fig. 1C) from which total decay-corrected activity, Aheart(t) (MBq), was measured for all mid-frame scan times, t (min). From the time–activity curve, the average decay-corrected activity in the late time frames, where tracer uptake had reached a stable maximum, was determined as the true reference value, Aref (MBq). Activity bias in each time frame was then calculated as:Eq. 1To compare dynamic range among scanners, Activity bias (t) was plotted as a function of total activity in the heart volume of interest, Adecay(t), where:Eq. 2For 82Rb and 13N-ammonia, the isotope decay constants are λ = ln(2)/1.27 and ln(2)/10 (min), respectively. The time, Tmax, of the earliest frame with 10% or less activity bias was identified (Fig. 1D), and total heart activity at t = Tmax was calculated as:Eq. 3Dead-time correction factors (DTFs) and prompt coincidences and/or singles counting rates associated with the highest activity, Amax, were tabulated as available in the image headers. The maximum weight–based activity recommended for patient studies was estimated as Amax divided by 50 kg, the representative body weight of the torso phantom, determined according to the attenuating cross-sectional area, which is approximately 2.3 times larger than the NEMA scatter phantom previously shown to represent a 21.5-kg patient (16). A repeated scan was obtained on the Discovery 690, 600, and the Biograph PET/CT-16 systems to assess reproducibility of injected activity/body weight values.
Scatter correction residual bias was estimated as the left ventricular (LV) cavity blood–to–myocardium ratio by plotting scatter bias(t) as a function of Adecay(t), where:Eq. 4Ccavity(t) represents average activity concentration in the heart cavity (Bq/cm3), and Cmyo(t) is the average concentration in the myocardial wall (Bq/cm3). Residual bias is an indicator of uncorrected scatter in the LV cavity and is important to measure because accurate scatter correction is required for quantitative MBF measurements using an image-derived input function. To extract myocardium and LV cavity blood time–activity curves, our in-house FlowQuant software was used (3). The blood-to-myocardium ratio was determined by taking the median of the cavity, base, and atrium time–activity curves and then dividing by the myocardium average time–activity curve.
Image quality was assessed as nonuniformity of the myocardium activity distribution, using the coefficient of variation (COV) of the LV polar map:Eq. 5where SDmyo(t) is the SD of the activity concentration in the myocardial wall polar map. Images were also inspected visually for count rate–dependent pile-up artifacts.
Patient Scans
Patient Population
Recommended weight-based activity and DTF limits defined by the phantom scans were validated using 82Rb PET images from 20 patients acquired on the Discovery 690 and 600 and the Scintron 3D PET cameras (Supplemental Table 2). All patients were referred for a clinically indicated myocardial perfusion scan for coronary artery disease diagnosis or risk stratification. The institutional review board (or equivalent) at each of the participating centers approved this study, and all subjects signed a written informed consent form.
Image Acquisition
On the Discovery cameras, the 82Rb rest scan was followed by a dipyridamole stress scan, whereas regadenoson stress was used on the Scintron. Injected activity of 10 MBq/kg of body weight was prescribed for patients scanned on the Discovery systems, and 8 MBq/kg on the Scintron, according to local clinical practice for MPI. At rest and stress, 6-min list-mode acquisitions were started at the time of injection to capture the first-pass transit of the tracer as required for MBF quantification (Supplemental Fig. 1).
Quantitative Analysis
DTFs were tabulated for each time frame to identify the peak counting rates and dead-time losses. Global LV MBF values were computed automatically using FlowQuant, as described for phantom scans. Blood and LV myocardium time–activity curves were used as input to a 1-tissue-compartment model with a constant distribution volume to estimate MBF (4,6).
Statistical Analysis
Values are presented as mean ± SD. Where applicable, means were compared via the Student t test or 1-way ANOVA using SPSS Statistics 23 (IBM). A P value of less than 0.05 was considered statistically significant.
RESULTS
Phantom Scans
Figure 1D shows the bias in measured activity as a function of time (t) and total activity in the heart phantom insert, Adecay(t), for a single scanner. The Amax with less than 10% bias was 325 MBq. Assuming a representative phantom mass of 50 kg, the highest recommended patient-equivalent injected activity/body weight was estimated as 6.5 MBq/kg. At this activity (frame 12 at 2.75 min = Tmax), the peak prompt and singles counting rates and DTF were 4.1 and 29 Mcps and 2.0, respectively (Table 1).
Across all investigated scanners, the maximum recommended injected activity/body weight, peak DTF, and counting rate for accurate dynamic, quantitative cardiac MBF imaging varied between 3 and 14.4 MBq/kg, 1.5–4.0 DTF, and 22–64 Mcps singles and 4–14 Mcps prompt counting rate, respectively (Table 1). As expected, scanners using optimized detector crystals (higher atomic number, shorter decay time, higher light output (17)) and/or improved processing electronics were found to accommodate higher injected activity/body weight while remaining quantitatively accurate (Fig. 2). Peak DTF values within the accurate range corresponded typically with peak coincidence dead-times of 50% or less. Peak counting rates varied considerably between scanners, and intercomparison was not possible in all cases, depending on the camera-specific information available. Repeated scans were within 4% ± 9% of the originally tabulated injected activity/body weight values (Discovery 690, 13.1 MBq/kg [+6.8%]; Discovery 600, 5.7 MBq/kg [−6.4%]; and Biograph PET/CT-16, 6.9 MBq/kg [+11.0%]), demonstrating good reproducibility of the proposed methodology.
Residual scatter bias varied from 2% to 10% within the accurate operating range (Table 1). Highly variable uncorrected scatter was observed for all scanners in early time frames when counts tend to pile up toward the center of the detector blocks (Fig. 3A). Within the accurate operating range only, scatter bias stabilized at a relatively constant level (Fig. 3B). This bias was found to be slightly higher on the lutetium oxyorthosilicate detector–based systems (7.8 ± 2.0) versus the other scanners (2.9 ± 1.1; P < 0.05), suggesting that the scatter-correction methods implemented on these 4 lutetium oxyorthosilicate–based scanners may benefit from further optimization to improve accuracy.
All phantom images showed high contrast and low noise over the entire range of activity. Assessment of the LV myocardium polar map nonuniformity demonstrated that COV was highest (lowest image quality) in early frames (t = 0-Tmax, Cmyo(t)COV = 10.2% ± 4.9%) (Fig. 4A). COV values stabilized within the recommended operating range (t = Tmax-7 min, Cmyo(t)COV = 8.9% ± 3.4%; P = not significant vs. early frames) (Fig. 4B), corresponding with the trend observed for residual scatter bias. These results suggest that despite high dead-time losses in the early time frames, image quality is not compromised and is not a limiting factor for quantitative accuracy. Visual image inspection confirmed the absence of any obvious counting rate–dependent pile-up artifacts for all PET systems.
Patient Scans
Delivered activity was 10.3 ± 0.3 and 9.9 ± 2.0 MBq/kg for patients imaged on the Discovery 690 and 600 systems. Peak DTF values were 1.5 ± 0.1 and 2.1 ± 0.2 (corresponding to 33% and 50% coincidence dead-time, respectively), similar to the recommended maxima suggested by the phantom scans (Table 1). For patients scanned on the Scintron with 8 MBq/kg, peak DTF was 1.6 ± 0.2 (38% coincidence dead-time), slightly lower than the phantom maximum value recommended to remain within the accurate dynamic operating range.
DISCUSSION
This study established methods to evaluate the accurate dynamic operating range of 3D PET systems for quantitative cardiac imaging with 82Rb. Decaying isotope phantom scans were obtained over 5 half-lives to determine the optimal operating range, defined by the maximum injected activity/body weight, and corresponding maximum singles, prompt coincidence counting rates, or peak dead-time factors. Patient scans were obtained near the suggested limits on 3 representative scanners and confirmed validity of the phantom scan recommendations. Evaluation of the scatter correction bias confirmed the effectiveness of manufacturer-implemented scatter corrections in 3D mode. Finally, LV polar map nonuniformity, and the absence of counting rate–dependent pile-up artifacts, was found to be adequate for diagnostic evaluation.
The results suggest that the evaluated PET scanners should be able to perform accurate quantitative imaging despite differences in manufacturing technology, including scintillation detectors, detector block size, coincidence processing hardware, and prompt-γ-correction availability. The most important factor to consider for quantitative imaging in patients is that peak dead-time, singles, and/or prompt coincidence counting rates remain below the threshold values determined from the phantom scans to obtain accurate images and prevent a biased MBF estimation. This technique allows for prospective determination of image accuracy, as opposed to retrospective evaluation of detector block saturation or other performance metrics after acquisition is completed (18). It can also be performed retrospectively as long as counting rate and/or dead-time parameters are contained in the image header files generated by the scanners. As opposed to the NEMA counting rate performance standard that is designed for whole-body oncology imaging, the proposed method measures myocardial activity using a more realistic cardiac imaging geometry, and the residual scatter fraction and myocardial image uniformity are measured at the highest counting rates typically encountered during the bolus first pass.
Scan header information obtained from most of the scanners did not include complete counting rate and dead-time information, which would make it impossible to retrospectively verify that patient scans were obtained within the accurate operating range, and therefore that quantitative MBF values were absolutely reliable. It would be beneficial if all manufacturers recorded this information in the scan headers to allow consistent evaluation of quantitative imaging performance.
Clinical Implications
An optimal PET imaging system should allow for conventional MPI and absolute MBF imaging with a single injection of tracer. For PET cameras with adequate dynamic range and for which the injected activity/body weight limit is high enough to obtain diagnostic-quality perfusion images, both static MPI and dynamic MBF images can be obtained with a single injection protocol. In our previous multicenter study, 10 MBq/kg was the benchmark established for diagnostic-quality MPI using 3D PET systems (15). If this threshold was applied to the cameras in the current study, a single injection protocol could be recommended on the first 4 systems listed in Table 1; the others would require a dual-injection protocol to first obtain diagnostic-quality perfusion images using a high-dose injection and then accurate images for MBF assessment with a lower-dose scan. However, all systems showed maintained image quality in the early frames outside of the accurate range, as measured by COV, suggesting that the dynamic range may be extended with vendor improvements in dead-time correction accuracy, potentially allowing for a single-injection protocol on additional systems.
If a single-injection protocol is used on a system without sufficient dynamic range, MBF values would be inaccurate because peak blood and myocardium activities would be underestimated due to high dead-time losses at activity values exceeding the dynamic range of the camera. A patient example demonstrating this effect is shown in Figure 5. A 170-cm, 100-kg female patient was scanned at rest on the Discovery 600 PET/CT system with an injected dose of 10 MBq/kg of 82Rb, above the maximum limit for accurate quantification determined by the phantom scan, and a 5-fold-lower dose of 2 MBq/kg, for comparison. At 10 MBq/kg, the peak DTF and prompt counting rates were 2.5 and 6.8 Mcps, exceeding the recommended limits of 2.0 and 4.1 Mcps, suggesting that the camera was not operating in the accurate dynamic range for that portion of the scan. As a result, the early peak values of the blood and LV myocardium time–activity curves are underestimated compared with the low-dose time–activity curves (Fig. 5A). Underestimation of the area under the blood curve causes overestimation of the MBF values, as shown in the LV myocardium polar maps (Fig. 5B).
Patient Scan Variability
To estimate the maximum weight–based activity to use for patient studies, Amax was divided by the representative weight of the phantom, estimated as 50 kg. A scanner-dependent estimate based on DTF values obtained using a particular activity/body weight protocol may be more appropriate. The patient data required to perform this estimate were available for the Discovery 690, 600, and Scintron 3D. The peak DTF values from the patient scans vary according to body weight (Fig. 6A), therefore for each scanner the phantom DTF curve was plotted as a function of the weight-based dose (MBq/kg), using representative phantom weights corresponding with (passing through) the mean and maximum patient DTF values (Fig. 6B). Adjustment of the representative phantom weight improves prediction of the recommended injected dose for a range of patient sizes and also accounts for differences in absolute calibration between the PET scanner and rubidium elution system (used to measure the phantom and patient injected activities, respectively). The difference in calibration was relatively small (−6% to +20%) between the RUBY-FILL (Jubilant DraxImage Inc.) elution system and Discovery scanners, whereas the difference was much larger (−68%) between the CardioGen-82 (Bracco Diagnostics Inc.) infuser and Scintron scanner. On the basis of the combined results, slightly lower injected activities of 9, 7, and 6 MBq/kg for the Discovery 690, 600, and Scintron would be required to keep all patient DTF values within the accurate operating range. For other scanners, the recommended injected activity/body weight values should be confirmed or adjusted using similar methods.
Comparison with Previous Studies
In previous work by Tout et al. (18), the dynamic range of the Biograph mCT for simultaneous 82Rb MPI and MBF assessment was investigated in patients. They determined that a dose of 1,110 MBq (30 mCi) resulted in lower detector block saturation (1%) versus the manufacturer-recommended dose of 1,480 MBq (40 mCi) (15% saturation). All patients were injected with the same activity rather than the weight-based dosing proposed here. When the weight-based value of 14.4 MBq/kg determined in the present study for this scanner model and the mean weight of 87 kg from the population studied in Tout et al. were used, an average injected activity of 1,250 MBq would be recommended, similar to their suggested value of 1,110 MBq. As we have previously presented (19), administered activity can be adjusted for patient weight to compensate for the tracer distribution volume in the body and increased attenuation. Diagnostic image quality would likely be improved with higher doses in larger patients, whereas smaller patients would benefit from lower radiation dose with maintained image quality. Most importantly in smaller patients, a standard dose may exceed the scanner dynamic range during the bolus first pass and prevent accurate MBF quantification. The method of Tout et al. relies on verifying detector block saturation post hoc, which is not possible in real-time during patient scanning and is a more complex procedure than simple observation of the peak DTF or system counting rates during the patient scan.
Recently, Kolthammer et al. (20) investigated the dynamic range of the Ingenuity TF, the successor to the Gemini TF evaluated in the present study. A cylinder phantom was infused with 4 separate doses of 82Rb, ranging from 370 to 1,480 MBq, with a 10-min PET acquisition starting simultaneously with the infusion. Dynamic images were reconstructed into 15-s time frames. From this experiment, they determined that 82Rb imaging was accurate up to a peak singles counting rate of 65 Mcps at an injected activity of 925 MBq. For an average-sized patient of 70 kg, 925 MBq corresponds to an injected activity/body weight of approximately 13 MBq/kg. For the Gemini TF investigated in the present study, we obtained a maximum recommended activity of approximately 5 MBq/kg, suggesting that the Ingenuity TF may accommodate higher injected activities due to improved detector crystals (lutetium-yttrium oxyorthosilicate) and electronics. Because the singles counting rates were not stored in the Gemini TF header files, these values could not be compared with the Ingenuity TF scanner.
In another recent study, O’Doherty et al. (21) investigated the effect of scanner dead-time on MBF values obtained from kinetic modeling of 13N-ammonia dynamic images acquired on the Discovery 710 PET/CT scanner. They showed that global LV MBF values in 4 patients were 8.9% ± 0.6% higher when the LV blood-pool input function was corrected for high dead-time losses in the early frames using the percentage difference between measured versus true activity obtained from phantom studies. These preliminary results again demonstrate the importance of using appropriate camera-specific maximum injected activity levels to obtain accurate MBF assessments.
CONCLUSION
Dynamic imaging to obtain accurate quantitative MBF measurements with 82Rb appears feasible on the 10 3D PET systems evaluated when the recommended peak dead-time, maximum counting rates, and injected activity limits are respected. Patient scans confirmed the validity of the injected activity/body weight recommendations to achieve accurate and reliable quantitative images.
DISCLOSURE
This study was funded by Canadian Institute of Health Research grant MIS-100935 (Rb-ARMI). Jennifer M. Renaud and Robert A. deKemp receive royalties from FlowQuant® sales. Robert A. deKemp receives royalties from rubidium PET technology licenses. Rob S. Beanlands and Robert A. deKemp received grant funding from a government/industry research program (Ontario Research Fund/Industry Partners: GE Healthcare, Nordion, Lantheus Medical Imaging, and Jubilant DraxImage Inc.). Jennifer M. Renaud, Rob S. Beanlands, and Robert A. deKemp are consultants for Jubilant DraxImage Inc. Rob S. Beanlands is a consultant for Lantheus Medical Imaging and a career scientist, supported by the Heart and Stroke Foundation of Ontario, Vered Chair of Cardiology and University of Ottawa Tier 1 Chair in Cardiovascular Imaging Research. Timothy G. Turkington is a consultant for Data Spectrum Corporation and has received grant funding from GE Healthcare. No other potential conflict of interest relevant to this article was reported.
Acknowledgments
We thank the technologists and research staff at the participating sites for their invaluable efforts in acquiring the data discussed herein. We acknowledge Jubilant DraxImage Inc. for provision of the 82Rb generators, elution system, training, and support.
Footnotes
Published online Aug. 18, 2016.
- © 2017 by the Society of Nuclear Medicine and Molecular Imaging.
REFERENCES
- Received for publication March 18, 2016.
- Accepted for publication July 17, 2016.