基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法技术

本发明专利技术公开了一种基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，包括：测量在舒张期心脏冠脉口的压力P

A method to calculate the ratio of instantaneous wave free and diastolic pressure based on the contrast image

【技术实现步骤摘要】
基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法
本专利技术涉及冠状动脉影像学评价领域，具体地涉及一种仅仅通过造影图像和主动脉压力确定造影瞬时无波型比率(caiFR)、造影舒张期压力比率(cadFR和caDFR)的方法。
技术介绍
血流储备分数(FFR)可表明冠脉狭窄病变对远端血流产生的影响，诊断心肌是否缺血，已经成为冠脉狭窄功能性评价的公认指标。FFR定义为狭窄冠状动脉提供给支配区域心肌的最大血流量与同一支冠状动脉正常时提供给心肌的最大血流量的比值。可以简化为心肌最大充血状态下的狭窄远端冠状动脉内平均压(Pd)与冠状动脉口部主动脉平均压(Pa)的比值，即FFR＝Pd/Pa。确定FFR时候，需要基于心肌最大充血状态下的血流速度和冠状动脉口部主动脉平均压，通过不同的手段获取狭窄远端冠状动脉内平均压来计算FFR。但心肌最大充血需要通过冠脉内或静脉内注射腺苷或ATP，注射腺苷或ATP会造成主动脉压力下降并有一定的副作用比如房室传导阻滞、窦缓、窦停等，禁忌症包括2度或3度房室传到阻滞、窦房结疾病、气管或支气管哮喘、对腺苷过敏。瞬时无波形比率(iFR)能提供和血流储备分数(FFR)类似的冠状动脉内压力测量方法。iFR不需要血管扩张剂、操作简单，将会更多的应用在冠状动脉介入治疗。ADVISE研究发现，当心脏舒张期的某段时间(称之为无波形期)，冠脉内微血管阻力相对是最稳定且是最低的，和腺苷等血管扩张药物所做成的冠脉充血期间达到的平均阻力相类似。即iFR＝PdWave-freeperiod/PaWave-freeperiod(PdWave-freeperiod：在无波形期间狭窄病变远端冠脉平均压。PaWave-freeperiod：在无波形期间主动脉平均压。瞬时无波形时期的运算时间：舒张期内无波形时期开始后25％的时间，到收缩期开始前5ms的时间停止计算)。顶级医学期刊NEJM上发表了一篇研究文章，在稳定型心绞痛或急性冠脉综合征患者中，IFR指导血运重建策略不劣于FFR指导的重建策略，在12个月内主要不良心脏事件发生率方面相似。静息态舒张期压力比率(dFR和DFR)可以表示为：dFR＝PdDiastolicperiod/PaDiastolicperiod(PdDiastolicperiod：在舒张态期间狭窄病变远端冠脉平均压。PaDiastolicperiod：在舒张态期间主动脉平均压)；DFR＝PdDiastolichyperemiafreeperiod/PaDiastolichyperemiafreeperiod(PdDiastolichyperemiafreeperiod：在主动脉压力小于主动脉平均压到主动脉压力最小的区间里，狭窄病变远端冠脉平均压。PaDiastolichyperemiafreeperiod：在主动脉压力小于主动脉平均压到主动脉压力最小的区间里，主动脉平均压)。进一步的研究表明静息态舒张期压力比率(dFR和DFR)跟瞬时无波形比率(iFR)基本上完全等价。因此，我们可以得到iFR≌DFR≌dFR＝PdDiastolicperiod/PaDiastolicperiod。目前，现有的瞬时无波型比率(iFR)和静息态舒张期压力比率(dFR和DFR)的测量方法主要为：压力导丝静息态下测量相应的舒张态区间来确定iFR、dFR、DFR。需要依靠压力导丝进行测量，压力导丝测量时需要介入血管末端，增加手术难度和风险，同时压力导丝昂贵的价格也限制其大规模应用。相对于静息态，在造影状态下可以导致一定程度的心肌充血，因此造影瞬时无波型比率(caiFR)、造影舒张期压力比率(cadFR和caDFR)应该更能展现心肌缺血状态。
技术实现思路
为了解决上述的技术问题，本专利技术目的是：提供了一种基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，针对冠心病患者通过常规冠状动脉造影手术检测心肌缺血情况，即无需使用血管扩张剂(即无需心肌最大充血状态且不使用腺苷或ATP)。通过常规造影图像、主动脉压和血流，来计算在造影状态下的瞬时无波型比率(caiFR)、舒张期压力比率(cadFR和caDFR)。本专利技术的技术方案是：一种基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，包括以下步骤：S01：通过血压传感器测量在舒张期心脏冠脉口的压力Pa；S02：通过造影图像获取血管的二维管径及长度，并通过两个呈30°以上夹角的造影图像生成三维血管网格模型并获得血管的三维管径及长度；S03：在心脏舒张期，测量包含造影剂的血液从一段指定血管的起始点到结束点所用的时间，并根据该时间和血管三维长度计算血流速度V1；S04：将计算得到的造影状态下的血流速度V1作为冠脉入口流速，计算冠脉入口到冠脉狭窄远端的压力降ΔP，在舒张期狭窄远端冠状动脉内平均压Pd＝Pa-ΔP，通过公式caiFR≌caDFR≌cadFR＝Pd/Pa计算得到造影瞬时无波型比率(caiFR)和造影静息态舒张期压力比率(cadFR和caDFR)。优选的技术方案中，所述步骤S01包括，使用血压传感器的压力管连接到多联三通，然后通过造影导管与心脏冠脉口部相连，在血压传感器的压力管内充满盐水，并保持血压传感器与心脏在同一水平位置，该血压传感器测量的压力波即为心脏冠脉口的压力波，在舒张期，瞬时压力的平均值即为Pa。优选的技术方案中，所述步骤S02中生成三维血管网格模型的方法包括以下步骤：S21：将两个呈30°以上夹角的X射线冠脉造影图像上，具有映射关系的两分段血管的2D结构数据进行三维重建，得到该分段血管的3D结构数据；S22：重复步骤S21，直到所有分段血管三维重建完成，再将重建后的分段血管合并，得到完整的三维血管网格模型。优选的技术方案中，步骤S03中计算血流速度V1的具体方法包括以下步骤：S31：获取指定的病人心率为H次/分钟，从造影图像信息中获取图像频率为S帧/秒，其帧数X的计算公式如下：X＝(1÷(H÷60))×S；S32：通过一个心跳周期舒张期内图像所走过的帧数，在二维起始帧和结束帧对应的图像上分别取得一个心跳周期舒张期的起始点和结束点，然后通过起始点和结束点在三维血管网格模型中截取一个心跳周期舒张期的血管长度；S33：通过公式V1＝L÷P，计算得到血液流动速度V1，L为血管长度，P为一个心跳周期舒张期所用的时间，P＝X÷S。优选的技术方案中，所述步骤S04中计算冠脉入口到冠脉狭窄远端的压力降ΔP的具体方法如下：S41：基于血流速度与三维血管网格模型，求解不可压缩流的基本公式，对三维血管网格模型进行求解，用数值法求解连续性和Navier-Stokes方程：其中P，ρ，μ分别为流速、压力、血流密度、血流粘性；入口边界条件为血流速度，出口边界条件为out-flow边界条件；S42：计算沿着血管中心线从入口到下游各点的压力降ΔP。与现有技术相比，本专利技术的优点是：针本文档来自技高网...

【技术保护点】
1.一种基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，其特征在于，包括以下步骤：/nS01：通过血压传感器测量在舒张期心脏冠脉口的压力Pa；/nS02：通过造影图像获取血管的二维管径及长度，并通过两个呈30°以上夹角的造影图像生成三维血管网格模型并获得血管的三维管径及长度；/nS03：在心脏舒张期，测量包含造影剂的血液从一段指定血管的起始点到结束点所用的时间，并根据该时间和血管三维长度计算血流速度V1；/nS04：将计算得到的造影状态下的血流速度V1作为冠脉入口流速，计算冠脉入口到冠脉狭窄远端的压力降ΔP，在舒张期狭窄远端冠状动脉内平均压Pd＝Pa-ΔP，通过公式caiFR≌caDFR≌cadFR＝Pd/Pa计算得到造影瞬时无波型比率(caiFR)和造影静息态舒张期压力比率(cadFR和caDFR)。/n

【技术特征摘要】
1.一种基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，其特征在于，包括以下步骤：
S01：通过血压传感器测量在舒张期心脏冠脉口的压力Pa；
S02：通过造影图像获取血管的二维管径及长度，并通过两个呈30°以上夹角的造影图像生成三维血管网格模型并获得血管的三维管径及长度；
S03：在心脏舒张期，测量包含造影剂的血液从一段指定血管的起始点到结束点所用的时间，并根据该时间和血管三维长度计算血流速度V1；
S04：将计算得到的造影状态下的血流速度V1作为冠脉入口流速，计算冠脉入口到冠脉狭窄远端的压力降ΔP，在舒张期狭窄远端冠状动脉内平均压Pd＝Pa-ΔP，通过公式caiFR≌caDFR≌cadFR＝Pd/Pa计算得到造影瞬时无波型比率(caiFR)和造影静息态舒张期压力比率(cadFR和caDFR)。


2.根据权利要求1所述的基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，其特征在于，所述步骤S01包括，使用血压传感器的压力管连接到多联三通，然后通过造影导管与心脏冠脉口部相连，在血压传感器的压力管内充满盐水，并保持血压传感器与心脏在同一水平位置，该血压传感器测量的压力波即为心脏冠脉口的压力波，在舒张期，瞬时压力的平均值即为Pa。


3.根据权利要求1所述的基于造影图像计算造影瞬时无波型比率和造影舒张期压力比率的方法，其特征在于，所述步骤S02中生成三维血管网格模型的方法包括以下步骤：
S21：将两个呈30°以上夹角的X射线冠脉造影图像上，具有映射关系的两分段血管的2D结构数据进行...

【专利技术属性】
技术研发人员：霍云飞，刘广志，吴星云，
申请(专利权)人：苏州润迈德医疗科技有限公司，
类型：发明
国别省市：江苏;32