一种基于图卷积空时特征融合选择的心律失常识别方法与流程

1.本发明涉及心律失常识别领域，具体而言，涉及一种基于图卷积空时特征融合选择的心律失常识别方法。


背景技术：

2.当前使用多导联ecg信号进行心律失常实时监测的方法主要包括三个方面：一是基于统计学规则的鉴别器方法；二是使用机器学习模型对ecg信号提取的特征信息作为输入进行分类；三是使用传统的基于欧氏距离的dnn模型对ecg信号进行分类。
3.在这类使用使用ecg信号进行心律失常类别的实时监测方法中，其描述ecg信号心律失常的信息量相对不够完整，在实际的监测过程中误判和漏判心律失常的概率较大；使用机器学习对ecg信号心律失常分类预测的过程中，得到具有较高准确性分类模型的关键是提取具有判别性的特征，但这种具有判别性的特征的提取需要有大量的先验知识，同时这样的特征工程十分费事费力；使用传统的基于欧氏距离的cnn模型对多导联ecg信号进行心律失常分类预测的过程中，由于cnn是通过欧氏距离对ecg信号的时间轴和导联的相关性进行描述，但对于使用多导联ecg信号进行心律失常的导联之间是一种不规则的关系，使用欧氏距离不能完全描述这种关联性，此时cnn在这类数据上无法保持位移不变性，也就导致实际的监测过程中对心律失常分类判别的准确性能不强，因此我们对此做出改进，提出一种基于图卷积空时特征融合选择的心律失常识别方法。


技术实现要素：

4.本发明的目的在于：针对目前存在的背景技术提出的问题。
5.为了实现上述发明目的，本发明提供了以下技术方案：一种基于图卷积空时特征融合选择的心律失常识别方法，包括，其特征在于，包括以下步骤：
6.步骤s1：将采集的历史数据标注样本后进行模型训练；
7.步骤s2：所述步骤s1进行模型训练后，使用滑动窗口法将实时采集的固定长度的多导心电信号，使用训练集共享的标准化参数进行标准化后作为输入；
8.步骤s21标准化的方法为，记变量序列为s1,
…
,sn，记其中平均值和标准差值分别为:μ和σ,则标准化后的序列为：
[0009][0010]
其中标准化后的序列值服从正态分布，μ为平均值，σ为标准差值，s为变量值，存储μ和σ的值；
[0011]
步骤s3：所述步骤s2标准化参数输入后，进行分类预测输出心律失常类别为标签。
[0012]
作为本技术优选的技术方案，所述步骤s1包括预处理单元、gcn空域特征提取单元、gru时域特征提取单元、特征融合单元，以及分类结果输出单元。
[0013]
作为本技术优选的技术方案，所述特征融合单元设置有注意力层，所述分类结果
输出单元设置有输出结果标签。
[0014]
作为本技术优选的技术方案，还包括构造邻接矩阵模块、gcn空域特征提取模块、swlstm时域特征提取模块、空时特征消除量纲影响模块、注意力机制融合特征模块和分类模块。
[0015]
作为本技术优选的技术方案，所述构造邻接矩阵模块，首先计算输入多导联ecg信号序列的协方差或相关矩阵，然后使用拉普拉斯laplacian归一化来估计邻接矩阵，具体过程如下：
[0016][0017]
其中in表示跟a相同阶数的单位矩阵，本发明中的n值为12。
[0018]
作为本技术优选的技术方案，所述gcn空域特征提取模块相较于cnn中的卷积运算，其过程为：首先对属性信息进行仿射变换，学习了属性特征之间的交互模式，然后迭代式的聚合邻居节点的特征，从而更新当前节点的特征，每多一层卷积运算，中心节点就能多融合进更外一圈的信息，可以更充分提取了多导联ecg信号中的特征信息并且筛选出具有区分度的信息，构建的gcn模型设置两层图卷积层，使用relu作为激活函数，通过全连接层将gcn提取的特征线性映射后输入特征融合模块。
[0019]
作为本技术优选的技术方案，所述swlstm时域特征提取模块，swlstm是将lstm中的输入门、遗忘门和输出门组合成一种称为共享门的新门，提取数据的时域特征，且在保障得到与lstm相近预测效果的同时大大降低了计算量并获得更快的收敛速度，其中构建的swlstm模型设置一层swlstm层，然后使用全连接层将swlstm提取的特征线性映射后输入特征融合模块。
[0020]
作为本技术优选的技术方案，所述空时特征消除量纲影响模块，将gcn提取的空域特征使用全连接层进行线性映射，同时将swlstm提取的时域特征使用相同单元的全连接层进行线性映射。
[0021]
作为本技术优选的技术方案，所述注意力机制融合特征模块，通过gcn提取的空域特征和swlstm提取的特征其在分类过程中的重要程度是不同的，因此，筛选出ecg信号经特征提取的有区分度的特征，注意力机制本质上为加权计算，通过训练过程中的权系数变化对空时域特征重要程度的筛选。
[0022]
作为本技术优选的技术方案，还包括建立多导ecg信号用以不同心律失常实时分类监测的gg-net模型，并对数据进行验证，包括正常n、房颤af、一级房室传导阻滞i-avb、左束支传导阻滞lbbb、右束支传导阻滞rbbb、房性早搏pac、室性早搏pvc、st段压低std和st段抬高ste九大类，使用准确性能指标：准确率、召回率、特异性和马太相关系数；
[0023]
其中准确率简称ac的计算公式为：
[0024]
其中召回率简称re的计算公式为：
[0025]
其中特异性简称sp的计算公式为：
[0026]
其中马太相关系数mc的计算公式为：
[0027]
其中tp、tn、fp和fn分别表示真阳性、真阴性、假阳性和假阴性。
[0028]
与现有技术相比，本发明的有益效果：
[0029]
在本技术的方案中：
[0030]
1)为充分提取多导ecg信号的特征信息，基于相关矩阵构造的图卷积神经网络和门控循环单元，提出了提取空时特征信息，然后通过相同全连接单元分别进行空时特征量纲消除，最后通过注意力机制进行特征融合的深度学习模型gtam；
[0031]
2)模型使用相关矩阵通过拉普拉斯laplacian归一化来估计邻接矩阵，作为图卷积的计算单元，充分表达出多导ecg信号各导联之间不规则的关系，获得了多导ecg信号用以不同心律失常实时分类监测的空域特征表示；
[0032]
3)利用swlstm在时序信号强大的处理能力，对多导ecg信号进行时域上的特征提取；
[0033]
4)通过注意力机制对提取的空域特征和时域特征融合，使得模型避免了人工特征所存在的有区分度特征表达不充分的问题，同时提升了模型的泛化性能，使得本发明提出的模型更适合基于多导ecg信号用以不同心律失常实时分类监测。
附图说明：
[0034]
图1为本技术提供的多导联ecg信号序列图；
[0035]
图2为本技术提供的相关矩阵图；
[0036]
图3为本技术提供的内部结构图；
[0037]
图4为本技术提供的流程图；
[0038]
图5为本技术提供的收敛曲线图。
具体实施方式
[0039]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图，对本发明实施例中的技术方案进行清楚、完整的描述。显然，所描述的实施例是本发明的一种具体实施方式，不限于全部的实施例。
[0040]
因此，以下对本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的部分实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0041]
如图1-图5，一种基于图卷积空时特征融合选择的心律失常识别方法，包括步骤s1：将采集的历史数据标注样本后进行模型训练；
[0042]
步骤s2：所述步骤s1进行模型训练后，使用滑动窗口法将实时采集的固定长度的
多导心电信号，使用训练集共享的标准化参数进行标准化后作为输入；
[0043]
步骤s21标准化的方法为，记变量序列为s1,
…
,sn，记其中平均值和标准差值分别为:μ和σ,则标准化后的序列为：
[0044][0045]
其中标准化后的序列值服从正态分布，μ为平均值，σ为标准差值，s为变量值，存储μ和σ的值。
[0046]
步骤s3：所述步骤s2标准化参数输入后，进行分类预测输出心律失常类别为标签。
[0047]
作为优选的实施方式，在上述方式的基础上，进一步的，所述步骤s1包括预处理单元、gcn空域特征提取单元、gru时域特征提取单元、特征融合单元，以及结果输出单元。
[0048]
作为优选的实施方式，在上述方式的基础上，进一步的，所述特征融合单元设置有注意力层，所述分类结果输出单元设置有输出结果标签。
[0049]
作为优选的实施方式，在上述方式的基础上，进一步的，还包括构造邻接矩阵模块、gcn空域特征提取模块、swlstm时域特征提取模块、空时特征消除量纲影响模块、注意力机制融合特征模块和分类模块。
[0050]
作为优选的实施方式，在上述方式的基础上，进一步的，所述构造邻接矩阵模块，首先计算输入多导联ecg信号序列的协方差或相关矩阵，然后使用拉普拉斯laplacian归一化来估计邻接矩阵，具体过程如下：
[0051][0052]
其中in表示跟a相同阶数的单位矩阵，本发明中的n值为12。
[0053]
作为优选的实施方式，在上述方式的基础上，进一步的，所述gcn空域特征提取模块相较于cnn中的卷积运算，其过程为：首先对属性信息进行仿射变换，学习了属性特征之间的交互模式，然后迭代式的聚合邻居节点的特征，从而更新当前节点的特征，每多一层卷积运算，中心节点就能多融合进更外一圈的信息，可以更充分提取了多导联ecg信号中的特征信息并且筛选出具有区分度的信息，构建的gcn模型设置两层图卷积层，使用relu作为激活函数，通过全连接层将gcn提取的特征线性映射后输入特征融合模块。
[0054]
作为优选的实施方式，在上述方式的基础上，进一步的，所述swlstm时域特征提取模块，swlstm是将lstm中的输入门、遗忘门和输出门组合成一种称为共享门的新门，提取数据的时域特征，且在保障得到与lstm相近预测效果的同时大大降低了计算量并获得更快的收敛速度，其中构建的swlstm模型设置一层swlstm层，然后使用全连接层将swlstm提取的特征线性映射后输入特征融合模块。
[0055]
作为优选的实施方式，在上述方式的基础上，进一步的，所述空时特征消除量纲影响模块，将gcn提取的空域特征使用全连接层进行线性映射，同时将swlstm提取的时域特征使用相同单元的全连接层进行线性映射。
[0056]
作为优选的实施方式，在上述方式的基础上，进一步的，所述注意力机制融合特征模块，通过gcn提取的空域特征和swlstm提取的特征其在分类过程中的重要程度是不同的，因此，筛选出ecg信号经特征提取的有区分度的特征，注意力机制本质上为加权计算，通过训练过程中的权系数变化对空时域特征重要程度的筛选。
[0057]
作为优选的实施方式，在上述方式的基础上，进一步的，还包括建立多导ecg信号用以不同心律失常实时分类监测的gg-net模型，并对数据进行验证，包括正常n、房颤af、一级房室传导阻滞i-avb、左束支传导阻滞lbbb、右束支传导阻滞rbbb、房性早搏pac、室性早搏pvc、st段压低std和st段抬高ste九大类，使用准确性能指标：准确率、召回率、特异性和马太相关系数；
[0058]
其中准确率简称ac的计算公式为：
[0059]
其中召回率简称re的计算公式为：
[0060]
其中特异性简称sp的计算公式为：
[0061]
其中马太相关系数mc的计算公式为：
[0062]
其中tp、tn、fp和fn分别表示真阳性、真阴性、假阳性和假阴性。
[0063]
工作模块：本发明在使用的过程中，步骤s1：将采集的历史数据标注样本后进行模型训练；步骤s2：所述步骤s1进行模型训练后，使用滑动窗口法将实时采集的固定长度的多导心电信号，使用训练集共享的标准化参数进行标准化后作为输入；步骤s21标准化的方法为，记变量序列为s1,
…
,sn，记其中平均值和标准差值分别为:μ和σ,则标准化后的序列为：其中标准化后的序列值服从正态分布，μ为平均值，σ为标准差值，s为变量值，存储μ和σ的值，步骤s3：所述步骤s2标准化参数输入后，进行分类预测输出心律失常类别为标签，还包括构造邻接矩阵模块、gcn空域特征提取模块、swlstm时域特征提取模块、空时特征消除量纲影响模块、注意力机制融合特征模块和分类模块，构造邻接矩阵模块，首先计算输入多导联ecg信号序列如图1所示的协方差或相关矩阵，本实施例以相关矩阵为例，然后对得到的相关矩阵如图2所示，然后使用拉普拉斯laplacian归一化来估计邻接矩阵，具体过程如下：
[0064][0065]
其中in表示跟a相同阶数的单位矩阵，本发明中的n值为12，所述gcn空域特征提取模块，相较于cnn中的卷积运算，gcn中的卷积运算最大的区别是具有显式的邻接矩阵表达。
其过程为：首先对属性信息进行仿射变换，学习了属性特征之间的交互模式，然后迭代式的聚合邻居节点的特征，从而更新当前节点的特征。其每多一层卷积运算，中心节点就能多融合进更外一圈的信息，可以更充分提取了多导联ecg信号中的特征信息并且筛选出具有区分度的信息，在提高了多导联ecg信号在不同心律失常实时分类监测方法的准确性能同时，也保证了分类系统的泛化性能。其中构建的gcn模型结果如图3所示，设置两层图卷积层，同时使用relu作为激活函数，然后使用全连接层将gcn提取的特征线性映射后输入特征融合模块，所述swlstm时域特征提取模块，swlstm是将lstm中的输入门、遗忘门和输出门组合成一种称为共享门的新门结构内部结构如图3所示，具有非线性拟合能力。能够有效提取数据的时域特征，且在保障得到与lstm相近预测效果的同时大大降低了计算量并获得更快的收敛速度。其中构建的swlstm模型结果如图3所示，设置一层swlstm层，然后使用全连接层将swlstm提取的特征线性映射后输入特征融合模块，所述空时特征消除量纲影响模块，将gcn提取的空域特征使用全连接层进行线性映射，同时将swlstm提取的时域特征使用相同单元的全连接层进行线性映射，这样就消除了不同域特征直接由于幅度不同导致对有区分度特征的提取不充分问题，其在模型训练过程中使用更加有效率的梯度下降以及反向传播，避免了梯度爆炸和梯度消失问题，同时简化计算过程，所述注意力机制融合特征模块，通过gcn提取的空域特征和swlstm提取的特征其在分类过程中的重要程度是不同的。因此，筛选出ecg信号经特征提取的有区分度的特征具有重要的价值，一个有效的方法是采用注意力机制。注意力机制本质上是一种加权计算，通过训练过程中的权系数变化实现空时域特征重要程度的筛选，还包括建立多导ecg信号用以不同心律失常实时分类监测的gg-net模型，并对数据进行验证，包括正常(n)、房颤af、一级房室传导阻滞i-avb、左束支传导阻滞lbbb、右束支传导阻滞rbbb、房性早搏pac、室性早搏pvc、st段压低std和st段抬高ste九大类，使用准确性能指标：准确率、召回率、特异性和马太相关系数。
[0066]
其中准确率简称ac的计算公式为：
[0067]
其中召回率简称re的计算公式为：
[0068]
其中特异性简称sp的计算公式为：
[0069]
其中马太相关系数mcc的计算公式为：
[0070]
其中tp、tn、fp和fn分别表示真阳性、真阴性、假阳性和假阴性。
[0071]
以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案，尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明，但本发明不局限于上述具体实施方式，因此任何对本发明进行修改或等同替换；而一切不脱离发明的精神和范围的技术方案及其改进，其均涵盖在本发明的权利要求范围当中。