【技术实现步骤摘要】
本专利技术涉及一种基于狭窄冠状动脉和修复冠状动脉之间的流量比计算ct-ffr的方法(ct-qffr)。
技术介绍
1、心外膜冠状动脉狭窄可能导致心肌缺血，而由此产生的阻塞性冠状动脉疾病是导致死亡的主要原因之一。血流储备分数(ffr)被定义为在冠状动脉存在狭窄病变的情况下，该血管所供心肌区域能获得的最大血流与同一区域理论上正常情况下所能获得的最大血流之比。然而，现有情况无法获得“同一区域理论上正常情况下所能获得的最大血流”，因此ffr被简化为狭窄远端平均压力与狭窄近端平均压力的比值。这样的简化，忽略了无狭窄时冠脉固有的压力损失，因此测量结果受所选择的测量位置影响，且被固有的压力损失干扰。
技术实现思路
1、本专利技术的目的是：基于ffr的定义，修复冠状动脉狭窄，计算狭窄冠状动脉和狭窄修复后冠状动脉之间的流量比，获得不依赖与测量位置的血流储备分数指标，即ct-qffr。
2、为了达到上述目的，本专利技术的技术方案是提供了一种基于狭窄冠状动脉和修复冠状动脉之间流量比计算ct-ffr的方法，其特征在于，包括以下步骤：
3、步骤1、基于冠脉cta图像，重建心外膜冠状动脉三维模型，其中，每个冠状动脉母血管在分叉过程中，分叉成两个子血管，分叉前冠状动脉母血管的直径和分叉后子血管的直径满足以下关系：
4、(d0)γ＝(d11)γ+(d12)γ＝2(d11)γ
5、其中，d0是分叉前冠状动脉母血管的直径，d11和d12是分叉后子血管的直径，γ是分叉幂律
6、步骤2、拟合分叉幂律指数γ和患者特异性心肌体积建立微循环模型并计算微循环阻力，设置为冠状动脉树出口阻力边界条件，其中，在建立微循环模型时，每条毛细血管的生长体积相同，第i个冠状动脉出口下游分配的生长区域的体积vi根据相应的截面等效直径di来分配，则有：
7、
8、其中，v代表所有冠状动脉生长区域的体积，n代表冠状动脉树几何模型上的冠状动脉出口总数；
9、在生长指数γ确定后，基于poiseuille定律求解各段血管的阻力，得到微循环阻力；
10、步骤3、患者肱动脉压设置为冠状动脉树入口压力边界条件；
11、步骤4、修复心外膜冠状动脉三维模型中狭窄的冠状动脉，进一步包括以下步骤：
12、步骤401、确定冠状动脉狭窄的位置，以及冠状动脉狭窄病灶的上游和下游非狭窄区域的横截面积；
13、步骤402、获得包含狭窄段的冠状动脉中心线，沿中心线内插值血管横截面以获得无狭窄的冠状动脉；
14、步骤403、将心外膜冠状动脉三维模型中的狭窄段替换为无狭窄的冠状动脉；
15、步骤5、利用cfd仿真原始的狭窄冠状动脉和通过步骤4修复后的非狭窄冠状动脉血流动力学环境；
16、步骤6、利用狭窄冠状动脉血管和非狭窄冠状动脉血管血流比例计算得到ct-qffr。
17、优选地，步骤2中拟合分叉幂律指数γ具体包括以下步骤：
18、步骤201、将生长指数γ设为初始值；
19、步骤202、在当前生长指数γ下计算冠状动脉各出口处的微血管阻力值；
20、步骤203、进行数值模拟，获得狭窄的冠状动脉的血流动力学参数；
21、步骤204、确定狭窄的冠状动脉下游ct-ifr与实际测量ifr之差：如果差值小于1％，则进入步骤206；如果差值大于1％，进入步骤205；
22、步骤205、根据牛顿迭代法确定新的生长指数γ，返回步骤202；
23、步骤206、计算多个病例的分叉动力法，并计算出所得到的所有生长指数γ的平均值，即为最终的生长指数γ。
24、优选地，步骤2中，基于poiseuille定律求解各段血管的阻力如下式所示：
25、
26、其中，r为血管阻力，l为血管长度，d为血管直径，μ为血液粘度。
27、优选地，所述血液粘度采用下式计算：
28、
29、
30、
31、其中，hd代表血细胞比容。
32、优选地，分叉前冠状动脉的生长区域相当于一个球体，在球面上取一点a，通过a点的平面将球体分为两个半球，两个半球的中心点分别为b和c，此时，线段ab和线段ac的长度分别为两个分叉的血管的长度，两个半球的体积是两个分叉后的血管生长区域的体积，相当于分叉后球体的体积。
33、本专利技术基于ffr的定义，即在冠状动脉存在狭窄病变的情况下，该血管所供心肌区域能获得的最大血流与同一区域理论上正常情况下所能获得的最大血流之比，提出计算狭窄冠状动脉和狭窄修复后冠状动脉之间的流量比，计算得到无创的基于流量比值的ct-ffr，即ct-qffr。本专利技术计算速度快，无创伤，腺苷非依赖，准确性高，可重复性高，诊断效能强，符合ffr原始定义，ct-qffr值不依赖测量点的选取。
本文档来自技高网...【技术保护点】
1.一种基于狭窄冠状动脉和修复冠状动脉之间流量比计算CT-FFR的方法，其特征在于，包括以下步骤：
2.如权利要求1所述的一种基于狭窄冠状动脉和修复冠状动脉之间流量比计算CT-FFR的方法，其特征在于，步骤2中拟合分叉幂律指数γ具体包括以下步骤：
3.如权利要求1所述的基于狭窄冠状动脉和修复冠状动脉之间流量比计算CT-FFR的方法，其特征在于，步骤2中，基于Poiseuille定律求解各段血管的阻力如下式所示：
4.如权利要求3所述的基于狭窄冠状动脉和修复冠状动脉之间流量比计算CT-FFR的方法，其特征在于，所述血液粘度采用下式计算：
5.如权利要求3所述的基于狭窄冠状动脉和修复冠状动脉之间流量比计算CT-FFR的方法，其特征在于，分叉前冠状动脉的生长区域相当于一个球体，在球面上取一点A，通过A点的平面将球体分为两个半球，两个半球的中心点分别为B和C，此时，线段AB和线段AC的长度分别为两个分叉的血管的长度，两个半球的体积是两个分叉后的血管生长区域的体积，相当于分叉后球体的体积。
【技术特征摘要】
1.一种基于狭窄冠状动脉和修复冠状动脉之间流量比计算ct-ffr的方法，其特征在于，包括以下步骤：
2.如权利要求1所述的一种基于狭窄冠状动脉和修复冠状动脉之间流量比计算ct-ffr的方法，其特征在于，步骤2中拟合分叉幂律指数γ具体包括以下步骤：
3.如权利要求1所述的基于狭窄冠状动脉和修复冠状动脉之间流量比计算ct-ffr的方法，其特征在于，步骤2中，基于poiseuille定律求解各段血管的阻力如下式所示：
4.如权利要求3所...
【专利技术属性】
技术研发人员：过伟锋，余龙，何玮，王盛章，曾蒙苏，沈雳，李晨光，
申请(专利权)人：复旦大学附属中山医院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。