|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Clinical Investigation |
1 The Weatherhead PET Center for Preventing and Reversing Atherosclerosis, University of Texas Medical School and Memorial Hermann Hospital, Houston, Texas; 2 Imaging Physics Department, M.D. Anderson Cancer Center, University of Texas, Houston, Texas; 3 The Weatherhead PET Center and Division of Cardiology, The Department of Medicine, University of Texas Medical School at Houston, Houston, Texas; 4 Cardiology Division, Department of Medicine, Northwestern Medical School, Chicago, Illinois; 5 Department of Medicine, University of Texas Medical School at Houston, Houston, Texas; and 6 The Weatherhead PET Center and Division of Cardiology, Department of Medicine, University of Texas Medical School at Houston, Houston, Texas
Correspondence: For correspondence or reprints contact: K. Lance Gould, MD, The Weatherhead PET Center, University of Texas Medical School, 6431 Fannin St., Room 4.256MSB, Houston, TX 77030. E-mail: K.Lance.Gould{at}uth.tmc.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: PET image registration heart attenuation artifact
| INTRODUCTION |
|---|
|
|
|---|
Although thoracic and cardiac PET/CT has proven feasible, misregistration of attenuation and emission images causes artifactual abnormalities on cardiac PET images that are false-positive defects (1624) or may artifactually worsen the severity of real perfusion defects. Misregistration of CT attenuation correction and PET emission data with associated artifactual PET defects is due to momentary helical CT "snapshots" at some point in the respiratory cycle of changing attenuating thoracicdiaphragmatic structures during breathing. The brief CT snapshots commonly do not match the actual average attenuation of the constantly changing thoracicdiaphragmatic structures over longer emission scans during normal breathing. This potentially major clinical limitation has not been definitively defined or resolved for current clinical protocols with commercially available cardiac PET/CT scanners and software. Thus, for PET/CT, the brief CT transmission data for attenuation correction are fundamentally different from earlier PET scanners using a rotating rod source for transmission data acquired slowly during normal breathing comparable to emission data acquisition.
Consequently, we examined the frequency, causes, consequences, and corrections of misregistration artifacts in current cardiac PET/CT protocols and software not previously reported, despite adverse effects on diagnostic accuracy.
| MATERIALS AND METHODS |
|---|
|
|
|---|
PET
Patients were instructed to fast for 4 h and abstain from caffeine, theophylline, and cigarettes for 24 h before the study. PET/CT was performed using a Discovery ST 16-slice PET/CT multislice bismuth germanate tomograph (GE Healthcare) in 2-dimensional mode with extended septa at a reconstructed in-plane resolution of 5.9-mm full width at half maximum (FWHM). Patients were positioned in the scanner using laser guides aligned to the base of the throat and confirmed by a CT scout scan. External body markers were used to check correct position throughout data acquisition. CT transmission scans were obtained as detailed below. Emission images were then obtained over 6 min after intravenous injection of 1,2951,850 MBq (3550 mCi) of generator-produced 82Rb and contained 2460 million total counts, of which 1230 million were true coincidence counts.
Immediately after completing the resting 82Rb scan, dipyridamole (0.142 mg/kg/min) was infused for 4 min. At 4 min after completion of the dipyridamole infusion, the same dose of 82Rb was given intravenously. Emission image acquisition was started at 70 s after the beginning of 82Rb infusion or 80 s after the beginning of infusion for patients with heart failure or heart rates below 55 beats per minute. After completion of stress emission imaging, a second helical CT transmission scan was obtained. For angina induced by dipyridamole, aminophylline (125 mg) was given intravenously.
The first 145 consecutive patients were imaged using helical CT attenuation correction of the PET data to determine the extent and causes of misregistration associated with artifactual PET defects as falsely positive studies. Helical CT transmission scans were obtained at 120 kV and 100 mA over 29 s during normal breathing in a subset of patients or over 4 s at end-expiratory breath-hold in another subset of patients. For the next 114 consecutive patients, in addition to the short end-expiratory helical CT, cine CT during normal breathing was added to the helical CT protocol. These subsets of patients served to compare quantitatively the extent of misregistration and associated artifactual defects between PET images reconstructed with helical CT attenuation correction data at end expiration and images reconstructed with cine CT averaged attenuation correction data during breathing in the same patient.
Imaging Protocols, Data Acquisition, and Attenuation Correction
After positioning in the gantry, patients had a CT scout scan using 120 kV and 10 mA in anteriorposterior format to bracket the heart. A cine CT scan for average CT attenuation correction was acquired during normal breathing at rest; this was followed by a helical CT scan at end expiration. Perfusion imaging of 82Rb was then obtained at rest and after dipyridamole stress; this was followed by another cine CT scan during normal breathing for stress average attenuation correction and another helical CT scan. The total time for a complete reststress study averaged 40 ± 6 min (mean ± 1.0 SD).
Cine-mode data acquisition is one of the 3 optional data acquisitions on the Discovery PET/CT scanner that include axial, helical, and cine modes. The helical CT acquisition used 120 kVp, x-ray collimation of 16 by 1.25 mm, gantry rotation of 0.5 s, and helical pitch of 1.75:1. The radiation exposure was 5.7 mGy for the helical CT scan.
The cine CT scan used 120 kVp, x-ray collimation of 8 by 2.5 mm, gantry rotation cycle of 0.5 s, cine duration of 10 s, and 14-cm coverage. Each reconstruction used 360° of CT data over 0.5 s. The cine data acquisition lasted for 78 s, of which 70 s was for 7 cine CT acquisitions of 10 s each and 8 s was for 6 table translations of 2 cm each between 2 cine CT acquisitions. Tube current for cine CT was 10 mA for patients weighing <100 kg, 15 mA for 100130 kg, and 20 mA for >130 kg. The radiation dose and effective dose equivalent for cine CT were 10 mGy and 2.36 mSv, respectively, for tube current of 10 mA. For each cine CT scan, there are 1,232 CT images reconstructed 2 min after the cine CT acquisition and 22 images per slice location during normal respiration before averaging into an average CT image over several breathing cycles. The processing time for an averaged cine CT attenuation scan was <1 min.
Interpolation of the average CT attenuation data was performed to make the average CT images align at the locations of the PET slices for attenuation correction of the PET data and for evaluating registration between the average CT image and the PET images corrected with the average CT attenuation data. In view of the inherent blurring introduced by the average CT but not by the helical CT data, we applied 4- and 10-mm FWHM gaussian filtering on the average CT and helical CT data, respectively, before CT attenuation correction of the PET data.
Software was also developed to shift the averaged cine CT data in horizontal and vertical directions in the transaxial plane and in the superiorinferior or z-axis direction of different planes to align the CT heart borders visually with the PET data to achieve good coregistration. The average cine CT slices shifted to match the PET slices are saved for repeat reconstruction of the PET images using the shifted average attenuation correction data. The cine CT PET has been previously described for cancer imaging (21), here adapted to cardiac imaging with the shift software.
Analysis of PET Images
Images were reconstructed using filtered backprojection with a Butterworth filter having a cutoff of 0.55, roll-off of 10, and pixel size of 3.27 x 3.27 mm. Fusion images were made by superimposing attenuation-corrected emission and CT transmission scans in horizontal, coronal, and sagittal views. Misregistration of CT and PET images was quantified in millimeters for comparison with quantitative size and severity of regional perfusion defects on PET images corresponding to the regions of misregistration on the PET/CT fusion display.
After attenuation correction, reconstructed PET emission images in DICOM (Digital Imaging and Communications in Medicine) format were exported to a Sun Ultra 60 workstation of a mPower-HZL PET scanner (Positron Corp.) for reorientation into long- and short-axis tomographic and topographic 3-dimensional (3D) displays using previously described automated quantitative software (110) for comparison with our earlier report on misregistration artifacts in PET using a rotating rod transmission source for attenuation correction (20). Quantitative data from 82Rb and 18F phantoms imaged on the GE Healthcare scanner are quantitatively and accurately exported to the Positron scanner and software for automated quantitative analysis and display of PET images.
A 3D restructuring algorithm generates true short- and long-axis views from reconstructed PET transaxial cardiac images, perpendicular to and parallel to the long axis of the left ventricle. From the tomographic data, circumferential profiles are used to reconstruct 3D topographic views of the left ventricle showing relative regional activity distribution divided into lateral, inferior, septal, and anterior quadrant views of the 3D topographic display corresponding to the coronary arteries as previously described (410,20).
Mean activity in each quadrant is normalized to the maximum 2% of pixels in the whole heart dataset. Regions of each quadrant are identified having values outside 2.0 SD and 2.5 SD (97.5% confidence intervals [CI]) outside normal values of 50 healthy volunteers with no risk factors by complete medical history (no hypercholesterolemia, hypertension, diabetes, smoking, obesity, drug or alcohol abuse, other systemic diseases, or family history of CAD). Percentages of circumferential profile units outside 2.0 SD and 2.5 SD were calculated automatically for each quadrant and the whole heart.
Severity of a perfusion defect is quantified as the average relative uptake in a defined quadrant or the lowest quadrant average relative activitythat is, the average relative activity for the quadrant having the lowest average activity of anterior, septal, lateral, and inferior quadrants for each subject, expressed as percentage of the highest 2% of activity in the image dataset. Size of perfusion defects is quantified as percentage of the whole cardiac image outside 2.0 SD or 2.5 SD or 97.5% CI of healthy control subjects in each quadrant and whole heart. Combined size and severity of perfusion defects is defined as percentage of the whole cardiac image with relative activity of <60% of maximum activity (100%), which is 3.0 SD below mean maximum activity of healthy control subjects.
Misregistration Artifacts
On fusion images of helical CT PET and separately of cine CT PET, electronic calipers were used to measure the distance in millimeters between the cardiac border on the PET image and the cardiac border on the CT image in the horizontal plane of maximal misregistration on transaxial and coronal views and in the vertical or z-axis of the coronal view. In addition, on the coronal view, the distance was measured from the zero reference position of helical and cine CT scans to the dome of the right diaphragm to determine systematically the effects of diaphragmatic position on misregistration artifacts.
On PET perfusion images, artifactual abnormalities due to misregistration were defined as criteria (i) plus (ii), (iii), or (iv) as follows:
Therefore, to be counted as having an abnormality due to PET/CT misregistration, every patient with abnormal helical CT PET at rest or stress had to have a stress perfusion study with no significant regional defect after correct PET/CT coregistration. Any patient with a defect that persisted on stress PET/CT images or after correct coregistration that was outside 2 SD of healthy subjects was classified as having a "real" defect, not an artifact, even if the defect improved after correct registration. Therefore, our observed prevalence of abnormalities due to PET/CT misregistration is a conservative underestimate.
For PET scans meeting the above criteria for misregistration artifacts, the severity, size, and combined sizeseverity of perfusion abnormalities in the same quadrant as the CT PET misregistration were objectively quantified by automated software. In addition to automated measures of defect size and severity, the perfusion abnormalities due to CT PET misregistration were visually graded according to the relative color scale bar by 2 independent readers, who were unaware of clinical information, as none (white or red), mild (yellow), moderate (green), or severe (blue) in the same quadrant for which misregistration was measured.
Statistical Methods
Statistical analyses were performed using SPSS version 11.5 (SPSS Inc.) for multivariate linear regression analysis, for multivariate logistic regression analysis, and for significance of differences in continuous variables among groups determined by a paired or unpaired 2-tailed t test with data reported as mean ± 1 SD. A Pearson
2 analysis was performed for discrete variables. A 2-tailed P value < 0.05 was considered statistically significant.
| RESULTS |
|---|
|
|
|---|
|
2 analysis in Table 1. With misregistration of >6 mm on the transaxial fusion images of helical CT and PET emission images, the PET perfusion scan showed artifactual defects that were significant diagnostic errors in 103 of 259 (40%) patients (P = 0.0000). Transaxial misregistration of >6 mm frequently caused artifactual defects with the proportion of moderate-to-severe artifactual defects progressively increasing for misregistration of 610 mm and for >10 mm.
|
|
Logistic multivariate regression analysis was done with the 3 misregistration measurements in millimeters (transaxial x-axis, coronal horizontal x-axis, and coronal vertical z-axis) as independent variables, with the dependent variable being separately either any artifact or only moderate-to-severe artifacts. Transaxial misregistration (x-axis) predicted misregistration artifactual defects so completely (coefficient of regression analysis, B = 0.94, P = 0.0000) that the other measurements added no more predictive power.
Similarly, multivariate regression analysis was done with the 3 misregistration measurements as independent variables and the dependent variable being separately the severity (mean quadrant activity), size (percentage outside 2.5 SD), or combined sizeseverity (percentage of the left ventricle outside 60% of maximum activity) of the abnormalities due to PET/CT misregistration. The coronal horizontal (x-axis) misregistration measurement was such a strong predictor of the 3 quantitative measures of defect severity (B = 2.23 for combined sizeseverity, B = 0.96 for severity, B = 0.59 for percentage outside 2.5 SD; for all, P = 0.0000) that other misregistration measurements added no predictive power. These results indicate that the 3 misregistration measurements are closely related such that the x-axis misregistration in either the transaxial or coronal view is a complete predictor of artifactual defects.
For patients with both helical and cine CT PET (Table 3), the artifactual defects due to misregistration of helical CT attenuation and PET emission data were markedly reduced by using the averaged cine CT data acquired during normal breathing that is comparable to the PET emission acquisition during normal breathing. Compared with all patients, this quantitative improvement was even greater for the subset of patients with artifactual defects due to PET/CT misregistration, excluding those with real defects as defined above. Transaxial PET/CT misregistration in millimeters was significantly less with cine CT PET compared with helical CT PET.
|
Most misregistration defects were anterior or lateral (76%) associated with corresponding anterior or lateral misregistration on the CT-PET fusion images. The remainder had inferior (8%) or apical (16%) abnormalities due to inferior or apical misregistration of helical CT PET or different diaphragmatic position during the helical CT compared with the diaphragm position on the cine CT average attenuation scan that resolved the artifactual defect. Both the artifactual defects and their associated misregistration on helical CT-PET fusion images resolved on cine CT PET in most cases.
However, of the 114 patients with both helical CT and cine CT PET, 22 (19%) had persisting artifactual defects with persisting PET/CT misregistration on the cine CT-PET fusion images. For these patients, the cine CT data were shifted to achieve proper visual coregistration with the PET images that were then reconstructed again with the shifted cine CT attenuation data. With this final coregistration, the artifactual defect disappeared and the images normalized.
Figure 2A illustrates severe anterior, apical, lateral, and basal inferior defects on the helical CT PET (top row) that remained but were less severe on cine CT PET (middle row). The anterior, apical, and lateral defects resolved only on the shifted cine PET/CT (bottom row) that achieved proper coregistration. However, the basal inferior defect persisted with correct registration, indicating right coronary stenosis, confirmed by coronary arteriography with no stenosis of the left coronary artery. Figures 2B2D show the corresponding misregistration on the helical CT PET and cine CT-PET fusion images but no misregistration on the shifted cine CT-PET fusion images.
|
|
|
This diaphragmheart mismatch on the stress fusion image caused attenuation overcorrection inferiorly and an associated relative anterior defect on the stress PET. In addition, the stomach air bubble is higher on the expiratory helical CT than the inferior heart border on the emission scan, thereby creating a small strip of inferior attenuation undercorrection despite the rest of the inferior half of the heart being overcorrected.
To resolve the question of severe worsening before cine CT was available for this patient, the PET study was repeated on the next day with the same restdipyridamole protocol, same medications, and same dose of 82Rb but on the Positron PET scanner with a rotating rod attenuation source. This repeat reststress study showed no misregistration and no significant stress-induced perfusion defects and was similar to the earlier PET scans, displayed for comparison with the helical CT PET of the prior day in Figure 3C. The absence of significant perfusion defects after correct coregistration indicates artifactual defects on the earlier PET/CT.
We also examined the frequency of artifactual defects due to misregistration in patients with slow 29-s helical CT scans during breathing compared with different patients with fast 4-s helical CT scans during end-expiratory breath-holding. Slow helical PET/CT during breathing had significantly fewer artifactual defects due to misregistration than fast helical PET/CT at end-expiratory breath-hold (Table 5). In addition, the artifactual defects with slow helical CT PET during breathing were of a different type associated with layering of the CT data and multiple diaphragm and heart border fragments, which caused jagged notched layers of CT attenuation correction data. Correspondingly, the artifactual defects consisted of jagged, irregular linear notches or bites out of the resting CT and PET perfusion images as illustrated in Figure 4.
|
|
| DISCUSSION |
|---|
|
|
|---|
Although slow helical CT scans during breathing have somewhat less-frequent artifactual defects than fast helical CT scans at end expiration, the layered attenuation data of the long helical CT during breathing distorts the PET images in a way that cannot be corrected by shifting the CT data to coregister with the PET data. Similarly, shifting the short helical CT data at end expiration to coregister with the PET data is also unsatisfactory because CT attenuation is greater during expiration and, therefore, causes attenuation overcorrection compared with the average attenuating structures during breathing as measured by cine CT during breathing.
Compared with our earlier study of PET attenuation misregistration using a rotating rod attenuation source during breathing, artifactual defects in helical CT PET are more frequent, more severe, and more complex. Preventing or correcting these artifactual defects in PET/CT requires cine CT and software for shifting the CT data to visually coregister with the PET data when needed.
Some early thoracic or cardiac studies on small numbers of patients reported satisfactory CT attenuation correction of PET emission data (16,19,22,23), whereas others reported problems of misregistration (17,18,24). A recent study on 28 patients reported substantial artifactual perfusion defects due to misregistration of PET/CT data in 21% of subjects (24). Visual realignment of PET/CT data to achieve coregistration eliminated the artifactual perfusion defects. However, our results indicate that a helical CT scan at one moment during the respiratory cycle may not fit the average attenuating structures during breathing of emission imaging even after visually matching the lateral heart borders.
Moreover, variable attenuation correction associated with respiratory motion causes marked heterogeneity in measurements of myocardial radionuclide uptake (22). Therefore, measurements of absolute myocardial perfusion in mL/min/g by PET/CT are appropriate only if the primary radionuclide update data are correctly obtained using cine or shifted cine CT attenuation correction that matches the average attenuation structures of emission data during breathing. We believe that the relatively poor specificity reported for PET/CT (13) is due to misregistration artifacts as no protocol for correcting misregistration of PET/CT data was described in that study.
One criticism of our study may be the lack of coronary arteriography to document the absence of stenosis. However, invasive coronary arteriography is not justified for a defect in a region of PET/CT misregistration that disappears after repeat reconstruction using correctly coregistrered CT attenuation data, particularly in view of well-documented inaccuracies of visually interpreted coronary arteriograms.
Our software implementation is applicable to all Discovery PET/CT scanners without any additional hardware. GE Healthcare provides an expensive hardware upgrade of their PET/CT console to the Dimension console that achieves some of the functions reported here. However, the Dimension console averages combined helical and cine CT data for attenuation correction, provides no average cine CT images for shifting to register with PET data, and is available only on the latest Discovery STE PET/CT. Other manufacturers do not have the interface for cine CT. Alternatively, a non-CT pure PET scanner with a 3-rod rotating attenuation source would acquire attenuation data very rapidly without the misregistration and relatively high radiation dose of PET/CT, thereby providing the basis for an inexpensive fast PET scanner for cardiac imaging.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| References |
|---|
|
|
|---|
Related articles in JNM:
This article has been cited by other articles:
![]() |
A. K. Buck, S. Nekolla, S. Ziegler, A. Beer, B. J. Krause, K. Herrmann, K. Scheidhauer, H.-J. Wester, E. J. Rummeny, M. Schwaiger, et al. SPECT/CT J. Nucl. Med., August 1, 2008; 49(8): 1305 - 1319. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. W. Townsend Dual-Modality Imaging: Combining Anatomy and Function J. Nucl. Med., June 1, 2008; 49(6): 938 - 955. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Gould, T. Pan, C. Loghin, N. P. Johnson, and S. Sdringola Reducing Radiation Dose in Rest-Stress Cardiac PET/CT by Single Poststress Cine CT for Attenuation Correction: Quantitative Validation J. Nucl. Med., May 1, 2008; 49(5): 738 - 745. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Gould Not All Randomized Trials Are Equal J. Am. Coll. Cardiol., November 13, 2007; 50(20): 2013 - 2015. [Full Text] [PDF] |
||||
![]() |
K. L. Gould and T. Pan Reply: Attenuation Correction for Stress and Rest PET 82Rb Myocardial Perfusion Images J. Nucl. Med., November 1, 2007; 48(11): 1913 - 1914. [Full Text] [PDF] |
||||
![]() |
R. L. Eisner and R. E. Patterson Attenuation Correction for Stress and Rest PET 82Rb Myocardial Perfusion Images J. Nucl. Med., November 1, 2007; 48(11): 1912 - 1913. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |