基于连续T波面积与深度特征的心肌缺血模型建模方法技术

本发明专利技术公开了基于连续T波面积与深度特征的心肌缺血模型建模方法，具体为：对原始ECG数据进行降噪后截取多个ECG记录片段；确定每个ECG记录片段中除最后一个T波外其余T波的起止点；根据每个ECG记录片段的连续T波面积变异系数得到一个变异系数特征向量；使用8层卷积模块提取ECG的深度特征；将两个方面的特征融合为混合特征，得到的混合特征向量输入到多层感知机中；根据以上步骤构建心肌缺血检测模型，并进行训练与测试。本发明专利技术关注心肌缺血发生时多个T波形态的变化，提出使用变异系数定量评估连续T波面积的波动程度，并结合ECG的深度特征，能够有效提高识别心肌缺血的准确率、高效性和实用性。性和实用性。性和实用性。

【技术实现步骤摘要】
基于连续T波面积与深度特征的心肌缺血模型建模方法


[0001]本专利技术属于计算机信号处理领域，具体涉及一种基于连续T波面积特征与深度特征的心肌缺血模型建模方法。

技术介绍

[0002]近些年来心血管疾病的发病率在不断攀升，在心血管疾病患者中，心肌缺血是主要的死亡原因之一。心肌缺血是由于冠状动脉闭塞导致心肌供血不足甚至心肌坏死引起的心血管疾病，在该病早期，患者常出现胸闷、胸痛等症状，但部分患者仍无明显症状，难以及时治疗，严重的将危及生命。因此，如何实现心肌缺血的早期诊断具有重要的临床价值，也是一项十分有意义的研究课题。
[0003]相比于医学影像检查，患者对于无创性心电检测技术更加青睐，其价格较低、实用性强且受众广泛。标准12导联的心电图(ECG)是心肌缺血的无创性常规检查方法之一，这些不同导联会反映心脏不同部位的电位活动，通过心电信号的某些特征能够区分健康人和心肌缺血患者，因此良好的检测结果有助于判断患者是否存在心肌缺血病症。
[0004]目前随着人工智能的发展，机器学习的方法大量应用在心电信号上。在使用心电信号检测患者心肌缺血时，以往的方法主要是提取心电信号的频域、时域、小波变换等特征。近年来随着深度学习的发展，神经网络在分类效果上有了较大幅度的提升。一维卷积神经网络可以实现对心电图特征的自动提取，在不考虑算法复杂度的情况下，一维卷积神经网络表现优于传统机器学习算法。在上述两种主流人工智能检测心肌缺血的方法中，前者能够提取的特征有限，后者虽检测能力上稍优于前者，但无法将更多的目光聚焦于心肌缺血发生时T波的形态改变上。

技术实现思路

[0005]本专利技术的目的是为克服上现有技术存在的缺陷，提出一种基于连续T波面积特征与深度特征的心肌缺血模型建模方法，该方法采用T波检测算法检测T波，提取心电信号的连续T波面积特征(变异系数)，然后使用卷积神经网络提取心电信号的深度特征，最后融合为混合特征后通过多层感知机检测分类为正常ECG或心肌缺血ECG，该方法在心电信号深度特征的基础上更多关注心肌缺血发生时T波的形态变化，具有良好的检测准确率、高效性和实用性。
[0006]本专利技术采用的技术方案是：
[0007]本专利技术基于连续T波面积特征与深度特征的心肌缺血模型建模方法，包括以下步骤：
[0008]1)对原始ECG数据进行预处理，预处理包含降噪和对降噪后心电信号截取多个ECG记录片段。
[0009]2)采用基于区域滑动窗口的方法确定每个ECG记录片段中除最后一个T波外其余T波的起止点。
[0010]3)对每个ECG记录片段中每根导联除最后一个心拍外其余每个心拍计算T波面积，然后计算每个ECG记录片段中每根导联的连续T波面积变异系数，最后根据每个ECG记录片段的所有变异系数得到一个变异系数特征向量。
[0011]4)首先，构建心肌缺血检测模型，具体过程如下：通过卷积神经网络模型提取每个ECG记录片段的深度特征向量并处理为一维特征向量，与该ECG记录片段对应的变异系数特征向量融合为混合特征向量，再对混合特征向量线性变换后使用Sigmoid激活函数转换为概率值，最后使用多层感知机根据概率值将每个ECG记录片段分类为正常ECG片段或心肌缺血ECG片段。构建心肌缺血检测模型后，选取ECG记录片段中的部分心肌缺血ECG片段和正常ECG片段作为训练集，剩余心肌缺血ECG片段和正常ECG片段作为测试集，对心肌缺血检测模型进行训练与测试，且测试时根据各ECG记录片段的模型分类结果，结合各ECG记录片段的真实分类，计算评价指标，对心肌缺血检测模型进行评价，若评价指标满足要求，得到心肌缺血检测模型中卷积神经网络模型和多层感知机的最优参数，否则调整卷积神经网络模型和多层感知机的参数，重新对心肌缺血检测模型进行训练与测试。
[0012]优选地，对原始ECG数据进行预处理的过程如下：
[0013]a1)选取Daubechies6小波函数作为母小波来过滤原始ECG数据的信号噪声，得到降噪后心电信号；其中，信号噪声包含肌电干扰、基线漂移和工频噪声；
[0014]a2)对降噪后心电信号，截取多个ECG记录片段，包括心肌缺血ECG片段和正常ECG片段；其中，心肌缺血ECG片段和正常ECG片段的数目差别小于200。
[0015]优选地，确定ECG记录片段中除最后一个T波外其余T波起止点的过程具体如下：
[0016]b1)采用Pan
‑
Tompkins算法定位ECG记录片段中各T波前的R峰位置；
[0017]b2)根据QT数据库中标注的QRS波群索引位置以及T波和P波的起止点索引位置，提取QT数据库中多个心跳信号的RR间期、RTon间期、RTend间期，对RR和RTon、RR和RTend做k
‑
means算法聚类分析，从而确定在不同RR间期下T波起点和终点的搜索范围[t1,t2]和[t3,t4]。
[0018]b3)首先根据各相邻两个R峰位置计算得到的RR间期确定ECG记录片段中除最后一个T波外其余T波的起点搜索范围，然后采用滑动窗口定位ECG记录片段中除最后一个T波外其余T波的起点；定位T波起点的过程具体如下：将滑动窗口放在搜索左边界t1处，然后将滑动窗口向右移动，当滑动窗口的积分面积取到最大时，滑动窗口左端点t的取值就是T波起点的位置。
[0019]b4)首先根据各相邻两个R峰位置计算得到的RR间期确定ECG记录片段中除最后一个T波外其余T波的终点搜索范围，然后采用滑动窗口定位ECG记录片段中除最后一个T波外其余T波的终点位置，定位T波终点的过程具体如下：将滑动窗口定义为[t
′‑
w,t
′
]，将滑动窗口从t3处开始向t4处滑动，当滑动窗口积分面积最大时，t
′
的取值就是T波终点的位置。
[0020]更优选地，定位T波起点的滑动窗口积分面积公式如下：
[0021][0022]其中s
j
表示滑动窗口中第j个采样点的信号幅值，表示滑动窗口左端点t的邻域p内信号幅值的平均值，的公式为：
[0023][0024]更优选地，定位T波终点的滑动窗口积分面积公式为：
[0025][0026]其中表示滑动窗口右端点t
′
的邻域p内信号幅值的平均值。
[0027]优选地，步骤3)具体包含以下步骤：
[0028]c1)对每个ECG记录片段中每根导联除最后一个心拍外其余每个心拍计算T波面积，T波面积用T波的积分面积表达；T波积分面积表示为：
[0029][0030]其中s
t
代表T波中第t个采样点的幅值大小，T
on
与T
off
分别表示T波起点与终点，表示T
on
的邻域p内信号幅值的平均值；
[0031]c2)计算每个ECG记录片段中每根导联的连续T波面积变异系数；变异系数表示为c
v
＝σ/μ，σ表示每根导联中本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.基于连续T波面积与深度特征的心肌缺血模型建模方法，其特征在于：包括以下步骤：1)对原始ECG数据进行预处理，预处理包含降噪和对降噪后心电信号截取多个ECG记录片段；2)采用基于区域滑动窗口的方法确定每个ECG记录片段中除最后一个T波外其余T波的起止点；3)对每个ECG记录片段中每根导联除最后一个心拍外其余每个心拍计算T波面积，然后计算每个ECG记录片段中每根导联的连续T波面积变异系数，最后根据每个ECG记录片段的所有变异系数得到一个变异系数特征向量；4)首先，构建心肌缺血检测模型，具体过程如下：通过卷积神经网络模型提取每个ECG记录片段的深度特征向量并处理为一维特征向量，与该ECG记录片段对应的变异系数特征向量融合为混合特征向量，再对混合特征向量线性变换后使用Sigmoid激活函数转换为概率值，最后使用多层感知机根据概率值将每个ECG记录片段分类为正常ECG片段或心肌缺血ECG片段；构建心肌缺血检测模型后，选取ECG记录片段中的部分心肌缺血ECG片段和正常ECG片段作为训练集，剩余心肌缺血ECG片段和正常ECG片段作为测试集，对心肌缺血检测模型进行训练与测试，且测试时根据各ECG记录片段的模型分类结果，结合各ECG记录片段的真实分类，计算评价指标，对心肌缺血检测模型进行评价，若评价指标满足要求，得到心肌缺血检测模型中卷积神经网络模型和多层感知机的最优参数，否则调整卷积神经网络模型和多层感知机的参数，重新对心肌缺血检测模型进行训练与测试。2.根据权利要求1所述基于连续T波面积与深度特征的心肌缺血模型建模方法，其特征在于：对原始ECG数据进行预处理的过程如下：a1)选取Daubechies6小波函数作为母小波来过滤原始ECG数据的信号噪声，得到降噪后心电信号；其中，信号噪声包含肌电干扰、基线漂移和工频噪声；a2)对降噪后心电信号，截取多个ECG记录片段，包括心肌缺血ECG片段和正常ECG片段；其中，心肌缺血ECG片段和正常ECG片段的数目差别小于200。3.根据权利要求1所述基于连续T波面积与深度特征的心肌缺血模型建模方法，其特征在于：确定ECG记录片段中除最后一个T波外其余T波起止点的过程具体如下：b1)采用Pan
‑
Tompkins算法定位ECG记录片段中各T波前的R峰位置；b2)根据QT数据库中标注的QRS波群索引位置以及T波和P波的起止点索引位置，提取QT数据库中多个心跳信号的RR间期、RTon间期、RTend间期，对RR和RTon、RR和RTend做k
‑
means算法聚类分析，从而确定在不同RR间期下T波起点和终点的搜索范围[t1,t2]和[t3,t4]；b3)首先根据各相邻两个R峰位置计算得到的RR间期确定ECG记录片段中除最后一个T波外其余T波的起点搜索范围，然后采用滑动窗口定位ECG记录片段中除最后一个T波外其余T波的起点；定位T波起点的过程具体如下：将滑动窗口放在搜索左边界t1处，然后将滑动窗口向右移动，当滑动窗口的积分面积取到最大时，滑动窗口左端点t的取值就是T波起点的位置；b4)首先根据各相邻两个R峰位置计算得到的RR间期确定ECG记录片段中除最后一个T波外其余T波的终点搜索范围，然后采用滑动窗口定位ECG记录片段中除最后一个T波外...

【专利技术属性】
技术研发人员：蒋明峰，卞飞彪，张鞠成，王志康，李杨，
申请(专利权)人：浙江理工大学，
类型：发明
国别省市：