基于造影图像修正血流速度和微循环参数的方法及装置与流程

本发明涉及冠状动脉医学技术领域，特别是涉及一种基于造影图像修正血流速度和微循环参数的方法、装置冠状动脉分析系统及计算机存储介质。

背景技术：

冠状动脉微循环功能异常对于心肌缺血的影响逐渐得到关注，冠脉系统由心外膜冠脉以及微循环组成。

一般来说，心外膜冠脉狭窄程度大于等于50％可导致心肌供血不足，临床诊断为冠心病。但是，临床研究显示，冠脉微循环异常也有可能导致心肌供血不足。

冠脉微循环是指微动脉和微静脉之间的血液循环，是血液与组织细胞进行物质交换的场所。研究表明尽管患者经皮冠状动脉介入治疗术成功后冠脉血流达timi3级，但仍有近30％的病人出现微血管功能异常，导致预后不佳。因此，随着研究的不断深入，人们逐渐认识到冠脉微血管功能异常是许多心脏疾病病理生理的一个重要机制，有必要准确评估冠脉微循环的功能状态。

冠脉微循环阻力指数(indexofmicrocirculatoryresistanceimr)是评估冠脉微循环功能状况的指标。

造影剂本身也对血管有扩张作用，早在1959年，研究人员就发现给狗冠脉内注射造影剂后，冠脉血流量将增加60％，提示造影剂可以诱发冠脉微循环的部分充血。1985年有研究发现，在有冠脉临界病变的患者冠脉内注射76％的泛影葡胺可明显增加跨病变压差；1995年，研究者进一步确定冠脉内注射8ml碘克沙醇270可以获得最大血流量的59％，而冠脉内注射腺苷200ug后也仅得到最大血流量的94％；到了2003年，基本明确了造影剂的微循环充血作用，但略弱于其它经典血管扩张剂的充血作用。后续研究发现造影剂的渗透压作用可促发血管内皮细胞钾离子通道开放，进而引起冠脉微循环的扩张。

基于造影剂的这些药理作用，临床专家探讨以造影剂代替腺苷诱导微循环充血，具体操作类似于腺苷等冠脉内途径给药过程。

目前多数研究采用的造影剂剂量为5～10ml。注射造影剂后冠脉微循环由充血状态恢复至基线状态的时间平均为12-30s。进行冠脉检查时，术者会针对检测血管从不同的体位角度进行造影。由于每次调节造影机的c臂到规定的角度的时间不同，以及开始造影的时间也有所不同，导致每次造影流速会受上次造影是否恢复到基线状态的影响。

技术实现要素：

本发明提供了一种基于造影图像修正血流速度和微循环参数的方法及装置，以降低现有技术中上次造影是否恢复到基线状态对造影流速产生的影响的问题。

为实现上述目的，第一方面，本申请提供了一种基于造影图像修正静息态下血流速度的方法，包括：

在造影态下，获取冠脉入口至冠脉狭窄远端的平均血流速度vh；

获取相邻两次推注造影剂时的开始时间之差δt；

根据时间差δt，获得修正系数k；

根据所述修正系数k、所述血流速度vh，获得静息态下血流速度vj。

可选地，上述的基于造影图像修正血流速度的方法，所述根据所述修正系数k、所述血流速度vh，获得静息态下血流速度vj的方法包括：

根据公式vj＝vh/k，获得静息态下血流速度vj。

可选地，上述的基于造影图像修正血流速度的方法，所述根据时间差δt，获得修正系数k的方法包括：；

如果δt≥30s，则k＝1；

如果20s≤δt＜30s，则1＜k≤1.5；

如果10s＜δt＜20s，则1.5＜k＜2.0；

如果δt≤10s，则k＝2。

可选地，上述的基于造影图像修正静息态下血流速度的方法，所述在造影态下，获取冠脉入口至冠脉狭窄远端的平均血流速度vh的方法包括：

获取所述心跳周期区域内包含的冠状动脉造影图像帧数；

其中，l表示心跳周期区域内造影剂流过的血管长度，n表示心跳周期区域包含的冠状动脉造影图像帧数，fps表示画面每秒传输帧数。

可选地，上述的基于造影图像修正静息态下血流速度的方法，l的取值范围为50～150mm；或l＝100mm。

可选地，上述的基于造影图像修正静息态下血流速度的方法，测量所述平均血流速度vh的方法包括：造影剂遍历距离算法、stewart—hamilton算法、first—pass分布分析法、光流法或者流体连续法。

第二方面，本申请提供了一种基于造影图像修正最大扩张血流速度的方法，包括：

上述的基于造影图像修正静息态下血流速度的方法；

根据静息态下血流速度vj获取最大扩张血流速度。

可选地，上述的基于造影图像修正最大扩张血流速度的方法，所述根据静息态下血流速度vj获取最大扩张血流速度的方法包括：

根据公式vmax＝avj+b；

其中，vmax表示最大扩张血流速度，a表示取值范围为1～3的常数、b表示取值范围为50～300的常数。

第三方面，本申请提供了一种基于造影图像修正冠状动脉微循环血管评定参数的方法，包括：

根据造影图像，获取心跳周期区域内的冠脉入口压力的平均值pa；

获取冠脉入口到冠脉狭窄远端的压力降δp；

根据上述的基于造影图像修正静息态下血流速度的方法获得的最大扩张血流速度vmax，以及δp、pa获得修正后的冠状动脉微循环血管评定参数。

第四方面，本申请提供了一种基于造影图像修正血流速度的装置，用于上述的基于造影图像修正静息态下血流速度的方法，包括：第一血流速度单元、时间差单元、修正系数单元和第二血流速度单元；所述第一血流速度单元与所述第二血流速度单元连接，所述修正系数单元分别与所述时间差单元所述第二血流速度单元连接；

所述第一血流速度单元，用于在造影态下，获取冠脉入口至冠脉狭窄远端的平均血流速度vh；

所述时间差单元，用于获取相邻两次推注造影剂时的开始时间之差δt；

所述修正系数单元，用于接收所述时间差单元传送的时间差δt，获得修正系数k；

所述第二血流速度单元，用于接收所述第一血流速度单元发送的造影态下的平均血流速度vh，以及接收所述修正系数单元发送的修正系数k，根据所述修正系数k、所述血流速度vh，获得静息态下血流速度vj。

第五方面，本申请提供了一种基于造影图像修正最大扩张血流速度的装置，用于上述的基于造影图像修正最大扩张血流速度的方法，包括：上述的基于造影图像修正血流速度的装置，以及与上述的基于造影图像修正血流速度的装置连接的第三血流速度单元；

所述第三血流速度单元，用于根据静息态下血流速度vj获取最大扩张血流速度。

第六方面，本申请提供了一种冠状动脉分析系统，包括：基体，以及均设置于所述基体上的血压采集装置和上述的基于造影图像修正最大扩张血流速度的装置。

第七方面，本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的基于造影图像修正静息态下血流速度的方法。

本申请实施例提供的方案带来的有益效果至少包括：

本申请提供了基于造影图像修正血流速度的方法，根据时间差δt，获得修正系数k；根据修正系数k、血流速度vh，获得静息态下血流速度vj，降低了现有技术中上次造影是否恢复到基线状态对造影流速产生的影响。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本发明的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1为本申请的基于造影图像修正血流速度的方法的实施例1的流程图；

图2为本申请的基于造影图像修正冠状动脉微循环血管评定参数的方法的流程图；

图3为参考图像；

图4为待分割的一目标图像；

图5为待分割的另一目标图像；

图6为增强后的导管图像；

图7为导管特征点的二值化图像；

图8为增强后的目标图像；

图9为冠状动脉所处位置的区域图像；

图10为结果图像；

图11为横切面截图；

图12为纵切面截图；

图13为两个体位造影图像；

图14的左图为血管长度与血管直径的曲线图；

图15为由图14结合体位角度以及冠脉中心线生成的冠状动脉三维结构图；

图16为分割图像的帧数个数图；

图17为冠脉入口压力测试图；

图18为imr测试图；

图19为基于造影图像修正血流速度的装置的结构框图；

图20为基于造影图像修正血流速度的装置的另一结构框图；

图21为三维建模装置的结构框图；

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明具体实施例及相应的附图对本发明技术方案进行清楚、完整地描述。显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

以下将以图式揭露本发明的多个实施方式，为明确说明起见，许多实务上的细节将在以下叙述中一并说明。然而，应了解到，这些实务上的细节不应用以限制本发明。也就是说，在本发明的部分实施方式中，这些实务上的细节是非必要的。此外，为简化图式起见，一些习知惯用的结构与组件在图式中将以简单的示意的方式绘示之。

目前多数研究采用的造影剂剂量为5～10ml。注射造影剂后冠脉微循环由充血状态恢复至基线状态的时间平均为12-30s。进行冠脉检查时，术者会针对检测血管从不同的体位角度进行造影。由于每次调节造影机的c臂到规定的角度的时间不同，以及开始造影的时间也有所不同，导致每次造影流速会受上次造影是否恢复到基线状态的影响。

实施例1：

为了解决每次造影流速会受上次造影是否恢复到基线状态的影响的问题，如图1所示，本申请提供了一种基于造影图像修正血流速度的方法，包括：

s100，在造影态下，获取冠脉入口至冠脉狭窄远端的平均血流速度vh；

s200，获取相邻两次推注造影剂时的开始时间之差δt；

s300，根据时间差δt获得修正系数k；

s400，根据修正系数k、血流速度vh，获得静息态下血流速度vj，具体公式为vj＝vh/k。

本申请的一个实施例中，δt分为4种情况对k产生影响，具体为：

(1)如果δt≥30s，则k＝1；

(2)如果20s≤δt＜30s，则1＜k≤1.5；

(3)如果10s＜δt＜20s，则1.5＜k＜2.0；

(4)如果δt≤10s，则k＝2。

本申请提供了基于造影图像修正血流速度的方法，根据时间差δt获得修正系数k；根据修正系数k、血流速度vh，获得静息态下血流速度vj，降低了现有技术中上次造影是否恢复到基线状态对造影流速产生的影响。

本申请的一个实施例中，s100的方法包括：

如果采用造影剂运输时间算法获取vh，则：获取心跳周期区域内包含的冠状动脉造影图像帧数，以及获取心跳周期区域内造影剂流过的血管长度；

根据公式计算vh；

其中，l表示心跳周期区域内造影剂流过的血管长度；n表示心跳周期区域包含的冠状动脉造影图像帧数，fps表示画面每秒传输帧数，优选地，fps＝15帧/秒；

本申请的一个实施例中，测量平均血流速度的方法包括：造影剂遍历距离算法、stewart—hamilton算法、first—pass分布分析法、光流法或者流体连续法。

本申请的一个实施例中，l的取值范围为50～150mm；或l＝100mm。

实施例2：

如图2所示，本申请提供了一种基于造影图像修正最大扩张血流速度的方法，包括：

上述的基于造影图像修正静息态下血流速度的方法；

s500，根据静息态下血流速度vj获取最大扩张血流速度，包括：根据公式vmax＝avj+b；其中，vmax表示最大扩张血流速度，a表示取值范围为1～3的常数、b表示取值范围为50～300的常数。

实施例2：

如图2所示，本申请提供了一种基于造影图像修正冠状动脉微循环血管评定参数的方法，包括：

s001，根据造影图像，获取心跳周期区域内的冠脉入口压力的平均值pa，具体方法为：通过血压采集装置实时测量pa；

s002，获取冠脉入口到冠脉狭窄远端的压力降δp，包括：

a，提取至少两个体位的冠脉造影图像；优选地，两个体位的拍摄角度＞30°；

b，对冠脉造影图像去噪，包括：静态噪声和动态噪声；

静态噪声为在时间中静止不变的噪声，如胸腔中的肋骨。

动态噪声为在时间中变化的噪声，如部分肺部组织、部分心脏组织，采用均值滤波，除去部分动态噪声；

以及包括：通过灰度直方图分析，利用阈值进一步去噪。

c，去除冠脉造影图像的干扰血管，得到如图10所示的结果图像，包括：

将有导管出现的第一帧分割图像定义为如图3所示的参考图像，将有完整冠状动脉出现的第k帧分割图像定义为如图4和图5所示的目标图像，k为大于1的正整数；

将如图3所示的参考图像减去如图4和图5所示的目标图像，提取导管的特征点o；对去噪后的图像进行图像增强；对如图6所示的增强后的导管图像进行二值化处理，得到如图3所示的具有一组导管特征点o的二值化图像；

将如图4和图5所示的目标图像减去如图3所示的参考图像；去噪，包括：静态噪声和动态噪声；采用多尺度海森矩阵对去噪后的图像进行图像增强；根据如图8所示的增强后的目标图像中各区域与导管特征点的位置关系，确定并提取冠状动脉的区域，即为如图9所示的冠状动脉所处位置的区域图像；

对如图9所示的冠状动脉所处位置的区域图像进行二值化处理，获得二值化冠状动脉图像；

对二值化冠状动脉图像进行形态学运算，以导管的特征点作为种子点，二值化冠状动脉图像依据种子点所处位置进行动态区域生长，获得如图10所示的结果图像；

d，沿着冠状动脉的延伸方向，提取每幅结果图像的冠脉中心线和直径；

e，将每根冠脉中心线和直径均投射于三维空间上进行三维建模，获得冠状动脉三维结构，包括：获取每幅冠脉造影图像的体位拍摄角度；将每根冠脉中心线结合体位拍摄角度以及血管长度l值和血管直径d值投射于三维空间上，生成冠状动脉三维结构；

f，对冠状动脉三维结构进行网格划分，如图11和图12所示，基于重构的冠状动脉三维结构，本申请的一个实施例采用标准的扫掠法进行网格划分，生成结构性三维六面体网格；进一步地，基于重构的冠脉三维模型，本申请也可以采用其它方法(例如：切分法、混合法)进行网格划分，生成结构性三维六面体网格；

g，以冠脉中心线作为纵轴，网格沿所述冠脉中心线划分为m个点，所述冠脉中心线每个点对应的横截面被划分为n个节点，δpi表示所述冠脉中心线上第i个点的横截面上所有节点的压力的平均值，即为冠脉入口到冠脉狭窄远端的压力降δp；

所述压力降δpi采用如下公式计算：

p1表示三维结构网格中第i个点的横截面上的第一个节点的压力值，p2表示三维结构网格中第i个点的横截面上的第二个节点的压力值，pn表示第i个点的横截面上的第n个节点的压力值，m、n均为正整数；pn的压力值采用纳维-斯托克斯方程计算得出；

s003，根据实施例1～3最大扩张血流速度vmax，以及δp、pa获得修正后的冠状动脉微循环血管评定参数。

如果冠状动脉微循环血管评定参数为微循环阻力指数imr，则imr＝(pa-δp)×l/vmax。

本申请通过最大扩张血流速度vmax获得的imr值更加准确，降低了上次造影时间以及推注造影剂时的推注压力对imr值计算准确性的影响。

实施例6：

如图19所示，本申请提供了一种基于造影图像修正血流速度的装置，包括：第一血流速度单元100、时间差单元200、修正系数单元400和第二血流速度单元500；第一血流速度单元100与第二血流速度单元500连接，修正系数单元400分别与时间差单元200、第二血流速度单元500连接；第一血流速度单元100用于在造影态下，获取冠脉入口至冠脉狭窄远端的平均血流速度vh；时间差单元200用于获取相邻两次推注造影剂时的开始时间之差δt；修正系数单元400用于接收时间差单元200传送的时间差δt，获得修正系数k；第二血流速度单元500用于接收第一血流速度单元100发送的平均血流速度vh，以及接收修正系数单元400发送的修正系数k，根据修正系数k、血流速度vh，获得静息态血流速度vj。

如图20所示，本申请的一个实施例中，还包括：与第一血流速度单元100连接的三维建模装置600，三维建模装置用于读取冠状动脉造影图像，选取冠状动脉造影图像的一个心跳周期区域，测量心跳周期区域内的血管的长度l，进行三维建模，获得冠状动脉三维结构。

如图21所示，本申请的一个实施例中，三维建模装置600包括图像读取模块610、分割模块620、血管长度测量模块630和三维建模模块640，分割模块620与图像读取模块610、血管长度测量模块630、三维建模模块640连接，血管长度测量模块630与第一血流速度单元100连接；图像读取模块610用于读取造影图像；分割模块620用于选取冠状动脉造影图像的一个心跳周期区域；血管长度测量模块630用于测量心跳周期区域内的血管的长度l，并将血管的长度l传递给第一血流速度单元10；三维建模模块640用于根据分割模块620选取的冠脉造影图像进行三维建模，获得冠状动脉三维结构。

实施例7：

如图21所示，本申请提供了一种基于造影图像修正最大扩张血流速度的装置，包括：实施例6中的基于造影图像修正血流速度的装置，以及与上述的基于造影图像修正血流速度的装置连接的第三血流速度单元700；第三血流速度单元400，用于根据静息态下血流速度vj获取最大扩张血流速度vmax。

本申请的一个实施例中，还包括：与所述第三血流速度单元400连接的冠状动脉微循环血管评定参数测量装置，以及与所述冠状动脉微循环血管评定参数测量装置连接的压力降测量模块。

下面结合具体样例对本申请进行具体阐述：

如图13所示，为一位患者拍摄的两个体位的冠脉造影图像；左图的体位角度为右前斜rao：25°和头位cra：23°；右图的体位角度为右前斜rao：3°和头位cra：30°；

如图14所示，冠状动脉三维结构的血管长度l值＝120mm；生成的冠状动脉三维结构如图15所示；

血管直径d值＝2～4mm；

如图16所示，

由于相邻两次推注造影剂的开始时间差为20s≤δt＜30s，因此，k取值1.1，因此，vj＝300/1.1＝272.7；

vmax＝272.7+295＝567.7

如图17所示，pa＝100mmhg；

如图18所示，δp＝7；因此imr＝(100-7)×120/567.7＝19.66；如果不修正，则计算得到的imr＝(100-7)×120/(300+295)＝18.75；

因此可知，通过系数k修正前与修正后的imr的测量结果之差为0.91，误差很大，因此采用系数对血流速度进行修正，获得更加准确的微循环血管评价参数是必要的，提高了测量结果的准确性。

本申请提供了一种冠状动脉分析系统，包括：基体，以及均设置于所述基体上的血压采集装置和上述的基于造影图像修正最大扩张血流速度的装置，压力降测量模块、血压采集装置、基于造影图像修正最大扩张血流速度的装置均与冠状动脉微循环血管评定参数测量装置连接。

本申请提供了一种计算机存储介质，计算机程序被处理器执行时实现上述的基于造影图像修正静息态下血流速度的方法。

所属技术领域的技术人员知道，本发明的各个方面可以实现为系统、方法或计算机程序产品。因此，本发明的各个方面可以具体实现为以下形式，即：完全的硬件实施方式、完全的软件实施方式(包括固件、驻留软件、微代码等)，或硬件和软件方面结合的实施方式，这里可以统称为“电路”、“模块”或“系统”。此外，在一些实施例中，本发明的各个方面还可以实现为在一个或多个计算机可读介质中的计算机程序产品的形式，该计算机可读介质中包含计算机可读的程序代码。本发明的实施例的方法和/或系统的实施方式可以涉及到手动地、自动地或以其组合的方式执行或完成所选任务。

例如，可以将用于执行根据本发明的实施例的所选任务的硬件实现为芯片或电路。作为软件，可以将根据本发明的实施例的所选任务实现为由计算机使用任何适当操作系统执行的多个软件指令。在本发明的示例性实施例中，由数据处理器来执行如本文的根据方法和/或系统的示例性实施例的一个或多个任务，诸如用于执行多个指令的计算平台。可选地，该数据处理器包括用于存储指令和/或数据的易失性储存器和/或用于存储指令和/或数据的非易失性储存器，例如，磁硬盘和/或可移动介质。可选地，也提供了一种网络连接。可选地也提供显示器和/或用户输入设备，诸如键盘或鼠标。

可利用一个或多个计算机可读的任何组合。计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举列表)将包括以下各项：

具有一个或多个导线的电连接、便携式计算机盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑盘只读存储器(cd-rom)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。

计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。

计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括(但不限于)无线、有线、光缆、rf等等，或者上述的任意合适的组合。

例如，可用一个或多个编程语言的任何组合来编写用于执行用于本发明的各方面的操作的计算机程序代码，包括诸如java、smalltalk、c++等面向对象编程语言和常规过程编程语言，诸如"c"编程语言或类似编程语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络--包括局域网(lan)或广域网(wan)-连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。

应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机程序指令实现。这些计算机程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些计算机程序指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。

也可以把这些计算机程序指令存储在计算机可读介质中，这些指令使得计算机、其它可编程数据处理装置、或其它设备以特定方式工作，从而，存储在计算机可读介质中的指令就产生出包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的指令的制造品(articleofmanufacture)。

还可将计算机程序指令加载到计算机(例如，冠状动脉分析系统)或其它可编程数据处理设备上以促使在计算机、其它可编程数据处理设备或其它设备上执行一系列操作步骤以产生计算机实现过程，使得在计算机、其它可编程装置或其它设备上执行的指令提供用于实现在流程图和/或一个或多个框图方框中指定的功能/动作的过程。

本发明的以上的具体实例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

起点商标作为专业知识产权交易平台，可以帮助大家解决很多问题，如果大家想要了解更多知产交易信息请点击 【在线咨询】或添加微信 【19522093243】与客服一对一沟通，为大家解决相关问题。

此文章来源于网络,如有侵权,请联系删除