CN106073894B - 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和*** - Google Patents
基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和*** Download PDFInfo
- Publication number
- CN106073894B CN106073894B CN201610379167.2A CN201610379167A CN106073894B CN 106073894 B CN106073894 B CN 106073894B CN 201610379167 A CN201610379167 A CN 201610379167A CN 106073894 B CN106073894 B CN 106073894B
- Authority
- CN
- China
- Prior art keywords
- geometric parameter
- implantation
- virtual bracket
- diameter
- vessel
- 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
Links
- 238000002513 implantation Methods 0.000 title claims abstract description 173
- 230000017531 blood circulation Effects 0.000 title claims abstract description 76
- 238000000034 method Methods 0.000 title claims abstract description 73
- 230000002792 vascular Effects 0.000 title claims description 27
- 208000031481 Pathologic Constriction Diseases 0.000 claims abstract description 90
- 230000036262 stenosis Effects 0.000 claims abstract description 90
- 208000037804 stenosis Diseases 0.000 claims abstract description 90
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 86
- 238000009795 derivation Methods 0.000 claims abstract description 5
- 210000004369 blood Anatomy 0.000 claims description 56
- 239000008280 blood Substances 0.000 claims description 56
- 238000002586 coronary angiography Methods 0.000 claims description 32
- 230000003902 lesion Effects 0.000 claims description 27
- 238000004364 calculation method Methods 0.000 claims description 19
- 239000007943 implant Substances 0.000 claims description 17
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 14
- 210000004351 coronary vessel Anatomy 0.000 claims description 12
- 230000004089 microcirculation Effects 0.000 claims description 11
- 238000010606 normalization Methods 0.000 claims description 10
- 239000002872 contrast media Substances 0.000 claims description 9
- 230000002107 myocardial effect Effects 0.000 claims description 9
- 230000000004 hemodynamic effect Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 239000011800 void material Substances 0.000 claims description 4
- 238000002601 radiography Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 239000003795 chemical substances by application Substances 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 96
- 238000000205 computational method Methods 0.000 description 7
- 230000002966 stenotic effect Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 241000196324 Embryophyta Species 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000002583 angiography Methods 0.000 description 2
- 238000000151 deposition Methods 0.000 description 2
- 230000035487 diastolic blood pressure Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000004393 prognosis Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002608 intravascular ultrasound Methods 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 230000035488 systolic blood pressure Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/04—Measuring blood pressure
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0891—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/101—Computer-aided simulation of surgical operations
- A61B2034/102—Modelling of surgical devices, implants or prosthesis
- A61B2034/104—Modelling the effect of the tool, e.g. the effect of an implanted prosthesis or for predicting the effect of ablation or burring
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B34/00—Computer-aided surgery; Manipulators or robots specially adapted for use in surgery
- A61B34/10—Computer-aided planning, simulation or modelling of surgical operations
- A61B2034/101—Computer-aided simulation of surgical operations
- A61B2034/105—Modelling of the patient, e.g. for ligaments or bones
-
- 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/481—Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Physiology (AREA)
- Vascular Medicine (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Hematology (AREA)
- Dentistry (AREA)
- Quality & Reliability (AREA)
- Robotics (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Media Introduction/Drainage Providing Device (AREA)
- Prostheses (AREA)
Abstract
本发明提供了一种植入虚拟支架后血流储备分数的评估方法及***,该方法包括:接收一段血管的几何参数,几何参数包括第一、第二、第三几何参数;以近端终点为参考点,基于几何参数和血管段上的点到参考点的距离,计算出参考管腔直径函数和虚拟支架植入前的直径狭窄率;接收虚拟支架植入位置;基于第三几何参数、虚拟支架植入位置和参考管腔直径函数,计算出植入虚拟支架后虚拟血管管腔几何参数、虚拟血管几何参数差异函数和直径狭窄率;在不同尺度下对虚拟支架植入后几何参数差异函数求导,再基于多尺度差值导数函数及植入支架后最大平均血流速度,获得血流储备分数。本发明通过植入虚拟支架后直径狭窄率和最大平均血流速度的改变,利用多尺度计算方法,评估植入虚拟支架后的血流储备分数。
Description
技术领域
本发明涉及应用于医疗领域,尤其涉及应用在基于影像植入虚拟支架后获取血管压力差以及血流储备分数(FFR)的计算方法及***。
背景技术
血流储备分数(Fractional Flow Reserve,FFR)是一种重要的诊断冠脉生理功能的技术。FFR通过测量冠脉在最大充血状态下,狭窄远端与近端的压力之比来反应狭窄对最大血流的限制,从而判断狭窄是否诱发缺血。到目前已经有大量临床试验提供了循证医学依据,提示FFR可以用于评估狭窄病变功能意义,尤其是对于临界狭窄病变治疗方案的选择具有重要的指导意义。
通过压力传感器对血管进行有创侵入性压力测量不仅测量时间长,操作复杂,费用昂贵,而且存在着损伤血管的风险。“血管压力差与血流储备分数的计算方法及***”(申请号:CN201510901329.X,申请日:2015.12.08)专利通过冠脉造影获得的血管几何参数和血流速度,通过多尺度计算方法,实现了在同一个血管存在不同程度病变情况下的血管压力差计算,实现了正确区分评估不同程度狭窄下几何参数改变对血流压力的不同影响。
对于存在不同程度狭窄的血管,支架植入的位置、支架的大小、数量以及植入顺序的选择都是临床的关注点。目前,没有方法可以快速实时预测支架植入对血流储备分数的影响,而快速准确的评估植入支架后的血流储备分数又具有十分重要的临床意义。
本发明在本发明人之前申请的“血管压力差与血流储备分数的计算方法及***”专利的基础上进行了针对上述技术问题的创新和改进,实现了对存在不同程度病变情况下的血管植入虚拟支架,并计算和评估植入虚拟支架后血管的压力差与血流储备分数,为临床医生对于多病变血管支架植入策略的选择、评估预后提供了帮助。
发明内容
本发明所要解决的技术问题是提供一种基于植入虚拟支架的血管血流储备分数的评估方法和***,具体方案包括:
一种基于植入虚拟支架后评估血管压力降数值的方法,包括:接收一段血管的几何参数,该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,即实际血管管腔,代表未植入支架前该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;基于所述第一几何参数、第二几何参数、第三几何参数以及第一位置的位置数据计算该血管第一位置处的参考管腔直径;接收虚拟支架植入位置,所述植入位置是虚拟支架近端和远端到血管近端终点的距离。基于第三几何参数、虚拟支架植入位置和第一位置处的参考管腔直径,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数,代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;接收该血管段植入支架前的平均血流速度;基于第三几何参数和第一位置处的参考管腔直径,计算植入支架前第一位置处的直径狭窄率;基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率和几何参数差异;基于植入支架前的平均血流速度、支架植入前后的直径狭窄率,查相对血流速度与直径狭窄率表,获得虚拟支架植入后所能获得的平均血流速度;基于植入虚拟支架后第一位置处的几何参数差异、平均血流速度V和平均血流速度的平方V2,获得所述血管植入虚拟支架后近端终点处的第一血流压力和第一位置处的第二血流压力之间的压力降数值ΔP。
优选的,所述方法进一步还包括:以近端终点为参考点,根据第一几何参数、第二几何参数,以及血管段上某点到参考点的距离x,计算出参考管腔直径函数;所述参考管腔直径函数用于代表该段血管不同位置上的参考管腔直径随着该位置到参考点的距离x变化的函数。
优选的,计算参考管腔直径函数包括对血管段的从血管段近端终点到远端终点范围内的位置参数进行线性归一化处理。
优选的,所述方法进一步还包括:以近端终点为参考点,根据虚拟支架植入位置,将虚拟支架位置处的第三几何参数用该位置处参考管腔直径代替,计算出植入虚拟支架后虚拟血管管腔几何参数,即第四几何参数。
优选的,所述方法进一步还包括:以近端终点为参考点,根据第四几何参数和参考管腔直径函数,计算出植入虚拟支架后血管段的几何参数差异函数;所述几何参数差异函数是参考管腔直径函数和虚拟血管管腔几何参数(即第四几何参数)的差异随着到参考点的距离x变化的函数。
优选的,所述方法进一步还包括:根据第三几何参数和参考管腔直径函数,计算出植入虚拟支架之前血管段的直径狭窄率DS%;根据第四几何参数和参考管腔直径函数,计算出植入虚拟支架之后血管段的直径狭窄率DS%;所述直径狭窄率是该血管段上最严重病变的狭窄程度。
优选的,所述方法进一步还包括:根据该血管段植入支架前的平均血流速度V,支架植入前后的直径狭窄率DS%、DS%′,通过查相对血流速度与直径狭窄率表,计算出支架植入后的平均血流速度V′。所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度值无关。
优选的,所述方法进一步还包括:计算出虚拟支架植入后血管段的几何参数差异函数在n个尺度下的差值导数函数,其中植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP是根据所述n个尺度下的差值导数计算求得;所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离,所述n个尺度为具有尺度各不相同的第一尺度、第二尺度、……、第n尺度;
所述第一尺度差值导数函数f1(x)用于检测第一种病变特征所引起的真实管腔直径和参考管腔直径之间的几何参数差异;
所述第二尺度差值导数函数f2(x)用于检测第二种病变特征所引起的实管腔直径和参考管腔直径之间的几何参数差异;
……
所述第n尺度差值导数函数fn(x)用于检测第n种病变特征所引起的实管腔直径和参考管腔直径之间的几何参数差异;其中,所述n为大于1的自然数值。
优选的,所述方法进一步还包括:通过上述n个尺度下差值导数函数f1(x)、…、fn(x)积分的加权以及植入支架后平均血流速度V和平均血流速度V的平方V2,可计算获得植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP。
优选的,所述植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP的计算公式为:
ΔP=α1[C1V+C2V2]*∫f1(x)dx+α2[C1V+C2V2]∫f2(x)dx+…+αn[C1V+C2V2]*∫fn(x)dx
其中,C1、C2分别代表支架植入后平均血流速度V和平均血流速度平方V2的参数系数,α1、α2、…、αn分别为不同尺度的差值导数函数f1(x)、…、fn(x)的加权系数。
优选的,所述第一位置的位置数据是第一位置到血管近端终点之间的距离;所述血管段的平均血流速度是血流从近端终点到远端终点之间的平均速度。
优选的,所述方法进一步还包括:接收某一个角度下的二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点。
优选的,所述方法进一步还包括:计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数。
优选的,所述方法进一步还包括:从灰度时间拟合函数中得到该血管段在冠脉造影过程中造影剂的平均流动速度。
优选的,该血管段植入虚拟支架之前的平均血流速度V近似等于从灰度时间拟合函数中计算的造影剂平均流动速度。
本发明还提供了一种基于植入虚拟支架的血管血流储备分数FFR的评估方法,包括:
获得病人未植入虚拟支架之前静息状态下的血管段平均血流速度V,可通过常规造影(无需微循环充分扩张)方法获取该血管段植入虚拟支架之前的平均血流速度V;基于该平均速度V获得病人植入虚拟支架之前心肌微循环充分扩张时的最大血流速度Vmax;基于植入虚拟支架之前最大血流速度Vmax,支架植入前后的直径狭窄率,通过查相对血流速度与直径狭窄率表,获得虚拟支架植入后最大血流速度V′max。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度无关。求解出上述虚拟支架植入后最大血流速度对应的压力降ΔPmax;
利用下面公式获得虚拟支架植入后FFR数值:FFR=(P1-ΔPmax)/P1;其中,P1代表某段血管的近端终点处的第一血流压力,可用心脏的舒张压与收缩压近似估计,或者通过造影导管进行准确测量获得所述P1的值。
优选的,该方法进一步包括:通过查对应表获得植入虚拟支架之前最大血流速度,所述对应表是病人在静息状态下冠脉的平均血流速度与心肌微循环充分扩张情况下对应的最大血流速度列表。
优选的,该方法进一步包括:采用前一技术方案中所述的方法获得某段血管的最大血流速度对应的压力差ΔPmax。
优选的,给定一个固定的最大血流速度Vmax,求解出FFR。
同时,本发明还提供了一种基于植入虚拟支架后血管压力降数值的评估***,所述***能够实现如上所述的方法,该***包括:几何参数数据获取模块,用于获得血管的几何参数,该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,即实际血管管腔,代表未植入虚拟支架前该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;位置数据获取模块,用于获取上述第一位置的位置数据;参考管腔直径计算模块,基于所述第一几何参数、第二几何参数、第三几何参数以及第一位置的位置数据计算该血管第一位置处的参考管腔直径;虚拟血管管腔计算模块,接收虚拟支架植入位置,并根据第三几何参数和参考管腔直径,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数,代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径。直径狭窄率计算模块,基于第三几何参数和第一位置处的参考管腔直径,计算未植入虚拟支架前第一位置处的直径狭窄率,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率;几何参数差异计算模块,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后虚拟血管第一位置处的几何参数差异;速度获取和计算模块,用于获得该血管段未植入虚拟支架之前的平均血流速度并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度;压力降计算模块,获得由几何参数差异计算模块输出的植入虚拟支架后第一位置处的几何参数差异数据、由速度获取模块输出的植入虚拟支架后平均血流速度以及平均血流速度平方值,计算出所述植入虚拟支架后血管近端终点处的第一血流压力和第一位置处的第二血流压力之间的压力降数值ΔP。
优选的,所述参考管腔直径计算模块,以近端终点为参考点,根据第一几何参数、第二几何参数,以及血管段上某点到参考点的距离x,计算出参考管腔直径函数;所述参考管腔直径函数用于代表该段血管不同位置上的参考管腔直径随着该位置到参考点的距离x变化的函数。
优选的,该***还包括归一化模块,用于在计算参考管腔直径函数过程中,对血管段的从血管段近端终点到远端终点范围内的位置参数进行线性归一化处理。
优选的,所述直径狭窄率是血管段上最严重病变的狭窄程度;
优选的,所述几何参数差异计算模块,以近端终点为参考点,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的几何参数差异。所述几何参数差异函数是参考管腔直径函数和虚拟血管管腔几何参数之间的差异随着到参考点的距离x变化的函数。
优选的,该***还包括多尺度差值导数计算模块,用于计算出植入虚拟支架后几何参数差异函数在不同尺度下的差值导数函数。
优选的,所述压力降计算模块,则通过多尺度差值导数计算模块获得的多个尺度下的差值导数函数积分的加权以及由速度获取模块输出的植入虚拟支架后平均血流速度V和平均血流速度V的平方V2,计算出植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP;其中,所述多个尺度包括两个或两个以上的尺度,所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离。
优选的,该***还包括二维冠脉造影模块,用于对血管的某一个角度下进行二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点。
优选的,所述速度获取与计算模块,基于上述二维冠脉造影模块输出的不同帧数下感兴趣区域进行图像配准数据,计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数,并从灰度时间拟合函数中得到该血管段在未植入虚拟支架前在冠脉造影过程中造影剂的平均流动速度。并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度。
优选的,所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表,其存储在平均血流速度获取和计算模块或者该***的另外一独立的存储模块中,本发明并不以此为限制。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度值无关。
此外,本发明还提供了一种基于植入虚拟支架的血管血流储备分数FFR评估***,包括:几何参数数据获取模块,用于获得该段血管的几何参数,其中该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,代表该未植入虚拟支架前血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;位置数据获取模块,用于获取上述第一位置的位置数据;参考管腔直径计算模块,以近端终点为参考点,基于所述几何参数和血管段上的点到参考点的距离,计算出参考管腔直径函数;虚拟血管管腔计算模块,接收虚拟支架植入位置,并根据第三几何参数和参考管腔直径函数,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数,代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径。直径狭窄率计算模块,基于第三几何参数和第一位置处的参考管腔直径,计算未植入虚拟支架前第一位置处的直径狭窄率,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率;几何参数差异计算模块,基于参考管腔直径函数和第四几何参数,计算出植入虚拟支架后虚拟血管几何参数差异函数;多尺度计算模块,用于在多个尺度下对植入虚拟支架后几何参数差异函数求导,得到多个尺度对应的差值导数函数;所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离;平均血流速度获取和计算模块,获取该血管段在常规冠脉造影的平均血流速度,并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度;最大平均血流速度计算模块,基于存储在该模块中的对应表通过查表方式获得该血管段植入虚拟支架后所能通过的最大平均血流速度;血流储备分数计算模块,基于植入虚拟支架后多尺度差值导数函数和最大平均血流速度,获得所述血管植入虚拟支架后第一位置处的第二血流压力与近端终点处的第一血流压力之间的比值,即血流储备分数。
优选的,所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表,其存储在平均血流速度获取和计算模块或者该***的另外一独立的存储模块中。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度值无关。
优选的,所述最大血流速度获取模块通过查对应表获得最大血流速度,所述对应表是病人在静息状态下冠脉的平均血流速度与心肌微循环充分扩张情况下对应的最大血流速度列表,其存储在该最大血流速度获取模块或者该***的另外一独立的存储模块中。
本发明的有益效果在于:在“血管压力差与血流储备分数的计算方法及***”的基础上,实现了对存在不同程度病变情况下的血管植入虚拟支架,并计算和评估植入虚拟支架后血管的压力差与血流储备分数,为临床医生对于多病变血管支架植入策略的选择、评估预后提供了帮助。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1A为本发明的未植入虚拟支架前实际血管管腔与参考血管管腔结构对比示意图;
图1B为本发明的植入虚拟支架后实际血管管腔与参考血管管腔结构对比示意图;
图2为本发明的植入虚拟支架后血管段参考管腔和实际管腔几何参数差异函数示意图;
图3A为本发明的第一尺度差值导数函数f1(x)示意图;
图3B为本发明的第二尺度差值导数函数f2(x)示意图;
图4为本发明的植入虚拟支架后血管压降***结构示意图。
具体实施方式
下面结合附图对本发明实施例一种应用程序推荐方法及装置进行详细描述。应当明确,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本领域技术人员应当知晓,下述具体实施例或具体实施方式,是本发明为进一步解释具体的发明内容而列举的一系列优化的设置方式,而该些设置方式之间均是可以相互结合或者相互关联使用的,除非在本发明明确提出了其中某些或某一具体实施例或实施方式无法与其他的实施例或实施方式进行关联设置或共同使用。同时,下述的具体实施例或实施方式仅作为最优化的设置方式,而不作为限定本发明的保护范围的理解。
实施例1:
本发明提供一种基于植入虚拟支架后评估血管压力降数值的方法,包括:接收一段血管的几何参数,该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,即实际血管管腔,代表未植入支架前该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;基于所述第一几何参数、第二几何参数、第三几何参数以及第一位置的位置数据计算该血管第一位置处的参考管腔直径;接收虚拟支架植入位置,所述植入位置是虚拟支架近端和远端到血管近端终点的距离。基于第三几何参数、虚拟支架植入位置和第一位置处的参考管腔直径,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数。所述第四几何参数代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;接收该血管段植入支架前的平均血流速度;
以近端终点为参考点,根据第一几何参数、第二几何参数,以及血管段上某点到参考点的距离x,计算出参考管腔直径函数;所述参考管腔直径函数用于代表该段血管不同位置上的参考管腔直径随着该位置到参考点的距离x变化的函数;
在一具体实施例中,计算参考管腔直径函数包括对血管段的从血管段近端终点到远端终点范围内的位置参数进行线性归一化处理。
在一具体实施例中,以近端终点为参考点,根据虚拟支架植入位置,将虚拟支架位置处的第三几何参数用该位置处参考管腔直径代替,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数。
在一具体实施例中,以近端终点为参考点,根据第四几何参数和参考管腔直径函数,计算出植入虚拟支架后血管段的几何参数差异函数;所述几何参数差异函数是参考管腔直径函数和虚拟血管管腔几何参数(即第四几何参数)的差异随着到参考点的距离x变化的函数。
在一具体实施例中,根据第三几何参数和参考管腔直径函数,计算出植入虚拟支架之前血管段的直径狭窄率DS%;根据第四几何参数和参考管腔直径函数,计算出植入虚拟支架之后血管段的直径狭窄率DS%′;所述直径狭窄率DS%和DS%′分别代表支架植入前和植入后该血管段上最严重病变的狭窄程度。
在一具体实施例中,根据该血管段植入支架前的平均血流速度V,支架植入前后的直径狭窄率DS%、DS%′,通过查相对血流速度与直径狭窄率表,查相对血流速度与直径狭窄率表,计算出支架植入后的平均血流速度V′;所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度值无关。
在一具体实施例中,计算出虚拟支架植入后血管段的几何参数差异函数在n个尺度下的差值导数函数,其中植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP是根据所述n个尺度下的差值导数计算求得;所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离,所述n个尺度为具有尺度各不相同的第一尺度、第二尺度、……、第n尺度;
所述第一尺度差值导数函数f1(x)用于检测第一种病变特征所引起的真实管腔直径和参考管腔直径之间的几何参数差异,忽略掉其他病变所引起的几何参数差异;
所述第二尺度差值导数函数f2(x)用于检测第二种病变特征所引起的实管腔直径和参考管腔直径之间的几何参数差异;
……
所述第n尺度差值导数函数fn(x)用于检测第n种病变特征所引起的实管腔直径和参考管腔直径之间的几何参数差异;其中,所述n为大于1的自然数值。
在一具体实施例中,通过上述n个尺度下差值导数函数f1(x)、…、fn(x)积分的加权以及植入支架后平均血流速度V和平均血流速度的平方V2,可计算获得植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP。
优选的,所述植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP的计算公式为:
ΔP=α1[C1V+C2V2]*∫f1(x)dx+α2[C1V+C2V2]∫f2(x)dx+…+αn[C1V+C2V2]*∫fn(x)dx
其中,C1、C2分别代表支架植入后平均血流速度V和平均血流速度平方V2的参数系数,α1、α2、…、αn分别为不同尺度的差值导数函数f1(x)、…、fn(x)的加权系数。
优选的,所述第一位置的位置数据是第一位置到血管近端终点之间的距离;所述血管段的平均血流速度是血流从近端终点到远端终点之间的平均速度。
在一具体实施例中,接收某一个角度下的二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点。
优选的,所述方法进一步还包括:计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数。
优选的,所述方法进一步还包括:从灰度时间拟合函数中得到该血管段在冠脉造影过程中造影剂的平均流动速度。
优选的,该血管段植入虚拟支架之前的平均血流速度V近似等于从灰度时间拟合函数中计算的造影剂平均流动速度。
下面将结合图1-3给予所述的计算方法进一步阐述,参见图1A,所述方法包括:接收一段血管的几何参数,包括:①血管近端终点P处的几何参数(横截面积或直径);②血管远端终点D处的几何参数(横截面积或直径);③并以P点为参考点,得到位于P和D之间的血管点的几何参数(横截面积或直径),以及该点距离参考点P的距离x。
其中,可以通过多种方法获得几何参数,包括二维或三维的冠脉造影,计算机断层扫描冠状动脉造影(CTA),IVUS或者OCT。血管的几何参数主要考虑的是血管的横截面积和直径。如果获得的是二维血管直径,可以假设血管横截面为圆形,计算得到血管横截面积。
根据以上三种几何数据①②③,可以得到未植入虚拟支架前的血管参考(血管无病变时)几何参数。并将它表示为到参考点P距离的线性函数。如图1所示,实线为实际血管管腔,虚线为参考血管管腔。
接收虚拟支架植入位置,将支架位置处几何数据③的几何参数(横截面积或直径)用参考血管管腔替代,得到第四几何数据,代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径。并将它表示为到参考点P距离的线性函数,如图1B实线代表植入虚拟支架后的血管管腔。
如图2所示,表示图1B中植入虚拟支架后血管段参考管腔和实际管腔几何参数差异函数F(x)。需要说明的,经过分析发现,当在单一尺度下计算该几何参数差异函数F(x)时,所求得的第一血流压力和第二血流压力之间的压力差在所述血管为单一性病变情况下的数值是准确的,但是,当血管为多种病变情况并存时,尤其是弥漫性病变,采用上述方法所计算的压力差则出现了较大的误差,当所述单一尺度较小时,血管严重狭窄处的差值导数函数值与轻度狭窄处的差值导数函数值近似,这样会低估严重狭窄对血流压力差的影响;当所述单一尺寸较大时,血管轻度狭窄处的差值导数函数值为零,即无法检测出轻度狭窄对血流压力差的影响。
因此,为解决上述现有单一尺度下计算方法带来的技术缺陷,当血管为多种病变情况并存时,优选采用计算出几何参数差异函数在n个尺度下的差值导数函数,并根据不同尺度下的差值导数计算求得第一血流压力和第二血流压力之间的压力差。
例如,优选的一个实施例中,采用计算几何参数差异函数在两个尺度下的差值导数函数,其中第一血流压力和第二血流压力之间的压力差是根据所述两个尺度(包括具有较大尺度的第一尺度和具有较小尺度的第二尺度)下的差值导数计算求得;所述第一尺度差值导数函数f1(x)用于检测大范围狭窄所引起的真实管腔直径和参考管腔直径之间的几何参数差异,忽略掉局限性狭窄病变所引起的几何参数差值;所述第二尺度差值导数函数f2(x)用于检测局部变化所引起的实管腔直径和参考管腔直径之间的几何参数差异。
对图2中的差值函数F(x)在不同尺度下进行求导,得到不同尺度的差值导数函数。
大尺度差值导数函数:
小尺度差值导数函数:
其中,Δh1>Δh2。
如图3A和3B所示,在大尺度Δh1下,局限性狭窄病变A、B处的F(X+Δh1)-F(X)约为零,所以f1(x)用于反应大范围狭窄所引起的真实管腔直径和参考管腔直径之间的几何参数差值,而忽略掉局限性狭窄病变所引起的几何参数差异。在小尺度Δh2下,可以同时找出局限性狭窄病变A、B处以及大范围狭窄C的几何参数差异。但在小尺度差值导数函数f2(x),不同程度狭窄所引起的差值导数在函数值上基本相同,因此不能从差值导数f2(x)的函数值上反映出不同程度狭窄对压力差的不同影响。如果对不同尺度差值导数函数f1(x)和f2(x)进行加权,则可以体现出不同程度狭窄对压力差的不同影响。
注:上述所涉及的多尺度下计算方法可参见“血管压力差与血流储备分数的计算方法及***”(申请号:CN201510901329.X,申请日:2015.12.08)专利中相关技术方案。
实施例2
本发明还提供了一种基于植入虚拟支架的血管血流储备分数FFR的评估方法,包括:获得病人未植入虚拟支架之前静息状态下的血管段平均血流速度V,可通过常规造影(无需微循环充分扩张)方法获取该血管段植入虚拟支架之前的平均血流速度V;基于该平均速度V获得病人植入虚拟支架之前心肌微循环充分扩张时的最大血流速度Vmax;基于植入虚拟支架之前最大血流速度Vmax,支架植入前后的直径狭窄率,通过查相对血流速度与直径狭窄率表,获得虚拟支架植入后最大血流速度V′max。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度无关。求解出上述虚拟支架植入后最大血流速度对应的压力降ΔPmax;利用下面公式获得虚拟支架植入后FFR数值:FFR=(P1-ΔPmax)/P1;其中,P1代表某段血管的近端终点处的第一血流压力,可用心脏的舒张压与收缩压近似估计,或者通过造影导管进行准确测量获得所述P1的值。
优选的,通过查对应表获得植入虚拟支架之前最大血流速度,所述对应表是病人在静息状态下冠脉的平均血流速度与心肌微循环充分扩张情况下对应的最大血流速度列表。
优选的,采用实施例1中所述的方法获得某段血管的最大血流速度对应的压力差ΔPmax。
优选的,给定一个固定的最大血流速度Vmax,进而求解出FFR。
实施例3
在一个具体的实施例中,本发明提供了一种基于植入虚拟支架后血管压力降数值的评估***,所述***能够实现如上述实施例中示出的计算一段植入虚拟支架血管压力差的方法,该***包括:几何参数数据获取模块,用于获得血管的几何参数,该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,即实际血管管腔,代表未植入虚拟支架前该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;位置数据获取模块,用于获取上述第一位置的位置数据;参考管腔直径计算模块,基于所述第一几何参数、第二几何参数、第三几何参数以及第一位置的位置数据计算该血管第一位置处的参考管腔直径;虚拟血管管腔计算模块,接收虚拟支架植入位置,并根据第三几何参数和参考管腔直径,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数,所述第四几何参数代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径。直径狭窄率计算模块,基于第三几何参数和第一位置处的参考管腔直径,计算未植入虚拟支架前第一位置处的直径狭窄率,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率;几何参数差异计算模块,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后虚拟血管第一位置处的几何参数差异;速度获取和计算模块,用于获得该血管段未植入虚拟支架之前的平均血流速度并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度;压力降计算模块,获得由几何参数差异计算模块输出的植入虚拟支架后第一位置处的几何参数差异数据、由速度获取模块输出的植入虚拟支架后平均血流速度以及平均血流速度平方值,计算出所述植入虚拟支架后血管近端终点处的第一血流压力和第一位置处的第二血流压力之间的压力降数值ΔP。
在一个具体的实施例中,所述参考管腔直径计算模块,以近端终点为参考点,根据第一几何参数、第二几何参数,以及血管段上某点到参考点的距离x,计算出参考管腔直径函数;所述参考管腔直径函数用于代表该段血管不同位置上的参考管腔直径随着该位置到参考点的距离x变化的函数;
在一个具体的实施例中,该***还包括归一化模块,用于在计算参考管腔直径函数过程中,对血管段的从血管段近端终点到远端终点范围内的位置参数进行线性归一化处理。
在一个具体的实施例中,直径狭窄率是血管段上最严重病变的狭窄程度;
在一个具体的实施例中,所述几何参数差异计算模块,以近端终点为参考点,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的几何参数差异。所述几何参数差异函数是参考管腔直径函数和虚拟血管管腔几何参数之间的差异随着到参考点的距离x变化的函数。
在一个具体的实施例中,该***还包括多尺度差值导数计算模块,用于计算出植入虚拟支架后几何参数差异函数在不同尺度下的差值导数函数。
在一个具体的实施例中,所述压力降计算模块,则通过多尺度差值导数计算模块获得的多个尺度下的差值导数函数积分的加权以及由速度获取模块输出的植入虚拟支架后平均血流速度V和平均血流速度V的平方V2,计算出植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP;其中,所述多个尺度包括两个或两个以上的尺度,所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离。
在一个具体的实施例中,该***还包括二维冠脉造影模块,用于对血管的某一个角度下进行二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点。
在一个具体的实施例中,所述速度获取与计算模块,基于上述二维冠脉造影模块输出的不同帧数下感兴趣区域进行图像配准数据,计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数,并从灰度时间拟合函数中得到该血管段在未植入虚拟支架前在冠脉造影过程中造影剂的平均流动速度。并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度。
在一个具体的实施例中,所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表,其存储在平均血流速度获取和计算模块或者该***的另外一独立的存储模块中。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度值无关。
实施例4
在又一个具体的实施例中,本发明提供了一种基于植入虚拟支架的血管血流储备分数FFR评估***,所述***能够实现如上述实施例中示出的计算一段植入虚拟支架血管血流储备分数FFR的方法,该***包括:几何参数数据获取模块,用于获得该段血管的几何参数,其中该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,代表该未植入虚拟支架前血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;位置数据获取模块,用于获取上述第一位置的位置数据;参考管腔直径计算模块,以近端终点为参考点,基于所述几何参数和血管段上的点到参考点的距离,计算出参考管腔直径函数;虚拟血管管腔计算模块,接收虚拟支架植入位置,并根据第三几何参数和参考管腔直径函数,计算出植入虚拟支架后虚拟血管管腔几何参数,记作第四几何参数,其代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径。直径狭窄率计算模块,基于第三几何参数和第一位置处的参考管腔直径,计算未植入虚拟支架前第一位置处的直径狭窄率,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率;几何参数差异计算模块,基于参考管腔直径函数和第四几何参数,计算出植入虚拟支架后虚拟血管几何参数差异函数;多尺度计算模块,用于在多个尺度下对植入虚拟支架后几何参数差异函数求导,得到多个尺度对应的差值导数函数;所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离;平均血流速度获取和计算模块,获取该血管段在常规冠脉造影的平均血流速度,并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度;最大平均血流速度计算模块,基于存储在该模块中的对应表通过查表方式获得该血管段植入虚拟支架后所能通过的最大平均血流速度;血流储备分数计算模块,基于植入虚拟支架后多尺度差值导数函数和最大平均血流速度,获得所述血管植入虚拟支架后第一位置处的第二血流压力与近端终点处的第一血流压力之间的比值,即血流储备分数。
在一个具体的实施例中,所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表,其存储在平均血流速度获取和计算模块或者该***的另外一独立的存储模块中。所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0,与病人具体血流速度值无关。
在一个具体的实施例中,所述最大血流速度获取模块通过查对应表获得最大血流速度,所述对应表是病人在静息状态下冠脉的平均血流速度与心肌微循环充分扩张情况下对应的最大血流速度列表,其存储在该最大血流速度获取模块或者该***的另外一独立的存储模块中。
需要指出的是,上述***及功能模块仅仅为示例性的给出了其实现该方案的基本结构,然而并非唯一的结构。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (31)
1.一种基于植入虚拟支架后血管压力降数值的评估方法,包括:
接收一段血管的几何参数,该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,即实际血管管腔,代表未植入支架前该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;
基于所述第一几何参数、第二几何参数、第三几何参数以及第一位置的位置数据,计算该血管第一位置处的参考管腔直径;
接收虚拟支架植入位置,所述植入位置是虚拟支架近端和远端到血管近端终点的距离;
基于第三几何参数、虚拟支架植入位置和第一位置处的参考管腔直径,计算植入虚拟支架后虚拟血管管腔的第四几何参数,所述第四几何参数代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;
接收该血管段植入支架前的平均血流速度;
基于第三几何参数和第一位置处的参考管腔直径,计算植入支架前第一位置处的直径狭窄率;基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率和几何参数差异;
基于植入支架前的平均血流速度、支架植入前后的直径狭窄率,查相对血流速度与直径狭窄率表,获得虚拟支架植入后的平均血流速度;
基于植入虚拟支架后第一位置处的几何参数差异、血管平均血流速度V和平均血流速度的平方V2,获得所述血管植入虚拟支架后近端终点处的第一血流压力和第一位置处的第二血流压力之间的压力降数值ΔP。
2.根据权利要求1所述的方法,其特征在于,所述方法进一步还包括:
以近端终点为参考点,根据第一几何参数、第二几何参数,以及血管段上某点到参考点的距离x,计算出参考管腔直径函数;所述参考管腔直径函数代表该段血管不同位置上的参考管腔直径随着该位置到参考点的距离x变化的函数。
3.根据权利要求2所述的方法,其特征在于:计算参考管腔直径函数,包括对血管段的从血管段近端终点到远端终点范围内的位置参数进行线性归一化处理。
4.根据权利要求2所述的方法,其特征在于,所述方法进一步还包括:
所述计算植入虚拟支架后虚拟血管管腔的第四几何参数,进一步包括,以近端终点为参考点,根据虚拟支架植入位置,将虚拟支架位置处的第三几何参数用该位置处参考管腔直径代替。
5.根据权利要求4所述的方法,其特征在于,所述方法进一步包括:
以近端终点为参考点,根据第四几何参数和参考管腔直径函数,计算出植入虚拟支架后血管段的几何参数差异函数;所述几何参数差异函数是参考管腔直径函数和第四几何参数的差异随着到参考点的距离x变化的函数。
6.根据权利要求4所述的方法,其特征在于,所述方法进一步包括:
根据第三几何参数和参考管腔直径函数,计算出植入虚拟支架之前血管段的直径狭窄率DS%;根据第四几何参数和参考管腔直径函数,计算出植入虚拟支架之后血管段的直径狭窄率DS%'。
7.根据权利要求6所述的方法,其特征在于,所述方法进一步包括:
根据该血管段植入支架前的平均血流速度V,支架植入前后的直径狭窄率DS%、DS%',通过查相对血流速度与直径狭窄率表,计算出支架植入后的平均血流速度V'。
8.根据权利要求7所述的方法,其特征在于:
所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表;所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0。
9.根据权利要求4所述的方法,其特征在于,所述方法还包括:
计算出虚拟支架植入后血管段的几何参数差异函数在n个尺度下的差值导数函数f1(x)、…、fn(x),其中植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP是根据所述n个尺度下的差值导数计算求得;
所述n为大于1的自然数值;所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离。
10.根据权利要求9所述的方法,其特征在于,所述方法进一步包括:
通过上述n个尺度下差值导数函数f1(x)、…、fn(x)积分的加权以及植入支架后平均血流速度V和平均血流速度的平方V2,计算植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP。
11.根据权利要求9或10所述的方法,其特征在于,所述方法进一步包括:
所述植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP的计算公式为:
ΔP=α1[C1V+C2V2]*∫f1(x)dx+α2[C1V+C2V2]∫f2(x)dx+…+αn[C1V+C2V2]*∫fn(x)dx
其中,C1、C2分别代表支架植入后平均血流速度V和平均血流速度平方V2的参数系数,α1、α2、…、αn分别为不同尺度的差值导数函数f1(x)、…、fn(x)的加权系数。
12.根据权利要求1-10中任一所述的方法,其特征在于:所述第一位置的位置数据是第一位置到血管近端终点之间的距离;所述血管段的平均血流速度是血流从近端终点到远端终点之间的平均速度。
13.根据权利要求1-10中任一所述的方法,其特征在于:该方法还包括接收某一个角度下的二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点。
14.根据权利要求13所述的方法,其特征在于:所述方法还包括计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数。
15.根据权利要求14所述的方法,其特征在于:该方法还包括从灰度时间拟合函数中得到该血管段在冠脉造影过程中造影剂的平均流动速度。
16.根据权利要求15所述的方法,其特征在于,该血管段植入虚拟支架之前的平均血流速度V近似等于从灰度时间拟合函数中计算的造影剂平均流动速度。
17.一种基于植入虚拟支架的血管血流储备分数FFR的评估方法,包括:
获得未植入虚拟支架前静息状态下的血管段平均血流速度V;
基于所述平均速度V获得植入虚拟支架之前心肌微循环充分扩张时的最大血流速度Vmax;
基于所述最大血流速度Vmax、支架植入前后的直径狭窄率,通过查相对血流速度与直径狭窄率表,获得虚拟支架植入后最大血流速度V′max;
计算虚拟支架植入后最大血流速度V'max对应的压力降ΔPmax;
利用下面公式计算虚拟支架植入后FFR数值:
FFR=(P1-ΔPmax)/P1;
其中,P1代表某段血管的近端终点处的第一血流压力。
18.根据权利要求17所述的方法,其特征在于,该方法进一步包括:
通过查对应表获得植入虚拟支架之前最大血流速度,所述对应表是病人在静息状态下冠脉的平均血流速度与心肌微循环充分扩张情况下对应的最大血流速度列表。
19.一种基于植入虚拟支架后血管压力降数值的评估***,所述***包括:
几何参数数据获取模块,用于获得血管的几何参数,所述血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,代表未植入虚拟支架前该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;
位置数据获取模块,用于获取所述第一位置的位置数据;
参考管腔直径计算模块,基于所述第一几何参数、第二几何参数、第三几何参数以及第一位置的位置数据,计算该血管第一位置处的参考管腔直径;
虚拟血管管腔计算模块,接收虚拟支架植入位置,并根据第三几何参数和参考管腔直径,计算出植入虚拟支架后虚拟血管管腔的第四几何参数,代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;
直径狭窄率计算模块,基于第三几何参数和第一位置处的参考管腔直径,计算未植入虚拟支架前第一位置处的直径狭窄率,以及基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率;
几何参数差异计算模块,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后虚拟血管第一位置处的几何参数差异;
速度获取和计算模块,用于获得该血管段未植入虚拟支架之前的平均血流速度,并通过查相对血流速度与直径狭窄率表获取植入虚拟支架之后的平均血流速度;
压力降计算模块,获得所述几何参数差异数据、由速度获取和计算模块输出的植入虚拟支架后平均血流速度V以及平均血流速度平方值V2,计算出所述植入虚拟支架后血管近端终点处的第一血流压力和第一位置处的第二血流压力之间的压力降数值ΔP。
20.根据权利要求19所述的***,其特征在于,所述参考管腔直径计算模块,以近端终点为参考点,根据第一几何参数、第二几何参数,以及血管段上某点到参考点的距离x,计算出参考管腔直径函数;所述参考管腔直径函数用于代表该段血管不同位置上的参考管腔直径随着该位置到参考点的距离x变化的函数。
21.根据权利要求20所述的***,其特征在于,所述***还包括归一化模块,用于在计算参考管腔直径函数过程中,对血管段的从血管段近端终点到远端终点范围内的位置参数进行线性归一化处理。
22.根据权利要求19所述的***,其特征在于,所述直径狭窄率是血管段上最严重病变的狭窄程度。
23.根据权利要求20所述的***,其特征在于,所述几何参数差异计算模块,以近端终点为参考点,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的几何参数差异;所述几何参数差异函数是参考管腔直径函数和虚拟血管管腔几何参数之间的差异随着到参考点的距离x变化的函数。
24.根据权利要求23所述的***,其特征在于,该***还包括多尺度差值导数计算模块,用于计算出植入虚拟支架后几何参数差异函数在不同尺度下的差值导数函数;
所述压力降计算模块,则通过多尺度差值导数计算模块获得的多个尺度下的差值导数函数积分的加权以及由速度获取和计算模块输出的植入虚拟支架后平均血流速度V和平均血流速度的平方V2,计算出植入虚拟支架后第一血流压力和第二血流压力之间的压力降ΔP;
其中,所述多个尺度包括两个或两个以上的尺度,所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离。
25.根据权利要求19所述的***,其特征在于,该***还包括二维冠脉造影模块,用于对血管的某一个角度下进行二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点。
26.根据权利要求25所述的***,其特征在于,所述速度获取和计算模块,基于所述二维冠脉造影模块输出的图像配准数据,计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数,并从灰度时间拟合函数中得到该血管段在未植入虚拟支架前在冠脉造影过程中造影剂的平均流动速度,并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度。
27.据权利要求26所述的***,其特征在于,所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表;所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0。
28.一种基于植入虚拟支架的血管血流储备分数FFR评估***,其特征在于,所述***包括:
几何参数数据获取模块,用于获得该段血管的几何参数,其中该血管包括近端终点和远端终点,所述几何参数包括第一几何参数,代表该血管段近端横截面的面积或直径;第二几何参数,代表该血管段远端横截面的面积或直径;第三几何参数,代表未植入虚拟支架前血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;
位置数据获取模块,用于获取所述第一位置的位置数据;
参考管腔直径计算模块,以近端终点为参考点,基于所述几何参数和血管段上的点到参考点的距离,计算参考管腔直径函数;
虚拟血管管腔计算模块,接收虚拟支架植入位置,并根据第三几何参数和参考管腔直径函数,计算出植入虚拟支架后虚拟血管的第四几何参数,代表植入虚拟支架后该血管段位于近端终点和远端终点之间的第一位置的横截面面积或直径;
直径狭窄率计算模块,基于第三几何参数和第一位置处的参考管腔直径,计算未植入虚拟支架前第一位置处的直径狭窄率,基于第四几何参数和第一位置处的参考管腔直径,计算植入虚拟支架后第一位置处的直径狭窄率;
几何参数差异计算模块,基于参考管腔直径函数和第四几何参数,计算植入虚拟支架后虚拟血管几何参数差异函数;
多尺度计算模块,用于在多个尺度下对植入虚拟支架后几何参数差异函数求导,得到多个尺度对应的差值导数函数;所述尺度指的是分辨率,即数值计算导数时相邻两点之间的距离;
平均血流速度获取和计算模块,获取该血管段在常规冠脉造影的平均血流速度;并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度;
最大平均血流速度计算模块,通过查表方式获得该血管段植入虚拟支架后所能通过的最大平均血流速度;
血流储备分数计算模块,基于所述多尺度差值导数函数和最大平均血流速度,获得所述血管植入虚拟支架后第一位置处的第二血流压力与近端终点处的第一血流压力之间的比值,即为血流储备分数。
29.根据权利要求28所述的***,其特征在于,该***还包括二维冠脉造影模块,用于执行上述常规冠脉造影,所述二维冠脉造影模块用于对血管的某一个角度下进行二维冠脉造影,并对不同帧数下感兴趣区域进行图像配准;该血管段冠脉造影感兴趣区域是从血管近端终点到远端终点;
所述速度获取和计算模块,基于上述二维冠脉造影模块输出的不同帧数下感兴趣区域进行图像配准数据,计算配准后的感兴趣区域内灰度直方图随心动周期变化的灰度时间拟合函数,并从灰度时间拟合函数中得到该血管段在未植入虚拟支架前在冠脉造影过程中造影剂的平均流动速度;并通过查相对血流速度与直径狭窄率表计算植入虚拟支架之后的平均血流速度。
30.据权利要求29所述的***,其特征在于,所述相对血流速度与直径狭窄率表是直径狭窄率与对应的相对血流速度列表,所述相对血流速度是对血流速度归一化后的值,最大值为1,最小值为0。
31.根据权利要求30所述的***,其特征在于,所述最大血流速度获取模块通过查对应表获得最大血流速度,所述对应表是病人在静息状态下冠脉的平均血流速度与心肌微循环充分扩张情况下对应的最大血流速度列表。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610379167.2A CN106073894B (zh) | 2016-05-31 | 2016-05-31 | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和*** |
JP2018548257A JP6723670B2 (ja) | 2016-05-31 | 2016-11-04 | バーチャルステント留置に基づくffrの評価方法およびシステム |
DE112016005442.4T DE112016005442B4 (de) | 2016-05-31 | 2016-11-04 | Verfahren und System zum Evaluieren der FFR auf Grundlage virtueller Stentimplantation |
PCT/CN2016/104685 WO2017206438A1 (zh) | 2016-05-31 | 2016-11-04 | 基于植入虚拟支架的ffr评估方法及*** |
US15/779,704 US10617473B2 (en) | 2016-05-31 | 2016-11-04 | Method and system for evaluating FFR on the basis of virtual stent implantation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610379167.2A CN106073894B (zh) | 2016-05-31 | 2016-05-31 | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106073894A CN106073894A (zh) | 2016-11-09 |
CN106073894B true CN106073894B (zh) | 2017-08-08 |
Family
ID=57230656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610379167.2A Active CN106073894B (zh) | 2016-05-31 | 2016-05-31 | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和*** |
Country Status (5)
Country | Link |
---|---|
US (1) | US10617473B2 (zh) |
JP (1) | JP6723670B2 (zh) |
CN (1) | CN106073894B (zh) |
DE (1) | DE112016005442B4 (zh) |
WO (1) | WO2017206438A1 (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11278206B2 (en) | 2015-04-16 | 2022-03-22 | Gentuity, Llc | Micro-optic probes for neurology |
JP6981967B2 (ja) | 2015-08-31 | 2021-12-17 | ジェンテュイティ・リミテッド・ライアビリティ・カンパニーGentuity, LLC | 撮像プローブおよびデリバリデバイスを含む撮像システム |
CN105326486B (zh) * | 2015-12-08 | 2017-08-25 | 博动医学影像科技(上海)有限公司 | 血管压力差与血流储备分数的计算方法及*** |
CN106539622B (zh) * | 2017-01-28 | 2019-04-05 | 北京欣方悦医疗科技有限公司 | 基于血流动力学分析的冠脉虚拟支架植入*** |
US11684242B2 (en) | 2017-11-28 | 2023-06-27 | Gentuity, Llc | Imaging system |
CN107978371B (zh) * | 2017-11-30 | 2021-04-02 | 博动医学影像科技(上海)有限公司 | 快速计算微循环阻力的方法及*** |
CN108022650B (zh) * | 2017-12-07 | 2020-11-17 | 博动医学影像科技(上海)有限公司 | 管腔建模方法及计算血管压力差的方法与*** |
US11311196B2 (en) * | 2018-02-23 | 2022-04-26 | Boston Scientific Scimed, Inc. | Methods for assessing a vessel with sequential physiological measurements |
CN108573488B (zh) * | 2018-03-27 | 2021-06-01 | 杭州脉流科技有限公司 | 一种计算瞬时无波形比率的装置 |
JP7314183B2 (ja) * | 2018-06-15 | 2023-07-25 | パイ メディカル イメージング ビー ヴイ | 定量的血行動態フロー分析のための方法および装置 |
CN108742667B (zh) * | 2018-06-20 | 2022-07-26 | 上海博动医疗科技股份有限公司 | 基于身体质量指数获取血流特征值的方法和装置 |
CN109044324B (zh) * | 2018-06-20 | 2021-11-19 | 博动医学影像科技(上海)有限公司 | 基于斑块位置修正血流特征值的方法及装置 |
CN109009061B (zh) * | 2018-06-20 | 2021-11-19 | 博动医学影像科技(上海)有限公司 | 基于血压修正获取血流特征值的计算方法及装置 |
CN110367965B (zh) * | 2018-09-19 | 2022-03-08 | 苏州润迈德医疗科技有限公司 | 便捷测量冠状动脉血管评定参数的方法、装置及*** |
CN109512450A (zh) * | 2018-10-18 | 2019-03-26 | 深圳市孙逸仙心血管医院(深圳市心血管病研究所) | 测量血管血流速度的方法 |
CN112022137B (zh) * | 2018-11-30 | 2021-07-13 | 博动医学影像科技(上海)有限公司 | 建立血管截面函数和血管应力的方法及装置 |
WO2020167631A1 (en) * | 2019-02-11 | 2020-08-20 | University Of Louisville Research Foundation, Inc. | System and method for determining a blood flow characteristic |
CN113827199B (zh) * | 2021-10-29 | 2024-01-23 | 苏州润迈德医疗科技有限公司 | 基于造影图像调节血管评定参数的方法、***及存储介质 |
CN113995388B (zh) * | 2021-12-13 | 2022-06-03 | 数坤(北京)网络科技股份有限公司 | 血流储备分数计算方法、装置、电子设备及可读存储介质 |
CN115272447B (zh) * | 2022-09-29 | 2022-12-20 | 全景恒升(北京)科学技术有限公司 | 基于多模态影像的血流储备分数计算方法、装置及设备 |
CN116211463B (zh) * | 2023-05-05 | 2023-08-01 | 杭州脉流科技有限公司 | 一种编织支架虚拟套叠装置 |
CN117064423B (zh) * | 2023-09-14 | 2024-06-04 | 柏意慧心(杭州)网络科技有限公司 | 用于确定支架过度扩张率的方法、计算设备和介质 |
CN116898472B (zh) * | 2023-09-14 | 2023-12-12 | 柏意慧心(杭州)网络科技有限公司 | 用于确定主动脉血管直径的方法、计算设备和介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103932694A (zh) * | 2014-05-07 | 2014-07-23 | 霍云龙 | 精确诊断心肌血流储备分数(ffr)的方法和设备 |
CN104768465A (zh) * | 2012-11-06 | 2015-07-08 | 皇家飞利浦有限公司 | 血流储备分数(ffr)指标 |
CN105249954A (zh) * | 2014-06-30 | 2016-01-20 | 西门子公司 | 预测支架术后血液动力学指标的方法和*** |
CN105326486A (zh) * | 2015-12-08 | 2016-02-17 | 上海交通大学 | 血管压力差与血流储备分数的计算方法及*** |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2817460B1 (fr) * | 2000-12-04 | 2003-09-05 | Ge Med Sys Global Tech Co Llc | Procede et systeme pour simuler l'elargissement de diametre d'une lesion d'un vaisseau sanguin, en particulier une stenose, a l'aide d'une prothese endovasculaire |
EP3363350B1 (en) * | 2009-09-23 | 2019-12-11 | Lightlab Imaging, Inc. | Lumen morphology and vascular resistance measurements data collection systems, apparatus and methods |
US8315812B2 (en) * | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
US8696584B2 (en) | 2011-10-05 | 2014-04-15 | 3Dt Holdings, Llc | Devices, systems and methods for determining fractional flow reserve in the presence of a catheter |
US10162932B2 (en) * | 2011-11-10 | 2018-12-25 | Siemens Healthcare Gmbh | Method and system for multi-scale anatomical and functional modeling of coronary circulation |
CN103300820A (zh) | 2012-03-13 | 2013-09-18 | 西门子公司 | 用于冠状动脉狭窄的非侵入性功能评估的方法和*** |
US9247918B2 (en) * | 2012-07-09 | 2016-02-02 | Siemens Aktiengesellschaft | Computation of hemodynamic quantities from angiographic data |
US10433740B2 (en) * | 2012-09-12 | 2019-10-08 | Heartflow, Inc. | Systems and methods for estimating ischemia and blood flow characteristics from vessel geometry and physiology |
US9858387B2 (en) * | 2013-01-15 | 2018-01-02 | CathWorks, LTD. | Vascular flow assessment |
NL2012459B1 (en) * | 2014-03-18 | 2016-01-08 | Medis Ass B V | Method and device for determining deviation in pressure in a blood vessel. |
-
2016
- 2016-05-31 CN CN201610379167.2A patent/CN106073894B/zh active Active
- 2016-11-04 DE DE112016005442.4T patent/DE112016005442B4/de active Active
- 2016-11-04 WO PCT/CN2016/104685 patent/WO2017206438A1/zh active Application Filing
- 2016-11-04 US US15/779,704 patent/US10617473B2/en active Active
- 2016-11-04 JP JP2018548257A patent/JP6723670B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104768465A (zh) * | 2012-11-06 | 2015-07-08 | 皇家飞利浦有限公司 | 血流储备分数(ffr)指标 |
CN103932694A (zh) * | 2014-05-07 | 2014-07-23 | 霍云龙 | 精确诊断心肌血流储备分数(ffr)的方法和设备 |
CN105249954A (zh) * | 2014-06-30 | 2016-01-20 | 西门子公司 | 预测支架术后血液动力学指标的方法和*** |
CN105326486A (zh) * | 2015-12-08 | 2016-02-17 | 上海交通大学 | 血管压力差与血流储备分数的计算方法及*** |
Also Published As
Publication number | Publication date |
---|---|
JP2019500986A (ja) | 2019-01-17 |
JP6723670B2 (ja) | 2020-07-15 |
DE112016005442B4 (de) | 2021-07-15 |
DE112016005442T5 (de) | 2018-08-09 |
US20180353241A1 (en) | 2018-12-13 |
WO2017206438A1 (zh) | 2017-12-07 |
CN106073894A (zh) | 2016-11-09 |
US10617473B2 (en) | 2020-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106073894B (zh) | 基于植入虚拟支架的血管压力降数值及血流储备分数的评估方法和*** | |
JP6792768B2 (ja) | 血流予備量比の計算方法およびシステム | |
CN107978371B (zh) | 快速计算微循环阻力的方法及*** | |
CN108186038B (zh) | 基于动脉造影影像计算冠脉血流储备分数的*** | |
US9977869B2 (en) | Vascular flow assessment | |
CN109065170B (zh) | 获取血管压力差的方法及装置 | |
CN106650029B (zh) | 基于cfd仿真的分叉血管压力差及ffr的快速计算方法及*** | |
Papafaklis et al. | Anatomically correct three-dimensional coronary artery reconstruction using frequency domain optical coherence tomographic and angiographic data: head-to-head comparison with intravascular ultrasound for endothelial shear stress assessment in humans | |
CN108742587B (zh) | 基于病史信息获取血流特征值的方法及装置 | |
CN109009061B (zh) | 基于血压修正获取血流特征值的计算方法及装置 | |
CN108717874B (zh) | 基于特定的生理参数获取血管压力值的方法及装置 | |
CN111134651A (zh) | 基于腔内影像计算血流储备分数的方法、装置、***以及计算机存储介质 | |
CN109620199A (zh) | 建立血管截面函数、血管压力差和血管应力的方法及装置 | |
CN106535746A (zh) | 用于脉管处置的设备、***和方法 | |
Govindaraju et al. | Evaluation of functional severity of coronary artery disease and fluid dynamics' influence on hemodynamic parameters: A review | |
CN108742570A (zh) | 基于冠脉优势类型获取血管压力差的方法及装置 | |
CN108742667B (zh) | 基于身体质量指数获取血流特征值的方法和装置 | |
CN109044324B (zh) | 基于斑块位置修正血流特征值的方法及装置 | |
CA2331106A1 (en) | Apparatus and method for identification and characterization of lesions and therapeutic success by flow disturbances analysis | |
CN108784676B (zh) | 基于年龄信息获取压力差的方法及装置 | |
CN108742562A (zh) | 基于高血脂信息获取血管压力差的方法及装置 | |
CN108777174A (zh) | 基于心梗史信息获取血管压力差的方法及装置 | |
CN108742586A (zh) | 基于糖尿病病史信息获取血流特征值的方法及装置 | |
Tar et al. | The effect of the sensor position of the pressure wire distal to a coronary stenosis on the calculated fractional flow reserve | |
Christiaens et al. | A simple determination system for optimal angiographic viewing angles and QCA parameters |
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 | ||
CP03 | Change of name, title or address |
Address after: 200233 16th floor, building 82, No. 1198, Qinzhou North Road, Xuhui District, Shanghai Patentee after: Shanghai Bodong Medical Technology Co.,Ltd. Address before: 200240 room g3054, floor 3, building 41, No. 398, Heqing Road, Minhang District, Shanghai Patentee before: PULSE MEDICAL IMAGING TECHNOLOGY (SHANGHAI) Co.,Ltd. |
|
CP03 | Change of name, title or address |