Abstract
Hybrid PET scanners offer the possibility of obtaining myocardial viability information from coincidence imaging of the positron emitter 18F-FDG and perfusion measurements from a single-photon tracer—potentially simultaneously. This new approach is less costly and more readily available than dedicated PET and offers potential for improved FDG resolution and sensitivity compared with SPECT with 511-keV collimators. Simultaneous imaging of the coincidence and single-photon events offers the further advantages of automatic image registration and reduced imaging time. However, the feasibility of simultaneous coincidence/single-photon imaging or even immediately sequential imaging is unknown. In this study, the potential of using standard low-energy high-resolution (LEHR) collimators with hybrid PET to obtain coincidence and SPECT data was assessed. Methods: Phantom and human studies were performed to investigate the effect of LEHR collimators on FDG coincidence imaging with a hybrid PET system, the effect of the presence of 99mTc during FDG coincidence imaging with LEHR collimators, and the effect of the presence of FDG during 99mTc SPECT imaging. Results: FDG images were somewhat degraded (a measure of myocardial nonuniformity increased 10%) with LEHR collimators. With 148 MBq (4 mCi) 99mTc present during FDG imaging of a phantom, image quality was maintained and the number of detected coincidences changed by <5%. With 99mTc/18F whole-body ratios of 7:1, crosstalk from 18F photons accounted for the majority of counts in the 99mTc SPECT images and resulted in severe artifacts. The artifacts were decreased with a simple crosstalk correction scheme but remained problematic. Conclusion: 99mTc/18F ratios of at least 9:1 and state-of-the-art reconstruction and crosstalk correction are likely to be required to perform immediately sequential coincidence/single-photon imaging of the myocardium with clinically useful results. Additional challenges remain before simultaneous imaging of coincidence events and single photons can be realized in practice.
Determination of myocardial perfusion and viability is important for diagnosing and ascertaining the prognosis of coronary artery disease. Relative myocardial perfusion can be measured with a variety of 99mTc-labeled or 201Tl single-photon emitters or with positron emitters. For the noninvasive assessment of myocardial viability, PET imaging of 18F-FDG is the current gold standard. FDG SPECT with 511-keV collimators has been shown to provide a clinically useful assessment of viability, though generally not equivalent to that provided by FDG PET (1,2). Compared with SPECT with 511-keV collimators, the new hybrid PET camera systems may offer benefits for cardiac viability imaging; coincidence detection provides better resolution and lacks the mammoth lead collimators that block most of the incident photons. Recently, some publications relevant to cardiac imaging with hybrid PET systems have appeared (3–5), although the lack of built-in attenuation correction methods has made cardiac imaging with hybrid PET impractical. This situation has changed lately, with several manufacturers now offering solutions.
Simultaneous measurement of myocardial perfusion and viability has several advantages over separate acquisitions, including intrinsic registration of the datasets, decreased motion artifacts, and increased throughput. For the case of SPECT with 511-keV collimators, energy discrimination permits the simultaneous imaging of FDG and a 99mTc-tagged pharmaceutical, although both agents are imaged with the 511-keV collimator. Such simultaneous viability and perfusion measurements using SPECT have been reported (6–8). A good review of FDG imaging with SPECT, including dual-isotope FDG/99mTc-sestamibi acquisitions, has been given by Sandler et al. (3,9).
Using coincidence-capable gamma cameras in fully 3-dimensional mode with low-energy high-resolution (LEHR) collimators mounted may enable simultaneous perfusion/viability imaging with improved resolution and a lowered dose of positron emitter, compared with 511-keV collimated imaging (10). Because the LEHR collimators are easily penetrated by the 511-keV photons, coincidence imaging is similar to imaging with “open-frame” collimators. However, several problems arise from the hybrid PET approach to simultaneous perfusion/viability imaging. One problem is that axial slat collimators cannot be used in conjunction with LEHR collimators. Axial slat collimators are narrow slats of lead septa running perpendicular to the z-axis of the scanner. The lead slats block interslice coincidences. Although some type of hybrid LEHR/lead slat collimator could be designed, the slats would block a significant portion of the low-energy data. Another potential problem is the impact of the presence of 99mTc on 18F imaging in terms of increasing dead time. And the final serious problem is the crosstalk from 511-keV photons scattering or depositing a fraction of their energy in the crystal and masquerading as 99mTc photons.
The goal of this study was to evaluate the feasibility of using a hybrid PET system to image FDG with LEHR collimators in place and 99mTc present and immediately thereafter acquiring 99mTc SPECT data. If such an approach to cardiac imaging is feasible, truly simultaneous acquisitions are then possible, although additional hardware may be necessary to minimize imaging time.
MATERIALS AND METHODS
The Effect of LEHR Collimators on Coincidence Imaging
Phantom Study.
A cardiac thorax phantom (Data Spectrum Corp., Hillsborough, NC) was filled with water, and a transmission scan using 2 scanning 133Ba point sources (365 keV) was obtained. 18F was added to give 10 MBq (0.27 mCi) in the myocardial insert (93.5 kBq/mL), 9.6 MBq (0.26 mCi) 18F in the liver (7.8 kBq/mL), and 37 MBq (1.0 mCi) 18F in the background (4.2 kBq/mL). The phantom was imaged with an IRIX hybrid PET scanner (Marconi Medical Systems, Cleveland, OH) with 2 coincidence detectors and used 60 gantry angles over 180°, with 20 s at each angle. Axial slat collimators were used. The axial slat collimators were then replaced with standard LEHR collimators (2.7 cm thick, 0.12-cm hole size), and 1 h later the coincidence acquisition was repeated with the phantom in the same position on the bed. The second acquisition used 60 gantry angles, with 30 s at each angle. At each gantry stop, list-mode coincidence data were acquired over a range of approximately ±20° (this is a function of the radius of the camera and the size of the detectors, as well as the presence of axial slat collimation). Both list-mode datasets were thus binned into 90 angles with single-slice rebinning (SSRB). SSRB is the simplest and most standard method for handling events in which 1 photon is detected in 1 transaxial slice and the corresponding photon is detected in another slice. SSRB simply places the event in the transaxial slice that bisects the 2 detected locations. The binned data were then corrected for “rotational weights” (11) and decay, multiplicatively corrected for attenuation effects, reconstructed with 8 iterations of ordered-subset expectation maximization (10 subsets), and postfiltered with a 3-dimensional Metz filter. No random or dead-time corrections were performed. The uniformity of the left ventricular wall distribution was assessed by calculating the SD of values from a maximum count circumferential profile program. Evaluation of contrast was made from manually drawn regions on the blood pool (BP) and the tissue of the left ventricle (LV) in a midventricular slice. The BP region comprised 4 voxels, and the LV region comprised 10 voxels. Reconstructed voxels for all phantom and patient studies were 0.47 cm on each side. Contrast was then calculated as (LV − BP)/(LV + BP).
A 22Na point source in air was also imaged with the axial slat collimators and with the LEHR collimators to measure resolution. The projections were binned and reconstructed in the same manner as the phantom data.
Patient Study.
A patient with a prior myocardial infarction was given an oral glucose load, and a transmission scan was acquired. Three external markers (IZI Medical Products Corp., Baltimore, MD) filled with 99mTc were affixed to the patient, and a brief SPECT scan was obtained to help register the scans. The LEHR collimators used with the transmission scan were then changed to axial slat collimators, and FDG was added to the markers. The patient then received a 255-MBq (6.9 mCi) injection of FDG. Imaging was started 55 min later (181 MBq [4.9 mCi] FDG at scan start) with axial slat collimators. The collimators were then changed to the LEHR type, and the patient was imaged a second time, starting 133 min after injection (111 MBq [3.0 mCi] FDG at scan start). Both scans acquired data at 30 gantry angles for 60 s at each angle. The list-mode data were processed and reconstructed as described for the phantom data. The emission datasets were registered with the attenuation map by shifting to visually align the markers using an overlay display. Short-axis and polar maps were created from the attenuation-uncorrected and -corrected reconstructed transaxial slices. The same measure of LV-to-BP contrast as described for the phantom study was used for the patient data. The BP region comprised 9 voxels from midventricular slices, and the LV region comprised 20 voxels from a normally appearing region of heart tissue.
The Effect of 99mTc on Coincidence Imaging
A total of 156 MBq (4.2 mCi) 99mTc was added to the 18F-loaded phantom to give a 7:1 99mTc/18F overall ratio at the time of imaging. Specifically, 10.3 MBq (0.28 mCi) 99mTc was added to the LV; a 3-cm3 defect with 18F present was left without 99mTc. Seventy megabecquerels (1.9 mCi) were added to the liver, and 75.5 MBq (2 mCi) 99mTc were added to the background. With LEHR collimators in place, the phantom was imaged over 60 gantry angles for 30 s at each angle. Sensitivity (total number of photopeak–photopeak coincidences) and image quality were compared with and without the 99mTc present.
The Effect of FDG on 99mTc SPECT Imaging
The same phantom loaded with 18F and 99mTc was imaged with SPECT with a 15% energy window centered at 140 keV. At the start of imaging, the 99mTc/18F ratio was 9:1 (3.5:1 ratio in the heart). A 3.2:1 ratio in the heart has been reported to be the minimum expected when a dose of 925 MBq (25 mCi) 99mTc-sestamibi and 370 MBq (10 mCi) FDG is used (3,6). Because much less 18F is sufficient with hybrid PET systems, this is a worst-case scenario. The acquisition was over 60 gantry angles for 20 s at each angle. Ten hours later, a 99mTc-alone scan with a 99mTc/18F overall ratio of 138:1 was obtained with 60 gantry angles and 60 s at each angle. All SPECT phantom datasets were reconstructed with ordered-subset expectation maximization using attenuation compensation and depth-dependent point response compensation.
Crosstalk Correction for 99mTc SPECT Imaging
A multiple–energy window method was evaluated for crosstalk correction. The method assumes that the crosstalk from 18F into the 99mTc photopeak window (p140(s,t,φ), where s and t denote the coordinates in a 2-dimensional projection image at angle φ) can be approximated as a weighted average of smoothed data from energy windows adjacent to the 140-keV photopeak:
Eq. 1
where ⊗ indicates convolution, g() is a 2-dimensional gaussian function g(σ)=
e(− s̃2−t̃2)/2σ2, and ku and kl are scaling factors. A similar method has been used for estimating the crosstalk in the 201Tl energy window when both 201Tl and 99mTc are present (12).
The parameters ku, σu, kl, and σl were determined from a separate thorax phantom experiment. In this experiment, the phantom had only 18F present but was imaged with SPECT with energy windows at 124 keV (7%), 140 keV (15%), and 162 keV (7%) to give plower, p140, and pupper. The 4 unknown parameters were estimated by minimizing:
Eq. 2
where
Eq. 3
A nonlinear optimization routine (MATLAB; The MathWorks, Inc., Natick, MA) constrained to give parameters in the range 0.1–10.0 was used.
The crosstalk correction was tested by imaging a healthy volunteer. A resting 99mTc-sestamibi study was performed first (68 angles over 204°, 17 min total time). SPECT data were acquired 50 min after the injection of 955 MBq (25.8 mCi) 99mTc-sestamibi. Then, 159 MBq (4.3 mCi) FDG were injected, and the SPECT scan was repeated 84 min later so that at scan start there was a 99mTc/18F whole-body ratio of 7.3:1. The energy windows were the same as those used for the 18F phantom study. The SPECT data were reconstructed after the scaled smoothed crosstalk estimates (from Eq. 2 and the parameters from the phantom study) were subtracted from the projection data. Reconstruction was also performed by adding the crosstalk estimate to the projections at each iteration. The crosstalk correction schemes were evaluated by a noise measure (percentage root-mean-square error in a region of interest) and by a contrast measure ([LV − BP]/[LV + BP]). The BP region comprised 12 voxels over 2 midventricular slices, and the LV region comprised 8 voxels from a normally appearing region of heart tissue. The reconstructions were filtered with a 3-dimensional Metz filter for display.
RESULTS
Effect of LEHR Collimators on Coincidence Imaging
Phantom Study.
Figure 1 compares 18F phantom data imaged with axial slat collimators and with LEHR collimators. After scaling for decay and acquisition time, 36% more events were acquired with the LEHR collimators than with the axial slat collimators. The measure of uniformity of the activity distribution in the heart wall indicated that the wall was slightly less uniform when LEHR collimators were used (0.10 vs. 0.11) and much less uniform when attenuation correction was not performed (0.22). Contrast was better with axial slat collimators (0.98—the BP region was close to zero) than with LEHR collimators (0.80). The LEHR data did appear with distortions, which are addressed in the Discussion.
18F-only short-axis images of cardiac thorax phantom obtained with axial slat collimators (uniformity, 0.10) (A), LEHR collimators (uniformity, 0.11) (B), and axial slat collimators but without attenuation correction (uniformity, 0.22) (C). Short-axis slices appear somewhat square in (B); this anisotropic resolution effect is detailed in Discussion.
The in-plane resolution of the scanner in coincidence mode was 5.8 mm full width at half maximum both with axial slat collimators and with LEHR collimators. The axial resolution was 9.3 mm at the center of the field of view when axial slat collimators were used. With LEHR collimators, the axial resolution increased to 13.4 mm. Away from the center, axial resolution degraded with either collimator, but a much greater impact was seen for LEHR collimators.
Patient Study.
The results from FDG imaging of the patient with axial slat and LEHR collimators are shown in Figure 2. Contrast was similar for the LEHR collimator images (0.45) and for the axial slat collimator images (0.47). Reduced contrast was expected from increased scatter in the LEHR images, but this increase might have been offset by additional clearing of the FDG from the BP and by FDG accumulation in tissues during the time between the 2 scans. Some extraneous activity was also apparent in the LEHR images.
Patient data shown by polar maps with attenuation correction (A), polar maps without attenuation correction (B), and short-axis slices without attenuation correction (C). For all images, left side is with axial slat collimators and right side is with LEHR collimators.
Effect of 99mTc on Coincidence Imaging
With the 7:1 99mTc/18F overall ratio used in the phantom, little change in coincidence imaging was noted. After scaling for scan time and decay, the total number of coincidenceevents recorded changed <5%, from 11.7 million without 99mTc present to 12.2 million with 99mTc present. Image quality was comparable between the 2 image sets, as seen in Figure 3.
Four noncontiguous short-axis slices of phantom. Top row shows 18F-only images obtained with LEHR collimators. Middle row shows images with 99mTc also present (99mTc/18F overall ratio, 7:1). Bottom row shows difference images.
Effect of FDG on SPECT Imaging
Figure 4 shows the SPECT reconstruction of the thorax phantom with 99mTc and 18F present. Severe artifacts are obvious because of crosstalk from the 18F. Approximately half the single photons detected in the 99mTc photopeak window were from 18F crosstalk. The 99mTc-only image is also shown. The true defect can be seen in the 99mTc-only image, and artifactual defects appear from 18F crosstalk.
Phantom study for crosstalk levels. Top row shows images with 99mTc and 18F present (99mTc/18F overall ratio, 9:1) and LEHR collimators. Bottom row shows 99mTc-only images (99mTc/18F overall ratio, 138:1). Arrows are pointing to true 99mTc defect.
The 18F-only phantom study to estimate the crosstalk correction parameters is shown in Figure 5. The windows adjacent to the 99mTc photopeak appear similar to the photopeak data. Subtraction of the estimated crosstalk data in the photopeak gives a noisy residual centered on zero (Fig. 5B). The volunteer study shows that the transaxial slices are improved visually with crosstalk correction (Fig. 6), although the correction fails to remove the artifactual inferolateral defect. With respect to noise and contrast measures, Table 1 shows that contrast is improved with crosstalk correction but at the cost of at least some increase in noise.
Phantom study for crosstalk compensation. (A) Reverse grayscale images. Top left shows 99mTc window. Top right shows 162-keV window. Bottom left shows smooth scatter estimate (from Eq. 1). Bottom right shows difference or residual between 99mTc window and smooth scatter estimate. (B) Horizontal profile through row 27 of images in A (162-keV window profile not shown). Scat est. = smooth scatter estimate.
Reconstructions from SPECT imaging of volunteer with only sestamibi present (A), with sestamibi and FDG present and crosstalk compensation (B), with sestamibi and FDG present and subtraction of crosstalk estimate (C), and with sestamibi and FDG present and inclusion of crosstalk estimate in iterative reconstruction algorithm (D). Even with correction, artifactual inferolateral defect remains.
Comparison of Reconstruction Methods
DISCUSSION
When investigating the feasibility of imaging coincidence events and single photons together using LEHR collimators mounted on a hybrid PET system, one must consider 3 major issues. The first regards the use of standard LEHR collimators for FDG coincidence imaging. Significant degradations to the images were expected because, as in 3-dimensional PET, a much larger scatter fraction is obtained than when slat collimators are used. Watson et al. (13) have reported that, for thorax imaging with 3-dimensional PET, >50% of detected coincidences have undergone scattering in the object, compared with 10%–15% in 2-dimensional imaging. In the studies performed here, contrast was indeed reduced for the thorax phantom. However, the patient study revealed similar contrast with axial slat and LEHR collimators, possibly because of increased physiologic contrast at the later imaging time of the LEHR scan.
Also important to consider is the large axial field of view of hybrid PET systems. When the thorax phantom with only FDG was imaged in coincidence mode with LEHR collimators, significant image distortions were evident when SSRB was used; axially oblique acceptance angles were not limited by collimation or by software. The distortions produce short-axis slices that appear somewhat square, likely because of the poor axial resolution inherent with SSRB. Thus, a cylinder with the long axis perpendicular to the z-axis of the scanner would appear somewhat square because the resolution is not isotropic. That is, the cylinder walls parallel to the z-axis will be blurred more than will the walls perpendicular to the z-axis. With axial slat collimators or with PET systems having a small axial extent, the oblique angles are likely small enough for the effect to be negligible. Others have shown squared-off (as in Figs. 1 and 3) phantom results with hybrid PET (5), though no comment was made. As a test of this anisotropic resolution hypothesis, the list-mode coincidence data were reconstructed with a fully 3-dimensional algorithm (3-dimensional tilt angle binning (14,15)) that provides more uniform spatial resolution than does SSRB. The myocardial short-axis slices became less square, as seen in Figure 7. Although the LEHR data did result in degraded images, even the SSRB images appeared to retain much of the clinical value (Fig. 2).
Comparison of reconstruction using 3-dimensional tilt angle binning and SSRB for 18F data acquired with LEHR collimators. Left image was obtained with axial slat collimators; middle image, with LEHR collimators and SSRB; and right image, with LEHR collimators and 3-dimensional tilt angle binning.
Although seemingly not an issue for the patient study shown here, activity outside the field of view may be a significant problem with LEHR collimators. Placing extra shields of lead on the ends or covering portions of the patient with curved lead slabs may improve the data by blocking these photons. This problem is found with 3-dimensional PET systems as well (16,17).
The second issue that was investigated was the effect of the presence of 99mTc while imaging FDG. This issue was found not to be critical when a threshold detector was used to discard events less than approximately 200 keV and when relatively low counting rates were used. Serious performance trade-offs are to be expected with other scenarios.
The third issue was to determine the degree of crosstalk into the 99mTc energy window and to test a simple correction scheme. Compared with dual-isotope imaging with 511-keV collimators, the situation is dramatically different with LEHR collimators. The scatter contribution has little structure because the collimator has little effect on the 511-keV photons. One group has found encouraging results after compensating for crosstalk (18,19). In 1 study (19), the group used a triple–energy window scatter subtraction technique with 18F/99mTc in a brain phantom with a hybrid PET system. Their method is equivalent to using Equation 1 without any smoothing and letting kl = ku = 1.5 and using 5% wide energy windows on either side of the 140-keV photopeak.
On the basis of our finding that approximately half the single photons detected in the 99mTc photopeak window were from 18F crosstalk, we determined that the 3.2:1 99mTc/18F ratio in the heart reported as useful for dual-isotope 511-keV SPECT is not high enough to perform dual-isotope coincidence/single-photon imaging. In addition, 18F not only adversely affects the imaging of 99mTc by crosstalk in the 99mTc photopeak, but the large number of 18F photons causing scintillations not in the 99mTc photopeak may affect the dead time of the camera so that fewer 99mTc photons are detected. A vital question, which was not addressed in this study, is how little FDG can be used without degrading clinical FDG viability images. A clinical report in which 19 women received simultaneous injections of 99mTc-sestamibi and FDG at stress (740 MBq [20 mCi] 99mTc and 74 MBq [2 mCi] 18F initially) and were imaged 50 min later with PET and then 4–6 h later with SPECT shows that minimizing the crosstalk from 18F may make such imaging possible (20,21). It should be noted that a standard SPECT camera with a thin crystal and therefore reduced crosstalk was used for imaging of 99mTc-sestamibi (20).
The crosstalk compensation method used here improved the images, but artifacts remained. As with the phantom images in Figure 4, using a similar 99mTc/18F ratio, artifactual defects result from the crosstalk of FDG. An alternative to Equation 1, using instead a very smooth version of the 511-keV coincidence projections (after binning with SSRB), was investigated both by itself and in conjunction with Equation 1. These alternative formulations did not match the measured p140 data as well. More sophisticated crosstalk compensation methods and the presence of less FDG should provide better results.
The advantage of sequential 99mTc/18F imaging with LEHR collimators is primarily that the data are likely to be registered, which can be important when assessing perfusion and glucose uptake matches and mismatches. As well, only a single transmission scan is needed to correct both sets of data. Throughput also increases relative to a sequential axial slat/LEHR collimator protocol because the collimators do not need to be changed. Imaging may not be possible with axial slat collimators and large amounts of 99mTc present because the hybrid PET systems have a poor counting-rate performance and the high number of 99mTc single photons may cause significant dead time. Finally, successful sequential imaging implies that truly simultaneous imaging will be possible.
CONCLUSION
If feasible, imaging 18F and 99mTc together with coincidence/single-photon imaging is an exciting way to extend the usefulness of the new hybrid PET systems. The results imply that with high 99mTc/18F whole-body ratios (>9:1), low FDG doses, and state-of-the-art reconstruction and crosstalk correction, performing immediately sequential coincidence/single-photon imaging with clinically useful results will likely be practical. Simultaneous coincidence/SPECT imaging may also be possible, although the counting-rate capability of the camera when both coincidence and single photons must be processed is of critical importance.
Acknowledgments
This study was supported by a grant from the Whitaker Foundation. Portions of this study were presented at the annual meeting of the Society of Nuclear Medicine, Los Angeles, CA, 1999, and at the midwinter meeting of the Society of Nuclear Medicine, New Orleans, LA, 2000.
Footnotes
Received Aug. 7, 2000; revision accepted Feb. 14, 2001.
For correspondence or reprints contact: Edward V.R. Di Bella, PhD, MIRL/Radiology, 729 Arapeen Dr., Salt Lake City, UT 84108.