一种基于深度学习预测冠状动脉狭窄阻力的系统与方法技术方案

一种基于深度学习预测冠状动脉狭窄阻力的系统与方法。属于人工智能领域，包括以下步骤：基于真实冠状动脉参数构建血管模型并记录血管模型参数(入口面积、狭窄长度等)；基于血管模型进行模型网格化预处理；基于几何多尺度血流动力学计算预处理模型并提取血管模型的狭窄阻力；基于数据提取与血管模型参数建立狭窄阻力训练集与预测集；基于BP神经网络建立神经网络框架；基于深度学习对狭窄阻力训练集进行训练与对预测集进行预测验证。

【技术实现步骤摘要】
一种基于深度学习预测冠状动脉狭窄阻力的系统与方法
本专利技术属于人工智能领域，是一种基于深度学习预测冠状动脉狭窄阻力的系统与方法。
技术介绍
FFR是一项直接描述功能性心肌缺血程度的可靠指标，可在相关研究中通过数值模拟获得，而冠状动脉狭窄阻力的确定对于数值计算FFR尤为重要。本项冠状动脉狭窄阻力确定方法拟为无创获取FFR(被定义为狭窄远端平均压力(Pd)与主动脉根部平均压力(Pa)的比值)值精确度提供计算依据。机器学习是一门多领域交叉学科，而深度学习是机器学习领域中一个新的研究方向，是学习样本数据的内在规律和表示层次，最终目的是让机器能够像人一样具有分析学习能力，能够识别文字，图像，形状等数据。本方法结合深度学习与数值模拟提出了一种基于深度学习预测冠状动脉狭窄阻力的算法实现了基于冠状动脉几何参数确定狭窄阻力，同时提高无创计算FFR数值模拟结果的准确性。对于个性化冠状动脉生理疾病的诊断、治疗策略研究具有一定理论价值。
技术实现思路
本专利技术提供一种基于深度学习预测冠状动脉狭窄阻力的算法，相较于数值模拟更加快速准确，普遍适用于确定血管狭窄产生的狭窄阻力计算。该项基于深度学习预测冠状动脉狭窄阻力的算法包括以下步骤：包括血管模型的构建、模型网格预处理、目标参数提取、建立数据库、建立神经网络、预测验证。为达到上述目的，本专利技术通过以下技术方案来实现：一种基于深度学习预测冠状动脉狭窄阻力的系统，其特征在于，是通过以下方法预测冠状动脉狭窄阻力的，包括以下步骤：步骤A1：基于真实冠脉解剖结构参数构建冠脉血管模型并记录冠脉血管模型参数(例：入口面积、狭窄面积、狭窄长度、狭窄面积等)；步骤A2：基于血管模型进行模型网格化预处理；步骤A3：基于几何多尺度血流动力学计算预处理模型并提取血管模型狭窄阻力；步骤A4:基于数据提取与血管模型参数建立狭窄阻力训练集与预测集；步骤A5：基于BP神经网络建立神经网络框架；步骤A6：基于深度学习对狭窄阻力训练集进行训练与对预测集进行预测验证。1.作为本专利技术的技术方案，步骤A1中所述特征，首先基于真实冠脉解剖结构参数构建冠脉血管模型,通过SoldWork软件构造冠状动脉血管模型。构建模型要求如下：入口面积、血管长度、狭窄率等形状参数应符合真实血管的参数。例如：血管直径D(2-6mm)、狭窄长度ls＝(2-40mm)、狭窄率DS％(20-90％)、狭窄直径Ds＝D*(1-DS％)、血管长度L＝ls+D*8+D*12)；作为本专利技术进一步的技术方案，步骤A2中所述特征，基于血管模型进行模型网格化预处理，计算软件支持网格文件的输入形式，因此对模型网格化前执行如下预处理步骤：步骤B1：通过SoldWorks软件对构建血管模型导出x_t格式；步骤B2：通过Ansysworkbench15.0网格划分子功能模块对构建的血管模型进行算前网格划分，该步骤要求对网格做敏感性分析，排除网格因子对计算结果的影响。作为本专利技术进一步的技术方案，步骤A3中所述特征，基于多尺度血流动力学计算预处理模型提取血管模型狭窄阻力。具体实施过程包括以下步骤：步骤C1：给予同静息流量确定对应零维后负荷模型，零维冠脉后负荷模型由基本电子元器件电感和电阻构成，血流动力学参数与物理电学参数等效关系如表1。表1血流动力学参数与物理电学参数等效关系可以假设一段血管长度为l，血管横截面积为A，一段血管的零维及三维对应参数关系见图5，则根据Frank建立的弹性腔理论，该段血管上的流量、压力及相关阻力之间的定量关系可以表示为：Δp＝Q*R′#(1)公式中ΔP表示在该段血管中产生的压降，Q表示该支血管的血流量，R’表示该段血管产生的粘性阻力；μ为正常生理下血液粘度值，其值为3.5e-3Pa·s；l表示该段血管长度，A表示该段血管横截面积。以电感的自感应性等效于血流的流动惯性，电阻的阻力效应等效于冠脉各分支微循环阻力。那么零维模型中电阻R及电感L值为：公式中P末表示该段血管末端压力，ρ表示血流密度值，定义其为1060kg/m3。步骤C2：将网格化模型导入Ansys-CFX流体计算软件中，赋予流体计算参数血流密度和血流动力粘度值(即正常生理下血液粘度值)；；步骤C3：对于边界条件的确定与三维模型的压力-流量-后负荷关系计算结果为依据，在血管入口给定平均动脉压，出口定义虚拟流量边界，并耦合零维后负荷模型步骤C4：数值模拟计算中，设置时间为1.2s，设置时间步长为0.01s，执行仿真计算过程。步骤C5：计算完成后在狭窄后30mm处提取血管压力，通过与血管入口压力的比较得出预处理模型的压差，通过压差除以流量得到狭窄阻力。并提取狭窄入口的血流量。作为本专利技术进一步的技术方案，步骤A4中所述特征，基于数据提取与血管模型参数建立狭窄阻力训练集与预测集；通过步骤A1、A2、A3对340根不同狭窄参数血管进行计算，提取的血管模型的狭窄阻力和血流量与血管模型的几何参数构建深度学习的数据集，将数据集分为训练集和测试集。作为本专利技术进一步的技术方案，步骤A5中所述特征，基于BP神经网络建立神经网络框架，具体实施包括以下步骤：步骤D1：建立深度学习隐藏层，设定隐藏层数、隐藏层节点数和激活函数等；步骤D2：选择合适的优化器、Loss函数和其他超参数。作为本专利技术进一步的技术方案，步骤A5中所述特征，基于深度学习对狭窄阻力训练集进行训练与对预测集进行预测验证。具体实施过程包括以下步骤：步骤E1：通过pandas模块将训练集和测试集导入程序中；步骤E2：通过相关系数矩阵分析血管几何参数对血管模型的狭窄阻力的影响，并分析避免变量之间的多重共线性，最终得到七个相关性较强的变量分别为入口面积、狭窄面积、狭窄入口长度、狭窄出口长度、最小狭窄长度、狭窄率、入口血流量。步骤E3：将处理后的数据导入神经网络框架中，进行超参数的调节。步骤E4：保存训练模型，执行预测血管模型的狭窄阻力。采用本专利技术的系统预测的冠状动脉狭窄阻力比较准确，误差小。附图说明：图1:本专利技术方法的流程图图2：血管模型图图3：零维三维耦合模型图图4：参数矩阵相关性分析图5：神经网络结构图6：预测数值与计算数值比对图具体实施方式下面结合实施例随本专利技术组进一步说明，但本专利技术并不限于以下是实施例。实施例1一种基于深度学习预测冠状动脉狭窄阻力的系统，是通过以下方法预测冠状动脉狭窄阻力的，包括以下步骤：步骤A1：基于真实冠脉解剖结构参数构建冠脉血管模型并记录冠脉血管模型参数(入口面积、狭窄面积、狭窄入口长度、狭窄出口长度、最小狭窄长度、狭窄率)；步骤A2：基于血管模型进行模型网格化预处理；步骤A3：基于几何多尺度血流动力学计算预处理模型并提取血管模型狭窄阻力；本文档来自技高网...

【技术保护点】
1.一种基于深度学习预测冠状动脉狭窄阻力的系统，其特征在于，是通过以下方法预测冠状动脉狭窄阻力的，包括以下步骤：/n步骤A1：基于真实冠脉解剖结构参数构建冠脉血管模型并记录冠脉血管模型参数(例：入口面积、狭窄面积、狭窄长度、狭窄面积等)；/n步骤A2：基于血管模型进行模型网格化预处理；/n步骤A3：基于几何多尺度血流动力学计算预处理模型并提取血管模型狭窄阻力；/n步骤A4：基于数据提取与血管模型参数建立狭窄阻力训练集与预测集；/n步骤A5：基于BP神经网络建立神经网络框架；/n步骤A6：基于深度学习对狭窄阻力训练集进行训练与对预测集进行预测验证。/n

【技术特征摘要】
1.一种基于深度学习预测冠状动脉狭窄阻力的系统，其特征在于，是通过以下方法预测冠状动脉狭窄阻力的，包括以下步骤：
步骤A1：基于真实冠脉解剖结构参数构建冠脉血管模型并记录冠脉血管模型参数(例：入口面积、狭窄面积、狭窄长度、狭窄面积等)；
步骤A2：基于血管模型进行模型网格化预处理；
步骤A3：基于几何多尺度血流动力学计算预处理模型并提取血管模型狭窄阻力；
步骤A4：基于数据提取与血管模型参数建立狭窄阻力训练集与预测集；
步骤A5：基于BP神经网络建立神经网络框架；
步骤A6：基于深度学习对狭窄阻力训练集进行训练与对预测集进行预测验证。


2.按照权利要求1所述的一种基于深度学习预测冠状动脉狭窄阻力的系统，其特征在于，步骤A1中所述特征，首先基于真实冠脉解剖结构参数构建冠脉血管模型，通过SoldWork软件构造冠状动脉血管模型，构建模型要求等形状参数应符合真实血管的参数。


3.按照权利要求1所述的一种基于深度学习预测冠状动脉狭窄阻力的系统，其特征在于，步骤A2中所述特征，基于血管模型进行模型网格化预处理，计算软件支持网格文件的输入形式，因此对模型网格化前执行如下预处理步骤：
步骤B1：通过SoldWorks软件对构建血管模型导出x_t格式；
步骤B2：通过Ansysworkbench15.0网格划分子功能模块对构建的血管模型进行算前网格划分，该步骤要求对网格做敏感性分析，排除网格因子对计算结果的影响。


4.按照权利要求1所述的一种基于深度学习预测冠状动脉狭窄阻力的系统，其特征在于，步骤A3中所述特征，基于多尺度血流动力学计算预处理模型提取血管模型狭窄阻力；具体实施过程包括以下步骤：
步骤C1：给予同静息流量确定对应零维后负荷模型，零维冠脉后负荷模型由基本电子元器件电感和电阻构成；
假设一段血管长度为l，血管横截面积为A，则根据Frank建立的弹性腔理论，该段血管上的流量、压力及相关阻力之间的定量关系可以表示为：
Δp＝Q*R′#(1)



公式中ΔP表示在该段血管中产生的压降，Q表示该支血管的血流量，R’表示该段血管产生的粘性阻力；μ为正常生理下血液粘度值，其值为3.5e-3Pa·s；l表示该段血管长度，A表示该段血管横截面积；
以电感的自感应性等效于血流的流动惯性，电阻的阻力效应等效于冠脉各分支微循环阻力。那么零维模型中电阻R及电感L值为：






公式中P末表示该段血管末端压力，ρ表示血流密度值，定义其为1060kg/m3。
步骤C2：将网格化模型导入Ansys-CFX流体计算软件中，赋予流体计算参数血流密度和血流动力粘度值(即正常生理下血液粘度值)；
步骤C3：对于边界条件的确定与三维模型的压力-流量-后负荷关系计算结果为依据，在血管入口给定平均动脉压，出口定义虚拟流量边界，并耦合零维后负荷模型；
步骤C4：数值模拟计算中，设置时间为1.2s...

【专利技术属性】
技术研发人员：刘有军，孙昊，刘金城，冯懿俐，席晓璐，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11