Multiple Linear Analysis for Generating Parametric Images of Irreversible Radiotracer

비가역 방사성추적자 파라메터 영상을 위한 다중선형분석법

  • Kim, Su-Jin (Departments of Nuclear Medicine, Seoul National University College of Medicine) ;
  • Lee, Jae-Sung (Departments of Nuclear Medicine, Seoul National University College of Medicine) ;
  • Lee, Won-Woo (Departments of Nuclear Medicine, Seoul National University College of Medicine) ;
  • Kim, Yu-Kyeong (Departments of Nuclear Medicine, Seoul National University College of Medicine) ;
  • Jang, Sung-June (Departments of Nuclear Medicine, Seoul National University College of Medicine) ;
  • Son, Kyu-Ri (Departments of Radiology, Seoul National University College of Medicine) ;
  • Kim, Hyo-Cheol (Departments of Radiology, Seoul National University College of Medicine) ;
  • Chung, Jin-Wook (Departments of Radiology, Seoul National University College of Medicine) ;
  • Lee, Dong-Soo (Departments of Nuclear Medicine, Seoul National University College of Medicine)
  • 김수진 (서울대학교 의과대학 핵의학교실) ;
  • 이재성 (서울대학교 의과대학 핵의학교실) ;
  • 이원우 (서울대학교 의과대학 핵의학교실) ;
  • 김유경 (서울대학교 의과대학 핵의학교실) ;
  • 장성준 (서울대학교 의과대학 핵의학교실) ;
  • 손규리 (서울대학교 의과대학 방사선과학교실) ;
  • 김효철 (서울대학교 의과대학 방사선과학교실) ;
  • 정진욱 (서울대학교 의과대학 방사선과학교실) ;
  • 이동수 (서울대학교 의과대학 핵의학교실)
  • Published : 2007.08.31

Abstract

Purpose: Biological parameters can be quantified using dynamic PET data with compartment modeling and Nonlinear Least Square (NLS) estimation. However, the generation of parametric images using the NLS is not appropriate because of the initial value problem and excessive computation time. In irreversible model, Patlak graphical analysis (PGA) has been commonly used as an alternative to the NLS method. In PGA, however, the start time ($t^*$, time where linear phase starts) has to be determined. In this study, we suggest a new Multiple Linear Analysis for irreversible radiotracer (MLAIR) to estimate fluoride bone influx rate (Ki). Methods: $[^{18}F]Fluoride$ dynamic PET scans was acquired for 60 min in three normal mini-pigs. The plasma input curve was derived using blood sampling from the femoral artery. Tissue time-activity curves were measured by drawing region of interests (ROls) on the femur head, vertebra, and muscle. Parametric images of Ki were generated using MLAIR and PGA methods. Result: In ROI analysis, estimated Ki values using MLAIR and PGA method was slightly higher than those of NLS, but the results of MLAIR and PGA were equivalent. Patlak slopes (Ki) were changed with different $t^*$ in low uptake region. Compared with PGA, the quality of parametric image was considerably improved using new method. Conclusion: The results showed that the MLAIR was efficient and robust method for the generation of Ki parametric image from $[^{18}F]Fluoride$ PET. It will be also a good alternative to PGA for the radiotracers with irreversible three compartment model.

동적 PET 데이터는 구획모델과 Nonlinear Least Squares(NLS)방법을 사용하여 분석함으로서 각종 생화학적 물질의 생체 대사율 등을 정량화 할 수 있다. 하지만 NLS방법은 변수의 초기 값이 적절하지 않을 경우 지역적인 최소점에 빠지거나 계산시간이 길어 동역학 변수들을 각 화소마다 구해야 하는 파라메터 영상(parametric image) 구성에는 효과적이지 않다. Patlak 도표분석법(PGA)은 비선형 방정식을 선형화하여 값을 추정함으로서 간단하면서 적은 계산량으로 인해 파라메터 영상을 구성하는데 많이 사용되고 있으나 잡음성분과 선형구간 선정에 따라 값이 영향을 받는 단점이 있다. 따라서 이 연구에서는 3구획 비가역 모델에 적합한 다중선형분석법(MLAIR)을 고안하였으며 3구획 비가역 모델의 대표적 예인 $[^{18}F]Fluoride$ PET을 이용하여 미니돼지에서의 뼈 섭취률을 계산하여 PGA방법과 비교 분석해 보았다. 대상 및 방법: 3마리의 미니돼지를 대상으로 100MBq의 $[^{18}F]Fluoride$를 대퇴부 정맥에 주사하면서 ECAT EXAET 47 PET 스캐너를 이용하여 60분간 PET 영상을 얻었다. 케타민과 자일라진을 이용하여 30분 간격으로 마취하였으며 실험동물을 진공 쿠션을 이용하여 반드시 누운 자세로 위치하도록 고정 시켰다. 입력함수인 혈장 내 농도곡선은 대퇴동맥으로부터 스캔 시작과 함께 혈액 채취를 통해 얻었다. ROI분석을 위해 대퇴골두, 척추 뼈, 근육에 ROI를 그려 조직 내 시간-방사능 곡선을 얻었다. $k_4$가 0인 3구획 비가역 모델로부터 MLAIR와 PGA방법을 사용하여 관심영역에서의 뼈 섭취률 $K_i$와 파라메터 영상을 구성하였다. PGA방법은 선형구간의 시작점인 $t^*$선택에 따른 영향을 보기 위해 분석구간을 변화시켜가며 분석하였다. 결과: ROI 분석결과 추정된 $K_i$값은 NLS방법에 비하여 MLAIR방법과 PGA방법 모두에서 과대 추정되었으나 두 분석방법 끼리는 비슷한결과를 보였다. Patlak 기울기는 $t^*$선택에 따라 값이 변하였으며 Patlak 상수는 Fluoride 섭취가 높은 대퇴골두나 척추뼈에서 30분이 지나서야 일정한 값으로 나타났고 섭취가 낮은 근육에서는 10분만에 일정해졌다. 파라메터 영상에서는 제안한 MLAIR방법이 PGA방법에 비해 영상의 질을 향상시킴을 알 수 있었다. 또한 PGA방법을 이용하여 구성한 파라메터 영상은 $t^*$값이 커질수록 급격히 영상의 질이 저하됨을 볼 수 있다. 특히 Fluoride 섭취가 높은 영역에서 Patlak 상수가 일정해지는 시간인 $t^*$값이 30분일 때 파라메터 영상은 MLAIR와 크게 차이가 났다. 결론 결론적으로 제안한 MLAIR방법은 선형구간을 정할 필요 없이 모든 데이터를 사용하는 이점이 있으며 선형적인 방법을 통해 $K_i$값을 얻을 수 있어 계산시간을 단축 시켜 줄 뿐 아니라 잡음성분에 강해 파라메터 영상의 질을 크게 향상 시켜 줌으로 비가역 3구획모델에서의 PGA방법을 대체할 새로운 파라메터 영상구성방법으로 적합할 것이다.

Keywords

References

  1. Feng D, Huang S-C, Wang Z, and Ho D. An unbiased parametric imaging algorithm for nonuniformly sampled biomedical system parameter estimation. IEEE Trans. Med Imaging 1996;15:512-18 https://doi.org/10.1109/42.511754
  2. Kimura Y, Hsu H, Toyama H, Senda M, and Alpert N. M. Improved signal-to-noise ratio in parametric images by cluster analysis Neuroimage 1999;9:554-61 https://doi.org/10.1006/nimg.1999.0430
  3. Lee JS, Lee DS, Ahn JY, Yeo JS, Cheon GJ, Kim SK, et al. Generation of parametric image of regional myocardial blood flow using $H_2^{15}O$dynamic PET and a linear least-squares method. J Nucl Med 2005;46:1687-95
  4. Lee JS, Lee DS, Ahn JY, Cheon GJ, Kim SK, Yeo JS, et al. Parametric image of myocardial blood flow generated from dynamic $H_2^{15}O$ PET using factor analysis and cluster analysis. Med Biol Eng Comput 2005;43:678-85 https://doi.org/10.1007/BF02351043
  5. Patlak CS, Blasberg RG, Fenstermacher JD. Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab 1983;3:1-7 https://doi.org/10.1038/jcbfm.1983.1
  6. Choi Y, Hawkins RA, Huang SC, Gambhir SS, Brunken RC, Phelps ME, el at. Parametric images of myocardial metabolic rate of glucose generated from dynamic cardiac PET and 2- [18F]fluoro-2-deoxy-d-glucose studies. J Nucl Med 1991;32:733-8
  7. Logan J, Fowler JS, Volkow ND, Wolf AP, Dewey SL, Schlyer DJ, et al. Graphical analysis of reversible radioligand binding from time activity measurements applied to [N-11C-methyl]-(-)-cocaine PET studies in human subjects. J Cereb Blood Flow Metab 1990;10:740-7 https://doi.org/10.1038/jcbfm.1990.127
  8. Hawkins RA, Choi Y, Huang SC, Hoh CK, Dahlbom M, Schiepers C, et al. Evaluation of the skeletal kinetics of fluorine-18-fluoride ion with PET. J Nucl Med 1992;33:633-42
  9. Cook GJ, Lodge MA, Blake GM, Marsden PK, Fogelman I. Differences in skeletal kinetics between vertebral and humoral bone measured by 18F-fluoride positron emission tomography in postmenopausal women. J Bone Miner Res 2000;15:763-9 https://doi.org/10.1359/jbmr.2000.15.4.763
  10. Schiepers C, Nuyts J, Bormans G, Dequeker J, Bouillon R, Mortelmans L, et al. Fluoride kinetics of the axial skeleton measured in vivo with fluorine-18-fluoride PET. J Nucl Med 1997;38:1970-6
  11. Brenner W, Vernon C, Muzi M, Mankoff DA, Link JM, Conrad EU, et al. Comparison of different quantitative approaches to $^{18}F-fluoride PET scans. J Nucl Med 2004;45:1493-500
  12. Piert M, Zittel TT, Becker GA, Jahn M, Stahlschmidt A, Maier G, et al. Assessment of porcine bone metabolism by dynamic [$^{18}F]fluoride ion PET: correlation with bone histomorphometry. J Nucl Med 2001;42:1091-100
  13. Brenner W, Kampen WU, Muller H, Scholz-Ahrens KE, Soti Z, Schrezenmeir J, et al. F-18 fluoride PET in minipigs with defined disorders of bone metabolism. J Nucl Med 2002;43:S37
  14. Kim JS, Lee JS, Lee DS, Chung JK, Lee MC. Performance Evaluation of Siemens CTI ECAT EXACT 47 Scanner Using NEMA NU2-2001. Korean J Nucl Med 2004;38:259-67
  15. Slifstein M, Laruelle M. Effects of statistical noise on graphic analysis of PET neuroreceptor studies. J Nucl Med 2000;41:2083-8
  16. Varga J, Szabo Z. Modified regression model for the Logan plot. J Cereb Blood Flow Metab 2002;22:240-4 https://doi.org/10.1097/00004647-200202000-00012
  17. Carson RE. PET parameter estimation using linear integration methods: bias and variability consideration. In: Quantification of brain function: tracer kinetics and image analysis in brain PET Amsterdam: Elsevier Science Publishers B.V., 1993. p499-507
  18. Evans AC. A double integral form of the three-compartment model, four rate-constant model for faster generation of parameter maps. J Cereb Blood Flow Metab 1987;7:S453