CN109288537B - 评估血流储备分数的系统、方法、设备及存储介质 - Google Patents

Abstract

本发明实施例公开了一种评估血流储备分数的系统、方法、设备及计算机可读存储介质，其中，所述系统包括：获取模块，用于获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；狭窄率确定模块，用于根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；模型确定模块，用于根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；FFR计算模块，用于依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

Description

评估血流储备分数的系统、方法、设备及存储介质

技术领域

本发明涉及生物医学工程技术，尤其涉及一种血流储备分数的系统、方法、设备及计算机可读存储介质。

背景技术

冠状血流储备分数(Fractional flow reserve，FFR)现已经成为冠脉狭窄功能性评价的公认指标，其最重要的功能是对一个未知影响的冠脉狭窄的功能后果进行准确评价。现有基于冠状动脉非创伤性血管成像(CT angiography,CTA)的无创评估是通过基于CTA图像重建技术和计算流体力学(Computational Fluid Dynamic,CFD)多尺度数值模拟计算FFR，然而在此评估过程中CFD计算模型复杂，往往需要计算几百万个网格，因此，所需时间较长，不适用于临床快速诊断，同时，由于计算量大，对于设备性能较高，需要在工作站或服务器上完成，导致成本上升。

发明内容

为解决现有存在的技术问题，本发明实施例提供一种评估血流储备分数FFR的系统、方法、设备及计算机可读存储介质，能够快速完成FFR计算。

为达到上述目的，本发明实施例的技术方案是这样实现的：

第一方面，本发明实施例提供一种评估血流储备分数FFR的系统，所述系统包括：

获取模块、狭窄率确定模块、模型确定模块以及FFR计算模块，其中，

所述获取模块，用于获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；

所述狭窄率确定模块，用于根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；

所述模型确定模块，用于根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；

所述FFR计算模块，用于依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

其中，所述狭窄率确定模块还用于所述获取冠状动脉的三维图像信息包括：获取冠状动脉的CTA图像信息，依据所述CTA图像信息，构建所述冠状动脉的三维图像信息；所述系统还包括图像处理模块，用于对所述冠状动脉的所述CTA图像信息和/或所述三维图像信息进行斑块处理。

其中，所述狭窄率确定模块包括截面积获取单元以及狭窄率确定单元，其中，所述截面积获取单元，用于依据龙骨节点的信息参数，获取所述龙骨节点的冠状动脉截面积；

所述狭窄率确定单元，用于根据全部所述龙骨节点的冠状动脉截面积，确定出冠状动脉的狭窄率。

其中，所述狭窄率确定单元，具体用于确定目标无狭窄处的龙骨节点的冠状动脉截面积；根据所有所述龙骨节点的冠状动脉截面积与所述与目标无狭窄处的龙骨节点的冠状动脉截面积的关系，确定所述各龙骨节点的狭窄率；比较所述各龙骨节点的狭窄率，选取所述最大龙骨节点的狭窄率为冠状动脉的狭窄率。

其中，所述模型确定模块，具体用于若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；若冠状动脉的狭窄率大于第一预设狭窄率，确定评估FFR的模型为第二模型。

其中，所述模型确定模块，具体用于若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；若冠状动脉的狭窄率不小于第二预设狭窄率，确定评估FFR的模型为第二模型；若冠状动脉的狭窄率大于第一预设狭窄率且小于第二预设狭窄率，确定评估FFR的模型为第三模型。

其中，所述评估FFR的模型为第一模型，所述FFR计算模块包括入口信息获取单元、压力分布获取单元以及第一计算单元，其中，

所述入口信息获取单元，用于根据所述冠状动脉的三维图像信息，获取各分支入口的信息；

所述压力分布获取单元，用于依据所述依据龙骨节点信息参数以及各分支入口的信息，获取龙骨节点的压力分布信息；

所述第一计算单元，用于根据龙骨节点的压力分布信息，计算出检测冠状动脉的FFR值。

其中，所述评估FFR的模型为第二模型，所述FFR计算模块包括流量信息获取单元、阻抗确定单元、压力确定单元、第二计算单元：

所述流量信息获取单元，用于根据从各出口龙骨节点的压力以及各龙骨节点信息参数，计算出相应分叉点龙骨节点的压力相等的关系，以及依据质量守恒定律，确定出所述冠状动脉入口龙骨节点流量与各个出口龙骨节点流量的关系，获取各个出口龙骨节点流量信息；

所述各个出口龙骨节点流量信息的表达式为

A_iQ_i ²+B_iQ_i+C_i＝0

其中，A_i、B_i、C_i都是与龙骨节点信息参数有关的常数，Q_i为第i个出口龙骨节点流量信息；

所述阻抗确定单元，用于根据所述冠状动脉的三维图像信息，获取冠状动脉入口信息以及各分支信息，确定各出口阻抗；

所述压力确定单元，用于依据各个出口龙骨节点流量信息和各出口阻抗，确定各个出口龙骨节点压力信息；

所述第二计算单元，用于根据各个出口龙骨节点压力信息和冠状动脉所有龙骨节点的信息参数，计算出待检测冠状动脉的FFR值。

其中，所述评估FFR的模型为第三模型，所述FFR计算模块包括第一FFR计算单元、第二FFR计算单元以及第三FFR计算单元，其中，

所述第一FFR计算单元，用于依据龙骨节点信息参数以及第一模型，计算出待检测冠状动脉的FFR值FFR₁；

所述第二FFR计算单元，用于依据龙骨节点信息参数以及第二模型，计算出待检测冠状动脉的FFR值FFR₂；

所述第三FFR计算单元，用于计算FFR₁和FFR₂加权平均值，确定依据所述第三模型计算出待检测冠状动脉的FFR值FFR₃；

所述FFR₁和FFR₂加权平均的计算公式为：

FFR₁(1-ε)+FFR₂ε＝FFR₃

其中，ε＝(θ-θ₁)/(θ₂-θ₁)，θ为冠状动脉的实际狭窄率，θ₁为第一预设狭窄率，θ₂为第二预设狭窄率。

第二方面，本发明实施例提供一种评估血流储备分数的方法，所述方法包括：

获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；

根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；

根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；

依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

第三方面，本发明实施例提供一种评估血流储备分数FFR的设备，所述设备包括本发明任一实施例所提供所述的评估FFR系统以及显示装置，其中，所述显示装置用于显示FFR信息。

第四方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现本发明任一实施例所述评估FFR方法。

本发明实施例所提供的评估血流储备分数的系统、方法、设备及计算机可读存储介质，该系统包括获取模块、狭窄率确定模块、模型确定模块以及FFR计算模块，其中，所述获取模块用于获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数，所述狭窄率确定模块用于根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率，所述模型确定模块用于根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型，所述FFR计算模块，用于依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。如此，采用龙骨节点信息来计算FFR，够简化计算FFR的模型，提高FFR计算速度，同时基于不同的狭窄率，采用不同的评估模型，能够提高评估FFR值的准确性。

附图说明

图1为本发明一实施例中评估血流储备分数FFR系统的结构示意图；

图2为本发明一实施例中冠状动脉的三维图像；

图3为本发明一实施例中冠状动脉龙骨节点示意图；

图4为本发明一实施例中确认FFR模型方法的流程示意图；

图5为本发明一实施例中基于第二模型的冠状动脉的简单示意图；

图6为本发明另一实施例中确认FFR模型方法的流程示意图；

图7为本发明一实施例中评估血流储备分数FFR方法的流程示意图；

图8为本发明一实施例中评估血流储备分数FFR设备的组成示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对发明的具体技术方案做进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。

在此需要说明的是，所述冠状动脉从主动脉根部发出，包绕整个心脏，主要分右冠状动脉和左冠状动脉两支。右冠状动脉由主干和末端分成后侧支和后降支的两个分支组成，左冠状动脉由主干、前降支和回旋支组成。本发明实施例中所述冠状动脉的龙骨节点是指冠状动脉血管的中心线上离散的空间分布点，在实际计算中，将冠状动脉血液流动等效成有限个管流元，每个管流元的长度为相邻龙骨节点的距离，以管流元的沿程压力损失h_l和局部阻力损失h_m为基础计算整体压力损失。其中，所述沿程阻力损失由流体在等效截面直管内的摩擦阻力引起的，所述局部阻力损失h_m是由于流道形状改变、流速受到扰动、流动方向变化引起的。通常认为，每种损失都能成分体现出来而且独立地不受其他损失的影响，因此，压力损失是可以叠加的。所述沿程阻力损失h_l可以表示为：

式中，△p、λ、g分别为沿程压力损失、沿程阻力系数和重力加速度，对于层流来说，λ＝64/Re，雷诺数

ρ和μ分别为血液的密度和粘度系数，U为流速，

Q_i为流量，A为管道截面积，l为管道长度，d为管道直径。

所述局部阻力损失h_m可以表示为：

式中，U为流速，

Q_i为流量，A为管道截面积，g为重力加速度，ζ为局部阻力损失系数，△p为局部压力损失。

本发明实施例一方面，提供一种评估血流储备分数FFR系统，请参阅图1，所述系统包括获取模块101、狭窄率确定模块102、模型确定模块103以及FFR计算模块104，其中，

所述获取模块101，用于获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；

其中，所述FFR系统获取冠状动脉的三维图像信息，可以是指基于冠状动脉非创伤性血管成像(CT angiography,CTA)图像数据，建立冠状动脉的三维图像模型，从而获取冠状动脉的三维图像信息。例如采用计算机断层扫描CT设备对心脏血管部位造影，生成心脏血管CTA图像数据，采用三维重建技术，建立心脏血管的三维图像数据，利用心脏各部分灰度差异、空间位置以及是否为连通域等信息，识别出冠状动脉的三维图像模型，请参阅图2，从而获取冠状动脉的图像信息。

所述FFR系统获取冠状动脉的三维图像信息，还可以是指通过特制的心导管进入冠状动脉，采用血管造影机进行冠状动脉造影，从而获取冠状动脉的三维图像信息。

其中，所述龙骨节点的信息参数至少包括龙骨节点索引、血流方向、空间分布参数、冠状动脉的等效直径、冠状动脉的曲率半径以及曲线切向方向等几何信息。

这里，所述依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数，具体是先基于冠状动脉的三维图像信息，基于图像形态学的细化算法，查找到冠脉血管的中心线，经过离散化计算，得到龙骨节点分布信息，请参阅图3，图中黑点代表冠状动脉的龙骨节点，然后根据龙骨节点分布信息，计算冠状动脉在垂直截面上的截面积和等效直径、曲率半径和曲线切向方向，以及血流方向等几何信息。

需要说明的是，若龙骨节点的总数是N，第i个节点的空间位置用p_i＝(x_i,y_i,z_i)，所述血流方向v_i可由p_i的一阶导数得到，所述截面积A_i指冠脉血管的垂直截面所包围的面积，所述周长C_i指冠脉血管的垂直截面的周长，所述等效直径，即水力直径d_i可用d_i＝4A_i/C_i，曲率半径r_i和曲线切向方向t_i可采用三维空间曲线拟合计算。

所述狭窄率确定模块102，用于根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；

这里，所述狭窄率确定模块包括截面积获取单元以及狭窄率确定单元，其中，所述截面积获取单元，用于依据龙骨节点的信息参数，获取所述龙骨节点的冠状动脉截面积，由于龙骨节点的信息参数中包括了冠状动脉截面积，因此可以直接获取冠状动脉截面积；

所述狭窄率确定单元，用于根据全部所述龙骨节点的冠状动脉截面积，确定出冠状动脉的狭窄率。具体地，确定目标无狭窄处的龙骨节点的冠状动脉截面积；根据所有所述龙骨节点的冠状动脉截面积与所述目标无狭窄处的龙骨节点的冠状动脉截面积的关系，确定所述各龙骨节点的狭窄率；比较所述各龙骨节点的狭窄率，选取所述最大龙骨节点的狭窄率为冠状动脉的狭窄率。

例如，当前龙骨节点的冠状动脉截面积为A_i，目标无狭窄处的龙骨节点的冠状动脉截面积A₀，所述当前龙骨节点的狭窄率θ_i为

按照此方法计算出所有龙骨节点的狭窄率，统计所有龙骨节点的狭窄率，选取最大的龙骨节点的狭窄率作为冠状动脉的狭窄率θ，其中，进行统计时可采取，当龙骨节点的狭窄率小于预设值时，不做统计。

这里，所述目标无狭窄处的龙骨节点是指当前龙骨节点血流反方向上的最近无狭窄处的龙骨节点。

在一个实施例中，所述目标无狭窄处的龙骨节点的方法包括：识别冠状动脉的狭窄区域和非狭窄区域，若识别当前龙骨节点位于狭窄区域，从当前龙骨节点沿血流反方向查找最近非狭窄区域的龙骨节点，该龙骨节点即为所述目标无狭窄处的龙骨节点；若识别当前龙骨节点位于非狭窄区域，从当前龙骨节点沿血流反方向查找上一个位于非狭窄区域龙骨节点，该龙骨节点即为所述目标无狭窄处的龙骨节点。

在另一个实施例中，所述确定目标无狭窄处的龙骨节点的方法还可以包括：比较分支各龙骨节点截面积与分支入口龙骨节点的截面积，若当前龙骨节点的截面积与入口龙骨节点的截面积的差值超出阈值，则沿血流反方向上查找出与入口龙骨节点的截面积的差值在阈值范围以内的最近龙骨节点，所述最近的龙骨节点即为目标无狭窄处的龙骨节点；若当前龙骨节点的截面积与入口龙骨节点的截面积的未超出阈值，则沿血流反方向上查找上一个龙骨节点，该龙骨节点即为目标无狭窄处的龙骨节点。

所述模型确定模块103，用于根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；

其中，在冠状动脉狭窄率较小时，可以通过心肌的调节功能来保证心肌的所需供血量一定，若冠状动脉狭窄率超出一定范围时，由于心肌的调节功能达到极限，从而无法保证心肌供血量不变，因此，评估FFR的模型时可采用不同的模型。

所述FFR计算模块104，用于依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

这里，根据狭窄率选择的相应的评估FFR模型，再依据冠状动脉的龙骨节点信息参数，计算各龙骨节点的压力值，依据FFR定义为最大充血状态下，存在狭窄病变时血管所能获得的最大血流量与正常状态下即无狭窄病变时血管所能获得的最大血流量的比值，可以等效成压力之比，因此，计算每个龙骨节点上的FFR_i，计算公式如下：

式中，p₀为冠状动脉入口的压力，p_i为第i个龙骨节点压力。

在本发明上述实施例中，该系统包括获取模块、狭窄率确定模块、模型确定模块以及FFR计算模块，其中，所述获取模块用于获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数，所述狭窄率确定模块用于根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率，所述模型确定模块用于根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型，所述FFR计算模块，用于依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。如此，采用龙骨节点信息来计算FFR，够简化计算FFR的模型，提高FFR计算速度，同时基于不同的狭窄率，采用不同的评估模型，能够提高评估FFR值的准确性。

在一个实施例中，所述预设狭窄率包括第一预设狭窄率，所述模型确定模块103，具体用于若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；若冠状动脉的狭窄率大于第一预设狭窄率，确定评估FFR的模型为第二模型。其中，所述第一模型可以为恒流模型，第二模型可以为恒阻模型，第一预设狭窄率可以是65％-75％，最优的是70％。请参阅图4，以评估FFR的模型包括恒流模型和恒阻模型、预设狭窄率为70％为例，所述评估血流储备分数FFR的系统进行血流储备分数FFR评估的方法包括：

步骤401、获取冠状动脉的三维图像信息；

步骤402、依据所述三维图像信息，确定冠状动脉的龙骨节点分布信息；

步骤403、根据冠状动脉的龙骨节点分布信息，获取冠状动脉的龙骨节点的信息参数；

步骤404、依据所述冠状动脉的龙骨节点的参数信息，确定冠状动脉的狭窄率θ；

步骤405、比较冠状动脉的狭窄率θ是否小于70％，若是，跳转406a，若不是，跳转到406b；

步骤406a、确定评估FFR模型为恒流模型；

步骤406b、确定评估FFR模型为恒阻模型。

在本发明上述实施例中，以心肌调节功能的能力为界，可将FFR模型分为第一模型和第二模型，依据冠状动脉的实际狭窄率与预设狭窄率的关系，确定具体采用哪种评估FFR模型，从而实现在心肌调节功能有效范围内，采用第一模型评估FFR，在心肌调节功能达到极限时，采用第二模型评估FFR，进而减少由于未考虑心肌调节功能的影响，造成评估FFR的误差，提高了FFR计算精度。

在一个实施例中，所述评估FFR的模型为第一模型，所述FFR计算模块104包括入口信息获取单元、压力分布获取单元以及第一计算单元，其中，

所述入口信息获取单元，用于根据所述冠状动脉的三维图像信息，获取各分支入口的信息；

这里，首先获取冠状动脉的总流量Q₁，Q₁可通过基于心肌模型的流量评估方法得到；

具体的是，通过冠状动脉的三维图像信息估计出左心室质量信息，则冠状动脉总流量信息为Q₁＝Q₀m^3/4，其中，Q₁为冠状动脉的血液总流量，Q₀为系数。

需要说明，冠状动脉的总流量也可以通过其他方式获取，例如磁共振实际测量冠状动脉左右分支的血液体积得到，也可设定冠状动脉入口总流量的平均值。

然后，对每条冠状动脉分支进行流量分配；

这里，可以采用进入分支的血液流量与血管直径的三次方成正比的关系，利用依据冠状动脉各分支入口直径的比例，计算分别进入冠状动脉各分支的血液流量。所述冠状动脉的入口直径d_j可根据冠状动脉的三维图像信息来确定，也就是说分支入口流量Q_j∝d_j ³，下标j为分支的索引。

所述压力分布获取单元，用于依据所述依据龙骨节点信息参数以及各分支入口的信息，获取龙骨节点的压力分布信息；

这里，根据分支入口流量及龙骨节点信息参数，可得到龙骨上每个节点的流量q_i，根据前述沿程阻力损失公式和局部阻力损失公式，计算当前节点的压力损失h_il和局部压力损失h_im，再根据冠状动脉入口压力对每个龙骨节点的阻力损失和局部阻力损失进行逐点积分，即可得到各龙骨节点的压力分布信息。

例如，当前子节点是离冠状动脉入口第二近龙骨节点，当前子节点流量为q₁，父节点流量为q_1p，流量

管道截面积

管道半径为

距离

可计算出当前子节点沿程阻力损失h_1l和局部阻力损失h_1m，进而求出沿程压力损失p_1l和局部压力损失p_1m，若入口压力为p₀，则当前子节点压力为p₁＝p₀-p_1l-p_1m，同理，若计算出第三龙骨节点的沿程压力损失p_2l和局部压力损失p_2m，则第三龙骨节点的压力为p₂＝p₁-p_2l-p_2m，依次类推，可得出各龙骨节点的压力分布信息。

需要说明的是，所述父节点为当前龙骨节点在血流反方向上距离该节点最近的龙骨节点，所述子节点为当前龙骨节点，入口节点是没有父节点，出口节点是没有子节点，采用按照血流反方向搜索的方法可以查找到各分支出口节点到入口节点的分支路径，例如冠状动脉有M个出口，则对应有M条分支，则每条分支L_j都是部分龙骨节点的集合，j最大为M。

所述第一计算单元，用于根据龙骨节点的压力分布信息，计算冠状动脉的FFR值。

这里，FFR定义为最大充血状态下，存在狭窄病变时血管所能获得的最大血流量与正常状态下即无狭窄病变时血管所能获得的最大血流量的比值，可以等效成压力之比，因此，计算每个龙骨节点上的FFR_i，计算公式如下：

式中，p₀为冠状动脉入口的压力，p_i为第i个龙骨节点压力。

在上述实施例中，以第一模型为恒流模型为例，所述模型采用从入口压力逐点递推的方式，确定各龙骨节点压力分布信息，进而确定FFR值，减少了由于未考虑心肌调节功能的影响，造成评估FFR的误差，同时相比现有采用CFD多尺度数值模拟计算FFR的方法来说，简化了计算FFR的方法，提高了计算效率和计算精度。

在一个实施例中，所述评估FFR的模型为第二模型，所述FFR计算模块包括流量信息获取单元、阻抗确定单元、压力确定单元、第二计算单元：

所述流量信息获取单元，用于根据从各出口龙骨节点的压力以及各龙骨节点信息参数，计算出相应分叉点龙骨节点的压力相等的关系，以及依据质量守恒定律，确定出所述冠状动脉入口龙骨节点流量与各个出口龙骨节点流量的关系，获取各个出口龙骨节点流量信息；

所述阻抗确定单元，用于根据所述冠状动脉的三维图像信息，获取冠状动脉入口信息以及各分支信息，确定各出口阻抗；

所述压力确定单元，用于依据各个出口龙骨节点流量信息和各出口阻抗，确定各个出口龙骨节点压力信息；

这里，首先根据冠状动脉入口总流量与总压力的关系计算出口总阻抗R，所述冠状动脉入口总流量与总压力的关系为：

P＝QR+P_v

式中，P为冠状动脉入口压力，Q为冠状动脉总流量，p_v为远端压力，R为总阻抗。

然后按照各分支出口直径三次方的倒数分配阻抗R_i，所述各分支出口直径的确定方法与各分支入口直径确定方法一致，这里不再细述。

由于从各个分支出口压力依次相应累加各分支中各个管元的沿程压力损失和局部压力损失，直至各自对应的分叉点，计算得到相应分支的分叉点应该压力相等，其中，所述分支出口压力的计算公式为：

P_io＝Q_ioR_i+P_v

式中，P_io为第i个分支出口压力，Q_io为第i个分支出口流量即各个出口龙骨节点流量信息，R_i为第i个出口阻抗，P_v为远端压力。

以及根据质量守恒定律，分支入口流量应该等于分叉点流量，以及等于各分叉点的分支出口的流量相加。

最后，根据整理出关于各个出口龙骨节点流量信息。

所述各个出口龙骨节点流量信息的表达式为

A_iQ_io ²+B_iQ_io+C_i＝0

其中，A_i、B_i、C_i都是与龙骨节点信息参数有关的常数，Q_io为第i个出口龙骨节点流量信息。

例如，将冠状动脉简化为只有一个分叉点和两个分支出口的模型，如图4所述，第1个点为冠状动脉入口龙骨节点，压力为P₁，流量为Q₁，第2个点为分叉点，第2个点为分叉龙骨节点，压力为P₂，流量为Q₂，第3个点为第1个出口龙骨节点，压力为P₃，流量为Q₃，出口阻抗R₃,第4个点为第2个出口龙骨节点，压力为P₄，出口阻抗R₄,流量为Q₄，其中，P₁和Q₁，R₃和R₄已知，由第3个点P₃依次累加沿程压力损失和局部压力损失，得到P₂：

式中，λ₂₃为沿程阻力系数，ζ₂₃为局部阻力系数，l₂₃为第2个龙骨节点到第3个龙骨节点的长度，d₂₃为平均直径，A₂₃为平均截面积，P_v为远端压力，ρ为血液密度。

由第4个点P₄依次累加沿程压力损失和局部压力损失，得到P₂：

式中，λ₂₄为沿程阻力系数，ζ₂₄为局部阻力系数，l₂₄为第2个龙骨节点到第4个龙骨节点的长度，d₂₄为平均直径，A₂₄为平均截面积，P_v为远端压力，ρ为血液密度。

由于两个出口龙骨节点求得的P₂应该相等，则有

根据质量守恒，有Q₁＝Q₂＝Q₃+Q₄,由分叉点龙骨节点压力P₂累加沿程压力损失和局部压力损失，可得P₁，即

上述两式子最后可改写为A Q ²+B Q+C＝0方程组形式，求解方程组即可求的Q₃和Q₄，按照公式P_io＝Q_ioR_i+P_v，可求出P₃＝Q₃R₃+P_v和P₄＝Q₄R₄+P_v，

按照公式(1)，可求出P₂。

需要说明的是，实际冠状动脉的形状比较复杂，但最终各出口龙骨节点流量都可以表示成

A_iQ_io ²+B_iQ_io+C_i＝0

其中，A_i、B_i、C_i都是与龙骨节点信息参数有关的常数，Q_io为第i个出口龙骨节点流量信息，各个分支的各个龙骨节点的压力值也是由各分支出口龙骨节点依次累加沿程压力损失局部压力损失求得。

所述第二计算单元，用于根据各个出口龙骨节点压力信息和冠状动脉所有龙骨节点的信息参数，计算冠状动脉的FFR值。

这里，FFR定义为最大充血状态下，存在狭窄病变时血管所能获得的最大血流量与正常状态下即无狭窄病变时血管所能获得的最大血流量的比值，可以等效成压力之比，因此，计算每个龙骨节点上的FFR_i，计算公式如下：

式中，p₀为冠状动脉入口的压力，p_i为第i个龙骨节点压力。

在本发明实施例中，以第二模型为恒阻模型为例，所述模型采用计算出各出口流量，从而确定各出口压力，从出口逐点递推回入口的方式，进而确定各龙骨节点的压力分布信息，最终确定FFR值，减少了由于心肌调节功能达到极限时仍采用第一模型造成的计算误差，同时，相比现有采用CFD多尺度数值模拟计算FFR的方法来说，简化了计算FFR的方法，提高了计算效率和计算精度提高了计算效率和计算精度。

在另一个实施例中，所述预设狭窄率包括第一预设狭窄率和第二预设狭窄率，所述模型确定模块103，若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；若冠状动脉的狭窄率不小于第二预设狭窄率，确定评估FFR的模型为第二模型；若冠状动脉的狭窄率大于第一预设狭窄率且小于第二预设狭窄率，确定评估FFR的模型为第三模型。

需要说明的是，在本发明实施例中第二预设狭窄率大于第一预设狭窄率，第一预设狭窄率为65％-70％，第二预设狭窄率为70％-75％，优选地，第一狭窄率为65％，第二狭窄率为75％，所述第三模型为混合模型，即通过其他模型计算出各自的FFR值，然后根据对各自模型的信任度，采用加权的方式计算出FFR值。请参阅图6，以评估FFR模型包括恒流模型、混合模型和恒阻模型为例、狭窄率为65％和75％为例，所述评估血流储备分数FFR的系统进行血流储备分数FFR评估的方法包括：

步骤601、获取冠状动脉的三维图像信息；

步骤602、依据所述三维图像信息，确定冠状动脉的龙骨节点分布信息；

步骤603、根据冠状动脉的龙骨节点分布信息，获取冠状动脉的龙骨节点的信息参数；

步骤604、依据所述冠状动脉的龙骨节点的参数信息，确定冠状动脉的狭窄率θ；

步骤605、比较冠状动脉的狭窄率θ与预设狭窄率的大小，所述预设狭窄率为65％和75％，若θ≤65％，跳转606a，若65％<θ<75％，跳转到606b，若θ≥75％，跳转到606c；

步骤606a、确认评估FFR模型为恒流模型；

步骤606b、确认评估FFR模型为混合模型；

步骤606c、确认评估FFR模型为恒阻模型。具体地，所述评估FFR的模型为第三模型，所述FFR计算模块包括第一FFR计算单元、第二FFR计算单元以及第三FFR计算单元，其中，

所述第一FFR计算单元，用于依据龙骨节点信息参数以及第一模型，计算出待检测冠状动脉的FFR值FFR₁；

其中，依据龙骨节点信息参数以及第一模型，计算出待检测冠状动脉的FFR值FFR₁的方法请参照前述发明实施例中依据第一模型计算FFR的方法，这里不再细述。

所述第二FFR计算单元，用于依据龙骨节点信息参数以及第二模型，计算出待检测冠状动脉的FFR值FFR₂；

其中，依据龙骨节点信息参数以及第二模型，计算出待检测冠状动脉的FFR值FFR₂前述发明实施例中依据第二模型计算FFR的方法，这里不再细述。

所述第三FFR计算单元，用于计算FFR₁和FFR₂加权平均值，确定依据所述第三模型计算冠状动脉的FFR值FFR₃；

所述FFR₁和FFR₂加权平均的计算公式为：

FFR₁(1-ε)+FFR₂ε＝FFR₃

其中，ε＝(θ-θ₁)/(θ₂-θ₁)，θ为冠状动脉的实际狭窄率，θ₁为第一预设狭窄率，θ₂为第二预设狭窄率。

在本发明实施例中，通过设置第一预设狭窄率和第二预设狭窄率，如此，细化了狭窄率对应的阈值区间，与此同时将评估FFR模型分为三类，其中，由于在第一预设狭窄率与第二预设狭窄率区间，心肌调节功能逐渐达到极限，采用加权平均的第三模型评估FFR值能减少用单一的模型计算FFR值的误差，进一步提高FFR计算精度。

在一个实施例中，所述狭窄率确定模块102还用于所述获取冠状动脉的三维图像信息包括：获取冠状动脉的CTA图像信息，依据所述CTA图像信息，构建所述冠状动脉的三维图像信息；所述系统还包括图像处理模块，用于对所述冠状动脉的所述CTA图像信息和/或所述三维图像信息进行斑块处理。

这里，所述冠状动脉的图像信息包括冠状动脉的三维图像信息和CTA图像信息。由于冠状动脉血管内存在脂肪类物质，当血管内皮破损，脂肪类物质进入血管壁，导致血管壁增厚，血管腔狭窄，其中，将使得血管内皮隆起的部分称为“斑块”，由于斑块的存在，通过对斑块的处理，可以进一步提升FFR计算的准确性，避免将狭窄率大的血管误认为狭窄率小的血管。

其中，所述斑块处理是根据斑块与血管内壁的灰度值不同，对冠状动脉的图像信息进行斑块处理，即识别出比邻近区域异常亮的区域，将该区域标记为斑块，删除标记为斑块的图像信息，从而确定出无斑块的冠状动脉的图像信息，再根据无斑块的冠状动脉的图像信息生成龙骨节点信息；也可以根据原始的冠状动脉的三维图像信息，确定了龙骨节点分布信息，获取了所述龙骨节点的信息参数，取其中一个龙骨节点的信息参数，确定由所述龙骨节点确定的截面上所在血管壁附近区域的图像数据，比较该区域的图像数据灰度值比邻近区域灰度值，若超出一定范围，则标记该区域为斑块，遍历所有龙骨节点，删除所有标记为斑块的三维图像数据，获取新的冠状动脉的三维图像信息，重新确定冠状动脉的龙骨节点，获取新的所述龙骨节点的信息参数。

在上述实施例中，通过对冠状动脉的图像信息进行斑块处理，能够减小由于斑块存在，造成重构冠状动脉的三维图像信息的误差，从而提高龙骨节点的信息参数的精度，进而提高FFR计算精度。

本发明实施例另一方面，提供一种评估血流储备分数FFR方法，请参照图7，所述方法包括：

步骤701，获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；

步骤702，根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；

包括：依据龙骨节点的信息参数，获取所述龙骨节点的冠状动脉截面积；

根据全部所述龙骨节点的冠状动脉截面积，确定出冠状动脉的狭窄率，

具体地，确定目标无狭窄处的龙骨节点的冠状动脉截面积；根据所有所述龙骨节点的冠状动脉截面积与所述目标无狭窄处的龙骨节点的冠状动脉截面积的关系，确定所述各龙骨节点的狭窄率；比较所述各龙骨节点的狭窄率，选取所述最大龙骨节点的狭窄率为冠状动脉的狭窄率。

步骤703，根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；

其中，若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；若冠状动脉的狭窄率大于第一预设狭窄率，确定评估FFR的模型为第二模型。

步骤704，依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

这里，所述评估FFR的模型为第一模型，步骤704包括：

根据所述冠状动脉的三维图像信息，获取各分支入口的信息；

依据所述依据龙骨节点信息参数以及各分支入口的信息，获取龙骨节点的压力分布信息；

根据龙骨节点的压力分布信息，计算冠状动脉的FFR值。

这里，所述评估FFR的模型为第二模型，步骤704包括：

根据从各出口龙骨节点的压力以及各龙骨节点信息参数，计算出相应分叉点龙骨节点的压力相等的关系，以及依据质量守恒定律，确定出所述冠状动脉入口龙骨节点流量与各个出口龙骨节点流量的关系，获取各个出口龙骨节点流量信息；

所述各个出口龙骨节点流量信息的表达式为

A_iQ_i ²+B_iQ_i+C_i＝0

其中，A_i、B_i、C_i都是与龙骨节点信息参数有关的常数，Q_i为第i个出口龙骨节点流量信息；

根据所述冠状动脉的三维图像信息，获取冠状动脉入口信息以及各分支信息，确定各出口阻抗；

依据各个出口龙骨节点流量信息和各出口阻抗，确定各个出口龙骨节点压力信息；

根据各个出口龙骨节点压力信息和冠状动脉所有龙骨节点的信息参数，计算冠状动脉的FFR值。

在本发明实施例中，在本发明实施例中，提出采用龙骨节点信息来计算FFR，够简化计算FFR的模型，提高FFR计算速度，同时基于不同的基于不同的狭窄率，采用不同的评估模型，能够减少由于心肌调节功能造成FFR值误差，能够提高评估FFR值的准确性。

在一个实施例中，步骤701，获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数，还包括：所述获取冠状动脉的三维图像信息包括：获取冠状动脉的CTA图像信息，依据所述CTA图像信息，构建所述冠状动脉的三维图像信息；所述系统还包括图像处理模块，用于对所述冠状动脉的所述CTA图像信息和/或所述三维图像信息进行斑块处理。

在本发明实施例中，通过对冠状动脉的图像信息进行斑块处理，能够减小由于斑块存在，造成重构冠状动脉的三维图像信息的误差，从而提高龙骨节点的信息参数的精度，进而提高FFR计算精度。

在一个实施例中，步骤703，所述根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型，包括：

若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；若冠状动脉的狭窄率不小于第二预设狭窄率，确定评估FFR的模型为第二模型；若冠状动脉的狭窄率大于第一预设狭窄率且小于第二预设狭窄率，确定评估FFR的模型为第三模型。

这里，所述评估FFR的模型为第三模型，在步骤704中，所述依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值，包括：

依据龙骨节点信息参数以及第一模型，计算出待检测冠状动脉的FFR值FFR₁；

依据龙骨节点信息参数以及第二模型，计算出待检测冠状动脉的FFR值FFR₂；

计算FFR₁和FFR₂加权平均值，确定依据所述第三模型计算冠状动脉的FFR值FFR₃；

所述FFR₁和FFR₂加权平均的计算公式为：

FFR₁(1-ε)+FFR₂ε＝FFR₃

其中，ε＝(θ-θ₁)/(θ₂-θ₁)，θ为冠状动脉的实际狭窄率，θ₁为第一预设狭窄率，θ₂为第二预设狭窄率。

在本发明实施例中，通过设置第一预设狭窄率和第二预设狭窄率，如此，细化了狭窄率对应的阈值区间，与此同时将评估FFR模型分为三类，其中，由于在第一预设狭窄率与第二预设狭窄率区间，心肌调节功能逐渐达到极限，采用加权平均的第三模型评估FFR值能减少用单一的模型计算FFR值的误差，进一步提高FFR计算精度。

本发明实施例又一方面，提供一种评估血流储备分数FFR的设备，请参照图8，所述设备包括本发明任一实施例中提供评估血流储备分数FFR系统801和显示装置802，其中，所述显示装置所述显示装置用于显示FFR信息。

本发明实施例再一方面，还提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被处理器执行时实现本发明任一实施例提供评估血流储备分数FFR方法的步骤。

除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。

本领域普通技术人员可以理解，上文中所公开方法中的全部或某些步骤、系统、装置中的功能模块/单元可以被实施为软件、固件、硬件及其适当的组合。在硬件实施方式中，在以上描述中提及的功能模块/单元之间的划分不一定对应于物理组件的划分；例如，一个物理组件可以具有多个功能，或者一个功能或步骤可以由若干物理组件合作执行。某些物理组件或所有物理组件可以被实施为由处理器，如中央处理器、数字信号处理器或微处理器执行的软件，或者被实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读介质上，计算机可读介质可以包括计算机可读存储介质(或非暂时性介质)和通信介质(或暂时性介质)。如本领域普通技术人员公知的，术语计算机可读存储介质包括在用于存储信息(诸如计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机可读存储介质包括但不限于RAM、ROM、EEPROM、闪存或其他存储器技术、CD-ROM、数字多功能盘(DVD)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。此外，本领域普通技术人员公知的是，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何信息递送介质。

以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。本发明的保护范围应以所述权利要求的保护范围以准。

Claims (11)

1.一种评估FFR的系统，其特征在于，所述系统包括获取模块、狭窄率确定模块、模型确定模块以及FFR计算模块，其中，

所述获取模块，用于获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；

所述狭窄率确定模块，用于根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；

所述模型确定模块，用于根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；

所述模型确定模块，具体用于若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；

所述FFR计算模块，用于依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

2.如权利要求1所述系统，其特征在于，所述获取冠状动脉的三维图像信息包括：获取冠状动脉的CTA图像信息，依据所述CTA图像信息，构建所述冠状动脉的三维图像信息；所述系统还包括图像处理模块，用于对所述冠状动脉的所述CTA图像信息和/或所述三维图像信息进行斑块处理。

3.如权利要求1所述系统，其特征在于，所述狭窄率确定模块包括截面积获取单元以及狭窄率确定单元，其中，

所述截面积获取单元，用于依据龙骨节点的信息参数，获取所述龙骨节点的冠状动脉截面积；

所述狭窄率确定单元，用于根据全部所述龙骨节点的冠状动脉截面积，确定出冠状动脉的狭窄率。

4.如权利要求3所述系统，其特征在于，所述狭窄率确定单元，具体用于确定目标无狭窄处的龙骨节点的冠状动脉截面积；根据所有所述龙骨节点的冠状动脉截面积与所述目标无狭窄处相邻的龙骨节点的冠状动脉截面积的关系，确定所述各龙骨节点的狭窄率；比较所述各龙骨节点的狭窄率，选取最大龙骨节点的狭窄率为冠状动脉的狭窄率。

5.如权利要求1所述系统，其特征在于，所述模型确定模块，具体还用于若冠状动脉的狭窄率大于第一预设狭窄率，确定评估FFR的模型为第二模型。

6.如权利要求1所述系统，其特征在于，所述模型确定模块，具体还用于若冠状动脉的狭窄率不小于第二预设狭窄率，确定评估FFR的模型为第二模型；若冠状动脉的狭窄率大于第一预设狭窄率且小于第二预设狭窄率，确定评估FFR的模型为第三模型。

7.如权利要求5或6所述系统，其特征在于，所述评估FFR的模型为第一模型，所述FFR计算模块包括入口信息获取单元、压力分布获取单元以及第一计算单元，其中，

所述入口信息获取单元，用于根据所述冠状动脉的三维图像信息，获取各分支入口的信息；

所述压力分布获取单元，用于依据所述依据龙骨节点信息参数以及各分支入口的信息，获取龙骨节点的压力分布信息；

所述第一计算单元，用于根据龙骨节点的压力分布信息，计算冠状动脉的FFR值。

8.如权利要求5或6所述系统，其特征在于，所述评估FFR的模型为第二模型，所述FFR计算模块包括流量信息获取单元、阻抗确定单元、压力确定单元、第二计算单元：

所述流量信息获取单元，用于根据从各出口龙骨节点的压力以及各龙骨节点信息参数，计算出相应分叉点龙骨节点的压力相等的关系，以及依据质量守恒定律，确定出所述冠状动脉入口龙骨节点流量与各个出口龙骨节点流量的关系，获取各个出口龙骨节点流量信息；

所述各个出口龙骨节点流量信息的表达式为

A_iQ_i ²+B_iQ_i+C_i＝0

其中，A_i、B_i、C_i都是与龙骨节点信息参数有关的常数，Q_i为第i个出口龙骨节点流量信息；

所述阻抗确定单元，用于根据所述冠状动脉的三维图像信息，获取冠状动脉入口信息以及各分支信息，确定各出口阻抗；

所述压力确定单元，用于依据各个出口龙骨节点流量信息和各出口阻抗，确定各个出口龙骨节点压力信息；

所述第二计算单元，用于根据各个出口龙骨节点压力信息和冠状动脉所有龙骨节点的信息参数，计算冠状动脉的FFR值。

9.如权利要求6所述系统，其特征在于，所述评估FFR的模型为第三模型，所述FFR计算模块包括第一FFR计算单元、第二FFR计算单元以及第三FFR计算单元，其中，

所述第一FFR计算单元，用于依据龙骨节点信息参数以及第一模型，计算出待检测冠状动脉的FFR值FFR₁；

所述第二FFR计算单元，用于依据龙骨节点信息参数以及第二模型，计算出待检测冠状动脉的FFR值FFR₂；

所述第三FFR计算单元，用于计算FFR₁和FFR₂加权平均值，确定依据所述第三模型计算冠状动脉的FFR值FFR₃；

所述FFR₁和FFR₂加权平均的计算公式为：

FFR₁(1-ε)+FFR₂ε＝FFR₃

其中，ε＝(θ-θ₁)/(θ₂-θ₁)，θ为冠状动脉的实际狭窄率，θ₁为第一预设狭窄率，θ₂为第二预设狭窄率。

10.一种评估FFR的设备，其特征在于，所述设备包括权利要求1至9中任一项所述的评估FFR的系统以及显示装置，其中，所述显示装置用于显示FFR信息。

11.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有评估FFR的程序，所述评估FFR的程序被处理器执行时实现如下步骤：

获取冠状动脉的三维图像信息，依据所述冠状动脉的三维图像信息，确定冠状动脉的龙骨节点的分布信息，获取所述龙骨节点的信息参数；其中，所述龙骨节点为冠状动脉血管的中心线上离散的空间分布点；

根据所述冠状动脉的龙骨节点的信息参数，确定冠状动脉的狭窄率；

根据冠状动脉的狭窄率与预设狭窄率的关系，确定评估FFR的模型；

若冠状动脉的狭窄率不大于第一预设狭窄率，确定评估FFR的模型为第一模型；

依据龙骨节点信息参数以及评估FFR模型，计算冠状动脉的FFR值。

Images