一种无创计算个性化冠状动脉分支血流储备能力的方法技术

一种无创计算个性化冠状动脉分支血流储备能力的方法，属于生物医学工程领域。运用心脏超声数据得到冠状动脉静息状态的流量。基于冠状动脉CT图像文件重建出冠状动脉三维模型并测量冠状动脉分支长度、体积和开口横截面积；然后根据冠状动脉分形系数和异速生长规律，建立基于冠状动脉分支开口横截面积的流量分配方法，结合冠状动脉狭窄段和无狭窄段对应阻力数学模型获得冠状动脉远端出口压力，通过流体动力学计算获得考虑狭窄作用的冠状动脉远端出口流量；最后依据冠状动脉最大充血状态流量是静息状态3倍的关系，获得最大充血状态冠状动脉远端出口流量，通过冠状动脉最大充血和静息状态血流量比值与阈值3比较，确定冠状动脉分支血流储备能力。

A noninvasive method for calculating individual coronary branch flow reserve capacity

【技术实现步骤摘要】
一种无创计算个性化冠状动脉分支血流储备能力的方法
本专利技术提出的一种无创计算个性化冠状动脉分支血流储备能力的方法，属于生物医学工程
，涉及一种从静息到最大充血状态无创计算冠状动脉分支血流储备能力的方法。
技术介绍
随着超声影像技术的发展，通过在静息和运动试验中心脏室壁运动情况对比，定性描述从静息到最大充血状态冠状动脉分支血流储备能力。常用的模拟冠状动脉最大充血状态方法包括运动负荷和药物负荷。(1)运动负荷是通过踏车或者活动平板等以增加心肌耗氧的方式来模拟冠状动脉最大充血状态，通过在静息与最大充血状态下心脏室壁运动情况对比，定性描述从静息到最大充血状态冠状动脉分支血流储备能力。虽然运动负荷可以模拟冠状动脉最大充血状态，但其往往需要通过运动的方式使冠状动脉达到最大充血状态，这对于不能运动的来说，运动负荷模拟冠状动脉最大充血状态具有一定的局限性。(2)药物负荷是通过注射血管扩张药等以扩张冠状动脉血管的方式来模拟冠状动脉最大充血状态。虽然对于不能运动的可通过药物负荷来模拟冠状动脉最大充血状态，但是其不能定量研究从静息到最大充血状态冠状动脉分支血流储备能力。虽然超声影像技术可无创评价冠状动脉分支血流储备能力，但是其不能定量研究从静息到最大充血状态冠状动脉分支血流储备能力。本专利技术从流体动力学方面研究冠状动脉血流情况，定量研究从静息到最大充血状态冠状动脉分支血流储备能力。
技术实现思路
本专利技术已知冠状动脉静息状态的流量，基于冠状动脉分支开口横截面积的流量分配方法，将计算得到的冠状动脉分支流量与阻力作乘积以便确定冠状动脉分支对于流量产生的压力差，并将其施加到边界条件中，通过对冠状动脉模型进行流体动力学计算，得到冠状动脉分支出口流量。通过调整冠状动脉入口流量来模拟最大充血状态，并分别比较两种状态冠状动脉分支出口流量。通过冠状动脉分支最大充血与静息状态流量的比值，从而完成冠状动脉分支由静息到最大充血状态血流储备能力的无创计算。本专利技术已知冠状动脉静息状态的流量，基于冠状动脉体积和长度参数，运用莫里原则和最小耗能定理，得到冠状动脉模型分形系数。在分形系数基础上，依据异速生长规律和冠状动脉左前降支LAD、左回旋支LCX、右冠状动脉RCA的流量，基于冠状动脉分支开口横截面积的流量分配方法，并确定冠状动脉分支流量，基本步骤包括：(1)模型的获取：基于冠状动脉CT图像的DICOM格式原始数据，重建冠状动脉模型。(2)分形系数的确定：依据冠状动脉长度和体积，可得到归一化冠状动脉长度和体积的数据点，在莫里原则和最小耗能定律基础上，以归一化冠状动脉长度为横坐标，归一化冠状动脉体积为纵坐标，通过拟合归一化长度和体积的指数关系曲线，得到关于分形系数的拟合方程并求得分形系数。(3)冠状动脉LAD、LCX、RCA流量的确定：已知冠状动脉的流量，依据冠状动脉长度和体积的指数关系及流经LAD、LCX、RCA的流量分别与它们的分支长度成正比，以冠状动脉分支长度为中间量，建立冠状动脉体积与LAD、LCX、RCA流量的关系，确定冠状动脉LAD、LCX、RCA的流量。(4)冠状动脉分支流量分配方法的确定：依据冠状动脉分形系数和异速生长规律，确定冠状动脉分支开口横截面积与流量的指数关系。(5)冠状动脉分支流量的确定：根据冠状动脉LAD、LCX、RCA的流量及基于冠状动脉分支开口横截面积的流量分配方法，确定冠状动脉分支流量。本专利技术将冠状动脉分支对于流量产生的阻力以压力差的方式加载到边界条件中，通过调整入口流量来模拟最大充血状态，并分别比较两种状态冠状动脉分支流量。通过最大充血与静息状态冠状动脉分支流量的比值，从而完成冠状动脉分支从静息到最大充血状态血流储备能力的无创计算，基本步骤包括：(1)基于泊肃叶定律、冠状动脉最狭窄处横截面积、狭窄长度及分支流量，确定冠状动脉分支对应阻力。(2)已知冠状动脉入口压力，通过冠状动脉分支流量与阻力的乘积，得到各分支对应压力差，沿着血流方向从冠状动脉近端到远端逐个分支计算并确定各分支出口压力。(3)以静息状态冠状动脉LAD、LCX、RCA流量为入口条件，通过逐个分支计算的出口压力为出口条件，将静息状态下冠状动脉LAD、LCX、RCA流量的3倍作为最大充血状态入口条件并计算出相应出口压力，通过对两种状态进行流体动力学计算，得到静息和最大充血状态对应出口流量分别为Ql和Qj(l或j取值为1-13)。(4)若冠状动脉分支没有狭窄，则在静息与最大充血状态下冠状动脉分支血流储备为Qj/Ql或3(l或j取值为1-13)；若冠状动脉分支有狭窄时，运用冠状动脉分支最大充血到静息状态出口流量的比值来计算冠状动脉分支血流储备能力。若比值大于或等于阈值3时，说明冠状动脉分支血流储备能力较强，若比值小于阈值3时，说明冠状动脉分支血流储备能力较弱，从而完成冠状动脉分支从静息到最大充血状态血流储备能力的无创计算。附图说明图1主动脉和冠状动脉模型俯视图，其中2-7代表冠状动脉LAD各分支，8-10代表冠状动脉LCX各分支，11-23代表冠状动脉RCA各分支，圆圈圈出的是冠状动脉狭窄部分，其位于冠状动脉LAD上。图2主动脉和冠状动脉模型主视图，其中outlet1-4代表冠状动脉LAD各分支出口，outlet5-6代表冠状动脉LCX各分支出口,outlet7-13代表冠状动脉RCA各分支出口。图3冠状动脉LAD模型，其中inlet为左冠状动脉入口，outlet1-4为冠状动脉LAD各分支出口。图4冠状动脉LAD狭窄部位局部放大图，outlet1-3为冠状动脉LAD狭窄上游分支，outlet4为冠状动脉LAD狭窄下游分支。图5从静息到最大充血状态冠状动脉分支血流储备能力的无创计算流程图具体实施方式下面参照本专利技术的流程示意图，对本专利技术的一个实施例进行详细阐述，以使本专利技术的优点和特征能更易于被本领域技术人员理解，从而对本专利技术的保护范围做出更为清楚明确的界定。实施例11.基于冠状动脉分支开口横截面积的流量分配方法。本专利技术利用冠状动脉分形系数和异速生长规律，得到冠状动脉分支开口横截面积与流量的指数关系式，沿着血流方向从冠状动脉近端到远端逐个分支计算得到冠状动脉分支流量占上一级流量的比例，在已知冠状动脉LAD、LCX、RCA的流量基础上，确定冠状动脉各分支流量，基本步骤包括：(1)模型获取：基于冠状动脉CT图像的DICOM格式原始数据，重建冠状动脉模型。(2)分形系数的确定：依据冠状动脉体积和长度，可得到归一化冠状动脉长度和体积的数据点，在莫里原则和最小耗能定律基础上，以归一化冠状动脉长度为横坐标，归一化冠状动脉体积为纵坐标，通过拟合归一化长度和体积的指数关系曲线，得到关于分形系数的拟合方程并求得分形系数。(3)冠状动脉分支流量分配方法的确定：依据冠状动脉分形系数和异速生长规律，确定冠状动脉分支开口横截面积与流量的关系，沿着血流方向从冠状动脉近端到远端逐个分支计算得到冠状动脉分支流量占上一级流量的比例，最后通过冠状动脉LAD、LCX、RCA的流量，确本文档来自技高网...

【技术保护点】
1.一种无创计算个性化冠状动脉分支血流储备能力的方法，其特征在于，包括以下步骤：/n步骤一，确定冠状动脉左前降支LAD、左回旋支LCX、右冠状动脉RCA的分形系数ε

【技术特征摘要】
1.一种无创计算个性化冠状动脉分支血流储备能力的方法，其特征在于，包括以下步骤：
步骤一，确定冠状动脉左前降支LAD、左回旋支LCX、右冠状动脉RCA的分形系数ε1，ε2，ε3；
依据CT图像重建冠状动脉模型，确定分支体积和长度；根据莫里原则和最小耗能定律，得到归一化长度与体积之间的指数关系，求得分形系数εi，这里i＝1,2,3；
步骤二，确定冠状动脉LAD、LCX、RCA流量：
运用心脏超声数据得到冠状动脉静息状态的流量，且流经LAD、LCX、RCA的流量分别与它们的分支长度成正比，根据步骤一中冠状动脉长度与体积的关系，以冠状动脉长度为中间量，建立冠状动脉LAD、LCX、RCA流量与体积之间的指数关系，确定LAD、LCX、RCA的流量；
步骤三，基于冠状动脉分支开口横截面积的流量分配方法：
在步骤二中得到的LAD、LCX、RCA流量基础上，依据冠状动脉分形系数和异速生长规律，得到冠状动脉分支开口横截面积与流量的关系，确定冠状动脉分支流量；
步骤四，确定冠状动脉分支对应阻力：
依据冠状动脉管径最小横截面积与参考横截面积比值来判断冠状动脉是否狭窄，当其比值小于1时，判定冠状动脉发生狭窄；若冠状动脉有狭窄段，则依据最狭窄处横截面积、狭窄长度、参考横截面积和分支流量，确定冠状动脉分支对应阻力；
若冠状动脉无狭窄段，则依据冠状动脉分支开口和出口横截面积的平均值和分支长度，并根据泊肃叶定律，确定冠状动脉分支对应阻力；
步骤五，确定冠状动脉分支出口压力并完成冠状动脉分支从静息到最大充血状态血流储备的无创计算：
通过冠状动脉分支流量与阻力的乘积得到冠状动脉分支对于流量产生的压力差；已知冠状动脉入口压力，将沿着血流方向从冠状动脉近端到远端逐个分支计算得到的静息状态下冠状动脉分支出口压力作为出口边界条件，通过步骤二确定的静息状态下LAD、LCX、RCA流量为入口条件，通过将其施加到冠状动脉模型的流体动力学计算中，得到冠状动脉分支出口流量Ql，这里l取值为1-13；
对于冠状动脉最大充血状态的模拟，将静息状态下LAD、LCX、RCA流量的3倍作为各自的入口条件，通过步骤三和四的计算得出最大充血状态下冠状动脉分支流量和阻力，并将分支流量和阻力做乘积得到冠状动脉分支对于流量产生的压力差，已知冠状动脉入口压力，将沿着血流方向从冠状动脉近端到远端逐个分支计算得到最大充血状态下冠状动脉分支出口压力作为边界条件进行流体动力学计算，随后确定冠状动脉分支出口流量Qj，并与Ql比较，这里l或j取值为1-13；
取冠状动脉分支最大充血与静息状态出口流量的比值来计算冠状动脉分支血流储备能力。


2.按照权利要求1所述的一种无创计算个性化冠状动脉分支血流储备能力的方法，其特征在于，依据莫里原则和最小耗能定律，确定冠状动脉分形系数：
依据CT图像获取冠状动脉分支长度和体积，并分别对其进行归一化处理：
冠状动脉体积的归一化：V/Vsum
冠状动脉长度的归一化：L/Lsum
其中V为冠状动脉分支体积，Vsum为冠状动脉分支及下游分支的体积之和，L为冠状动脉分支长度，Lsum为冠状动脉分支及下游分支的长度之和；
根据莫里原则和最小耗能定律，以L/Lsum为横坐标，V/Vsum为纵坐标，通过拟合归一化长度和体积的指数关系曲线，得到关于εi的拟合方程并求得分形系数εi：



其中i取1,2,3,分别对应LAD、LCX、RCA。


3.按照权利要求1所述的一种无创计算个性化冠状动脉分支血流储备能力的方法，其特征在于，依据归一化冠状动脉长度和体积的指数关系，确定冠状动脉LAD、LCX、RCA流量：
因流经LAD、LCX、RCA的流量分别与它们的分支长度成正比：
Q0＝KLh(2)
其中Q0为冠状动脉流量，K为比例常数，由冠状动脉入口和末端压力差决定，h取LAD，LCX，RCA；
依据归一化冠状动脉长度和体积的指数关系，以冠状动脉长度为中间量，确定冠状动脉LAD、LCX、RCA的流量：



其中VLAD、VLCX和VRCA分别为LAD、LCX和RCA体积，Vsum为冠状动脉总体积。


4.按照权利要求1所述的一种无创计算个性...

【专利技术属性】
技术研发人员：乔爱科，张宏辉，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11