CN105748093B - 大脑灰质作参考区的人脑局部水分布容积测定方法 - Google Patents

大脑灰质作参考区的人脑局部水分布容积测定方法 Download PDF

Info

Publication number
CN105748093B
CN105748093B CN201511009429.8A CN201511009429A CN105748093B CN 105748093 B CN105748093 B CN 105748093B CN 201511009429 A CN201511009429 A CN 201511009429A CN 105748093 B CN105748093 B CN 105748093B
Authority
CN
China
Prior art keywords
human brain
water distribution
gray matter
distribution volume
reference area
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.)
Active
Application number
CN201511009429.8A
Other languages
English (en)
Other versions
CN105748093A (zh
Inventor
吴义根
郑改革
赵德林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Lansheng brain medical research Co., Ltd
Original Assignee
Nanjing University of Information Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201511009429.8A priority Critical patent/CN105748093B/zh
Publication of CN105748093A publication Critical patent/CN105748093A/zh
Application granted granted Critical
Publication of CN105748093B publication Critical patent/CN105748093B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/40Arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4057Arrangements for generating radiation specially adapted for radiation diagnosis by using radiation sources located in the interior of the body
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/103Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/107Measuring physical dimensions, e.g. size of the entire body or parts thereof
    • A61B5/1073Measuring volume, e.g. of limbs
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/501Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the head, e.g. neuroimaging or craniography

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Neurology (AREA)
  • Neurosurgery (AREA)
  • Image Processing (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明公开了一种大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,包括如下步骤:步骤一:对人脑做PET/H2 15O动态扫描;步骤二:将人脑的大脑灰质作为参考区,计算出人脑中每个待求部位的相对水分布容积值P1;步骤三:求出人脑每个部位的相对水分布容积值P1后,根据需要作人脑的归一化处理;所得数值与该部位的绝对水分布容积仅相差一比例常数,用该所得数值作为绝对水分布容积的指标。本发明既不需要动态抽取动脉血也不需要可能失真的模版或模型,更不需要做复杂的修正等人工干预,因而便于实现计算机的自动化处理。同时测定的结果与精确局部水分布容积定量测定方法没有偏差,具有很好的应用情景。

Description

大脑灰质作参考区的人脑局部水分布容积测定方法
技术领域
本发明涉及一种人脑局部水分布容积测定方法,具体涉及一种大脑灰质作参考区的人脑局部水分布容积测定方法。本发明属于生物及医学成像领域。
背景技术
利用正电子放射性核素15O标记的水(即H2 15O)作显像剂的正电子发射断层显像(Positron Emission Tomography,即PET)方法(简记为PET/H2 15O方法),是对静脉注射H215O后进行PET动态扫描,而后用Kety-Schmidt单室模型(文献见Kety SS.(1951)Thetheory and applications of the exchange of inert gas at the lungs andtissues.Pharmacol Rev 3:1-41.)等方法作局部水分布容积逐像素定量计算,这样可以生成局部水分布容积的图像。
早期的定量计算中输入函数需要通过动态抽取动脉血或者动脉化的静脉血的方法来测得。尽管这种方法很准确和有效,但它是有损伤的,因此这种方法日常很少应用。为了减少或者消除这种抽血方法的影响,一些研究者发明了一种称为源自图像的(image-derived)的方法(包括两份文献,文献1见Weinberg I,Huang SC,Hoffman E,et al.(1988)Validation of PET-acquired input functions for cardiac studies.J Nucl Med 29:241-247,文献2见Chen K,Bandy D,Reiman E,et al.(1998)Noninvasive Quantificationof the Cerebral Metabolic Rate for Glucose Using Positron EmissionTomography,18F-Fluoro-2-Deoxyglucose,the Patlak Method,and an Image-DerivedInput Function.J Cereb Blood Flow Metab 18:716–723)。将PET图像中主要反映血浆特性的动脉血管或左心室的PET数据作为输入函数。这些方法虽然成功地用来代替动态抽血,但不准确或者需要许多复杂的手工修正。有鉴于此,一种从一群人血样数据中计算出的输入函数模板或对此血样数据进行拟合给出的输入函数解析模型(文献见Feng D,Huang SC,Wang X(1993)Models for computer simulation studies of input functions fortracer kinetic modeling with positron emission tomography.Int J Biomed Comput32: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)将大脑灰质作为参考区(以下简称大脑灰质),然后直接采用下列公式:
计算出人脑 中每个待求部位的相对水分布容积值P1,式中cto(t)为t时刻待求部位的单位体积的H2 15O放 射性活度,ctr(t)为t时刻大脑灰质的单位体积的H2 15O放射性活度,P2、P3是其它两个拟合参 数;通过动态测量采样,拟合得到P1、P2、P3
(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 (4)

1.大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,包括如下步骤:
步骤一:对人脑做PET/H2 15O动态扫描;
步骤二:将人脑的大脑灰质作为参考区,计算出人脑中每个待求部位的相对水分布容积值P1
步骤二中,采用公式:
计算人脑中每个待求部位的相对水分布容积值P1;其中:cto(t)为t时刻待求部位的单位体积的H2 15O放射性活度,ctr(t)为t时刻大脑灰质的单位体积的H2 15O放射性活度,P2、P3是其它两个拟合参数;通过动态测量采样,拟合得到P1、P2、P3
步骤三:求出人脑每个部位的相对水分布容积值P1后,根据需要作人脑的归一化处理;所得数值与该部位的绝对水分布容积仅相差一比例常数,用该所得数值作为绝对水分布容积的指标;
步骤三中,在需要测定绝对水分布容积值时,将测定出的相对水分布容积值P1乘以大脑灰质的绝对水分布容积值即可。
2.根据权利要求1所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,所述的PET/H2 15O动态扫描为用H2 15O作显像剂的PET动态扫描方法。
3.根据权利要求1所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,基于像素进行定量测定。
4.根据权利要求3所述的大脑灰质作参考区的人脑局部水分布容积测定方法,其特征在于,归一化后的图像至少能够用于单个人体不同状态或不同人体间的比较研究。
CN201511009429.8A 2015-12-28 2015-12-28 大脑灰质作参考区的人脑局部水分布容积测定方法 Active CN105748093B (zh)

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 CN105748093A (zh) 2016-07-13
CN105748093B true 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)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110530902B (zh) * 2019-08-26 2022-04-12 中国科学院地质与地球物理研究所 一种天然气水合物组分运移pet成像方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101172038A (zh) * 2007-09-27 2008-05-07 南京信息工程大学 用于人脑等器官局部葡萄糖代谢率的无损伤定量计算方法
CN101320468A (zh) * 2007-06-07 2008-12-10 株式会社东芝 数据处理装置、医用诊断装置、数据处理方法以及医用诊断方法
CN102068269A (zh) * 2011-01-21 2011-05-25 南京信息工程大学 用于人体器官局部血流的无损伤定量测定方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5433240B2 (ja) * 2009-01-21 2014-03-05 株式会社東芝 超音波診断装置及び画像表示装置
WO2012113069A1 (en) * 2011-02-24 2012-08-30 Dog Microsystems Inc. Method and apparatus for isolating a potential anomaly in imaging data and its application to medical imagery
JP6222795B2 (ja) * 2012-06-20 2017-11-01 東芝メディカルシステムズ株式会社 画像診断装置及びその制御方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101320468A (zh) * 2007-06-07 2008-12-10 株式会社东芝 数据处理装置、医用诊断装置、数据处理方法以及医用诊断方法
CN101172038A (zh) * 2007-09-27 2008-05-07 南京信息工程大学 用于人脑等器官局部葡萄糖代谢率的无损伤定量计算方法
CN102068269A (zh) * 2011-01-21 2011-05-25 南京信息工程大学 用于人体器官局部血流的无损伤定量测定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
脑局部葡萄糖代谢率PET/FDG定量计算的误差分析;吴义根等;《自然科学进展》;20020228;第12卷(第2期);全文 *

Also Published As

Publication number Publication date
CN105748093A (zh) 2016-07-13

Similar Documents

Publication Publication Date Title
US10580526B2 (en) System and method for calculating vessel flow parameters based on angiography
JP5953438B2 (ja) 冠状動脈を通る血流量を決定する方法
Feng et al. A technique for extracting physiological parameters and the required input function simultaneously from PET image measurements: Theory and simulation study
Boellaard et al. Evaluation of basis function and linear least squares methods for generating parametric blood flow images using 15 O-water and positron emission tomography
So et al. Quantitative myocardial perfusion measurement using CT perfusion: a validation study in a porcine model of reperfused acute myocardial infarction
Islam et al. Estimation of arterial input by a noninvasive image derived method in brain H215O PET study: confirmation of arterial location using MR angiography
GB2491942A (en) Measuring Activity of a Tracer in Medical Imaging
CN102068269B (zh) 用于人体器官局部血流的无损伤定量测定方法
CN105748093B (zh) 大脑灰质作参考区的人脑局部水分布容积测定方法
US10867383B2 (en) Determination of a clinical characteristic using a combination of different recording modalities
Su et al. Quantification method in [18F] fluorodeoxyglucose brain positron emission tomography using independent component analysis
CN105686844B (zh) 前a%区域作参考区的人脑局部水分布容积测定方法
Kudomi et al. Reconstruction of an input function from a dynamic PET water image using multiple tissue curves
Zhang et al. Quantification of effects of mean blood pressure and left ventricular mass on noninvasive fast fractional flow reserve
CN105997120A (zh) 前a%区域作参考区的人脑局部水分布容积测定方法
Ahmed et al. New approach for simultaneous respiratory and cardiac motion correction in cardiac PET (NAMC-CPET)
Li et al. Reproducible computer-assisted quantification of myocardial perfusion with contrast-enhanced ultrasound
JPWO2016009957A1 (ja) コンピュータプログラム、画像処理装置及び方法
Wang et al. Non-invasive quantification of the brain [18 F] FDG-PET using inferred blood input function learned from total-body data with physical constraint
Barna et al. A direct method for reconstructing dynamic SPECT images
JP2020094905A (ja) 心筋核医学画像データのスコアリング
Cui et al. Comparison of parameter estimations using dual-Input and arterial-input in liver kinetic studies of FDG metabolism
Slomka et al. Quantification of myocardial perfusion
Evans et al. Comparison of first pass bolus AIFs extracted from sequential 18F-FDG PET and DSC-MRI of mice
Mateos-Pérez et al. Automatic TAC extraction from dynamic cardiac PET imaging using iterative correlation from a population template

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
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