CN101172038A - 用于人脑等器官局部葡萄糖代谢率的无损伤定量计算方法 - Google Patents
用于人脑等器官局部葡萄糖代谢率的无损伤定量计算方法 Download PDFInfo
- Publication number
- CN101172038A CN101172038A CNA2007101332998A CN200710133299A CN101172038A CN 101172038 A CN101172038 A CN 101172038A CN A2007101332998 A CNA2007101332998 A CN A2007101332998A CN 200710133299 A CN200710133299 A CN 200710133299A CN 101172038 A CN101172038 A CN 101172038A
- Authority
- CN
- China
- Prior art keywords
- metabolic rate
- glucose metabolic
- organs
- brain
- human brain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Abstract
本发明涉及一种用于人脑及其他器官局部葡萄糖代谢率定量计算方法,其特征在于:对人体做PET/FDG动态扫描,扫描部位根据情况而定;根据需要选定大脑或其他器官的全部作为参考区,也可以选择他们的任意部分区域作为参考区,然后采用公式:计算出大脑或其他器官中每个部位的相对葡萄糖代谢率K,计算中只要cio *(t)与cir *(t)涵盖的区域大部分不重叠就可以,完全重叠时,K将自动取1;当求出选定大脑或其他器官每个部位的相对葡萄糖代谢率K值后,作全脑或全器官的归一化处理。本发明可以对所分析的器官逐像素进行定量计算,进而生成葡萄糖代谢率的功能图像,无需输入函数,因而既不需要抽血,又不需要做复杂的修正,准确性高,便于实际应用。
Description
技术领域
本发明涉及人脑及其他器官局部葡萄糖代谢率PET/FDG定量计算方法,该方法可以对所分析的器官逐像素进行定量计算,进而生成葡萄糖代谢率的功能图像,属于生物及医学成像领域。
背景技术
利用正电子放射性核素18F标记的脱氧葡萄糖(英文名:18F-fluoro-2-deoxy-D-glucose,即FDG)作显像剂的正电子发射断层显像(Positron Emission Tomography,即PET)方法(简记为PET/FDG方法),是对动脉注射FDG后进行PET动态扫描,而后用传统的Patlak方法或三室模型等方法作葡萄糖代谢率逐像素定量计算,这样可以生成葡萄糖代谢率的图像。
以往的定量计算方法中需要输入函数,早期的研究中输入函数是通过动态抽取动脉血或者动脉化的静脉血的方法来测得。尽管这种方法很准确和有效,但它是有损伤的,容易带来不舒适,甚至会引起血管方面的其他危险,并有血样放射物处理以及对人的干扰等麻烦,因此这种方法日常很少应用。
为了减少或者消除这种抽血方法的影响,一些研究者发明了一种称为源自图像的(image-derived)的方法,将PET图像中主要反映血浆特性的特定感兴趣区(ROI)的PET数据作为输入函数。这方面的研究有用左心室的血池作为ROI,但只在心脏的PET研究中适用;在肝、肾的PET研究中也有用腹部大动脉来做的;脑部的PET研究中用颈动脉来做。但这些方法存在这样的问题:a)需要快速的图像采集方案来反映输入函数的时间曲线特征,特别是注射FDG后早期的迅速变化;b)需要有效解决局部体积效应;c)需要修正从围绕所选ROI的脑区到这些ROI的外流效应(spillover)。当示踪剂在围绕所选ROI的脑区中开始积聚时,这种组织到血浆的外流将变得很严重。因而这些方法虽然成功地用来代替动态抽血,但不准确或者需要许多复杂的手工修正。
有鉴于此,研究者又提出一些近似计算输入函数的方法,他们通过将输入函数模型化来减少抽血的次数,这样的模型包括生理学模型以及用已有的成组数据来代替个体的血样数据。采用了数学模型后,输入函数将通过拟合少量抽血数据获得模型参数后来确定。但由于输入函数的形状强烈依赖于个体的生理状态(如心脏输出等)以及实际操作因素(如FDG的注射速度等),这样的模型不一定能很好地拟合各类实际数据,因而难以应用。
由于上述困难,目前通行的所谓的定量方法是标准摄取值(Standard Uptakevalue即SUV)方法。但它存在很多争议,此方法只是放射性活度的变形,仍然含有一定量的未被代谢FDG的成分,与真正的葡萄糖代谢率存在一定的偏差(仅凭PET数据无法校正),因而严格来讲并不是真正的葡萄糖代谢率定量计算方法,常发生与已有的定量方法相背离的情况,因而有学者发表论文戏称它为“可笑无用值”(Silly Useless Value)。实践上之所以采用它,只是因为它不像已有的葡萄糖代谢率定量计算方法那样需要抽血,且比之单纯的放射性活度,它又能用于不同人体间的简单比较。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种有助于提高计算准确性的用于人脑及其他器官局部葡萄糖代谢率PET/FDG定量计算方法。
为实现上述发明目的,本发明是采取以下的技术方案来实现的:
一种用于人脑等器官局部葡萄糖代谢率PET/FDG无损伤定量计算方法,其特征在于:
(1)对人体做PET/FDG动态扫描,扫描的部位根据情况而定;
(2)、根据需要选定大脑或者其他器官的全部作为参考区,也可以选择他们的任意部分区域作为参考区,然后采用下列公式:
计算出大脑或其他器官中每个部位的相对葡萄糖代谢率K,式中t0是注射FDG后PET起始扫描时间,cio *(t)为t时刻待求部位的单位体积的FDG放射性活度,cir *(t)为所选定参考区的相应值,W1、W2是两个拟合参数,通过动态测量采样,拟合得到K、W1、W2,计算中只要cio *(t)与cir *(t)涵盖的区域有大部分不重叠就可以了;
(3)当求出选定大脑或其他器官每个部位的相对葡萄糖代谢率K值后,作全脑或其他全器官的归一化处理。所得数值与该部位的绝对葡萄糖代谢率仅相差一比例常数,可以用它作绝对葡萄糖代谢率的指标。这个归一化后的图像可以用于单个人体的定量计算或不同人体间的比较研究等。
前述的用于人脑等器官局部葡萄糖代谢率定量计算方法,其特征在于计算的区域与所选参考区完全重叠时,K取为1;
前述的用于人脑等器官局部葡萄糖代谢率定量计算方法,其特征在于所述的部位可以是像素。
前述的用于人脑等器官局部葡萄糖代谢率PET/FDG无损伤定量计算方法,其特征在于所述的注射FDG后PET起始扫描时间t0为5分钟以上。
本发明的有益效果是:本发明可以对所分析的器官逐像素进行葡萄糖代谢率的定量计算,进而生成葡萄糖代谢率的功能图像,无需输入函数,因而既不需要抽血,又不需要做复杂的修正,准确性高,便于实际应用。这种完全无损伤的定量方法具有很大的应用情景。
附图说明
图1是噪声水平为1的情况下本发明与Patlak方法的回归比较示意图;
图2是本发明中任选两个不同参考区在噪声水平高达2的情况下所得结果的回归比较示意图。
具体实施方式
下面结合附图介绍本发明与国际公认的金标准Patlak方法计算结果的比较:
图1是噪声水平为1的情况下本发明与Patlak方法的回归比较。图中蓝色散点对应每个像素的结果,绿色直线为两种方法的回归直线,其斜率为0.9909,非常接近于1;截距为0.0814,接近于0。这表明本发明的结果与Patlak方法几乎是正比的,因而二者相等价。同时我们直接比较了本发明与Patlak方法的差异,在噪声水平直至4的情况下,二者相差小于5%,而本发明受噪声的影响不差于Patlak方法。这些结果显示本发明和Patlak方法的计算结果近似,但不像传统的方法那样需要输入函数,从而免去了前述背景技术中所描述的那些困惑。
图2是本发明中任选两个不同参考区在噪声水平高达2的情况所得结果的回归比较。图上蓝色散点对应每个像素的结果,绿色直线为用两个参考区所得结果的回归直线,其斜率为0.9982,非常接近于1;截距为0.0006,非常接近于0。两个参考区所得结果几乎一致,表明本发明方法的计算结果不依赖于参考区的选取。同时本发明与Patlak方法一样都采用线性参数拟合,说明本发明的计算相当稳定和快捷,结果很可靠,非常适合实践需要。
由上所述可见,本发明选用非特定参考区后,利用前述我们建立的无损伤葡萄糖代谢率定量计算公式,输入函数将消失。而由于不像前述源自图像的方法那样选特定感兴趣区(比如反映血浆特性的ROI),因而不需对它的特征有所要求,进而不用像源自图像方法那样需要复杂的修正。由于不需要输入函数,从而既不需要动态抽取动脉血也不需要做复杂的修正等人工干预,因而便于实现计算机的自动化处理。同时我们的结果与精确葡萄糖代谢率定量计算方法没有偏差,因而也避免了目前通行的SUV方法常遇到的困境,本发明有望代替SUV方法。
上述实施例不以任何形式限制本发明,凡采取等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (4)
1.用于人脑等器官局部葡萄糖代谢率定量计算方法,其特征在于:
(1)、对人体做PET/FDG动态扫描,扫描的部位根据情况而定;
(2)、根据需要选定大脑或者其他器官的全部作为参考区,也可以选择他们的任意部分区域作为参考区,然后采用下列公式:
计算出大脑或其他器官中每个部位的相对葡萄糖代谢率K,式中t0是注射FDG后PET起始扫描时间,cio *(t)为t时刻待求部位的单位体积的FDG放射性活度,cir *(t)为所选定参考区的相应值,W1、W2是两个拟合参数,通过动态测量采样,拟合得到K、W1、W2,计算中只要cio *(t)与cir *(t)涵盖的区域有大部分不重叠就可以了;
(3)、当求出选定大脑或其他器官每个部位的相对葡萄糖代谢率K值后,作全脑或其他全器官的归一化处理。
2.根据权利要求1所述的用于人脑等器官局部葡萄糖代谢率定量计算方法,其特征在于计算的区域与所选参考区完全重叠时,K取为1;
3.根据权利要求1所述的用于人脑等器官局部葡萄糖代谢率定量计算方法,其特征在于所述的部位可以是像素。
4.根据权利要求1所述的用于人脑等器官局部葡萄糖代谢率定量计算方法,其特征在于所述的注射FDG后PET起始扫描时间t0为5分钟以上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101332998A CN100569183C (zh) | 2007-09-27 | 2007-09-27 | 用于人体器官局部葡萄糖代谢率的无损伤定量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101332998A CN100569183C (zh) | 2007-09-27 | 2007-09-27 | 用于人体器官局部葡萄糖代谢率的无损伤定量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101172038A true CN101172038A (zh) | 2008-05-07 |
CN100569183C CN100569183C (zh) | 2009-12-16 |
Family
ID=39420955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007101332998A Expired - Fee Related CN100569183C (zh) | 2007-09-27 | 2007-09-27 | 用于人体器官局部葡萄糖代谢率的无损伤定量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100569183C (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102068269A (zh) * | 2011-01-21 | 2011-05-25 | 南京信息工程大学 | 用于人体器官局部血流的无损伤定量测定方法 |
CN104737016A (zh) * | 2012-08-20 | 2015-06-24 | 大塚制药株式会社 | 测定糖代谢能力的方法及用于该方法的组合物 |
CN105686844A (zh) * | 2015-12-28 | 2016-06-22 | 南京信息工程大学 | 前a%区域作参考区的人脑局部水分布容积测定方法 |
CN105748093A (zh) * | 2015-12-28 | 2016-07-13 | 南京信息工程大学 | 大脑灰质作参考区的人脑局部水分布容积测定方法 |
CN106361305A (zh) * | 2016-09-19 | 2017-02-01 | 爱国者电子科技有限公司 | 糖代谢率的测量方法和装置 |
CN108272440A (zh) * | 2018-01-24 | 2018-07-13 | 上海联影医疗科技有限公司 | 新陈代谢速率显示方法、装置、终端及医疗设备 |
CN114936989A (zh) * | 2022-07-22 | 2022-08-23 | 首都医科大学宣武医院 | 基于医学图像的葡萄糖代谢率定量分析方法及系统 |
-
2007
- 2007-09-27 CN CNB2007101332998A patent/CN100569183C/zh not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102068269A (zh) * | 2011-01-21 | 2011-05-25 | 南京信息工程大学 | 用于人体器官局部血流的无损伤定量测定方法 |
CN104737016A (zh) * | 2012-08-20 | 2015-06-24 | 大塚制药株式会社 | 测定糖代谢能力的方法及用于该方法的组合物 |
CN105686844A (zh) * | 2015-12-28 | 2016-06-22 | 南京信息工程大学 | 前a%区域作参考区的人脑局部水分布容积测定方法 |
CN105748093A (zh) * | 2015-12-28 | 2016-07-13 | 南京信息工程大学 | 大脑灰质作参考区的人脑局部水分布容积测定方法 |
CN105686844B (zh) * | 2015-12-28 | 2018-11-16 | 南京信息工程大学 | 前a%区域作参考区的人脑局部水分布容积测定方法 |
CN105748093B (zh) * | 2015-12-28 | 2018-11-16 | 南京信息工程大学 | 大脑灰质作参考区的人脑局部水分布容积测定方法 |
CN106361305A (zh) * | 2016-09-19 | 2017-02-01 | 爱国者电子科技有限公司 | 糖代谢率的测量方法和装置 |
CN108272440A (zh) * | 2018-01-24 | 2018-07-13 | 上海联影医疗科技有限公司 | 新陈代谢速率显示方法、装置、终端及医疗设备 |
CN114936989A (zh) * | 2022-07-22 | 2022-08-23 | 首都医科大学宣武医院 | 基于医学图像的葡萄糖代谢率定量分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN100569183C (zh) | 2009-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100569183C (zh) | 用于人体器官局部葡萄糖代谢率的无损伤定量计算方法 | |
Ritter et al. | Medical image analysis | |
CN102422307B (zh) | 用于交互式肝脏血管和胆管系统评估的方法、系统、装置和计算机程序产品 | |
Sari et al. | Estimation of an image derived input function with MR-defined carotid arteries in FDG-PET human studies using a novel partial volume correction method | |
CN108257134A (zh) | 基于深度学习的鼻咽癌病灶自动分割方法和系统 | |
US20090097726A1 (en) | Systems, methods and apparatus automatic segmentation of liver in multiphase contrast-enhanced medical images | |
RU2455689C2 (ru) | Создание изображения на основе ограниченного набора данных | |
Veni et al. | Echocardiography segmentation based on a shape-guided deformable model driven by a fully convolutional network prior | |
CN102301394B (zh) | 透壁灌注梯度图像分析 | |
CN102068269B (zh) | 用于人体器官局部血流的无损伤定量测定方法 | |
CN102419864B (zh) | 一种提取脑部ct图像骨骼方法及装置 | |
GB2461991A (en) | Positron emission tomography data with time activity derived framing intervals. | |
CN113889261B (zh) | 基于病理特征辅助的pet/ct自动肺癌诊断分类模型训练方法 | |
US20150049935A1 (en) | Region extraction system | |
Wang et al. | An interactive software module for visualizing coronary arteries in CT angiography | |
CN113850794A (zh) | 一种图像处理方法及装置 | |
Truong et al. | Medical images sequence normalization and augmentation: improve liver tumor segmentation from small dataset | |
CN105686844B (zh) | 前a%区域作参考区的人脑局部水分布容积测定方法 | |
Yang et al. | Development and Usability Testing of a User-Centered 3D Virtual Liver Surgery Planning System | |
CN105997120A (zh) | 前a%区域作参考区的人脑局部水分布容积测定方法 | |
Suetens et al. | Medical image analysis: what and why | |
CN105748093A (zh) | 大脑灰质作参考区的人脑局部水分布容积测定方法 | |
CN104574385A (zh) | 一种基于动态pet影像的心肌脂肪酸代谢量化检测方法 | |
Bai et al. | Progressive medical image annotation with convolutional neural network-based interactive segmentation method | |
Skurski et al. | Image processing methods for diagnostic and simulation applications in cardiology |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091216 Termination date: 20130927 |