阅读说明：本技术 归一化冠状动脉微循环阻力指数计算方法 (Normalized coronary artery microcirculation resistance index calculation method ) 是由 谢辛舟 廉坤 陈焱 李成祥 张瑞晨 赵帅 郭晓雄 李玮昊 谢松云 于 2021-07-23 设计创作，主要内容包括：本发明提出了归一化冠状动脉微循环阻力指数计算方法,包括：测量感兴趣血管最大充血状态下的远、近端压力值；基于造影图像重建感兴趣血管三维模型,同时提取血管远端截面面积；仿真计算感兴趣血管血流阻力模型参数；结合前述压力值和血流阻力模型参数计算对应血流量；基于前述压力值、血流量和血管远端截面面积计算归一化冠状动脉微循环阻力指数。本发明基于造影图像和侵入式测量的压力值,提出归一化冠状动脉微循环阻力指数计算方法,以使得可在冠状动脉不同部位、不同大小的血管间建立统一的微循环阻力计算指标。(The invention provides a method for calculating a normalized coronary artery microcirculation resistance index, which comprises the following steps: measuring the far and near end pressure values of the blood vessel of interest in the maximal hyperemia state; reconstructing a three-dimensional model of the interested blood vessel based on the contrast image, and extracting the cross section area of the far end of the blood vessel; simulating and calculating blood flow resistance model parameters of the interested blood vessel; calculating corresponding blood flow by combining the pressure value and the blood flow resistance model parameter; and calculating a normalized coronary artery microcirculation resistance index based on the pressure value, the blood flow and the cross-sectional area of the far end of the blood vessel. The invention provides a normalized coronary artery microcirculation resistance index calculation method based on a contrast image and a pressure value measured in an invasive manner, so that a unified microcirculation resistance calculation index can be established among blood vessels with different sizes and different parts of a coronary artery.)

归一化冠状动脉微循环阻力指数计算方法

技术领域

本发明涉及计算冠状动脉生理学领域，尤其涉及应用血流动力学仿真计算方法计算冠状动脉微循环阻力的方法。

背景技术

本发明是基于本人已授权专利(“基于冠状动脉CT血管造影的狭窄病灶血流阻力计算方法”，ZL201710505979.1)的专利发明的发明。

冠状动脉微循环阻力(microcirculatory resistance,MR)是反映冠状动脉微循环功能的有效指标，其定义为冠脉远端动脉压力(P_d)除以最大充血状态下的冠脉血流量(Q_max)。目前可通过以下方法计算：

(1)第一类方法通过集成了压力传感器和血流量传感器的导丝置入冠状动脉远端，在药物诱导最大充血状态的条件下同时测量压力和血流量，进而计算MR(Williams,R.P.,et al.(2018)."Doppler Versus Thermodilution-Derived CoronaryMicrovascular Resistance to Predict Coronary Microvascular Dysfunction inPatients With Acute Myocardial Infarction or Stable Angina Pectoris."Am JCardiol 121(1):1-8.)；

(2)通过序列造影图像估算Q_max，基于血流动力学模型仿真计算最大充血状态下远、近端压力差(ΔP)，基于测量的静息态冠脉入口压力和深度学习模型估算最大充血状态下冠脉入口压力(P_a)(“快速计算微循环阻力的方法与系统”，中国201711258493.9[P])或直接测量最大充血状态下P_a(“基于造影图像和流体力学模型的微循环阻力指数计算方法”，中国201810413391.8[P])；

(3)结合侵入式压力导丝测量的充血态冠状动脉压力值和基于多角度冠脉造影图像重建的冠状动脉三维模型，采用计算流体力学方法仿真计算最大充血状态下的冠脉血流量(Q_max)(Morris PD,Gosling R,Zwierzak I,Evans H,Aubiniere-Robb L,Czechowicz K,et al.A novel method for measuring absolute coronary blood flow andmicrovascular resistance in patients with ischaemic heart disease.CardiovascRes.2021；117:1567-77.)。

但上述各种方法计算的微循环阻力均为绝对微循环阻力值，其在不同部位、不同大小的血管间数值分布差异巨大，难以直接为临床诊断冠状动脉微循环障碍提供精确的参考。

发明内容

有鉴于此，本发明提供一种归一化冠状动脉微循环阻力指数计算方法，该方法将侵入式压力导丝测量的远、近端压力值作为已知条件，基于多角度冠脉造影图像重建冠状动脉三维模型，利用血流阻力模型仿真计算冠脉血流量，并进一步基于压力值和重建三维模型中获取的血管远端截面面积，计算归一化微循环阻力指数。所采取的技术方案如下：

1.侵入式压力导丝测量感兴趣血管的远、近端压力值；

2.基于造影图像重建感兴趣血管三维模型，提取血管远端截面面积A；

3.通过数值方法计算感兴趣血管段血流阻力(压降-血流曲线)模型参数；

4.结合测量压力值和血流阻力模型参数计算冠脉血流量(Q_max)；

5.基于血管远端截面面积A和计算得到的Q_max，计算归一化微循环阻力指数MR_no_malized。

附图说明

图1.归一化冠状动脉微循环阻力指数计算方法的整体流程图。

具体实施方式

下面结合附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

归一化冠状动脉微循环阻力指数计算方法的整体流程图如图1所示。以下将结合图1对具体实施方式进行详细说明。

1.测量感兴趣血管最大充血状态下的远、近端压力值：通过药物诱导最大充血状态，利用压力导丝测量感兴趣血管远、近端压力值(P_a、P_d)；

2.基于造影图像重建感兴趣血管三维模型：多角度拍摄的冠脉造影图像进行分割，分别得到冠脉中心线及直径；通过三维空间投影计算生成冠脉三维模型，同时提取血管远端截面面积A；

3.基于重建的冠脉三维模型，采用本人已授权的专利发明方法(“基于冠状动脉CT血管造影的狭窄病灶血流阻力计算方法”，ZL201710505979.1)计算血流阻力(压降-血流曲线)模型参数f和s；

4.基于f、s和前述测量压力值P_a、P_d，求解如下一元二次方程，得到Q_max：

P_a-P_d＝f·Q_max+s·Q_max ² (1)

5.基于如下公式计算归一化冠状动脉微循环阻力指数：

其中MR_nomalized为归一化冠状动脉微循环阻力指数，P_d为前述1步骤中测量的远端压力P_d，Q_max为上述4步骤中计算的最大充血血流量，k为固定的模型参数，典型取值在2到3之间，r_equ为等效血管半径，其基于步骤2中提取的血管远端截面面积A按如下公式计算得到：

其中π为圆周率。

本发明的有益效果在于：

基于造影图像和侵入式测量的压力值，提供了归一化的微循环阻力指数计算方法，以使得可在冠状动脉不同部位、不同大小的血管间建立统一的微循环阻力计算指标。