CN105748093A - 大脑灰质作参考区的人脑局部水分布容积测定方法 - Google Patents
大脑灰质作参考区的人脑局部水分布容积测定方法 Download PDFInfo
- Publication number
- CN105748093A CN105748093A CN201511009429.8A CN201511009429A CN105748093A CN 105748093 A CN105748093 A CN 105748093A CN 201511009429 A CN201511009429 A CN 201511009429A CN 105748093 A CN105748093 A CN 105748093A
- Authority
- CN
- China
- Prior art keywords
- water distribution
- human brain
- distribution volume
- gray matter
- cerebral gray
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/40—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
- A61B6/4057—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis by using radiation sources located in the interior of the body
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/107—Measuring physical dimensions, e.g. size of the entire body or parts thereof
- A61B5/1073—Measuring volume, e.g. of limbs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/50—Clinical applications
- A61B6/501—Clinical applications involving diagnosis of head, e.g. neuroimaging, craniography
Abstract
本发明公开了一种大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,包括如下步骤:步骤一:对人脑做PET/H2 15O动态扫描;步骤二:将人脑的大脑灰质作为参考区,计算出人脑中每个待求部位的相对水分布容积值P1;步骤三:求出人脑每个部位的相对水分布容积值P1后,根据需要作人脑的归一化处理;所得数值与该部位的绝对水分布容积仅相差一比例常数,用该所得数值作为绝对水分布容积的指标。本发明既不需要动态抽取动脉血也不需要可能失真的模版或模型,更不需要做复杂的修正等人工干预,因而便于实现计算机的自动化处理。同时测定的结果与精确局部水分布容积定量测定方法没有偏差,具有很好的应用情景。
Description
技术领域
本发明涉及一种人脑局部水分布容积测定方法,具体涉及一种大脑灰质作参考区的人脑局部水分布容积测定方法。本发明属于生物及医学成像领域。
背景技术
利用正电子放射性核素15O标记的水(即H2 15O)作显像剂的正电子发射断层显像(PositronEmissionTomography,即PET)方法(简记为PET/H2 15O方法),是对静脉注射H215O后进行PET动态扫描,而后用Kety-Schmidt单室模型(文献见KetySS.(1951)Thetheoryandapplicationsoftheexchangeofinertgasatthelungsandtissues.PharmacolRev3:1-41.)等方法作局部水分布容积逐像素定量计算,这样可以生成局部水分布容积的图像。
早期的定量计算中输入函数需要通过动态抽取动脉血或者动脉化的静脉血的方法来测得。尽管这种方法很准确和有效,但它是有损伤的,因此这种方法日常很少应用。为了减少或者消除这种抽血方法的影响,一些研究者发明了一种称为源自图像的(image-derived)的方法(包括两份文献,文献1见WeinbergI,HuangSC,HoffmanE,etal.(1988)ValidationofPET-acquiredinputfunctionsforcardiacstudies.JNuclMed29:241-247,文献2见ChenK,BandyD,ReimanE,etal.(1998)NoninvasiveQuantificationoftheCerebralMetabolicRateforGlucoseUsingPositronEmissionTomography,18F-Fluoro-2-Deoxyglucose,thePatlakMethod,andanImage-DerivedInputFunction.JCerebBloodFlowMetab18:716–723)。将PET图像中主要反映血浆特性的动脉血管或左心室的PET数据作为输入函数。这些方法虽然成功地用来代替动态抽血,但不准确或者需要许多复杂的手工修正。有鉴于此,一种从一群人血样数据中计算出的输入函数模板或对此血样数据进行拟合给出的输入函数解析模型(文献见FengD,HuangSC,WangX(1993)Modelsforcomputersimulationstudiesofinputfunctionsfortracerkineticmodelingwithpositronemissiontomography.IntJBiomedComput32:95-110.)被引进,但这些模版或模型不能很好地符合个别被试的特殊血样数据,这样的输入函数无法简单应用于水分布容积的定量计算。
正如前述,现有技术存在下述问题:
1)通过动态抽取动脉血或者动脉化的静脉血获取输入函数的方法,尽管很准确和有效,但它是有损伤的,容易带来不舒适,甚至会引起血管方面的其他危险,并有血样放射物处理以及对人的干扰等麻烦,因此这种方法日常很少应用。
2)源自图像的(image-derived)的方法,存在这样的问题:a)需要快速的图像采集方案来反映输入函数的时间曲线特征,特别是注射FDG后早期的迅速变化;b)需要有效解决局部体积效应;c)需要修正从围绕所选ROI的脑区到这些ROI的外流效应(spillover)。当示踪剂在围绕所选ROI的脑区中开始积聚时,这种组织到血浆的外流将变得很严重。因而这些方法虽然成功地用来代替动态抽血,但不准确或者需要许多复杂的手工修正。
3)对于输入函数模版或模型方法,由于每个个体的特殊血样数据不一定具有相同的特定值,如果统一选用这样的模版或模型将导致给出的输入函数失真,无法简单应用于水分布容积的定量计算。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种大脑灰质作参考区的人脑局部水分布容积测定方法,以解决现有技术在测定人脑局部水分布容积时准确性不足的技术问题。
为了实现上述目标,本发明采用如下的技术方案:
大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,包括如下步骤:
步骤一:对人脑做PET/H2 15O动态扫描;
步骤二:将人脑的大脑灰质作为参考区,计算出人脑中每个待求部位的相对水分布容积值P1;
步骤三:求出人脑每个部位的相对水分布容积值P1后,根据需要作人脑的归一化处理;所得数值与该部位的绝对水分布容积仅相差一比例常数,用该所得数值作为绝对水分布容积的指标。
前述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,所述的PET/H2 15O动态扫描为用H2 15O作显像剂的PET动态扫描方法。
前述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,步骤二中,采用公式: 计算人脑中每个待求部位的相对水分布容积值P1;其中:cto(t)为t时刻待求部位的单位体积的H2 15O放射性活度,ctr(t)为t时刻大脑灰质的单位体积的H2 15O放射性活度,P2、P3是其它两个拟合参数;通过动态测量采样,拟合得到P1、P2、P3。
前述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,步骤三中,在需要测定绝对水分布容积值时,将测定出的相对水分布容积值P1乘以大脑灰质的绝对水分布容积值即可。
前述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,可以基于像素进行定量测定。
前述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,归一化后的图像至少能够用于单个人体不同状态或不同人体间的比较研究。
本发明的有益之处在于:本发明可以对所分析的器官逐像素进行局部水分布容积的定量测定,进而生成人脑局部水分布容积的功能图像,无需输入函数,从而既不需要动态抽取动脉血也不需要可能失真的模版或模型,更不需要做复杂的修正等人工干预,因而便于实现计算机的自动化处理。同时测定的结果与精确局部水分布容积定量测定方法没有偏差,这种完全无损伤的定量测定方法具有很好的应用情景。
附图说明
图1是本发明大脑灰质作参考区的人脑局部水分布容积测定方法的流程示意图;
图2是噪声水平为1的情况下本发明与Kety-Schmidt单室模型方法测定的相对水分布容积的回归比较示意图;
图3是是噪声水平为1的情况下本发明与Kety-Schmidt单室模型方法测定的绝对水分布容积的回归比较示意图。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
参照图1所示,本发明大脑灰质作参考区的人脑局部水分布容积测定方法,为一种人脑局部水分布容积的PET/H2 15O无损伤定量测定方法,该方法可以对人脑逐像素进行定量测定,进而生成人脑局部水分布容积局的功能图像。为实现上述发明目的,本发明是采取以下的技术方案来实现的:
大脑灰质作参考区的人脑局部水分布容积定量测定方法,其特征在于:
(1)对人脑做PET/H2 15O动态扫描;PET/H2 15O动态扫描为用H2 15O作显像剂的PET动态扫描方法;
(2)将大脑灰质作为参考区(以下简称大脑灰质),然后直接采用下列公式:
(3)当求出人脑每个部位的相对水分布容积值P1后,根据需要可作人脑的归一化处理。所得数值与该部位的绝对水分布容积仅相差一比例常数,可以用它作绝对水分布容积的指标。这个归一化后的图像可以用于单个人体不同状态或不同人体间的比较研究等。
当需要测定绝对水分布容积值时,将测定出的相对水分布容积值P1乘以大脑灰质的绝对水分布容积值即得。本发明方法测定的绝对水分布容积值全域分布的特征与真实情况相比没有偏差。本发明可以基于像素进行定量测定。
下面再结合图2、图3介绍用模拟数据(根据已有文献报道的实验数据参数用计算机随机产生的PET/H2 15O动态扫描数据)进行的本发明与国际公认的金标准Kety-Schmidt单室模型方法测定结果的比较:
图2、图3分别是选用大脑灰质作参考区(绝对血流值取为0.8毫升/毫升/分钟,绝对水分布容积值取1.02毫升/毫升,此值取之于已有文献报道的实验结果,此值为文献中已有实验所有被测者的平均值。)在噪声水平为1的情况下大脑逐像素本发明与Kety-Schmidt单室模型方法两种测定出的相对(图2)和绝对(图3)水分布容积值的回归比较。图中散点对应每个像素的结果,直线为两种方法的回归直线,其斜率都为1.02,很接近于1;截距都为0.005,非常接近于0。这表明本发明的结果与Kety-Schmidt单室模型方法几乎是成正比的,因而二者相等价。这些结果显示本发明和Kety-Schmidt单室模型方法的计算结果近似,但不像传统的方法那样需要输入函数,从而免去了前述背景技术中所描述的那些困惑。
由上所述可见,本发明选用大脑灰质作参考区后,直接利用前述的无损伤局部水分布容积定量计算公式,输入函数将消失。当需要测定绝对水分布容积时,将测定出的相对水分布容积值乘以大脑灰质的绝对水分布容积值即可。这样计算出每个像素或区域的绝对水分布容积值与真值相比仅相差一小的比例常数,甚至计算出每个像素或区域的绝对水分布容积值与真值相等。也就是所计算出绝对水分布容积值的全域分布特征与真实情况相比不会出现偏差。
由于本发明方法不像源自图像方法那样需要复杂的修正,且不需要输入函数,从而既不需要动态抽取动脉血也不需要做复杂的修正等人工干预,因而便于实现计算机的自动化处理。同时我们的结果与精确局部水分布容积定量测定方法没有偏差,本发明有望应用于日常需要。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (6)
1.大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,包括如下步骤:
步骤一:对人脑做PET/H2 15O动态扫描;
步骤二:将人脑的大脑灰质作为参考区,计算出人脑中每个待求部位的相对水分布容积值P1;步骤三:求出人脑每个部位的相对水分布容积值P1后,根据需要作人脑的归一化处理;所得数值与该部位的绝对水分布容积仅相差一比例常数,用该所得数值作为绝对水分布容积的指标。
2.根据权利要求1所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,所述的PET/H2 15O动态扫描为用H2 15O作显像剂的PET动态扫描方法。
3.根据权利要求2所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,步骤二中,采用公式:
4.根据权利要求3所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,步骤三中,在需要测定绝对水分布容积值时,将测定出的相对水分布容积值P1乘以大脑灰质的绝对水分布容积值即可。
5.根据权利要求4所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,基于像素进行定量测定。
6.根据权利要求5所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,归一化后的图像至少能够用于单个人体不同状态或不同人体间的比较研究。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009429.8A CN105748093B (zh) | 2015-12-28 | 2015-12-28 | 大脑灰质作参考区的人脑局部水分布容积测定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511009429.8A CN105748093B (zh) | 2015-12-28 | 2015-12-28 | 大脑灰质作参考区的人脑局部水分布容积测定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105748093A true CN105748093A (zh) | 2016-07-13 |
CN105748093B CN105748093B (zh) | 2018-11-16 |
Family
ID=56342216
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511009429.8A Active CN105748093B (zh) | 2015-12-28 | 2015-12-28 | 大脑灰质作参考区的人脑局部水分布容积测定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105748093B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110530902A (zh) * | 2019-08-26 | 2019-12-03 | 中国科学院地质与地球物理研究所 | 一种天然气水合物组分运移pet成像方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101172038A (zh) * | 2007-09-27 | 2008-05-07 | 南京信息工程大学 | 用于人脑等器官局部葡萄糖代谢率的无损伤定量计算方法 |
CN101320468A (zh) * | 2007-06-07 | 2008-12-10 | 株式会社东芝 | 数据处理装置、医用诊断装置、数据处理方法以及医用诊断方法 |
US20100185094A1 (en) * | 2009-01-21 | 2010-07-22 | Kabushiki Kaisha Toshiba | Ultrasound diagnosis apparatus, medical image display apparatus and medical image displaying method |
CN102068269A (zh) * | 2011-01-21 | 2011-05-25 | 南京信息工程大学 | 用于人体器官局部血流的无损伤定量测定方法 |
US20140010430A1 (en) * | 2011-02-24 | 2014-01-09 | Dog Microsystems Inc. | Method and apparatus for isolating a potential anomaly in imaging data and its application to medical imagery |
US20140177934A1 (en) * | 2012-06-20 | 2014-06-26 | Toshiba Medical Systems Corporation | Image diagnosis device and control method thereof |
-
2015
- 2015-12-28 CN CN201511009429.8A patent/CN105748093B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320468A (zh) * | 2007-06-07 | 2008-12-10 | 株式会社东芝 | 数据处理装置、医用诊断装置、数据处理方法以及医用诊断方法 |
CN101172038A (zh) * | 2007-09-27 | 2008-05-07 | 南京信息工程大学 | 用于人脑等器官局部葡萄糖代谢率的无损伤定量计算方法 |
US20100185094A1 (en) * | 2009-01-21 | 2010-07-22 | Kabushiki Kaisha Toshiba | Ultrasound diagnosis apparatus, medical image display apparatus and medical image displaying method |
CN102068269A (zh) * | 2011-01-21 | 2011-05-25 | 南京信息工程大学 | 用于人体器官局部血流的无损伤定量测定方法 |
US20140010430A1 (en) * | 2011-02-24 | 2014-01-09 | Dog Microsystems Inc. | Method and apparatus for isolating a potential anomaly in imaging data and its application to medical imagery |
US20140177934A1 (en) * | 2012-06-20 | 2014-06-26 | Toshiba Medical Systems Corporation | Image diagnosis device and control method thereof |
Non-Patent Citations (1)
Title |
---|
吴义根等: "脑局部葡萄糖代谢率PET/FDG定量计算的误差分析", 《自然科学进展》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110530902A (zh) * | 2019-08-26 | 2019-12-03 | 中国科学院地质与地球物理研究所 | 一种天然气水合物组分运移pet成像方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105748093B (zh) | 2018-11-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sullivan et al. | Metrology standards for quantitative imaging biomarkers | |
Kim et al. | Histogram analysis of intravoxel incoherent motion for differentiating recurrent tumor from treatment effect in patients with glioblastoma: initial clinical experience | |
US11769594B2 (en) | Deep learning model learning device and method for cancer region | |
Boes et al. | Image registration for quantitative parametric response mapping of cancer treatment response | |
Krishna Kumar et al. | Accurate lumen diameter measurement in curved vessels in carotid ultrasound: an iterative scale-space and spatial transformation approach | |
CN107945878A (zh) | 一种基于放射组学的肝静脉压力梯度计算模型的构建方法 | |
CN102068269B (zh) | 用于人体器官局部血流的无损伤定量测定方法 | |
GB2491942A (en) | Measuring Activity of a Tracer in Medical Imaging | |
O'Brien et al. | Automated left ventricle ischemic scar detection in CT using deep neural networks | |
Qian et al. | Segmentation of myocardium from cardiac MR images using a novel dynamic programming based segmentation method | |
CN100569183C (zh) | 用于人体器官局部葡萄糖代谢率的无损伤定量计算方法 | |
Hameeteman et al. | Carotid wall volume quantification from magnetic resonance images using deformable model fitting and learning-based correction of systematic errors | |
Chao et al. | TumourMetrics: a comprehensive clinical solution for the standardization of DCE-MRI analysis in research and routine use | |
Beache et al. | Fully automated framework for the analysis of myocardial first‐pass perfusion MR images | |
Veeraraghavan et al. | Simultaneous segmentation and iterative registration method for computing ADC with reduced artifacts from DW‐MRI | |
Teo et al. | Regional ejection fraction and regional area strain for left ventricular function assessment in male patients after first-time myocardial infarction | |
Corsi et al. | Improved quantification of left ventricular volumes and mass based on endocardial and epicardial surface detection from cardiac MR images using level set models | |
US10052076B2 (en) | Diagnostic brain imaging | |
Tan et al. | Automatic localization of the left ventricular blood pool centroid in short axis cardiac cine MR images | |
WO2021067853A1 (en) | Characterization of lesions via determination of vascular metrics using mri data | |
Kadir et al. | LV wall segmentation using the variational level set method (LSM) with additional shape constraint for oedema quantification | |
CN105748093A (zh) | 大脑灰质作参考区的人脑局部水分布容积测定方法 | |
CN105686844A (zh) | 前a%区域作参考区的人脑局部水分布容积测定方法 | |
Caiani et al. | Automated interpretation of regional left ventricular wall motion from cardiac magnetic resonance images | |
US10417764B2 (en) | System and methods for diagnostic image analysis and image quality assessment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200831 Address after: Room 158, 14 / F, building 2, 588 Zixing Road, Minhang District, Shanghai, 201100 Patentee after: Shanghai Lansheng brain medical research Co., Ltd Address before: 210016 No. 219, Ning six road, Nanjing, Jiangsu Patentee before: NANJING University OF INFORMATION SCIENCE & TECHNOLOGY |
|
TR01 | Transfer of patent right |