CN108742667A - 基于身体质量指数获取血流特征值的方法和装置 - Google Patents
基于身体质量指数获取血流特征值的方法和装置 Download PDFInfo
- Publication number
- CN108742667A CN108742667A CN201810637733.4A CN201810637733A CN108742667A CN 108742667 A CN108742667 A CN 108742667A CN 201810637733 A CN201810637733 A CN 201810637733A CN 108742667 A CN108742667 A CN 108742667A
- Authority
- CN
- China
- Prior art keywords
- interest
- area
- model
- mass index
- cross
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 230000017531 blood circulation Effects 0.000 claims abstract description 79
- 230000000877 morphologic effect Effects 0.000 claims abstract description 49
- 239000008280 blood Substances 0.000 claims abstract description 37
- 210000004369 blood Anatomy 0.000 claims abstract description 36
- 238000006243 chemical reaction Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 230000008859 change Effects 0.000 claims description 19
- 230000002792 vascular Effects 0.000 claims description 19
- 230000003902 lesion Effects 0.000 claims description 16
- 239000000203 mixture Substances 0.000 claims description 12
- 210000004204 blood vessel Anatomy 0.000 claims description 6
- 239000012530 fluid Substances 0.000 claims description 6
- 230000037396 body weight Effects 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000001435 haemodynamic effect Effects 0.000 claims description 3
- 230000010412 perfusion Effects 0.000 claims description 3
- 210000005166 vasculature Anatomy 0.000 abstract 1
- 210000004351 coronary vessel Anatomy 0.000 description 15
- 238000010586 diagram Methods 0.000 description 9
- 230000006872 improvement Effects 0.000 description 9
- 238000005259 measurement Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- OIRDTQYFTABQOQ-KQYNXXCUSA-N adenosine Chemical compound C1=NC=2C(N)=NC=NC=2N1[C@@H]1O[C@H](CO)[C@@H](O)[C@H]1O OIRDTQYFTABQOQ-KQYNXXCUSA-N 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 201000000057 Coronary Stenosis Diseases 0.000 description 3
- 208000029078 coronary artery disease Diseases 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000000004 hemodynamic effect Effects 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 230000002107 myocardial effect Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000002126 C01EB10 - Adenosine Substances 0.000 description 2
- 208000031481 Pathologic Constriction Diseases 0.000 description 2
- 229960005305 adenosine Drugs 0.000 description 2
- 210000000709 aorta Anatomy 0.000 description 2
- 238000005452 bending Methods 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000035487 diastolic blood pressure Effects 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 238000002608 intravascular ultrasound Methods 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 230000036262 stenosis Effects 0.000 description 2
- 208000037804 stenosis Diseases 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 238000002586 coronary angiography Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 125000003473 lipid group Chemical group 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000001404 mediated effect Effects 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 230000004218 vascular function Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/504—Apparatus 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 blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/46—Arrangements for interfacing with the operator or the patient
- A61B6/467—Arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/469—Arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis 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 or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus 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/507—Apparatus 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 determination of haemodynamic parameters, e.g. perfusion CT
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Surgery (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Computer Graphics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Vascular Medicine (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Architecture (AREA)
- Human Computer Interaction (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供了一种基于身体质量指数获取血流特征值的方法和装置,所述基于身体质量指数获取血流特征值的方法包括:基于个体的身体质量指数调节X射线发出装置的工作参数,以获取个体感兴趣区域特定的几何参数,并建立几何模型;基于所述几何模型建立感兴趣区域的形态差异函数f(x);获取个体冠脉系统在感兴趣区域的血流模型,并结合个体在感兴趣区域的几何模型,获取感兴趣区域的血流速度V、血流压力Pa和压力差数值ΔP等反应血流特征的数值。本发明提供的基于身体质量指数获取血流特征值的方法和装置,通过引入形态学的概念,并结合个体的身体质量指数,可保证几何模型建立的准确性,进一步达到提高血流特征值计算结果准确性的目的。
Description
技术领域
本发明涉及一种基于身体质量指数获取血流特征值的方法及装置,属于医疗技术领域。
背景技术
人体血液中的脂类及糖类物质在血管壁上的沉积将在血管壁上形成斑块,继而导致血管狭窄;特别是发生在心脏冠脉附近的血管狭窄将导致心肌供血不足,诱发冠心病、心绞痛等病症,对人类的健康造成严重威胁。据统计,我国现有冠心病患者约1100万人,心血管介入手术治疗患者数量每年增长大于10%。
冠脉造影CAG、计算机断层扫描CT等常规医用检测手段虽然可以显示心脏冠脉血管狭窄的严重程度,但是并不能准确评价冠脉的缺血情况。为提高冠脉血管功能评价的准确性,1993年Pijls提出了通过压力测定推算冠脉血管功能的新指标——血流储备分数(Fractional Flow Reserve,FFR),经过长期的基础与临床研究,FFR已成为冠脉狭窄功能性评价的金标准。
血流储备分数(FFR)通常是指心肌血流储备分数,定义为病变冠脉能为心肌提供的最大血流与该冠脉完全正常时最大供血流量之比,研究表明,在冠脉最大充血状态下,血流量的比值可以用压力值来代替。即FFR值的测量可在冠脉最大充血状态下,通过压力传感器对冠脉远端狭窄处的压力和冠脉狭窄近端压力进行测定继而计算得出。近年来,基于压力导丝测量FFR值的方法逐渐进入临床应用,成为冠心病患者获得精准诊断的有效方法;然而,由于压力导丝在介入过程中易对病人的血管造成损伤;同时,通过压力导丝对FFR值进行测定需要注射腺苷/ATP等药物保证冠脉达到最大充血状态,部分病人会因药物的注射感到不适,使得基于压力导丝测量FFR值的方法存在较大的局限性。此外,虽然基于压力导丝引导的FFR的测定是冠脉狭窄血液动力学的重要指标,但是由于压力导丝的造价高,介入血管过程操作困难,因此严重限制了基于压力导丝测量FFR值的方法的推广及使用。
随着CT与三维造影重建技术的发展及3D冠状动脉几何重建技术在血液力学研究领域的推广应用,同时,为减少FFR值测量过程中对人体带来的伤害及测量成本,基于医疗影像学的FFR计算技术已成为研究重点。
现有技术中,Taylor等人将计算机流体力学应用于计算机断层扫描冠状动脉造影(CTA)中,利用CTA得到冠脉解剖数据,包括血管供应心肌的体积和质量等,估算出最大冠脉血流量,模拟出血管下游微循环阻力,作为计算流体力学仿真的边界条件进行流体方程求解,得到计算FFR的非侵入式方法FFRCT。
事实上,现有技术虽然从不同角度、不同方法中给出了确定血流储备分数(FFR)的方法,但其实质均是通过感兴趣区域近端终点处的血流压力Pa和感兴趣区域近端终点处和远端终点处的血流压力的差值ΔP来计算FFR。而在血液流动的实际过程中,即血流压力的差值ΔP的实际计算过程中,病变的位置、大小和类型等因素均会对血流压力的差值ΔP的计算产生影响;现有技术中,病变的位置、大小和类型可通过建立感兴趣区域的几何模型获得,特别地,在使用X射线装置获取感兴趣区域病变的几何参数时,几何参数的准确性与个体的身体质量指数相关,即当使用相同的工作参数通过X射线装置获取个体的几何参数时,肥胖个体获取的几何参数偏差大,无法准确反映个体在几何形态,致使通过几何参数建立的几何模型的偏差大进一步使得血流特征值的计算出现偏差。
有鉴于此,确有必要提供一种新的基于身体质量指数获取血流特征值的方法,以解决上述问题。
发明内容
本发明的目的在于提供一种基于身体质量指数获取血流特征值的方法和装置,以至少解决现有技术中存在的技术问题之一。本发明提供的基于身体质量指数获取血流特征值的方法,通过引入形态学的概念,并结合个体的身体质量指数,可准确的获得个体感兴趣区域的几何参数,保证几何模型建立的准确性,进一步达到提高血流特征值计算结果准确性的目的。
为实现上述发明目的,本发明提供了一种基于身体质量指数获取血流特征值的方法,包括以下步骤:
获取个体的特异性身体质量指数,根据所述身体质量指数调节X射线发出装置的工作参数,以获取个体感兴趣区域特定的几何参数;
根据所述几何参数建立感兴趣区域的几何模型,根据所述几何模型建立感兴趣区域在近端终点和远端终点之间各个位置处的横截面形态模型;并依据所述横截面形态模型建立感兴趣区域的形态差异函数f(x);
获取个体在感兴趣区域特定的血流参数,并结合个体在感兴趣区域的几何模型,获取感兴趣区域的血流模型,并根据所述血流模型获取感兴趣区域近端终点处的血流压力Pa;及相应状态下感兴趣区域近端终点处与远端终点处的压力差数值ΔP;
所述血流特征值包括所述血流压力Pa、压力差数值ΔP以及基于所述血流压力Pa、压力差数值ΔP和形态差异函数f(x)计算的反应血流特征的数值。
作为本发明的进一步改进,所述工作参数包括X射线发出装置中X射线管的管电压和/或管电流。
作为本发明的进一步改进,所述身体质量指数为个体体重与身高的平方的比值;当所述个体的身体质量指数<18.5时,所述X射线管的管电压为 80~120kV;当所述个体的身体质量指数在18.5~23.9之间时,所述X射线管的管电压为120kV;当所述个体的身体质量指数>24时,所述X射线管的管电压为120~140kV。
作为本发明的进一步改进,所述横截面形态模型的建立包括:
S1、定义感兴趣区域近端终点处的横截面为参考面,通过中心线提取与建立方法获得的所述几何模型的中心径线;
S2、以所述参考面的中心点为原点建立坐标系,沿垂直所述中心径线的方向对所述感兴趣区域进行分割,将各横截面内外边缘投影在所述坐标系中,以获取感兴趣区域在各个位置处管腔横截面的平面几何图像,横截面形态模型建立结束;
其中,所述横截面形态模型包括各横截面上斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化。
作为本发明的进一步改进,所述形态差异函数f(x)用于表示感兴趣区域不同位置处的横截面形态变化随着该位置到参考点的距离x变化的函数;所述形态差异函数f(x)的获取包括:
基于所述横截面形态模型,建立各横截面的形态函数;
对相邻两横截面的形态函数进行拟合,并获取相邻两横截面在不同尺度下的差异变化函数;
以感兴趣区域的近端终点为参考点,根据差异变化函数获取管腔形态随着到参考点的距离x的变化率,对感兴趣区域从近端终点到远端终点范围内的位置参数进行归一化处理,获取形态差异函数f(x)。
作为本发明的进一步改进,所述压力差数值ΔP通过感兴趣区域血管管腔在不同尺度下的形态差异函数f(x)和血流速度V计算获得,所述ΔP在不同尺度下的计算公式为:
ΔP=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
其中,V为血流速度,为通过所述血流模型直接/间接获得,且所述血流速度V可以为常数;
c1、c2、…、cm分别代表血流速度V的参数系数;
α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数f1(x)、 f2(x)、…、fn(x)的加权系数;
m为大于等于1的自然数;
n为尺度为大于等于1的自然数。
优选的,所述不同尺度包括第一尺度、第二尺度、……、第n尺度;
所述第一尺度形态差异函数f1(x)用于检测第一种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;
所述第二尺度形态差异函数f2(x)用于检测第二种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;
……
所述第n尺度形态差异函数fn(x)用于检测第n种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;其中,所述n为大于等于1的自然数。
为实现上述发明目的,本发明还提供了一种基于身体质量指数获取血流特征值的装置,所述基于身体质量指数获取血流特征值的装置包括:
图像采集器,所述图像采集器用于接收个体的身体质量指数,并根据所述身体质量指数调节图像采集器的工作参数,获取个体感兴趣区域特定的几何参数;
数据采集器,所述数据采集器用于接收、传递及存储图像采集器传递的几何参数;
血流特征处理器,所述血流特征处理器用于获取感兴趣区域的血流模型,和基于所述几何参数建立对应感兴趣区域的几何模型;
所述血流特征处理器还用于通过所述血流模型获取感兴趣区域近端终点处的血流压力Pa=,并基于所述几何模型和所述血流模型,建立横截面形态模型以获取形态差异函数f(x);并根据所述形态差异函数f(x)、血流压力Pa、血流速度V和血流动力学,获取所述感兴趣区域的血流特征值。
作为本发明的进一步改进,所述图像采集器包括X射线发出装置,所述图像采集器可用于X射线发出装置工作参数的调节,且所述工作参数为X射线发出装置中X射线管的管电压和/或管电流。
作为本发明的进一步改进,所述身体质量指数为个体体重与身高的平方的比值;当所述个体的身体质量指数<18.5时,所述图像采集器调节所述X 射线管的管电压为80~120kV;当所述个体的身体质量指数在18.5~23.9之间时,所述图像采集器调节所述X射线管的管电压为120kV;当所述个体的身体质量指数>24时,所述图像采集器调节所述X射线管的管电压为 120~140kV。
作为本发明的进一步改进,所述几何模型为所述血流特征处理器通过接收所述数据采集器传递的个体感兴趣区域特定的几何参数,对所述几何参数进行测算,并拟合校准获得;所述横截面形态模型为通过所述几何模型直接 /间接获得;其中,所述横截面形态模型包括由所述血流特征处理器基于所述数据采集器获取板块信息,所述斑块信息各横截面上斑块的有无、斑块的位置、斑块的大小、斑块的组成、斑块组成的变化、斑块的形状及斑块形状的变化。
作为本发明的进一步改进,所述基于身体质量指数获取血流特征值的装置还包括速度采集器,所述速度采集器用于获取感兴趣区域的血流速度V;所述速度采集器包括速度计算模块及速度提取模块;所述速度提取模块可通过所述数据采集器直接获取血流速度V的信息,也可通过所述血流模型直接提取血流速度V;所述速度计算器包括速度转换模块及速度测算模块,所述血流速度V可通过血管中流体充盈的速度经所述速度转换模块转换获得,还可通过几何模型中血管树的形态经速度测算模块计算获得;优选的,所述血管树的形态至少包括所述血管树的面积、体积和中血管段的管腔直径中的一种或几种;所述血流速度V通过所述血管树的形态计算获得时,所述几何参数还包括所述血管树中血管段的长度、灌注面积及分支角度中的一种或几种。
本发明的有益效果是:本发明的基于身体质量指数获取血流特征值的方法和装置通过建立横截面形态模型,并通过对不同位置处的横截面形态模型进行拟合建立形态差异函数f(x);进一步在几何模型获取的过程中引入身体质量指数对几何模型的准确性进行调解,综合考虑了血管管腔中斑块的位置、形状及身体质量指数对血流特征值计算的影响;使得通过本发明的基于身体质量指数获取血流特征值的方法计算得到的血流特征值更加准确,可准确反映感兴趣区域的血流特征变化,且结果准确可靠。
附图说明
图1是本发明感兴趣区域的一种形态下的几何模型的示意图。
图2是图1中D1位置处横截面形态模型的结构示意图。
图3是图1中D2位置处横截面形态模型的结构示意图。
图4是图2和图3中D1和D2位置处横截面形态模型拟合后的结构示意图。
图5是本发明感兴趣区域的另一种形态下的几何模型的示意图。
图6是图5中D1位置处横截面形态模型的结构示意图。
图7是图5中D2位置处横截面形态模型的结构示意图。
图8是图6和图7中D1和D2位置处横截面形态模型拟合后的结构示意图。
图9是本发明基于身体质量指数获取血流特征值的装置的结构框图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图和具体实施例对本发明进行详细描述。
本发明提供了一种基于身体质量指数获取血流特征值的方法,所述基于身体质量指数获取血流特征值的方法包括以下步骤:
获取个体的特异性身体质量指数,根据所述身体质量指数调节X射线发出装置的工作参数,以获取个体感兴趣区域特定的几何参数;
根据所述几何参数建立感兴趣区域的几何模型,根据所述几何模型建立感兴趣区域在近端终点和远端终点之间各个位置处的横截面形态模型;并依据所述横截面形态模型建立感兴趣区域的形态差异函数f(x),所述尺度为计算形态差异函数f(x)时相邻两横截面之间的距离;
获取个体在感兴趣区域特定的血流参数,并结合个体在感兴趣区域的几何模型,获取感兴趣区域的血流模型,并根据所述血流模型获取=感兴趣区域近端终点处的血流压力Pa;及相应状态下感兴趣区域近端终点处与远端终点处的压力差数值ΔP;
所述血流特征值包括所述血流压力Pa、压力差数值ΔP以及基于所述血流压力Pa、压力差数值ΔP和形态差异函数f(x)计算的反应血流特征的数值。
以下说明书部分将对所述基于身体质量指数获取血流特征值的方法进行详细阐述。
所述身体质量指数为个体体重(kg)与身高(m)的平方的比值;在使用 X射线装置获取个体的几何参数时,相同的放射剂量作用在肥胖的个体时,将导致X射线装置获取的影像质量差,边缘位置影像噪声大,进一步致使获取的几何参数存在误差;因此,当个体的身体质量指数不同时,应调整X射线装置中X射线管的工作参数,以提高几何参数的准确性。
具体来讲,当所述个体的身体质量指数<18.5时,所述X射线管的管电压为80~120kV;当所述个体的身体质量指数在18.5~23.9之间时,所述X射线管的管电压为120kV;当所述个体的身体质量指数>24时,所述X射线管的管电压为120~140kV;如此设置,可在保证个体放射剂摄入量不变时,获得的几何参数准确,进一步提高几何模型建立的准确性。
进一步的,当个体的身体质量指数在24~27之间时,所述X射线管的管电压为120~130kV;当个体的身体质量指数在28~32之间时,所述X射线管的管电压为130~135kV;当个体的身体质量指数>32时,所述X射线管的管电压为135~140kV。
所述横截面形态模型为通过所述几何模型直接/间接获得,且在本发明中所述几何模型至少包括所述感兴趣区域的形状、直径和面积等几何参数,进一步的,所述几何参数还包括血管段的弯曲角度等可以反映感兴趣区域实际形态的参数。具体来讲,所述横截面形态模型的建立包括以下步骤:
S1、定义感兴趣区域近端终点处的横截面为参考面,通过中心线提取与建立方法,获得所述几何模型的中心径线;
S2、以所述参考面的中心点为原点建立坐标系,沿垂直所述中心径线的方向对所述感兴趣区域进行分割,将各横截面内外边缘投影在所述坐标系中,以获取感兴趣区域在各个位置处管腔横截面的平面几何图像,横截面形态模型建立结束。
其中,所述横截面形态模型包括各横截面位置处的斑块信息,所述斑块信息即为感兴趣区域的病变信息,且大量数据表明:当斑块(即为病变)的长度>20mm时,将导致感兴趣区域压力差数值ΔP的升高,进一步导致血流特征值如血流储备分数FFR的计算出现误差;而当同一横截面处斑块的组成复杂或尺寸过大致使感兴趣区域的狭窄率高,则会进一步导致感兴趣区域压力差数值ΔP的升高;同时当所述斑块处于不同的位置处时,目标血管所供应的心肌体积区域不同,将导致病变位置与非病变位置处的比例发生变化,进一步影响血流速度V,从而导致感兴趣区域血流特征值的计算发生偏差。
因此,在建立所述横截面形态模型时,所述斑块信息还需包括斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化,且在本发明中,各个位置处的管腔横截面的平面几何图像均需以步骤S2中建立的坐标系为参考,明确各横截面上斑块的位置,以方便横截面形态模型的后续拟合。
需要说明的是,在所述横截面形态模型的建立过程中,当所述解剖模型参数为采用CT、OCT、IVUS等检测手段获取时,所述横截面形态模型可通过所述几何模型直接获取,只需保证每个所述横截面形态模型的原点及坐标方向一致即可;当所述解剖模型参数为采用造影设备获取时,由于所述几何模型为沿血流方向延伸的立体模型,则在通过所述几何模型建立所述横截面形态模型时,需对所述几何模型进行坐标转换,以准确反应各个横截面的截面形态。
所述获取血管压力差的方法还包括对不同尺度下的所述横截面形态模型进行拟合,计算感兴趣区域血管管腔的形态差异函数f(x)。其中,所述形态差异函数f(x)用于表示感兴趣区域不同位置处的横截面形态变化随着该位置到参考点的距离x变化的函数;且所述形态差异函数f(x)的获取包括:
基于横截面形态模型,建立各横截面的形态函数;
对相邻两横截面的形态函数进行拟合,并获取相邻两横截面在不同尺度下的差异变化函数;
以感兴趣区域的近端终点为参考点,根据差异变化函数获取管腔形态随着到参考点的距离x的变化率,对感兴趣区域从近端终点到远端终点范围内的位置参数进行归一化处理,以最终获取形态差异函数f(x)。
所述形态函数包括面积函数、直径函数或边缘距离函数,即在本发明中可通过各横截面面积、直径或边缘距离函数之间的拟合,获取相邻两横截面在不同尺度下的差异变化函数;进一步的,通过差异变化函数获取管腔形态随着到参考点的距离x的变化率,获得形态差异函数f(x)。
具体来讲,当所述形态函数为面积函数时,如图1至图4,对D1和D2位置处的两横截面形态模型进行拟合,D1、D2位置处的横截面形态模型拟合后,有血管管腔斑块增加的区域为A1,对应的面积S1;血管管腔减少的区域为A2,对应的面积S2。由于所述D1和D2位置处的血管管腔(斑块)不重叠,因此当血流经D1处流向D2处时,血流压力将随之发生变化;此时,差异变化函数即为血管管腔中非重叠区域(S1、S2)与重叠区域之间面积(S3)的比值,或者为非重叠区域的面积(S1、S2)与总面积(S1、S2、S3)的比值;且此时,所述形态差异函数f(x)>0,即横截面D1和D2之间存在压力差。进一步的,当所述D1和D2位置处的血管管腔(斑块)完全重叠时,如图5至图8,所述区域A1与A2完全重叠,即非重叠区域A1与A2的面积S1=S2=0,此时,差异变化函数为0,即所述形态差异函数f(x)=0,此时,横截面D1和D2之间不存在压力差。
当所述形态函数为距离函数时,此时,确立选取的第一管腔边界上每个点与第二管腔边界上每个点的对应关系,然后求出第一管腔边界上的每个点与第二管腔边界上的每个点所对应的距离,减去沿着血管中心径线的距离,并获取所有点的距离之和或者是平均距离。具体来讲,若第一管腔边界与第二管腔边界的对应点到中心经线的距离均为y,则第一管腔与第二管腔的形态完全一致,即所述形态差异函数f(x)=0;若第一管腔边界与第二管腔边界的对应点到中心经线的距离不同,则第一管腔与第二管腔的形态不完全一致,即所述形态差异函数f(x)>0。
所述血流模型包括固定血流模型及个性化血流模型;其中所述固定血流模型即为经验值血流模型,为根据临床实际经验,通过大数据采集及模拟的方法直接建立;所述个性化血流模型包括静息态血流模型和负荷态血流模型。
具体来讲,在本发明中,所述血流模型中还包括感兴趣区域的血流速度V,且所述血流速度V和所述第一血流压力P1均通过所述血流模型直接/间接获取。进一步的,当所述血流模型为固定血流模型或静息态血流模型时,所述血流压力 Pa可通过获取的个体的收缩压混合舒张压计算获得,且此时所述血流压力 Pa=1/3收缩压+2/3舒张压;当所述血流模型为负荷态血流模型,此时所述血流压力Pa可通过所述负荷态血流模型直接测量获得。
在所述血流速度V的获取过程中,当所述血流模型为固定血流模型时,所述血流速度V可从固定血流模型中直接获取;当所述血流模型为静息态血流模型时,所述血流速度V可通过流体充盈的速度计算获得;当所述血流模型为负荷态血流模型,此时所述血流速度V为注射腺苷血管充分扩张后的血流速度V,且此时,所述血流速度V为最大血流速度Vmax。
特别地,在本发明中所述血流速度V包括感兴趣区域处于最大充血状态下的血流速度Vmax和静息状态下的血流速度Vqc,当感兴趣区域位于冠脉区域时,所述血流速度V为最大充血状态下的血流速度Vmax,进一步的血流速度Vmax可直接通过血流模型获取,或通过血流模型计算的血流速度V转换获得;当感兴趣区域位于外周血管系统时,所述血流速度V为静息态下的血流速度Vqc。
在本明中所述血流特征值中压力差数值ΔP为通过不同尺度下的形态差异函数f(x)和感兴趣区域的血流速度V计算获得,且所述压力差数值ΔP在不同尺度下的计算公式为:
ΔP=(c1V+c2V2+…+cmVm)
*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
其中,V为血流速度,为通过所述血流模型直接/间接获得,且所述血流速度V可以为常数;
c1、c2、…、cm分别代表血流速度的参数系数,所述参数系数包括血液粘度影响因数、血液湍流影响因数及粘滞系数等多个参数系数;进一步的, m为大于等于1的自然数,以分别代表不同参数系数对血流速度V的影响,以对压力差数值ΔP进行修正,保证压力差数值ΔP计算的准确性。优选的,在本发明中所述m的取值为2,且当所述m为2时,c1为因血液流动摩擦产生的参数系数,c2为血液湍流产生的参数系数。
所述α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数 f1(x)、f2(x)、…、fn(x)的加权系数,其中,n为尺度为大于等于1的自然数;进一步的,所述加权系数的增加可进一步对形态差异函数f(x)进行修正,保证两横截面之间形态差异拟合计算的准确性。
具体来讲,所述不同尺度包括第一尺度、第二尺度、……、第n尺度;
所述第一尺度形态差异函数f1(x)用于检测第一种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;
所述第二尺度形态差异函数f2(x)用于检测第二种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;
……
所述第n尺度形态差异函数fn(x)用于检测第n种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异。
请参参阅图9所示,本发明还提供了一种基于身体质量指数获取血流特征值的装置,所述基于身体质量指数获取血流特征值的装置包括:
图像采集器,所述图像采集器用于接收个体的身体质量指数,并根据所述身体质量指数调节图像采集器的工作参数,获取个体感兴趣区域特定的几何参数;
数据采集器,所述数据采集器用于接收、传递及存储图像采集器传递的几何参数;
血流特征处理器,所述血流特征处理器用于获取感兴趣区域的血流模型,和基于所述几何参数建立对应感兴趣区域的几何模型;
所述血流特征处理器还用于通过所述血流模型获取感兴趣区域近端终点处的血流压力Pa,并基于所述几何模型和所述血流模型,建立横截面形态模型以获取形态差异函数f(x);并根据所述形态差异函数f(x)、血流压力Pa和血流动力学,获取所述感兴趣区域的血流特征值。
在本发明中,所述图像采集器包括X射线发出装置,所述图像采集器可用于X射线发出装置工作参数的调节,且所述工作参数为X射线发出装置中X射线管的管电压和/或管电流。
进一步的,在本发明中,当所述个体的身体质量指数<18.5时,所述图像采集器可调节所述X射线管的管电压为80~120kV;当所述个体的身体质量指数在18.5~23.9之间时,所述图像采集器可调节所述X射线管的管电压为120kV;当所述个体的身体质量指数>24时,所述图像采集器可调节所述 X射线管的管电压为120~140kV。如此设置,可在保证个体放射剂摄入量不变时,获得的几何参数准确,进一步提高几何模型建立的准确性。
所述几何模型为所述血流特征处理器通过接收所述数据采集器传递的个体感兴趣区域特定的几何参数,对所述几何参数进行测算,并拟合校准获得;所述横截面形态模型为通过所述几何模型直接/间接获得。具体来讲,所述血流特征处理器获取的几何模型至少包括所述感兴趣区域的形状、直径和面积等几何参数,所述几何参数还包括血管段的弯曲角度等可以反映感兴趣区域实际形态的参数;即在本发明中,所述几何模型既可为单一血管段也可为血管树,且所述血管树包括至少一段主动脉或是至少一段主动脉和由所述主动脉发出的多个冠状动脉。
具体来讲,当所述几何参数的图像数据为通过CT、OCT和IVUS等设备获取时,所述数据采集器可直接对所述图像数据进行收集,并传递至所述血流特征处理器进行拟合建立几何模型;而当所述几何参数的图像数据为通过造影的方法获取时,所述数据采集器在对所述图像数据进行采集时,所述图像数据不少于两组,任意两组所述图像数据之间存在采集角度差,且所述采集角度差不小于20 度,如此设置,所述血流特征处理器获取的几何模型时,可保证几何模型的建立准确。
所述横截面形态模型为通过所述几何模型直接/间接获得,所述横截面形态模型包括各横截面上斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化。
进一步的,所述形态差异函数f(x)为所述血流特征处理器通过对任意两位置处的所述横截面形态模型进行拟合获得,且所述形态差异函数f(x)的获取包括:
通过所述血流特征处理器基于所述横截面形态模型,建立各横截面的形态函数;
对相邻两横截面的形态函数进行拟合,并获取相邻两横截面在不同尺度下的差异变化函数;
以感兴趣区域的近端终点为参考点,根据差异变化函数获取管腔形态随着到参考点的距离x的变化率,对感兴趣区域从近端终点到远端终点范围内的位置参数进行归一化处理,获取形态差异函数f(x)。
所述基于身体质量指数获取血流特征值的装置还包括速度采集器,所述速度采集器用于获取感兴趣区域的血流速度V,并将所述血流速度V传递至所述压力差处理器用以推算所述感兴趣区域近端终点与远端终点之间的任意两位置处之间的压力差数值ΔP。
在本发明中,所述速度采集器包括速度计算模块及速度提取模块;所述速度提取模块可通过所述数据采集器直接获取血流速度V的信息,也可通过所述血流模型直接提取血流速度V。
所述速度计算器包括速度转换模块及速度测算模块,所述血流速度V可通过血管中流体充盈的速度经所述速度转换模块转换获得,还可通过几何模型中血管树的形态经速度测算模块计算获得。
在本发明中,所述血管树的形态至少包括所述血管树的面积、体积和血管树中血管段的管腔直径中的一种或几种;所述血流速度V通过所述血管树的形态计算获得时,所述几何参数还包括所述血管树中血管段的长度、灌注面积及分支角度中的一种或几种。
在本发明的一个实施例中,所述血流特征值还包括所述感兴趣区域近端终点位置处和远端终点位置处之间的压力差数值ΔP,所述压力差数值ΔP通过感兴趣区域血管管腔在不同尺度下的形态差异函数f(x)和血流速度V计算获得;优选的,所述压力差数值ΔP通过如下公式计算获得:
ΔP=(c1V+c2V2+c3V3+…+cmVm)*(α1*∫f1(x)dx+α2*∫f2(x)dx
+…+αn*∫fn(x)dx)
其中,V为血流速度,为通过所述血流模型直接/间接获得,且所述血流速度V可以为常数;c1、c2、c3、…、cm为血流速度V的参数系数,此处的参数系数包括血液粘度影响因素、血液湍流影响因素及粘滞系数等多个参数系数。m 为大于等于1的自然数,以对压力差数值ΔP进行修正,保证压力差数值ΔP计算的准确性;优选的,本发明中m的取值为2,且当m为2时,c1为因血液流动摩擦产生的参数系数,c2为血液湍流产生的参数系数。
α1、α2...αn为不同尺度下形态差异函数f1(x),f2(x)…fn(x)的加权系数;所述加权系数的增加可进一步对形态差异函数f(x)进行修正,保证两横截面之间形态差异拟合计算的准确性。n为尺度为大于等于1的自然数。
进一步的,在本发明的另一实施例中,所述血流特征值的计算还可与所述感兴趣区域的血流速度V无关,以所述血流特征值为压力差数值ΔP为例,当所述血流特征值的计算与所述血流速度V无关时,所述ΔP在不同尺度下的计算公式为:
ΔP=k*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
其中,k为修正参数,且k为常数;进一步的,所述修正参数k为基于个体信息直接/间接获取的数值。
进一步的,所述血流特征值还包括所述感兴趣区域的血流储备分数FFR,所述血流储备分数通过感兴趣区域血管管腔在不同尺度下的形态差异函数f(x)、血流速度V和感兴趣区域近端终点位置处的血流压力Pa计算获得,且所述血流储备分数FFR通过如下公式计算获得:
需要指出的是,上述装置及功能模块仅仅为示例性的给出实现该技术方案的基本结构,而非唯一结构。
综上所述,本发明的基于身体质量指数获取血流特征值的方法和装置通过建立横截面形态模型,并通过对不同位置处的横截面形态模型进行拟合建立形态差异函数;进一步在几何模型获取的过程中引入身体质量指数对几何模型的准确性进行调解,综合考虑了血管管腔中斑块的位置、形状及身体质量指数对血流特征值计算的影响;使得通过本发明的基于身体质量指数获取血流特征值的方法计算得到的血流特征值更加准确,可准确反映感兴趣区域的血流特征变化,且结果准确可靠。
以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围。
Claims (11)
1.一种基于身体质量指数获取血流特征值的方法,其特征在于,包括以下步骤:
获取个体的特异性身体质量指数,根据所述身体质量指数调节X射线发出装置的工作参数,以获取个体感兴趣区域特定的几何参数;
根据所述几何参数建立感兴趣区域的几何模型,根据所述几何模型建立感兴趣区域在近端终点和远端终点之间各个位置处的横截面形态模型;并依据所述横截面形态模型建立感兴趣区域的形态差异函数f(x);
获取个体在感兴趣区域特定的血流参数,并结合个体在感兴趣区域的几何模型,获取感兴趣区域的血流模型,并根据所述血流模型获取感兴趣区域近端终点处的血流压力Pa;及相应状态下感兴趣区域近端终点处与远端终点处的压力差数值ΔP;
所述血流特征值包括所述血流压力Pa、压力差数值ΔP以及基于所述血流压力Pa、压力差数值ΔP和形态差异函数f(x)计算的反应血流特征的数值。
2.根据权利要求1所述的基于身体质量指数获取血流特征值的方法,其特征在于:所述工作参数包括X射线发出装置中X射线管的管电压和/或管电流。
3.根据权利要求2所述的基于身体质量指数获取血流特征值的方法,其特征在于:所述身体质量指数为个体体重与身高的平方的比值;当所述个体的身体质量指数<18.5时,所述X射线管的管电压为80~120kV;当所述个体的身体质量指数在18.5~23.9之间时,所述X射线管的管电压为120kV;当所述个体的身体质量指数>24时,所述X射线管的管电压为120~140kV。
4.根据权利要求1所述的基于身体质量指数获取血流特征值的方法,其特征在于:所述横截面形态模型的建立包括:
S1、定义感兴趣区域近端终点处的横截面为参考面,通过中心线提取与建立方法获得的所述几何模型的中心径线;
S2、以所述参考面的中心点为原点建立坐标系,沿垂直所述中心径线的方向对所述感兴趣区域进行分割,将各横截面内外边缘投影在所述坐标系中,以获取感兴趣区域在各个位置处管腔横截面的平面几何图像,横截面形态模型建立结束;
其中,所述横截面形态模型包括各横截面上斑块的有无、斑块的位置、斑块的大小、斑块形成的角度、斑块的组成及斑块组成的变化、斑块的形状及斑块形状的变化。
5.根据权利要求1所述的基于身体质量指数获取血流特征值的方法,其特征在于:所述形态差异函数f(x)用于表示感兴趣区域不同位置处的横截面形态变化随着该位置到参考点的距离x变化的函数;所述形态差异函数f(x)的获取包括:
基于所述横截面形态模型,建立各横截面的形态函数;
对相邻两横截面的形态函数进行拟合,并获取相邻两横截面在不同尺度下的差异变化函数;
以感兴趣区域的近端终点为参考点,根据差异变化函数获取管腔形态随着到参考点的距离x的变化率,对感兴趣区域从近端终点到远端终点范围内的位置参数进行归一化处理,获取形态差异函数f(x)。
6.根据权利要求1所述的基于身体质量指数获取血流特征值的方法,其特征在于:所述压力差数值ΔP通过感兴趣区域血管管腔在不同尺度下的形态差异函数f(x)和血流速度V计算获得,所述ΔP在不同尺度下的计算公式为:
ΔP=(c1V+c2V2+…+cmVm)*[α1*∫f1(x)dx+α2*∫f2(x)dx+…+αn*∫fn(x)dx]
其中,V为血流速度,为通过所述血流模型直接/间接获得,且所述血流速度V可以为常数;
c1、c2、…、cm分别代表血流速度V的参数系数;
α1、α2、…、αn分别为不同尺度下血管管腔的形态差异函数f1(x)、f2(x)、...、fn(x)的加权系数;
m为大于等于1的自然数;
n为尺度为大于等于1的自然数。
优选的,所述不同尺度包括第一尺度、第二尺度、......、第n尺度;
所述第一尺度形态差异函数f1(x)用于检测第一种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;
所述第二尺度形态差异函数f2(x)用于检测第二种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;
……
所述第n尺度形态差异函数fn(x)用于检测第n种病变特征所引起的相邻两横截面形态模型所对应的几何形态差异;其中,所述n为大于等于1的自然数。
7.一种基于身体质量指数获取血流特征值的装置,其特征在于,包括:
图像采集器,所述图像采集器用于接收个体的身体质量指数,并根据所述身体质量指数调节图像采集器的工作参数,获取个体感兴趣区域特定的几何参数;
数据采集器,所述数据采集器用于接收、传递及存储图像采集器传递的几何参数;
血流特征处理器,所述血流特征处理器用于获取感兴趣区域的血流模型,和基于所述几何参数建立对应感兴趣区域的几何模型;
所述血流特征处理器还用于通过所述血流模型获取感兴趣区域近端终点处的血流压力Pa,并基于所述几何模型和所述血流模型,建立横截面形态模型以获取形态差异函数f(x);并根据所述形态差异函数f(x)、血流压力Pa和血流动力学,获取所述感兴趣区域的血流特征值。
8.根据权利要求7所述的基于身体质量指数获取血流特征值的装置,其特征在于:所述图像采集器包括X射线发出装置,所述图像采集器可用于X射线发出装置工作参数的调节,且所述工作参数为X射线发出装置中X射线管的管电压和/或管电流。
9.根据权利要求8所述的基于身体质量指数获取血流特征值的装置,其特征在于:所述身体质量指数为个体体重与身高的平方的比值;当所述个体的身体质量指数<18.5时,所述图像采集器调节所述X射线管的管电压为80~120kV;当所述个体的身体质量指数在18.5~23.9之间时,所述图像采集器调节所述X射线管的管电压为120kV;当所述个体的身体质量指数>24时,所述图像采集器调节所述X射线管的管电压为120~140kV。
10.根据权利要求7所述的基于身体质量指数获取血流特征值的装置,其特征在于:所述几何模型为所述血流特征处理器通过接收所述数据采集器传递的个体感兴趣区域特定的几何参数,对所述几何参数进行测算,并拟合校准获得;所述横截面形态模型为通过所述几何模型直接/间接获得;其中,所述横截面形态模型包括由所述血流特征处理器基于所述数据采集器获取板块信息,所述斑块信息各横截面上斑块的有无、斑块的位置、斑块的大小、斑块的组成、斑块组成的变化、斑块的形状及斑块形状的变化。
11.根据权利要求10所述的基于身体质量指数获取血流特征值的装置,其特征在于:所述基于身体质量指数获取血流特征值的装置还包括速度采集器,所述速度采集器用于获取感兴趣区域的血流速度V;所述速度采集器包括速度计算模块及速度提取模块;所述速度提取模块可通过所述数据采集器直接获取血流速度V的信息,也可通过所述血流模型直接提取血流速度V;所述速度计算器包括速度转换模块及速度测算模块,所述血流速度V可通过血管中流体充盈的速度经所述速度转换模块转换获得,还可通过几何模型中血管树的形态经速度测算模块计算获得;优选的,所述血管树的形态至少包括所述血管树的面积、体积和中血管段的管腔直径中的一种或几种;所述血流速度V通过所述血管树的形态计算获得时,所述几何参数还包括所述血管树中血管段的长度、灌注面积及分支角度中的一种或几种。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810637733.4A CN108742667B (zh) | 2018-06-20 | 2018-06-20 | 基于身体质量指数获取血流特征值的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810637733.4A CN108742667B (zh) | 2018-06-20 | 2018-06-20 | 基于身体质量指数获取血流特征值的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108742667A true CN108742667A (zh) | 2018-11-06 |
CN108742667B CN108742667B (zh) | 2022-07-26 |
Family
ID=63979455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810637733.4A Active CN108742667B (zh) | 2018-06-20 | 2018-06-20 | 基于身体质量指数获取血流特征值的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108742667B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019242160A1 (zh) * | 2018-06-20 | 2019-12-26 | 博动医学影像科技(上海)有限公司 | 基于特定的生理参数获取血管压力值的方法及装置 |
CN113288087A (zh) * | 2021-06-25 | 2021-08-24 | 成都泰盟软件有限公司 | 一种基于生理信号的虚实联动实验系统 |
CN113855431A (zh) * | 2021-11-02 | 2021-12-31 | 毕威泰克(浙江)医疗器械有限公司 | 防褥疮调节监测系统及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104068886A (zh) * | 2013-03-28 | 2014-10-01 | 上海联影医疗科技有限公司 | 医学摄影系统及其摄影方法 |
CN105326486A (zh) * | 2015-12-08 | 2016-02-17 | 上海交通大学 | 血管压力差与血流储备分数的计算方法及系统 |
CN105518684A (zh) * | 2013-08-27 | 2016-04-20 | 哈特弗罗公司 | 用于预测冠状动脉病变的位置、开始、和/或变化的系统和方法 |
CN106073894A (zh) * | 2016-05-31 | 2016-11-09 | 博动医学影像科技(上海)有限公司 | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和系统 |
CN108742562A (zh) * | 2018-06-20 | 2018-11-06 | 博动医学影像科技(上海)有限公司 | 基于高血脂信息获取血管压力差的方法及装置 |
CN109009061A (zh) * | 2018-06-20 | 2018-12-18 | 博动医学影像科技(上海)有限公司 | 基于血压修正获取血流特征值的计算方法及装置 |
-
2018
- 2018-06-20 CN CN201810637733.4A patent/CN108742667B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104068886A (zh) * | 2013-03-28 | 2014-10-01 | 上海联影医疗科技有限公司 | 医学摄影系统及其摄影方法 |
CN105518684A (zh) * | 2013-08-27 | 2016-04-20 | 哈特弗罗公司 | 用于预测冠状动脉病变的位置、开始、和/或变化的系统和方法 |
CN105326486A (zh) * | 2015-12-08 | 2016-02-17 | 上海交通大学 | 血管压力差与血流储备分数的计算方法及系统 |
CN106073894A (zh) * | 2016-05-31 | 2016-11-09 | 博动医学影像科技(上海)有限公司 | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和系统 |
CN108742562A (zh) * | 2018-06-20 | 2018-11-06 | 博动医学影像科技(上海)有限公司 | 基于高血脂信息获取血管压力差的方法及装置 |
CN109009061A (zh) * | 2018-06-20 | 2018-12-18 | 博动医学影像科技(上海)有限公司 | 基于血压修正获取血流特征值的计算方法及装置 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019242160A1 (zh) * | 2018-06-20 | 2019-12-26 | 博动医学影像科技(上海)有限公司 | 基于特定的生理参数获取血管压力值的方法及装置 |
CN113288087A (zh) * | 2021-06-25 | 2021-08-24 | 成都泰盟软件有限公司 | 一种基于生理信号的虚实联动实验系统 |
CN113855431A (zh) * | 2021-11-02 | 2021-12-31 | 毕威泰克(浙江)医疗器械有限公司 | 防褥疮调节监测系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108742667B (zh) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105326486B (zh) | 血管压力差与血流储备分数的计算方法及系统 | |
CN109065170B (zh) | 获取血管压力差的方法及装置 | |
CN108186038B (zh) | 基于动脉造影影像计算冠脉血流储备分数的系统 | |
CN106073894B (zh) | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和系统 | |
CN106650029B (zh) | 基于cfd仿真的分叉血管压力差及ffr的快速计算方法及系统 | |
CN108742587B (zh) | 基于病史信息获取血流特征值的方法及装置 | |
CN111134651B (zh) | 基于腔内影像计算血流储备分数的方法、装置、系统以及计算机存储介质 | |
CN107978371A (zh) | 快速计算微循环阻力的方法及系统 | |
CN108717874A (zh) | 基于特定的生理参数获取血管压力值的方法及装置 | |
CN109009061A (zh) | 基于血压修正获取血流特征值的计算方法及装置 | |
WO2019169779A1 (zh) | 一种无需血管扩张剂测量血流储备分数的方法 | |
CN103932694A (zh) | 精确诊断心肌血流储备分数(ffr)的方法和设备 | |
CN108742667A (zh) | 基于身体质量指数获取血流特征值的方法和装置 | |
CN108742570A (zh) | 基于冠脉优势类型获取血管压力差的方法及装置 | |
US11039804B2 (en) | Apparatus and method for determining a fractional flow reserve | |
CN109620199A (zh) | 建立血管截面函数、血管压力差和血管应力的方法及装置 | |
CN109009037A (zh) | 基于个体性别获取血管压力差的计算方法及装置 | |
CN111091913A (zh) | 基于血流储备分数和冠状动脉ct造影图像的微循环阻力计算方法 | |
CN109044324A (zh) | 基于斑块位置修正血流特征值的方法及装置 | |
CN108784676A (zh) | 基于年龄信息获取压力差的方法及装置 | |
CN108777174A (zh) | 基于心梗史信息获取血管压力差的方法及装置 | |
CN108742586B (zh) | 基于糖尿病病史信息获取血流特征值的方法及装置 | |
CN108742562A (zh) | 基于高血脂信息获取血管压力差的方法及装置 | |
CN108742547B (zh) | 基于吸烟史信息获取压力差的方法及装置 | |
CN109662703A (zh) | 带有mems质量传感器的医用介入器械和血流参数测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: No. 11982, Qinhui North Road, Shanghai Applicant after: Shanghai Bodong Medical Technology Co.,Ltd. Address before: Room 402, building 2, 481 Guiping Road, Xuhui District, Shanghai 200233 Applicant before: PULSE MEDICAL IMAGING TECHNOLOGY (SHANGHAI) Co.,Ltd. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |