(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210491604.5 (22)申请日 2022.04.29 (71)申请人 西安理工大 学 地址 710048 陕西省西安市碑林区金花 南 路5号 (72)发明人 谢国　叶欣　穆凌霞　李艳恺　 王承兰　王博　李佳婕　金永泽　 上官安琪　 (74)专利代理 机构 西安弘理专利事务所 61214 专利代理师 王敏强 (51)Int.Cl. G16H 50/50(2018.01) G16H 30/40(2018.01) G06F 30/23(2020.01) G06F 30/28(2020.01)G06V 10/44(2022.01) G06V 10/82(2022.01) G06N 3/04(2006.01) G06N 3/08(2006.01) G06T 5/00(2006.01) G06T 7/00(2017.01) G06T 7/11(2017.01) G06T 7/136(2017.01) G06T 17/20(2006.01) G06F 111/10(2020.01) G06F 113/08(2020.01) G06F 119/14(2020.01) (54)发明名称 一种结合多运动强度下的狭窄冠脉血运变 化的分析方法 (57)摘要 本发明公开了一种结合多运动强度下的狭 窄冠脉血运变化的分析方法， 首先对冠脉狭 窄CT 图构建图像样本库， 进行标签类别处理； 对由心 脏CTA提取出的CT原始序列图采用小波相关方法 进行去噪、 增强预处理； 对CT图像采用深度学习 方法中卷积神经网络模型CNN进行特征提取， 并 采用Otsu方法进行分割， 最后进行三维重建； 将 三维重建后的图像进行有限元分析； 最后将测试 患者在静息状态、 不同运动强度的情况下， 病人 的壁面剪切力、 血流速度、 壁面压力以及壁面形 变变化情况。 本发明测量冠脉粥样硬化病人在静 息及三个不同等级的运动强度下的血流动力学 分析， 实现在控制病情上的准确判断， 免去有创 手术给患者带来的伤害。 权利要求书4页 说明书9页 附图3页 CN 114864095 A 2022.08.05 CN 114864095 A 1.一种结合多运动强度下的狭窄冠脉血运变化的分析方法， 其特征在于， 具体按照以 下步骤实施： 步骤1、 对冠 脉狭窄CT图构建图像样本库， 进行 标签类别处理； 步骤2、 对步骤1中由心脏CTA提取出的CT原始序列图采用小波相关方法进行去噪、 增强 预处理； 步骤3、 对CT图像采用深度学习方法中卷积神经网络模型CNN进行特征提取， 并采用 Otsu方法进行分割， 最后进行三维重建； 步骤4、 按照7： 3的原则， 采用随机索引法， 将入组的M例样本分为训练组和测试组， 将步 骤3中三维重建后的图像进行有限元分析； 步骤5、 将步骤4中筛选后的病人统一集中， 保证一切正常情况下， 测试患者在静息状 态、 不同运动强度的情况下， 病 人的壁面剪切力、 血流速度、 壁面压力以及壁面形变变化情 况。 2.根据权利要求1所述的一种结合多运动强度下的狭窄冠脉血运变化的分析方法， 其 特征在于， 所述 步骤1具体按照以下步骤实施： 在医院数据系 统中收集近一年来同时做过心脏CTA和冠脉造影DSA检查的患者信息和 图像， 即能够将 CT图像和冠脉狭 窄金指标数据进 行对应， 隐去图像中患者的基本信息之后， 进行标签类别标注， 筛选出符合影像质量、 年龄在20～50岁、 且狭窄范围在25％～49％的 400例患者， 分别进 行静息以及不同强度运动等级下的冠脉CT图像作为选定的输入样本， 其 中， 运动强度等级依据最大耗氧量VO2max分为 五级， 一级： 最大摄氧量小于45％； 二级： 最大 摄氧量为55％～65％； 三级： 最大摄氧量为75％～8 5％； 四级： 最大摄氧量为90％～95％； 五 级： 最大摄氧量大于95％。 3.根据权利要求2所述的一种结合多运动强度下的狭窄冠脉血运变化的分析方法， 其 特征在于， 所述 步骤2具体按照以下步骤实施： 步骤2.1、 首先， 进行图像的去噪处理， 将步骤1中作为选定的输入样本作归一化处理， 图片尺寸统一成像素为512*512， 将样本中的原始信号和含噪信号首先采用Mallat算法进 行小波分解， 分解为多层的低频信号小波系 数和多层的细节信号小波系 数， 然后采用小波 局部阈值的方法进行去噪处 理， 最后小 波重构， 得到去噪信号的样本； 其中， 采用小 波局部阈值的方法具体如下： 阈值的选取如下： 在对角方向上 取： 其中， λ1为对角方向的阈值， σ 为噪声标准方差， N 为信号的尺度或长度； 在水平H、 垂直V方向上 取： 其中， λ2为水平H、 垂直V方向上的阈值， σ 为噪声标准方差， N 为信号的尺度或长度； σ 估计为： 权　利　要　求　书 1/4 页 2 CN 114864095 A 2其中， M表示小 波变换的第1层 细节信号小 波系数中值的绝对值； 阈值函数的选取： 阈值函数是修正小波系数的规则， 阈值 函数为： 其中， wj、k是二维小 波系数， 阈值 λ取值 为 λ1、 λ2； 步骤2.2、 对去噪后的图像作增强处理， 采用小波图像子带增强算法， 第一步进行2层小 波分解， 分别提取低频子图像以及低频子图像的小波系数； 第二步确定高频小波阈值； 第三 步计算各层的增益系 数； 第四步去 噪并增强按照图像增强的要求分别对低频、 高频子图像 的小波系数乘以不同的权数对小 波系数进行处 理。 4.根据权利要求3所述的一种结合多运动强度下的狭窄冠脉血运变化的分析方法， 其 特征在于， 所述步骤2.2中小波图像子带增强算法 中， 低频子图像的增强权值为 高频子 图像采用公式进行处 理： 其中， 选取图像增强的阈值l＝0.1， 垂直边缘子图HL、 水平边缘子图LH的增强系数k＝ 20， 对角线子图HH的增强系数k＝23， 最后对增强后的小波系数进行小波逆变换得到增强后 的图像。 5.根据权利要求4所述的一种结合多运动强度下的狭窄冠脉血运变化的分析方法， 其 特征在于， 所述 步骤3具体如下： 步骤3.1、 对步骤2中经过去噪和增强处理后得到的CT图像， 通过卷积－池化层提冠状 动脉的特征， 并生 成冠状动脉的特征图； 然后在生 成特征图后使用候选区域网络RPN提取疑 似冠状动脉 的候选区域， 图像的特征提取和候选框提取是共享卷积的， 图像的特征提取网 络使用的是5 0层卷积的残差网络， 提取 出冠状动脉 灰度区域图； 步骤3.2、 将步骤3.1中提取出的冠状动脉灰度区域图， 采用改进的Otsu算法对冠状动 脉灰度区域图进行分割提取， 选定阈值T， 遍历图像的像素， 将灰度值小于T的像素分为背 景， 将灰度值大于T的像素分为 目标区域， 对待分割 灰度图像作归一化处理， 阈值选择准则 函数改进为： 其中， σω(t)代表带指数参数的类内方差， μB(t)和 μo(t)的代表背景和目标的类内灰度 均值， pi表示不同灰度级出现的概 率， t表示时间， i表示灰度值的不同取值； 分割后的图像进行连通 区域的计算； 对落入各连通 区域中的被认为冠状动脉的像素点 进行计数； 选取计数值最大的判断为冠状动脉所在的区域， 计数值小的区域则去除， 从而完 成冠状动脉区域的分割提取；权　利　要　求　书 2/4 页 3 CN 114864095 A 3