Abstract
The aim of this work was to compare the quantitative accuracy of iteratively reconstructed cardiac 18F-FDG PET with that of filtered backprojection for both 2-dimensional (2D) and 3-dimensional (3D) acquisitions and to establish an optimal procedure for imaging myocardial viability with 18F-FDG PET. Methods: Eight patients underwent dynamic cardiac 18F-FDG PET using an interleaved 2D/3D scan protocol, enabling comparison of 2D and 3D acquisitions within the same patient and study. A 10-min transmission scan was followed by a 10-min, 25-frame dynamic 3D scan and then by a series of 10 alternating 5-min 3D and 2D scans. Images were reconstructed with filtered backprojection (FBP) or attenuation-weighted ordered-subsets expectation maximization (OSEM), combined with Fourier rebinning (FORE) for 3D acquisitions, applying all usual corrections. Regions of interest (ROIs) were drawn in the myocardium, left ventricle, and ascending aorta, with the last 2 being used to define image-derived input functions (IDIFs). Patlak graphical analysis was used to compare net 18F-FDG uptake in the myocardium, calculated from either 2D or 3D data, after reconstruction with FBP or OSEM. Either IDIFs or arterial sampling was used as the input function. The same analysis was performed on parametric images. Results: A good correlation (r2 > 0.99) was found between net 18F-FDG uptake values for a myocardium ROI determined using each acquisition and reconstruction method and blood-sampling input functions. A similar result was found for parametric images. The ascending aorta was the best choice for IDIF definition. Conclusion: Good correlation and no bias of net 18F-FDG uptake in relation to that based on FBP images, combined with less image noise, make 3D acquisition with FORE plus attenuation-weighted OSEM reconstruction the preferred choice for cardiac 18F-FDG PET studies.
The use of iterative reconstruction methods, such as ordered-subsets expectation maximization (OSEM), rather than filtered backprojection (FBP) has been investigated in several oncologic PET studies. Apart from assessing lesion detectability (1–4), some of these studies have also investigated the effects of iterative reconstruction on quantitative accuracy (5–7), not only in oncology but also in quantitative brain studies. In contrast, little is known about the use of iterative reconstruction in cardiac PET (8). The same is true for the comparison of PET studies with (2-dimensional [2D]) and without (3-dimensional [3D]) interplane septa. Removal of the interplane septa in 3D PET significantly increases not only sensitivity but also the detection of scattered radiation, random coincidences, and detector dead time—all of which require appropriate correction methods. The quantitative accuracy of 3D brain studies has been well documented, and the use of 3D PET in lesion detection in oncology has also been studied (9–14). The use of 3D PET for quantitative cardiac studies has been investigated less extensively (15–17), although the scatter correction method used on many clinical scanners (18,19) was developed with quantitative cardiac 3D imaging in mind. Furthermore, many commercially available PET or PET/CT systems are no longer equipped with septa and are thus capable only of 3D acquisitions.
Myocardial viability studies that use 18F-FDG with the aim of quantifying glucose metabolism require the measurement of plasma radioactivity. Instead of arterial blood sampling, an image-derived input function (IDIF) can be defined using the radioactivity concentration in the left ventricle, left atrium, ascending aorta, or descending aorta (20–25). In a previous study, the ascending aorta was found to be the best choice for defining IDIFs (24), using 2D acquisitions and reconstruction by filtered backprojection.
The purpose of the present study was to assess the quantitative accuracy of both iterative reconstruction and 3D acquisition for cardiac PET studies. To this end, an interleaved 2D/3D scan protocol was used to enable comparison of 2D and 3D imaging within the same patient and study, and data were reconstructed with both FBP and OSEM. In addition, both arterial blood sampling and IDIFs were used as plasma input curves.
MATERIALS AND METHODS
PET Scanner
All measurements were performed using an ECAT Exact HR+ scanner (CTI/Siemens). Standard settings for data compression (maximum ring difference, mash, and span) were used, together with standard corrections for normalization, random coincidences, and dead time (26,27). Data were corrected for attenuation using a 10-min windowed transmission scan acquired with 3 rotating 68Ge/68Ga rod sources. Scatter correction was performed using the methods included in the scanner software. This correction involves a stationary convolution subtraction method in 2D mode (28) and a simulation-based method in 3D mode (29). The scaling factor in the 3D scatter-correction algorithm was modified to ensure accurate cross-calibration of 2D and 3D acquisitions (30).
Image Reconstruction
The 3D data were converted to 2D sinograms using Fourier rebinning (FORE) (31), the standard procedure in the ECAT 7.2 software, which does not include full 3D reconstruction algorithms. Images were reconstructed using FBP or attenuation-weighted OSEM (32) (6 iterations, 16 subsets) as implemented in the ECAT 7.2 software, applying a Hanning filter of 5 mm in full width at half maximum to ensure resolution matching between reconstruction methods (5–7), resulting in a transverse image resolution of around 7 mm in the center of the FOV.
Patients and Scan Procedure
Eight patients (all male; mean age, 65 y; range, 56–76 y), admitted for assessment of myocardial viability, were included. After a 10-min transmission scan, approximately 370 MBq of 18F-FDG dissolved in 5 mL of saline were administered as a bolus and a 10-min 3D dynamic emission scan was acquired: 6 × 5 s, 6 × 10 s, 3 × 20 s, 5 × 30 s, and 5 × 60 s. This dynamic scan was then followed by a series of ten 5-min scans alternating between 3D and 2D modes, separated by the time necessary to retract or insert the septa.
Blood Sampling
The radioactivity concentration of arterial blood was measured online during the course of the scan (33). At set times, additional discrete blood samples were collected and the arterial line was flushed to avoid clotting. The arterial blood curve was corrected for delay relative to the measured total trues counting rates, by matching times at which 50% of the peak counting rate was reached in the up-slopes of trues and blood curves. The whole-blood and plasma radioactivity concentrations of the discrete samples were measured in a well counter, cross-calibrated against the PET scanner. The plasma input function was then derived by multiplying the measured whole-blood curve with a linear fit to the plasma-to-blood ratios obtained from the discrete samples.
ROI Analysis
A region of interest (ROI) was drawn around the contour of the myocardium (6 image planes, total volume of approximately 40 cm3) in transaxial images reconstructed using either FBP or OSEM and acquired in either 2D or 3D mode. 18F-FDG influx rate constants (Ki) for the whole myocardium were calculated with graphical analysis according to the Patlak method (34), using only the last 5 frames in either 2D or 3D mode for determination of the Patlak slope. Plasma radioactivity concentration determined by arterial blood sampling was used as the input function. Because net 18F-FDG uptake, described by Ki, is directly proportional to glucose metabolic rate by multiplication with the plasma glucose concentration divided by the lumped constant (35,36), the Ki results of the present study are valid for glucose metabolic rates as well.
Parametric Image Analysis
All images were resampled to 1-cm3 voxels. Ki images were calculated using the Patlak method, again for each combination of acquisition and reconstruction methods and with blood-sampling data as input function. Scatterplots were generated, comparing measured Ki in all voxels within a 10 × 10 × 15 cm box around the heart for the various combinations of acquisition and reconstruction methods. Correlation coefficients and slopes of linear regressions were calculated. The data were also analyzed by Bland-Altman plots (37), assessing the differences between voxel values for various combinations of reconstruction and acquisition methods.
IDIFs
ROIs were drawn in the left ventricle (diameter, 15 mm; length, 6 image planes) and ascending aorta (diameter, 13 mm; length, 6 image planes) for definition of IDIFs. The same ROIs were used for images obtained using all combinations of acquisition and reconstruction methods. For definition of the 3D IDIFs, all frames acquired in 3D mode were used. For the 2D IDIFs, the first 10 min were assumed to be equal to the 3D IDIFs, after which the 2D frames were added. These IDIFs were then interpolated to the same time points as used for the continuous sampling curves by linear interpolation, and the plasma radioactivity concentration was calculated using the same plasma-to-blood data as used to correct the sampling data. Mean values of the IDIFs in the last 5 frames in either 2D or 3D mode were calculated and compared with the mean value of the sampling plasma curve at the corresponding midpoint times. The same comparison (3D only) was made for the areas under the first 90 s of the IDIF and the blood-sampling plasma curve. Subsequently, the IDIFs were used as the input function for calculation of myocardial Ki.
RESULTS
ROI Analysis
Figure 1 shows scatterplots of Ki based on myocardium ROIs in images acquired in either 2D or 3D mode and reconstructed using FBP or OSEM. A clear linear correlation without bias exists between the Ki values for all methods, with the square of the Pearson correlation coefficient being at least 0.99 for all scatterplots.
Parametric Image Analysis
18F-FDG Ki parametric images of a typical patient are shown in Figure 2. Figure 3 shows scatterplots for all combinations of reconstruction and acquisition methods for another patient. Table 1 summarizes the results, giving mean values of the scatterplot slopes and intercepts, the Pearson correlation coefficients, and the absolute differences between the results obtained using the different methods. The results shown are for all voxels above a threshold at the mean of all voxel values within the image (a threshold chosen to include mainly voxels inside the myocardium).
IDIFs
The mean ratios of the IDIF relative to the blood-sampling curve during the last 10 frames and the area under the first 90s of the curve are given in Table 2. The area under the peak of the IDIF is significantly lower than the area under the peak of the blood-sampling curve, especially for images reconstructed by OSEM. The 2D IDIFs generally overestimate the input functions at later times, especially the left ventricle IDIFs. For images acquired in 3D mode, the left ventricle IDIFs overestimate the input function whereas the ascending aorta IDIFs underestimate the input function.
Figure 4 shows scatterplots of Ki calculations based on IDIFs versus those based on blood-sampling input functions. With 8 data points and a significance of 0.05, an r2 value below 0.45 is required to reject the hypothesis of a linear correlation with zero intercept, which is the case only for 2D FBP left ventricle IDIFs. Table 3 shows the ratios of Ki based on IDIFs and sampling data. Figure 4 and Table 3 clearly show that ascending aorta IDIFs resemble blood-sampling input functions better than do left ventricle IDIFs and that there is a larger spread in left ventricle data than in ascending aorta data. In addition, no significant differences were found between Ki values calculated using FBP and OSEM images with similar IDIFs and acquisition modes, and 3D images led to significantly higher Ki values than did 2D images with similar IDIFs.
DISCUSSION
Using an interleaved 2D/3D scanning protocol, the effects of both 2D and 3D acquisition protocols and analytic and iterative reconstruction techniques on quantitative cardiac 18F-FDG PET studies were investigated. First, the effects of acquisition and reconstruction methods on Ki, calculated with blood-sampling data as the input function, were studied using either ROI or parametric analyses. Next, the quantitative effects of using an IDIF were assessed.
Good correlations were observed between Ki obtained from various combinations of acquisition and reconstruction methods, using continuous blood sampling as the input function and using either ROI analysis or parametric analysis. Neither 3D acquisition nor OSEM reconstruction led to a bias in Ki. No significant quantitative differences between cardiac studies reconstructed with OSEM and FBP and acquired in 2D or 3D mode were found.
3D IDIFs underestimate the blood radioactivity concentration during the first 90 s after administration, with a larger underestimation for images reconstructed by OSEM. During the late part of the input curve, all IDIFs except those using the ascending aorta in a 3D image overestimate the input function, with the 2D values being significantly higher than the corresponding 3D values. This overestimation is largest for 2D left ventricle IDIFs, and FBP and OSEM did not significantly differ in this case. The apparent overestimation of left ventricle IDIFs is likely caused by spillover of radioactivity from the myocardium. Because this effect is the same in 2D and 3D modes, the difference in overestimation between 2D and 3D left ventricle IDIFs may be caused by the different scatter-correction methods. Underestimation of the 3D ascending aorta IDIFs and the early part of the left ventricle IDIFs may be caused by the limited resolution of the scanner, resulting in a recovery of 80%–90% for a structure the size of the ascending aorta (38). In addition, the dead-time correction of the scanner may be inaccurate for a situation in which a large amount of radioactivity is concentrated in a small volume, such as the left ventricle and ascending aorta during the start of this scan, and especially in 3D acquisition mode because of the higher counting rates. The underestimation in the early part of left ventricle IDIFs, as well as the overestimation during the late part of left ventricle IDIFs and better agreement between ascending aorta IDIFs and online input function, was also seen for 2D FBP images in an earlier publication (24).
The ratios between Ki for the whole myocardium calculated using ascending aorta IDIFs and blood sampling in 2D mode are consistent with those previously published (24), leading to the conclusion that the ascending aorta is the preferred choice for IDIF definition in 2D mode. This also seems to be the case in 3D mode, in which ascending aorta IDIFs lead to a minor overestimation and left ventricle IDIFs lead to a significant underestimation of Ki, compared with the blood-sampling input function.
Partial-volume correction of the ascending aorta IDIFs may be possible, but knowledge of the exact shape and size of the ascending aorta is needed for accurate correction. In the late part of the plasma curve, mainly determining the slope of the Patlak plot, the limited recovery in the myocardium may cancel out the limited recovery in the ascending aorta and thus not influence the Patlak slope. Partial-volume correction of the IDIF may lead to a smaller Patlak slope, but this effect is probably small. The limited recovery could also be improved by using less smoothing in OSEM reconstructions of 3D scans, possible because of the better statistics of 3D data. It should be noted that the 2D IDIFs contained a peak measured in 3D mode, which, as for the late part of the IDIFs, might have been higher had it been acquired in 2D mode. This higher peak, in turn, would have led to slightly smaller Patlak slopes for the 2D IDIFs. The linear correlation shown in Figure 4 leads to the conclusion that ascending aorta IDIFs can be used for Patlak analysis of cardiac 18F-FDG scans without further partial-volume correction of the ascending aorta data.
The use of IDIFs enables modeling and quantitative analysis of tracer uptake without the need for arterial blood sampling. Earlier studies showed that the ascending aorta is the preferred choice for IDIF definition in 2D mode (24). This result was confirmed by the current study. Ki values based on IDIFs were significantly higher for 3D acquisitions than for 2D acquisitions. Further investigations—of, for example, the effect of the different scatter-correction algorithms in 2D and 3D PET studies—are required. Ascending aorta IDIFs are preferred to left ventricle IDIFs because they have a smaller deviation from the online plasma input function and are more reproducible, as shown by the lower SDs in Table 3. For nonlinear analysis using compartmental modeling, continuous blood sampling probably remains the method of choice, but for Patlak linear analysis of cardiac 18F-FDG PET scans, ascending aorta IDIFs may be a sufficiently accurate alternative.
The use of 3D acquisition results in higher-quality images with the same amount of radioactivity because of the increase in scanner sensitivity. Alternatively, compared with 2D acquisition, 3D acquisition enables the administration of less radioactivity without loss of image quality. This aspect is especially helpful in research studies on healthy volunteers, in whom the relatively high absorbed dose of 18F prohibits repeated injections. Furthermore, PET/CT scanners from most manufacturers are not equipped with septa, thus dictating 3D scanning. Therefore, it would be advantageous if cardiac PET could be performed in 3D mode without loss of quantitative accuracy.
This study showed a good correlation between Ki determined using 2D and 3D images when blood-sampling input functions were used. In addition, no significant effects of reconstruction method on quantitative results were observed. When blood sampling is used for determination of the input function, 3D acquisition and OSEM reconstruction are quantitatively as accurate as 2D acquisition and FBP reconstruction, provided the scaling factor in the 3D scatter-correction algorithm is modified to ensure accurate cross-calibration of 2D and 3D acquisitions (30).
CONCLUSION
The use of 3D acquisition, combined with blood sampling for determination of the input function, is a good alternative for quantification of cardiac 18F-FDG PET studies. The ascending aorta is the preferred choice for definition of IDIFs for 2D and 3D acquisitions reconstructed with FBP as well as OSEM, resulting in Ki values that show a linear correlation and no bias compared with those obtained using a blood-sampling input function. Iterative reconstruction using attenuation-weighted OSEM is a quantitatively accurate alternative to FBP reconstruction for cardiac 18F-FDG PET. Furthermore, the use of 3D acquisition and OSEM reconstruction may enable less smoothing of the images, leading to a better spatial resolution, which is of importance, for example, in measurements of transmural differences in glucose metabolism. Good correlation and no bias of Ki in relation to Ki values based on FBP images, combined with less image noise, make 3D acquisition with attenuation-weighted OSEM reconstruction the preferred choice for cardiac 18F-FDG PET studies.
Acknowledgments
The authors thank the staff at the PET Centre of VU University Medical Centre for assistance during scanning and Dr. Paul Knaapen for assistance with ROI definition.
Footnotes
Received Apr. 22, 2004; revision accepted Jul. 15, 2004.
For correspondence or reprints contact: Mark Lubberink, PhD, Department of Nuclear Medicine and PET Research, VU University Medical Centre, P.O. Box 7057, 1007 MB Amsterdam, The Netherlands.
E-mail: mark.lubberink{at}vumc.nl