CN108922580A - 一种获取血流储备分数的方法、装置、系统和计算机存储介质 - Google Patents
一种获取血流储备分数的方法、装置、系统和计算机存储介质 Download PDFInfo
- Publication number
- CN108922580A CN108922580A CN201810517185.1A CN201810517185A CN108922580A CN 108922580 A CN108922580 A CN 108922580A CN 201810517185 A CN201810517185 A CN 201810517185A CN 108922580 A CN108922580 A CN 108922580A
- Authority
- CN
- China
- Prior art keywords
- blood flow
- flow reserve
- reserve score
- model
- coronary artery
- 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.)
- Pending
Links
- 230000017531 blood circulation Effects 0.000 title claims abstract description 88
- 238000000034 method Methods 0.000 title claims abstract description 45
- 210000004351 coronary vessel Anatomy 0.000 claims abstract description 52
- 230000000004 hemodynamic effect Effects 0.000 claims abstract description 33
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 18
- 230000002068 genetic effect Effects 0.000 claims abstract description 15
- 210000000709 aorta Anatomy 0.000 claims description 32
- 230000036772 blood pressure Effects 0.000 claims description 27
- 230000000747 cardiac effect Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 16
- 210000001367 artery Anatomy 0.000 claims description 14
- 238000004088 simulation Methods 0.000 claims description 12
- 230000035487 diastolic blood pressure Effects 0.000 claims description 11
- 238000007634 remodeling Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000035488 systolic blood pressure Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 206010020565 Hyperaemia Diseases 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 230000004087 circulation Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 4
- 210000004204 blood vessel Anatomy 0.000 description 20
- 239000008280 blood Substances 0.000 description 17
- 238000010586 diagram Methods 0.000 description 17
- 210000004369 blood Anatomy 0.000 description 12
- 230000002792 vascular Effects 0.000 description 7
- 238000002583 angiography Methods 0.000 description 6
- 230000003902 lesion Effects 0.000 description 6
- 239000003990 capacitor Substances 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 5
- 238000009396 hybridization Methods 0.000 description 5
- 210000003462 vein Anatomy 0.000 description 5
- 238000003745 diagnosis Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 210000002837 heart atrium Anatomy 0.000 description 4
- 210000005240 left ventricle Anatomy 0.000 description 4
- 201000000057 Coronary Stenosis Diseases 0.000 description 3
- 206010011089 Coronary artery stenosis Diseases 0.000 description 3
- 230000035772 mutation Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000002601 radiography Methods 0.000 description 3
- 210000005241 right ventricle Anatomy 0.000 description 3
- 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 2
- 210000002376 aorta thoracic Anatomy 0.000 description 2
- 210000001772 blood platelet Anatomy 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000013160 medical therapy Methods 0.000 description 2
- 230000036651 mood Effects 0.000 description 2
- 230000002107 myocardial effect Effects 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 230000010412 perfusion Effects 0.000 description 2
- 230000035790 physiological processes and functions Effects 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 210000001013 sinoatrial node Anatomy 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000002560 therapeutic procedure Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- 239000002126 C01EB10 - Adenosine Substances 0.000 description 1
- 206010064571 Gene mutation Diseases 0.000 description 1
- 102000001554 Hemoglobins Human genes 0.000 description 1
- 108010054147 Hemoglobins Proteins 0.000 description 1
- 206010020751 Hypersensitivity Diseases 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 229960005305 adenosine Drugs 0.000 description 1
- 208000026935 allergic disease Diseases 0.000 description 1
- 230000007815 allergy Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 208000029078 coronary artery disease Diseases 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001595 flow curve Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000001435 haemodynamic effect Effects 0.000 description 1
- 238000005534 hematocrit Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 210000004115 mitral valve Anatomy 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000036262 stenosis Effects 0.000 description 1
- 208000037804 stenosis Diseases 0.000 description 1
- 230000002966 stenotic effect Effects 0.000 description 1
- 238000011144 upstream manufacturing Methods 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Classifications
-
- 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
- G16H10/00—ICT specially adapted for the handling or processing of patient-related medical or healthcare data
- G16H10/40—ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Primary Health Care (AREA)
- Public Health (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种获取血流储备分数的方法、装置、系统和计算机存储介质,其中获取血流储备分数的方法,包括步骤S1,获取与冠状动脉以及心脏相关的图像数据和生理参数,通过对图像数据的处理构建相应的三维模型;步骤S2,根据所述三维模型求解血液动力学控制方程,以获得三维模型所表达区域的冠状动脉的血液动力学参数分布;其中三维模型边界部分的血液动力学参数是利用集总模型且通过遗传算法求解获得;步骤S3,根据步骤S2获得的相关参数,计算血管狭窄部位的血流储备分数;步骤S4,基于步骤S3的结果生成相应的血流储备分数报告。本发明可精确计算获取患者冠状动脉FFR。
Description
技术领域
本发明涉及医学影像的数据处理,尤其涉及一种基于医学影像数据获取血流储备分数(FFR)的方法、装置、系统和计算机存储介质。
背景技术
冠状动脉血管造影曾被认为是诊断冠心病的“金标准”,但它只能定性地评价病变狭窄程度,却无法定量地评价病变狭窄对冠状动脉生理功能的影响,因此可能高估或低估病变的严重程度,导致需要治疗的病变没有处理或过度处理不需治疗的病变。1993年NicoPijls等提出了通过血压测定推算冠状动脉血流的新指标—血流储备分数(FractionalFlow Reserve,FFR)。经过长期的基础与临床研究,FFR已经成为冠状动脉狭窄功能性评价的公认指标。
血流储备分数(FFR)是一种用于医疗诊断冠状动脉生理功能的参数,指在冠状动脉存在狭窄病变的情况下,该血管所供心肌区域能获得的最大血流量与同一区域理论上正常情况下所能获得的最大血流量之比,其等效压强比定义如下:冠状动脉在最大充血状态下,狭窄近端和心脏主动脉的压强的比值即为血流储备分数。
可通过侵入式测试获取FFR,例如诊断性心导管插入术,其可包括进行常规冠状动脉血管造影(CCA)以可使冠状动脉病变可视化,同时在静脉内施用腺苷诱导的条件下(冠状动脉处于最大充血状态),计算由压强传感器获取的冠状动脉狭窄近端压强和心脏主动脉压强的比值从而得到FFR。侵入式测试的缺点在于会对患者造成增加的风险和更多的费用。
因此需要一种以非侵入式获取FFR的方法,降低诊断的风险和费用,此方法和系统需要提供患者冠状动脉特异性医学影像和患者生理信息,影像可由计算机断层成像造影(CTA)、旋转造影(RA)、核磁共振成像造影(MRA)和数字减影造影(DSA)等方法获取。另外此方法和系统还可在不能直接测量的条件下(例如运动,体质过敏不适)获取FFR,预测冠状动脉血流和心机灌注的医疗治疗、介入治疗和手术治疗的结果。
发明内容
本发明提供一种在通过医学造影术获取患者特异性冠状动脉医学影像和患者生理信息后,利用血流动力学建模获取血流储备分数(Fractional Flow Reverse,FFR)的方法。
一种获取血流储备分数的方法,包括:
步骤S1,获取与冠状动脉以及心脏相关的图像数据和生理参数,通过对图像数据的处理构建相应的三维模型;
步骤S2,根据所述三维模型求解血液动力学控制方程,以获得三维模型所表达区域的冠状动脉的血液动力学参数分布;
其中三维模型边界部分的血液动力学参数是利用集总模型且通过遗传算法求解获得;
步骤S3,根据步骤S2获得的相关参数,计算血管狭窄部位的血流储备分数;
步骤S4,基于步骤S3的结果生成相应的血流储备分数报告。
本发明采用非侵入式方法,降低诊断的风险和费用,快速便捷地获取FFR,可预测冠状动脉血流和心机灌注的医疗治疗、介入治疗和手术治疗的结果。本发明方法和系统综合了CFD和基于病人特异性边界条件仿真的优点,可以通过CFD精确计算冠状动脉正常和狭窄区域的FFR,也可以利用集总模型涵盖了非CFD仿真几何区域,如微小血管和毛细血管,为CFD模型提供精确计算的边界条件。
可选的,构建三维模型时,包括:
获取冠状动脉图像数据,对所获得的冠状动脉图像数据进行重构,得到三维形式的冠状动脉重构模型;
获取心脏图像数据,对所获得的心脏图像数据进行重构,得到三维形式的心脏重构模型;
根据所述冠状动脉重构模型建立冠状动脉的几何模型和网格模型。
可选的,所述三维模型中,边界部分至少包括冠状动脉各分支血管出口边界,血液动力学控制方程与集总模型在冠状动脉的各分支血管出口边界耦合,耦合设置条件为血液流量或血压。
可选的,所述集总模型内的目标参数通过遗传算法进行求解,包括:
依照预设的边界条件,将所述生理参数中相关的实测数据作为目标曲线;
集总模型的目标参数作为遗传算法的种群的个体,在种群逐代进化时,每代均获得对应的进化后目标参数,依据该进化后目标参数拟合获得相应的模拟曲线并与目标曲线相比;
循环生成模拟曲线,直至模拟曲线与目标曲线的误差值小于设定值,停止进化,获得目标参数的最优解,即获得冠状动脉在三维模型边界处的血液动力学参数。
可选的,所述目标曲线为主动脉出口血压的实测波形曲线,或主动脉出口血液流量的实测波形曲线。
可选的,所述误差值为yi是目标曲线,y′i是模拟波形曲线。
可选的,步骤S3中,血流储备分数计算公式为:
Pdistal为血管狭窄位置处下游2~3cm的血压;
Pa取近似值,等于主动脉运动参考平均压强
可选的,所述的计算方式为:
Psys-hyperemia=Psys-rest×(1+γsys)
Pdias-hyperemia=Pdias-rest×(1+γdias)
Pdias-hyperemia为主动脉运动舒张压;
Psys-hyperemia为主动脉运动收缩压;
Pdias-rest为静止状态舒张压;
Psys-rest为静止状态收缩压;
Sd为心舒张期占比;
St为心收缩期占比;
γsys为收缩压变化率;
γdias为舒张压变化率。
可选的,所述血流储备分数报告中,在所述三维模型中模拟显示各部位的血流储备分数。
本发明还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被计算机处理器执行时实现所述的获取血流储备分数的方法。
本发明还提供一种获取血流储备分数的装置,包括计算机存储器、计算机处理器、以及存储在所述计算机存储器中并可在所述计算机处理器上执行的计算机程序,所述计算机处理器执行所述计算机程序时,实现所述的获取血流储备分数的方法。
本发明还提供一种获取血流储备分数的系统,包括终端以及服务器,所述服务器包括计算机存储器、计算机处理器、以及存储在所述计算机存储器中并可在所述计算机处理器上执行的计算机程序,所述服务器从终端获取与冠状动脉以及心脏相关的图像数据和生理参数;
所述计算机处理器执行所述计算机程序时,实现所述的获取血流储备分数的方法。
本发明通过获取病人特异性冠状动脉三维几何结构,患者特异性生理和临床检测参数,依据这些患者特异性参数信息,医学统计大数据信息,或者心脏和冠状动脉的集总模型,确定并优化边界条件和血液特征参数,进而求解血液动力学控制方程,通过计算流体力学(Computational Fluid Dynamics,CFD)精确计算获取患者冠状动脉FFR云图分布(包括狭窄区域的FFR值)。
附图说明
图1为冠状动脉图像数据的示意图;
图2为冠状动脉重构模型的示意图;
图3为冠状动脉的几何模型示意图;
图4为冠状动脉的网格模型图示意图;
图5为图4中狭窄部位冠状动脉的的示意图;
图6为集总模式的原理示意图;
图7为图6的集总模型中心脏模型的示意图;
图8为三元Windkessel模型的示意图;
图9为二元Windkessel模型的示意图;
图10为四元Windkessel模型的示意图;
图11为双腔Windkessel模型的示意图;
图12为遗传算法流程原理示意图;
图13为主动脉出口血液流量的模拟曲线与目标曲线的变化趋势示意图,图中横坐标为时间(s),纵坐标为血液流量(ml/s);
图14为冠状动脉的某一出口血液流量的模拟曲线与目标曲线的变化趋势示意图,图中横坐标为时间(s),纵坐标为血液流量(ml/s);
图15为三维模型表达区域中,冠状动脉的血压分布示意图;
图16为FFR数值的分布示意图;
图17为图16中冠状动脉狭窄区域的局部图;
图18为本发明流程示意图。
具体实施方式
本发明一种获取血流储备分数的方法,包括:
步骤S1,获取与冠状动脉以及心脏相关的图像数据和生理参数,通过对图像数据的处理构建相应的三维模型;
步骤S2,根据三维模型求解血液动力学控制方程,以获得三维模型所表达区域的冠状动脉的血液动力学参数分布;
其中三维模型边界处的血液动力学参数是利用集总模型且通过遗传算法求解获得;
步骤S3,根据步骤S2获得的相关参数,计算血流储备分数。
步骤S4,基于步骤S3的结果生成相应的血流储备分数报告。
以下结合相应附图对各步骤详细说明。
步骤S1中,生理参数包括:身高、体重、血液密度、血液粘度、血细胞比容、血小板、肱动脉收缩压和舒张压曲线、心率、主动脉血压波形图、通过多普勒超声测得的主动脉和主要冠状动脉流量以及波形图、病史等。
步骤S1中,构建三维模型时,包括:
获取冠状动脉图像数据,对所获得的冠状动脉图像数据进行重构,得到三维形式的冠状动脉重构模型;
获取心脏图像数据,对所获得的心脏图像数据进行重构,得到三维形式的心脏重构模型;
根据冠状动脉重构模型建立冠状动脉的几何模型和网格模型。
步骤S1中,与冠状动脉以及心脏相关的图像数据包括但不限于计算机断层成像(CT)、动态(CT)、核磁共振、脉管造影术、超声、单光子发射计算断层成像(SPECT)和任意其他类型的医学成像模态。图像数据可以是2D、3D、4D,可以从一个或多个图像采集设备,诸如CT扫描仪、MR扫描仪、脉管造影术扫描仪、超声设备等接收,或者患者有记录的医学图像数据。还可根据需要采集包含其他部位的影像数据,例如脊椎、肺、前胸等。
例如图1中冠状动脉图像数据为CTA方式获取的dicom格式。
本发明中,冠状动脉图像数据以及相应的三维模型中除冠状动脉外还包括部分主动脉,部分主动脉至少由左心室出口延伸至冠状动脉入口。
部分主动脉一般包括升主动脉和降主动脉,或至少包括由左心室至冠状动脉入口部位的主动脉。
冠状动脉图像数据重构时,需要计算冠状动脉相关的结构和形态参数,如:
升主动脉,降主动脉;
左冠状动脉(LCA):左主干,左前降支(左室前支,右室前支,前间隔支),左回旋支(窦房结支,左房支,左室前支,钝缘支,左室后支);
右冠状动脉(RCA):右主干,右圆锥支,窦房结支,右室前支,锐缘支,后降支,左室后支,右心房支。
因CTA空间分辨率的局限性,图1和图2示意的部分只是重构出主动脉和直径大于0.5mm的冠状动脉血管。冠状动脉的微小动脉血管和毛细血管不在重构模型内。也正基于此,本发明的后续的步骤用利用集总模型对冠状动脉的微小动脉血管和毛细血管进行模拟,获取重构模型中处在血管边界部位的血液动力学参数,例如血液流量或血压,以提高FFR计算的精确性。
根据冠状动脉图像数据分辨率,三维模型优选包括主动脉和直径大于0.2mm的冠状动脉。
心脏图像数据重构时,至少包括与主动脉相连的左心室部分,还可以根据需要重构右心室,左心房或右心房等。
冠状动脉、心脏图像数据重构可采用现有算法,例如重构算法选自Level Sets,Graph Cuts,CNN神经网络模型,深度学习模型或边缘学习模型等。重构模型例如为stl格式的几何模型,可参见图2中冠状动脉重构模型的示意图。
步骤S1中,根据冠状动脉重构模型建立冠状动脉的几何模型时,需要标注所有血管的血管壁,进口和出口。
本实施例中冠状动脉的网格模型包含一个主动脉进口,一个主动脉出口,和多个分支血管出口,例如图3中所标识的主动脉进口A、主动脉出口B、冠状动脉的分支血管出口C、主动脉血管壁D、冠状动脉管壁E。
因为冠状动脉的血管系统具有非常复杂的三维拓扑结构,依据血管几何尺寸和弯曲度,网格在血管各个细分区域做自适应划分,最大尺寸为3mm,最小尺寸为0.1mm,效果示意可参见图4。
步骤S2中,血液动力学控制方程为:
与血液动力学控制方程相应的边界控制方程为:
AP=RQ
Q=Qo(D/Do)a
R=Po/(Qo(D/Do)a)
其中:
△P为一段血管的压力降,(例如上游血管或下游血管);
Q为边界处血流量;
R为血流阻力;
D为血管直径;
Do为某参照处血管的血管直径;
Po为与某参照处血管相对应的血压值;
Qo为与某参照处血管相对应的流量;
ρ为血液密度;
μ为血液粘度;
F为体积力,如重力;
t为时间;
V为血流速度;
P为血压。
边界控制方程根据Poiseuille Law(血管阻力计算方程)推导。
Poiseuille Law:
R为血流阻力;
r为血管半径;
L为血管长度;
μ为血液粘度;
Q为血液流量;
P为血压;
hem为血红细胞比容;
c为血管弹性。
步骤S2中,血液动力学参数分布为冠状动脉的血流速度或血压分布。基于冠状动脉的网格模型以及预设的边界条件,求解血液动力学控制方程可以获得冠状动脉的血流速度或血压分布。
由于冠状动脉的三维模型仅表达了局部区域,对于未重构出的冠状动脉的微小动脉血管和毛细血管,需用通过集总模型(结构可参照图6)做解析计算。
求解血液动力学控制方程可使用CFD计算方法,具体可包括但不限于有限元法(FEM),有限体积法(FVM),有限差分法(FDM),边界元法(BEM),浸入边界法(IBM),格子玻尔兹曼法(LBM),光滑粒子法(SPH),半隐式运动粒子法(MPS),有限体积粒子法(FVP)等。
求解血液动力学控制方程以及集中模型所需的血流动力学参数可预先通过生理参数、三维模型、历史数据等方式获得。
血流动力学参数可以包括患者身高,体重,心率,收缩期和舒张期血压,血细胞比容,血液粘度,血液密度,血流量等。
血流动力学参数还可以包括生理医学测试结果(心动周期,血压,血流,血红蛋白,血小板,心电图,基因,家族病史等等),图像数据/分割数据/重构几何数据(心脏尺寸,冠状动脉分支和拓扑结构,狭窄位置,狭窄长度,狭窄截面,钙化斑块,等等)。
冠状动脉的三维模型中,边界部分包括血管壁边界,血管进口边界(主动脉进口边界),血管出口边界(主动脉出口边界,冠状动脉各分支血管出口边界);
针对血管进口边界,血液流量或血压可来自步骤S1中针对患者测得的生理参数。
主动脉运动血流量QCOhyperemia为主动脉进口边界条件的一个选项,对应边界设置。QCOhyperemia为主动脉运动流量,计算方式为:
QCOhyperemia=QCOrest×χ
QCOrest=Qcoronary/β
mmyo=ρmyo×VLV
其中α,β,χ,为参数系数;通过医学统计大数据,依据患者个体特征作相应调整,其中为心肌单位质量血流消耗率;
QCOrest为主动脉静息流量;即静息状态下主动脉血流量;
Qcoronary为冠状动脉总流量;可通过医学经验公式得出;
mmyo为心肌质量;
ρmyo为心肌密度;一般可以根据经验值确定;
VLV为左心室体积;可从心脏重构模型中获得。
冠状动脉各分支血管出口边界的血液动力学参数是利用集总模型且通过遗传算法求解获得;
血液动力学控制方程与集总模型在冠状动脉的各分支血管出口耦合,耦合设置条件为血液流量或血压。
例如采用血压作为边界条件,通过求解血液动力学控制方程得到血管出口边界的血液流量;再将血液流量作为输入,通过集总模型求解得到血管出口边界的血压;再将血压作为冠状动脉在三维模型边界处的血液动力学参数。
就集中模型本身而言,可采用现有技术构建,本发明重点之一在于利用遗传算法求解集中模型。在集总模型中,电阻表征血流在流动过程中受到的阻力影响,电容表征血管或心肌的体积变化对血液流动的影响,电感表征血液流动时的惯性作用。
集总模型包括但不限于二元Windkessel模型,三元Windkessel模型,四元Windkessel模型,双腔Windkessel模型,心脏集总模型,冠状动脉-心脏耦合集总模型,以及其他高阶集总模型。
一般的对应关系为电流-血流量,电压-血压,电阻-血流阻力,电感-血流惯性,电容-血流顺应性。
参见图7,心脏部分的模型包括左心房压强P(t),二尖瓣(与P(t)部分连接的二极管),心房-心室瓣膜电阻Rav,心房-心室瓣膜电感Lav,主动脉瓣(与Lav连接的二极管),心室-动脉瓣膜电阻RV-ar1,心室-动脉瓣膜电感Lv-ar和时变倒电容Ev(t)。
时变倒电容值可以由以下公式表示,
Ev(t)=(Emax-Emin)·En(t)+Emin
Emax是心室在收缩末期压强-体积比,Emin是心室在舒张末期的压强-体积比,En(t)可以表示为双希尔方程:
其中tn=t/Tmax,Tmax通过心动周期tc计算:
Tmax=0.2+0.15tc
参见图8,以集总模型中三元Windkessel模型为例(其他同理),使用RCR电路模拟微小动脉和毛细血管网的血压,血液流量和血管阻力,具体为:
电阻R对应血管阻力,电容C对应血管弹性,电压对应血压P,电流对应血管流量Q,其电路方程为三元Windkessel方程:
集总模型中,心脏模型内有4个变量即目标参数(两个电阻+两个电感),主动脉三元Windkessel模型有3个变量(一个电容+两个电阻),每一支冠状动脉出口双腔Windkessel模型有5个变量(两个电容+三个电阻)。电阻的取值范围设为(0,500),电容的取值范围取为(0,1)。
遗传算法种群个体由这些变量组成的二进制串,初始种群由若干个个体组成,且由算法随机产生。种群通过一代代自然选择,杂交及基因变异向最优解进化。遗传算法的原理可参见图12。
遗传算法借鉴了进化生物学中的自然选择,杂交,基因突变等现象。在集总模型的参数优化问题中,未知量的解都由二进制表示(即0和1的字符组成),随机产生一定个体数目的种群,进化从这些随机个体开始。在每一代进化中,通过评价该种群的适应度,基于适应度从该种群中随机选择若干个个体,通过自然选择,杂交及变异产生新的种群,成为下一次迭代的初始种群,促使种群向最优解进化。自然选择从种群中选取若干较优个体并将基因传递到下一代,自然选择的方法包括但不限于轮盘赌选择策略,竞标赛选择策略等。杂交操作会根据交叉概率将两个个体中的部分基因按某种方式交换得出新的个体,杂交的策略包括但不限于单点交换,两点交换,三点交换等。变异操作会根据变异概率随机改变某点基因,变异的策略包括但不限于内变异,外变异,单点变异和多点变异等。
集总模型内的目标参数通过遗传算法进行求解,包括:
依照预设的边界条件(例如为血压或血液流量),将生理参数中相关的实测数据作为目标曲线;
集总模型的目标参数作为遗传算法的种群的个体,在种群逐代进化时,每代均获得对应的进化后目标参数,依据该进化后目标参数拟合获得相应的模拟曲线并与目标曲线相比;
循环生成模拟曲线,直至模拟曲线与目标曲线的误差值小于设定值,停止进化,获得目标参数的最优解,即获得冠状动脉在三维模型边界处的血液动力学参数。
目标曲线为主动脉出口血压的实测波形曲线,或主动脉出口血液流量的实测波形曲线。
为了进一步提高精度,还可以比较多条目标曲线,例如目标曲线还包括冠状动脉某一出口血压的实测波形曲线,或冠状动脉某一出口血液流量的实测波形曲线。
误差值为yi是目标曲线,y'i是模拟波形曲线,整个遗传算法就是在求解该误差函数的最小值。当误差值小于0.5%,计算所得的集总模型中的各目标参数值即可认定为最优解。模拟曲线与目标曲线的变化趋势可参见图13和图14。
通过集中模型与血液动力学控制方程的耦合,可以获得三维模型表达区域中,冠状动脉的血压分布,参见图15。
步骤S3中,血流储备分数(FFR)计算公式为:
Pdistal为血管狭窄位置处下游2~3cm的血压;
血管狭窄位置一方面可以通过步骤S1的图像数据和生理参数获知,还可以通过步骤S1构建的三维模型(例如重构模型)测量获得;
Pa取近似值,等于主动脉运动参考平均压强
计算方式为:
Psys-hyperemia=Psys-rest×(1+γsys)
Pdias-hyperemia=Pdias-rest×(1+γdias)
Pdias-hyperemia为主动脉运动舒张压;
Psys-hyperemia为主动脉运动收缩压;
Pdias-rest为静止状态舒张压;
Psys-rest为静止状态收缩压;
Sd为心舒张期占比;
St为心收缩期占比;
γsys为收缩压变化率;
γdias为舒张压变化率;
Sd、St、γsys、γdias均为医学统计大数据,可依据患者个体特征作相应调整。
步骤S4中基于步骤S3的结果生成相应的血流储备分数报告,血流储备分数报告的形式可以是文本、图片或三维影像等。
作为血流储备分数报告的形式之一,作为优选,血流储备分数报告中,在三维模型中模拟显示各部位的血流储备分数。模拟显示效果参见图16和图17。
本实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被计算机处理器执行时实现上述获取血流储备分数的方法。
本实施例还提供一种获取血流储备分数的装置,包括计算机存储器、计算机处理器、以及存储在所述计算机存储器中并可在所述计算机处理器上执行的计算机程序,所述计算机处理器执行所述计算机程序时,实现上述获取血流储备分数的方法。
本实施例装置可配置的远端,通过与之相连的远程终端获得与冠状动脉以及心脏相关的图像数据和生理参数,还可以是本实施例装置本身就配置在终端,直接通过医学影像设备或人工输入的方式获得与冠状动脉以及心脏相关的图像数据和生理参数。
本实施例还提供一种获取血流储备分数的系统,包括终端以及服务器,所述服务器包括计算机存储器、计算机处理器、以及存储在所述计算机存储器中并可在所述计算机处理器上执行的计算机程序,所述服务器从终端获取与冠状动脉以及心脏相关的图像数据和生理参数;
所述计算机处理器执行所述计算机程序时,实现上述获取血流储备分数的方法。
本实施例系统中的终端以及云端的服务器上安装有分别安装有软件系统,服务器系统通过C,C++,Java,Python,HTM5等编程语言开发。终端软件系统包括但不限于Mac OS版,Windows版,Unix/Linux版,Androids版或Apple版。
以上公开的仅为本发明的具体实施例,但是本发明并非局限于此,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。显然这些改动和变型均应属于本发明要求的保护范围保护内。此外,尽管本说明书中使用了一些特定的术语,但这些术语只是为了方便说明,并不对本发明构成任何特殊限制。
Claims (12)
1.一种获取血流储备分数的方法,其特征在于,包括:
步骤S1,获取与冠状动脉以及心脏相关的图像数据和生理参数,通过对图像数据的处理构建相应的三维模型;
步骤S2,根据所述三维模型求解血液动力学控制方程,以获得三维模型所表达区域的冠状动脉的血液动力学参数分布;
其中三维模型边界部分的血液动力学参数是利用集总模型且通过遗传算法求解获得;
步骤S3,根据步骤S2获得的相关参数,计算血管狭窄部位的血流储备分数;
步骤S4,基于步骤S3的结果生成相应的血流储备分数报告。
2.如权利要求1所述的获取血流储备分数的方法,其特征在于,步骤S1中,构建三维模型时,包括:
获取冠状动脉图像数据,对所获得的冠状动脉图像数据进行重构,得到三维形式的冠状动脉重构模型;
获取心脏图像数据,对所获得的心脏图像数据进行重构,得到三维形式的心脏重构模型;
根据所述冠状动脉重构模型建立冠状动脉的几何模型和网格模型。
3.如权利要求1所述的获取血流储备分数的方法,其特征在于,所述三维模型中,边界部分至少包括冠状动脉各分支血管出口边界,血液动力学控制方程与集总模型在冠状动脉的各分支血管出口边界耦合,耦合设置条件为血液流量或血压。
4.如权利要求3所述的获取血流储备分数的方法,其特征在于,所述集总模型内的目标参数通过遗传算法进行求解,包括:
依照预设的边界条件,将所述生理参数中相关的实测数据作为目标曲线;
集总模型的目标参数作为遗传算法的种群的个体,在种群逐代进化时,每代均获得对应的进化后目标参数,依据该进化后目标参数拟合获得相应的模拟曲线并与目标曲线相比;
循环生成模拟曲线,直至模拟曲线与目标曲线的误差值小于设定值,停止进化,获得目标参数的最优解,即获得冠状动脉在三维模型边界处的血液动力学参数。
5.如权利要求4所述的获取血流储备分数的方法,其特征在于,所述目标曲线为主动脉出口血压的实测波形曲线,或主动脉出口血液流量的实测波形曲线。
6.如权利要求4所述的获取血流储备分数的方法,其特征在于,所述误差值为yi是目标曲线,y′i是模拟波形曲线。
7.如权利要求1所述的获取血流储备分数的方法,其特征在于,步骤S3中,血流储备分数计算公式为:
Pdistal为血管狭窄位置处下游2~3cm的血压;
Pa取近似值,等于主动脉运动参考平均压强
8.如权利要求7所述的获取血流储备分数的方法,其特征在于,所述的计算方式为:
Psys-hyperemia=Psys-rest×(1+γsys)
Pdias-hyperemia=Pdias-rest×(1+γdias)
Pdias-hyperemia为主动脉运动舒张压;
Psys-hyperemia为主动脉运动收缩压;
Pdias-rest为静止状态舒张压;
Psys-rest为静止状态收缩压;
Sd为心舒张期占比;
St为心收缩期占比;
γsys为收缩压变化率;
γdias为舒张压变化率。
9.如权利要求1所述的获取血流储备分数的方法,其特征在于,所述血流储备分数报告中,在所述三维模型中模拟显示各部位的血流储备分数。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被计算机处理器执行时实现如权利要求1~9任一项所述的获取血流储备分数的方法。
11.一种获取血流储备分数的装置,包括计算机存储器、计算机处理器、以及存储在所述计算机存储器中并可在所述计算机处理器上执行的计算机程序,其特征在于,所述计算机处理器执行所述计算机程序时,实现如权利要求1~9任一项所述的获取血流储备分数的方法。
12.一种获取血流储备分数的系统,包括终端以及服务器,所述服务器包括计算机存储器、计算机处理器、以及存储在所述计算机存储器中并可在所述计算机处理器上执行的计算机程序,其特征在于,所述服务器从终端获取与冠状动脉以及心脏相关的图像数据和生理参数;所述计算机处理器执行所述计算机程序时,实现如权利要求1~9任一项所述的获取血流储备分数的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810517185.1A CN108922580A (zh) | 2018-05-25 | 2018-05-25 | 一种获取血流储备分数的方法、装置、系统和计算机存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810517185.1A CN108922580A (zh) | 2018-05-25 | 2018-05-25 | 一种获取血流储备分数的方法、装置、系统和计算机存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108922580A true CN108922580A (zh) | 2018-11-30 |
Family
ID=64410640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810517185.1A Pending CN108922580A (zh) | 2018-05-25 | 2018-05-25 | 一种获取血流储备分数的方法、装置、系统和计算机存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108922580A (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109567776A (zh) * | 2018-12-31 | 2019-04-05 | 深圳北芯生命科技有限公司 | 用于测试ffr主机系统的导管模拟器 |
CN109602401A (zh) * | 2019-01-28 | 2019-04-12 | 徐州市心血管病研究所 | 一种微血管血流动力学参数分析仪及分析方法 |
CN109686450A (zh) * | 2018-12-22 | 2019-04-26 | 北京工业大学 | 一种基于超声和ct成像技术的冠脉血流储备分数计算方法 |
CN110501278A (zh) * | 2019-07-10 | 2019-11-26 | 同济大学 | 一种基于YOLOv3和密度估计的细胞计数方法 |
CN110742688A (zh) * | 2019-10-30 | 2020-02-04 | 北京理工大学 | 血管模型建立方法、装置及可读取存储介质 |
CN111067494A (zh) * | 2019-12-27 | 2020-04-28 | 西北工业大学 | 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法 |
CN111241759A (zh) * | 2020-01-13 | 2020-06-05 | 北京工业大学 | 一种基于零维血流动力学模型的ffr快速计算系统模型 |
CN111312375A (zh) * | 2020-03-12 | 2020-06-19 | 上海杏脉信息科技有限公司 | 虚拟冠脉手术的搭桥分析方法、系统、介质及设备 |
CN111652881A (zh) * | 2020-07-01 | 2020-09-11 | 杭州脉流科技有限公司 | 基于深度学习的冠脉重构和血流储备分数计算方法、装置、设备以及可读存储介质 |
WO2020211051A1 (zh) * | 2019-04-18 | 2020-10-22 | 麦层移动健康管理有限公司 | 一种测量心肌组织运动特征的非侵入性方法及系统 |
CN111932652A (zh) * | 2020-10-12 | 2020-11-13 | 杭州阿特瑞科技有限公司 | 一种基于血流量参数的冠状动脉建模动态调整方法 |
CN112071419A (zh) * | 2019-06-11 | 2020-12-11 | 西门子医疗有限公司 | 使用递归神经网络对血管的血液动力学分析 |
CN112690814A (zh) * | 2020-11-06 | 2021-04-23 | 杭州阿特瑞科技有限公司 | 一种低误差的冠状动脉血流储备分数测量方法 |
CN112711881A (zh) * | 2020-12-29 | 2021-04-27 | 深圳技术大学 | 生理动态血流模拟方法、装置、计算机设备及存储介质 |
CN113545846A (zh) * | 2021-07-22 | 2021-10-26 | 强联智创(北京)科技有限公司 | 一种血流动力学仿真模拟方法、装置以及设备 |
CN113940650A (zh) * | 2021-10-19 | 2022-01-18 | 福州市第二医院(福建省福州中西医结合医院、福州市职业病医院) | 一种检测按摩推拿对血流动力学影响的方法和存储设备 |
CN114287954A (zh) * | 2021-12-13 | 2022-04-08 | 杭州脉流科技有限公司 | 颅内血管狭窄的功能学评估方法、系统、电子设备及介质 |
CN116564525A (zh) * | 2023-07-10 | 2023-08-08 | 天津恒宇医疗科技有限公司 | 基于冠脉血流量分配的血流储备分数预测方法和系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103300820A (zh) * | 2012-03-13 | 2013-09-18 | 西门子公司 | 用于冠状动脉狭窄的非侵入性功能评估的方法和系统 |
CN105877767A (zh) * | 2016-03-31 | 2016-08-24 | 北京思创贯宇科技开发有限公司 | 一种冠脉图像造影方法和装置 |
US20170132788A1 (en) * | 2015-11-10 | 2017-05-11 | General Electric Company | System and method for estimating arterial pulse wave velocity |
CN108511075A (zh) * | 2018-03-29 | 2018-09-07 | 向建平 | 一种非侵入式获取血流储备分数的方法和系统 |
-
2018
- 2018-05-25 CN CN201810517185.1A patent/CN108922580A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103300820A (zh) * | 2012-03-13 | 2013-09-18 | 西门子公司 | 用于冠状动脉狭窄的非侵入性功能评估的方法和系统 |
US20170132788A1 (en) * | 2015-11-10 | 2017-05-11 | General Electric Company | System and method for estimating arterial pulse wave velocity |
CN108471970A (zh) * | 2015-11-10 | 2018-08-31 | 通用电气公司 | 用于估计动脉脉搏波速度的系统和方法 |
CN105877767A (zh) * | 2016-03-31 | 2016-08-24 | 北京思创贯宇科技开发有限公司 | 一种冠脉图像造影方法和装置 |
CN108511075A (zh) * | 2018-03-29 | 2018-09-07 | 向建平 | 一种非侵入式获取血流储备分数的方法和系统 |
Non-Patent Citations (3)
Title |
---|
吴再彬著: "《血液动力学》", 31 January 1981 * |
闫剀: ""冠状动脉集中参数模型个性化研究"", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
陈晓阳: ""利用冠状动脉CT血管造影计算血流储备分数"", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109686450A (zh) * | 2018-12-22 | 2019-04-26 | 北京工业大学 | 一种基于超声和ct成像技术的冠脉血流储备分数计算方法 |
CN109567776A (zh) * | 2018-12-31 | 2019-04-05 | 深圳北芯生命科技有限公司 | 用于测试ffr主机系统的导管模拟器 |
CN109602401A (zh) * | 2019-01-28 | 2019-04-12 | 徐州市心血管病研究所 | 一种微血管血流动力学参数分析仪及分析方法 |
CN109602401B (zh) * | 2019-01-28 | 2021-08-27 | 徐州市心血管病研究所 | 一种微血管血流动力学参数分析仪及分析方法 |
WO2020211051A1 (zh) * | 2019-04-18 | 2020-10-22 | 麦层移动健康管理有限公司 | 一种测量心肌组织运动特征的非侵入性方法及系统 |
CN112071419A (zh) * | 2019-06-11 | 2020-12-11 | 西门子医疗有限公司 | 使用递归神经网络对血管的血液动力学分析 |
CN110501278A (zh) * | 2019-07-10 | 2019-11-26 | 同济大学 | 一种基于YOLOv3和密度估计的细胞计数方法 |
CN110501278B (zh) * | 2019-07-10 | 2021-04-30 | 同济大学 | 一种基于YOLOv3和密度估计的细胞计数方法 |
CN110742688A (zh) * | 2019-10-30 | 2020-02-04 | 北京理工大学 | 血管模型建立方法、装置及可读取存储介质 |
CN110742688B (zh) * | 2019-10-30 | 2021-05-25 | 北京理工大学 | 血管模型建立方法、装置及可读取存储介质 |
CN111067494B (zh) * | 2019-12-27 | 2022-04-26 | 西北工业大学 | 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法 |
CN111067494A (zh) * | 2019-12-27 | 2020-04-28 | 西北工业大学 | 基于血流储备分数和血流阻力模型的微循环阻力快速计算方法 |
CN111241759B (zh) * | 2020-01-13 | 2024-03-26 | 北京工业大学 | 一种基于零维血流动力学模型的ffr快速计算方法 |
CN111241759A (zh) * | 2020-01-13 | 2020-06-05 | 北京工业大学 | 一种基于零维血流动力学模型的ffr快速计算系统模型 |
CN111312375A (zh) * | 2020-03-12 | 2020-06-19 | 上海杏脉信息科技有限公司 | 虚拟冠脉手术的搭桥分析方法、系统、介质及设备 |
CN111652881A (zh) * | 2020-07-01 | 2020-09-11 | 杭州脉流科技有限公司 | 基于深度学习的冠脉重构和血流储备分数计算方法、装置、设备以及可读存储介质 |
CN111932652B (zh) * | 2020-10-12 | 2021-01-08 | 杭州阿特瑞科技有限公司 | 一种基于血流量参数的冠状动脉建模动态调整方法 |
CN111932652A (zh) * | 2020-10-12 | 2020-11-13 | 杭州阿特瑞科技有限公司 | 一种基于血流量参数的冠状动脉建模动态调整方法 |
CN112690814A (zh) * | 2020-11-06 | 2021-04-23 | 杭州阿特瑞科技有限公司 | 一种低误差的冠状动脉血流储备分数测量方法 |
CN112711881A (zh) * | 2020-12-29 | 2021-04-27 | 深圳技术大学 | 生理动态血流模拟方法、装置、计算机设备及存储介质 |
CN112711881B (zh) * | 2020-12-29 | 2023-03-31 | 深圳技术大学 | 生理动态血流模拟方法、装置、计算机设备及存储介质 |
CN113545846A (zh) * | 2021-07-22 | 2021-10-26 | 强联智创(北京)科技有限公司 | 一种血流动力学仿真模拟方法、装置以及设备 |
CN113940650A (zh) * | 2021-10-19 | 2022-01-18 | 福州市第二医院(福建省福州中西医结合医院、福州市职业病医院) | 一种检测按摩推拿对血流动力学影响的方法和存储设备 |
CN113940650B (zh) * | 2021-10-19 | 2023-04-11 | 福州市第二医院(福建省福州中西医结合医院、福州市职业病医院) | 一种检测按摩推拿对血流动力学影响的方法和存储设备 |
CN114287954A (zh) * | 2021-12-13 | 2022-04-08 | 杭州脉流科技有限公司 | 颅内血管狭窄的功能学评估方法、系统、电子设备及介质 |
CN116564525A (zh) * | 2023-07-10 | 2023-08-08 | 天津恒宇医疗科技有限公司 | 基于冠脉血流量分配的血流储备分数预测方法和系统 |
CN116564525B (zh) * | 2023-07-10 | 2023-09-26 | 天津恒宇医疗科技有限公司 | 基于冠脉血流量分配的血流储备分数预测方法和系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108922580A (zh) | 一种获取血流储备分数的方法、装置、系统和计算机存储介质 | |
JP6409114B2 (ja) | 患者固有の幾何学的形状モデルを変更することによって治療を決定する方法及びシステム | |
CN105380598B (zh) | 用于针对动脉狭窄的自动治疗规划的方法和系统 | |
KR101818645B1 (ko) | 혈류 특성 모델링 감도 분석 방법 및 시스템 | |
CN108511075B (zh) | 一种非侵入式获取血流储备分数的方法和系统 | |
CN105249954B (zh) | 预测支架术后血液动力学指标的方法和系统 | |
KR102404554B1 (ko) | 환자-특정적 해부학적 데이터로부터 개인화된 맥관 이식물을 식별하기 위한 시스템 및 방법 | |
US10354744B2 (en) | Non-invasive functional assessment of coronary artery stenosis including simulation of hyperemia by changing resting microvascular resistance | |
KR102103126B1 (ko) | 환자별 혈류 모델링 방법 및 시스템 | |
US20240081910A1 (en) | Patient-specific modeling of hemodynamic parameters in coronary arteries | |
JP2018509971A (ja) | 定量的フロー分析のための方法および装置 | |
CN108992059A (zh) | 用于根据血管几何形状和生理学估计血流特性的系统和方法 | |
KR102620470B1 (ko) | 관상 동맥의 혈류역학적 파라미터의 환자별 모델링 | |
CN106537392A (zh) | 用于冠状动脉中的血液动力学计算的方法和系统 | |
JP6987282B1 (ja) | 対象者固有の血管に対する血流シミュレーション方法及び装置 | |
CN114052764B (zh) | 获取血流储备分数的方法、装置、系统和计算机存储介质 | |
CN114947910A (zh) | 冠脉末端微血管阻力的计算方法及ffr的计算方法、系统 | |
US20210290076A1 (en) | System and a method for determining a significance of a stenosis | |
WO2020083390A1 (zh) | 获取心表大动脉的血流量的方法、装置、系统及存储介质 | |
EA042942B1 (ru) | Специфическое для пациента моделирование гемодинамических параметров в коронарных артериях | |
CN116313101A (zh) | 一种冠状动脉血流储备分数确定方法、系统、设备及介质 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181130 |
|
RJ01 | Rejection of invention patent application after publication |