|
|
||||||||
BASIC SCIENCE INVESTIGATIONS |
Department of Nuclear Medicine, Seoul National University College of Medicine, Seoul; Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul; and Department of Physics, Ewha Womans University, Seoul, Korea
| ABSTRACT |
|---|
|
|
|---|
Key Words: H215O PET factor analysis regional myocardial blood flow
| INTRODUCTION |
|---|
|
|
|---|
Since its first introduction in nuclear medicine, factor analysis has been widely used to extract pure physiologic timeactivity curves (TACs) and their factor images from dynamic nuclear images (59). The purpose of factor analysis is to separate a few underlying but unobservable quantities, or factors, using the covariance relationship among the observed data. By singular value decomposition, the covariancevariance matrix of the observed data can be divided into the eigenvalue and eigenvector matrices. Thus, the observed data can be represented with a few sets of eigenvalues and eigenvectors that compose the greatest proportion of the total eigenvalues and eigenvectors.
Hermansen et al. (9) acquired myocardial factor images from H215O myocardial PET using a fixed-flow model and applied the myocardial factor images to define tissue ROIs without C15O blood-pool PET. However, these authors employed a modeling approach that used the lung TAC to obtain the left and right ventricular TACs, with consideration of the time delay between the lungs and the ventricles, and the computation time was long.
In this study, we measured rMBF using two-step factor analysis with dynamic H215O PET images without an additional C15O blood-pool PET scan. The purpose of dividing the procedure into two steps was to obtain the pure regional-tissue TACs by factor analysis. In the first step of the factor analysis, we obtained the input function and the spatial information for each component (myocardium and left and right ventricles) so that we could determine the myocardial regions (septal, inferior, lateral, and anterior). In the second step of the factor analysis, by drawing an ROI on a quarter of the myocardium, we obtained regional-tissue TACs. After validation of input functionderived factor analysis by that derived from arterial sampling, rMBF measured by factor analysis was validated by microsphere studies performed simultaneously with H215O PET.
| MATERIALS AND METHODS |
|---|
|
|
|---|
A transmission scan was acquired using triple 67Ge rod sources. After a bolus injection of 518666 MBq H215O, dynamic H215O PET images were acquired at rest for 4 min (5 s x 12, 10 s x 9, and 30 s x 3). Adenosine was infused at a continuous velocity (0.14 mg/kg/min) for 7 min. Three minutes after the start of the adenosine infusion, dynamic H215O PET images were acquired for 4 min while the adenosine was infused.
Transaxial images were reconstructed using a filtered backprojection method with a Shepp-Logan filter having a cutoff frequency of 0.3 cycles per pixel for 128 x 128 x 47 matrices and a pixel size of 2.5 x 2.5 x 3.3 mm. The dynamic images were reoriented to the short axis using static images already reoriented to 128 x 128 x 64 matrices and a pixel size of 2.5 x 2.5 x 2.5 mm. Twenty-four reoriented slices that included myocardium were selected and summed to become six slices to increase the signal-to-noise ratio. The resulting images were six short-axis slices with a per-slice thickness of 1 cm, which was same as the width of myocardium cut after microsphere studies.
Microsphere Studies
Three kinds of microspheres (15 µm; Dupont, Billerica, MA)those labeled with 85Sr, those labeled with 46Sc, and those labeled with 113Snwere used for the rest and stress study. The number of microspheres was approximately 2.8 x 106 per unit volume for 85Sr and 46Sc and 2.4 x 106 for 113Sn. Thirty-seven megabecquerels of microspheres were injected into the left ventricle through a catheter that had already been inserted into the left ventricular cavity through the femoral artery. The microspheres were injected simultaneously with injection of H215O through a forelimb vein.
After the dogs were killed, the heart was extracted and ventricular blood was washed gently with normal saline. The myocardium was cut along the long axis at 1-cm intervals. Six slices were obtained. Every slice was cut into four (anterior, lateral, septal, and posterior) regions, and then each region was cut into pieces and counted (counts per gram). By comparing the tissue counts with the blood sample counts, blood flow per unit mass (gram) of tissue was calculated (10).
Factor Analysis
The reoriented images were masked (32 x 32) to include the cardiac region using a midventricular slice. The same mask was applied to all six slices. Only the first 2 min (18 frames) of the H215O dynamic data were chosen for factor analysis. The masked images (32 x 32 x 18 x 6) were summed to be 8 x 8 x 18 x 6 matrices to reduce statistical fluctuation. We used the factor analysis algorithm based on the methods of Di Paola et al. (11) and Barber (12). In brief, the normalized TAC from each pixel, called a dixel, was assumed as a linear combination of the physiologic TACs (myocardium, left and right ventricle, and the other regions). After principal component analysis with the normalized dixels, a few eigenvectors were obtained and a new affine space was defined from these vectors. Then, using the positive constraint, this vector space was rotated until all factors had positive values, because all factors should be the pure physiologic data.
The procedure of our factor analysis was as follows. First, all six slices were used to extract factor images for the left and right ventricles and the myocardial tissues. The left ventricular TAC was obtained as an input function. Second, ROIs were drawn on the myocardial factor images to include the whole left ventricular cavity and a quarter (anterior, lateral, septal, and posterior) of a slice of myocardial tissue. The ROIs were projected onto the dynamic image, and then the second factor analysis was performed. Because the extracted TACs by factor analysis were in normalized units, the scaling factor to recover the true scale of the TACs was calculated using factor images (5,6). The average over a threshold (70%) of the maximum value of a factor image was used as the scaling factor. Thus, the scaling factor for input function was calculated with the left ventricular factor image from the first step, and the scaling factor for tissue TAC was calculated with the tissue factor image from the second step. Multiplying the normalized TACs by the scaling factor restored the true scale of the input function and tissue TACs. The whole procedure is shown as a schematic flowchart in Figure 1.
|
Calculation of rMBF by Factor Analysis and ROI Method
rMBF was calculated using the following single-compartment model (4).
![]() |
is partition coefficient. Flow, tissue recovery coefficient, and spillover fraction were estimated using a nonlinear least squares fitting method. The partition coefficient was set to 0.91 (3). rMBF was estimated with the regional-tissue TACs by factor analysis and the conventional ROI method as follows. In the second step of the factor analysis, ROIs were drawn over four regions (anterior, posterior, septal, and lateral) on the short-axis myocardial factor images including the left ventricle to extract left ventricular TAC simultaneously with tissue TACs. Tissue TACs were obtained for these four quarters of a myocardial slice and were used to estimate rMBF.
In the ROI method, tissue ROIs over the same four regions of a myocardial factor image were drawn narrowly to minimize spillover and partial-volume effects, and these ROIs were projected onto the reoriented dynamic H215O PET images. For an input function, an ROI was drawn on the left ventricle of the midventricular level of the myocardial factor image and projected onto the reoriented dynamic H215O PET image.
rMBFs calculated by factor analysis and the ROI method were compared with rMBF obtained by microsphere studies. All calculations, factor analyses, and rMBF estimations were performed using MATLAB (The MathWorks, Inc., Natick, MA) on an AlphaStation 600 workstation (Digital Equipment Corp., Maynard, MA).
Data Analysis
Input functions derived from factor analysis and arterial sampling were compared with each other by calculating the relative error of AUC ([AUCfactor analysis - AUCsampling]/AUCsampling). The correlation between rMBF by factor analysis, rMBF by the ROI method, and rMBF by microsphere studies was analyzed.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
The main finding of this study is that rMBF can be calculated directly from dynamic H215O PET images through two-step factor analysis without the aid of C15O blood-pool scans. Usually, H215O PET has required an additional C15O blood-pool scan to determine ROIs because of the poor quality of dynamic H215O images. This study showed that C15O scanning is no longer needed for myocardial H215O PET. Therefore, the additional cost and radiation dose from C15O scanning can be reduced by virtue of factor analysis, and the method can easily be applied to patients undergoing H215O PET for rMBF estimation even in PET centers without instruments for gas studies.
In the first step, we used all the slices as input data for the factor analysis; thus, three factor images (right and left ventricles and myocardium) for all slices and input function were obtained. The relative error of AUC between the input functions obtained by factor analysis and those obtained by arterial sampling was approximately 24%, mainly because of fluctuations in the tails of the curves. The shape of both curves, however, was nearly the same except for the expected delay and dispersion of radial artery sampling. The slight difference in the tails of both curves may have been caused by counting errors in the arterial sampling. Because the shape of the initial part of the input function is mainly what affects rMBF, fluctuations in the tails of the curves are not a problem in flow estimation.
The second step of the factor analysis enabled extraction of pure regional-myocardium TACs by factor analysis without careful drawing of subtraction-defined ROIs on the myocardium (13). In the conventional ROI method, the regional-tissue TACs are user dependent; however, using factor analysis we could get pure regional-tissue TACs. Therefore, rMBF could be calculated robustly with the single-compartment model using the tissue TACs obtained from factor analysis. Although the estimated rMBFs obtained by factor analysis were similar to those obtained by the conventional ROI method, factor analysis is preferable to the ROI method because of the subjectivity of the latter.
Factor analysis has been applied successfully to images with a high signal-to-noise ratio (SNR), such as dynamic 13NH3 or 18F-FDG PET images (5,6). Because factor analysis for dynamic images with a low SNR was considered difficult, a modeling method was proposed for dynamic H215O PET (9,15). In this study, we reduced statistical fluctuation by aggregating into 6 planes the 24 planes that were selected to include only the cardiac region. The pixels were also aggregated into an 8 x 8 region before factor analysis was applied to the dynamic images. Moreover, in the first step of the factor analysis, the inherent difficulty of applying factor analysis to dynamic H215O PET could be overcome using the multiplane dynamic data.
The low SNR of dynamic H215O PET images might have influenced factor analysis, causing inaccurate oblique rotation during apex seeking (16,17). Inaccuracy in apex seeking affects the quantitative results for rMBF by using the TACs misrepresented by factor analysis. A trend was seen toward slightly higher rMBFs by factor analysis than by the ROI method and microsphere studies. This bias might have been caused by a low SNR for the dynamic H215O PET images in the ROI chosen for the second step of the factor analysis. Contrary to the first step of the factor analysis, which used all six planes, the second step of the factor analysis used only a quarter of a slice for extraction of the regional-tissue TAC. Although several methods have been proposed to improve the efficiency of apex seeking by overcoming this shortcoming, these methods were not applied in this study.
Spillover fraction was 0.1 ± 0.1 by factor analysis and smaller (0.2 ± 0.1) by the ROI method, but the difference was not significant. Although we expected that use of TACs obtained from the second step of the factor analysis would remove the spillover effect from rMBF estimation, a fraction of spillover remained, indicating that our factor analysis did not perfectly separate tissue TACs from blood TACs. Nonetheless, we speculate that factor analysis was better than the subjective ROI drawing method at least for obtaining regional-tissue TACs. Further analysis is warranted to find the differences between blood-pool factor images and C15O blood-pool images.
Activity was found in the lowest planes outside the myocardium. We suspect that this activity was from the liver. Although the earlier part of a TAC for liver would be different from that for myocardium because the liver is perfused mostly from the portal vein, the liver could have mimicked the TAC of myocardium. Whether similar activity would show up on human myocardial dynamic H215O PET images is unknown. If such adjacent-organ activity were to appear, further studies would be required to exclude it.
Hermansen et al. (9) reported that the SNR of a bolus injection was higher than that of a slow infusion because the shapes of the blood and myocardium TACs differ more from each other with a bolus injection. If the tracer is infused slowly, the left ventricular and myocardial factor TACs may be less distinctive on factor analysis. If one wishes to separate the TAC of the injected H215O in the myocardium from the TACs of radioactivity in the ventricles, rapid injection of the tracer by bolus may be better than slow infusion. In contrast, bolus injection with a high dose of H215O may induce a dead-time effect, underestimate the input function, and thus overestimate rMBF. Therefore, injection time should be determined carefully after considering both the dead-time effect of a bolus injection and the lower probability that a slow infusion will produce clearly distinctive factor images.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
For correspondence or reprints contact: Dong Soo Lee, MD, PhD, Department of Nuclear Medicine, Seoul National University College of Medicine, 28 Yungundong Chongnogu, Seoul 110-744, Korea.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. S. Lee, D. S. Lee, J. Y. Ahn, J. S. Yeo, G. J. Cheon, S.-K. Kim, K. S. Park, J.-K. Chung, and M. C. Lee Generation of Parametric Image of Regional Myocardial Blood Flow Using H215O Dynamic PET and a Linear Least-Squares Method J. Nucl. Med., October 1, 2005; 46(10): 1687 - 1695. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. A. Kaufmann and P. G. Camici Myocardial Blood Flow Measurement by PET: Technical Aspects and Clinical Applications J. Nucl. Med., January 1, 2005; 46(1): 75 - 88. [Full Text] [PDF] |
||||
![]() |
C. Katoh, K. Morita, T. Shiga, N. Kubo, K. Nakada, and N. Tamaki Improvement of Algorithm for Quantification of Regional Myocardial Blood Flow Using 15O-Water with PET J. Nucl. Med., November 1, 2004; 45(11): 1908 - 1916. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. N. Mazzadi, M. F. Janier, B. Brossier, X. Andre-Fouet, D. Revel, and P. Croisille Tagged MRI and PET in severe CAD: discrepancy between preoperative inotropic reserve and intramyocardial functional outcome after revascularization Am J Physiol Heart Circ Physiol, November 1, 2004; 287(5): H2226 - H2233. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. N. Mazzadi, M. F. Janier, B. Brossier, X. Andre-Fouet, E. McFadden, D. Revel, and P. Croisille Dobutamine-tagged MRI for inotropic reserve assessment in severe CAD: relationship with PET findings Am J Physiol Heart Circ Physiol, May 1, 2004; 286(5): H1946 - H1953. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Koepfli, T. F. Hany, C. A. Wyss, M. Namdar, C. Burger, A. V. Konstantinidis, T. Berthold, G. K. von Schulthess, and P. A. Kaufmann CT Attenuation Correction for Myocardial Perfusion Quantification Using a PET/CT Hybrid Scanner J. Nucl. Med., April 1, 2004; 45(4): 537 - 542. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. N. Mazzadi, X. Andre-Fouet, J. Duisit, V. Gebuhrer, N. Costes, P. Chevalier, C. Rodriguez, J.-J. Schott, H. Le Marec, P. Guicheney, et al. Cardiac retention of [11C]HED in genotyped long QT patients: a potential amplifier role for severity of the disease Am J Physiol Heart Circ Physiol, August 7, 2003; 285(3): H1286 - H1293. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. A. Wyss, P. Koepfli, K. Mikolajczyk, C. Burger, G. K. von Schulthess, and P. A. Kaufmann Bicycle Exercise Stress in PET for Assessment of Coronary Flow Reserve: Repeatability and Comparison with Adenosine Stress J. Nucl. Med., February 1, 2003; 44(2): 146 - 154. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| JOURNAL OF NUCLEAR MEDICINE TECHNOLOGY | THE JOURNAL OF NUCLEAR MEDICINE |