一种基于谱空间表示学习的ECG波形聚类方法及系统技术方案

本发明专利技术提出了一种基于谱空间表示学习的ECG波形聚类方法及系统，涉及心电波形聚类技术领域，获取多导联心电图数据，构建由多导联心电图数据及心搏的分导联谱空间相对特征组成的训练数据集；基于训练数据集，对构建的特征提取模型进行训练，学习谱空间表示；将待聚类的多导联心电图数据输入到训练好的特征提取模型中，得到心搏的分导联绝对特征；对心搏的分导联绝对特征进行聚类，得到心搏的聚类结果；本发明专利技术利用神经网络学习谱空间中心搏的相对特征，对波形进行聚类，注重特征表达质量方面的提升，避免传统波形聚类方法直接提取特征后续特征处理的繁琐问题，解决动态心电实际应用分析过程中的数据量庞大、计算复杂度高的问题。题。题。

【技术实现步骤摘要】
一种基于谱空间表示学习的ECG波形聚类方法及系统


[0001]本专利技术属于心电波形聚类
，尤其涉及一种基于谱空间表示学习的ECG波形聚类方法及系统。

技术介绍

[0002]本部分的陈述仅仅是提供了与本专利技术相关的
技术介绍
信息，不必然构成在先技术。
[0003]心电图(Electrocardiogram,ECG)作为能够直观展现心脏兴奋活动过程的技术，自诞生以来，被广泛应用于心血管等领域的临床健康检查；而为了能够更直观、实时观察监测患者一个生理周期24小时的心脏工作情况，动态心电图技术诞生，经过几十年的发展，其在无创心电辅助诊断
应用非常广泛，成为临床上最为常用的心脏检查和医疗监护工具之一。
[0004]动态心电图通常采集时间为24小时，这使得心电图数据一般会包括8万左右数量的心搏，对于如此大规模的心电数据，医生必须依赖可靠的心电分析软件对患者心电图进行分析诊断。
[0005]在动态心电分析中，波形数据的叠加处理分析是必不可缺的，多家心电方面的厂商将心电分析的叠加处理分析部分称为模板分析，即将叠加相似度高的心搏归为一个模板(类别)。叠加相似度使用交并比(IOU)衡量更为准确，而在通过对两两相似度进行聚类的方法中，谱聚类是一个较直观有效的方法；若使用谱聚类方法，动态心电图通常包括8万左右的心搏，需要处理维度约为(8万
×
8万)的矩阵，因此，现有动态心电实际应用分析过程存在数据量庞大、计算复杂度高的问题。
[0006]还存在一种基于深度学习的波形聚类方法，采用深度学习提取特征，同样因为数据量庞大，存在特征丢失及重复问题，而且提取特征后，还需要大量的初级特征处理操作，所以特征后续处理繁琐，因此整体效率不高。

技术实现思路

[0007]为克服上述现有技术的不足，本专利技术提供了一种基于谱空间表示学习的ECG波形聚类方法及系统，利用神经网络学习谱空间中心搏的相对特征，对波形进行聚类，注重特征表达质量方面的提升，避免传统波形聚类方法直接提取特征后续特征处理的繁琐问题，解决动态心电实际应用分析过程中的数据量庞大、计算复杂度高的问题。
[0008]为实现上述目的，本专利技术的一个或多个实施例提供了如下技术方案：
[0009]本专利技术第一方面提供了一种基于谱空间表示学习的ECG波形聚类方法；
[0010]一种基于谱空间表示学习的ECG波形聚类方法，包括：
[0011]获取多导联心电图数据，构建由多导联心电图数据及心搏的分导联谱空间相对特征组成的训练数据集；
[0012]基于训练数据集，以多导联心电图数据为输入、心搏的分导联绝对特征为输出，对
构建的特征提取模型进行训练，学习谱空间表示，其中，基于心搏的分导联谱空间相对特征与分导联绝对特征的欧式距离构建损失函数，对特征提取模型的训练进行约束；
[0013]将待聚类的多导联心电图数据输入到训练好的特征提取模型中，得到心搏的分导联绝对特征；
[0014]对得到的心搏的分导联绝对特征进行聚类，得到心搏的聚类结果。
[0015]进一步的，所述训练数据集的构建步骤为：
[0016]对多导联心电图数据中的心搏位置进行人工识别，按照时间顺序排列识别出来的心搏位置，构建心搏位置数据集合；
[0017]基于心搏位置数据集合，对每个导联的心搏两两做IOU计算，构建每个导联的邻接矩阵；
[0018]根据每个导联的邻接矩阵，计算心搏的分导联谱空间相对特征；
[0019]多导联心电图数据与最终计算得到的心搏的分导联谱空间相对特征，组成训练数据集。
[0020]进一步的，所述对每个导联的心搏两两做IOU计算，具体为：
[0021]对每个导联上的任意两个心搏数据进行补位操作；
[0022]通过计算补位后的两个心搏数据之间的IOU值，构建每个导联的连接矩阵，IOU值的具体计算方式为：
[0023][0024][0025][0026]其中，pre、pos表示两个心搏数据，pre_beat、pos_beat分别表示补位后的两个心搏数据，pre_up、pos_up分别表示pre_beat、pos_beat小于零的部分置为0的波形，pre_down、pos_down分别表示pre_beat、
[0027]pos_beat大于零的部分置为0的波形。
[0028]进一步的，所述计算心搏的分导联谱空间相对特征，具体为：
[0029]根据每个导联的邻接矩阵，计算度矩阵，由邻接矩阵矩阵和度矩阵，计算得到拉普拉斯矩阵；
[0030]对拉普拉斯矩阵进行标准化，并对标准化后的拉普拉斯矩阵特征分解,得到多个特征值及对应的初级特征向量；
[0031]将多个特征值进行升序排列，并根据特征值的顺序对特征向量进行排序，利用二
阶差分方法寻找升序排列的拐点对应的初级特征值，记拐点对应的初级特征值在特征值升序排列中的序号为K；
[0032]取排序后的前K个特征向量，作为心搏的相对特征，即得到当前导联中心搏的分导联谱空间相对特征。
[0033]进一步的，所述分导联绝对特征的计算方式为：
[0034]通过卷积神经网络提取多导联心电图数据中的初级特征；
[0035]依据导联及心搏位置，得到每个导联上每个心搏对应的初级特征；
[0036]将每个导联上的心搏特征取平均值，得到每个导联上的心搏绝对特征，构成心搏的分导联绝对特征。
[0037]进一步的，所述基于心搏的分导联谱空间相对特征与分导联绝对特征的欧式距离构建损失函数，具体为：
[0038]计算分导联谱空间相对特征的欧式距离矩阵和分导联绝对特征的欧式距离矩阵
[0039][0040][0041]其中，E
re,g
表示第g个导联的心搏谱空间相对特征，E
re,g
(i)为E
re,g
的第i行数据，[E
re,g
(i)]1为E
re,g
的第i行第l列的数据，E
ab,g
表示第g个导联的心搏绝对特征，K、E_S分别为心搏谱空间相对特征和心搏绝对特征的维数；
[0042]所述损失函数，两个欧式距离矩阵做均方误差mse，具体为：
[0043][0044]其中，A1
g
、A2
g
均为(N,N)维矩阵。
[0045]进一步的，所述特征提取模型的损失函数，还包括重合心搏的均方误差mse，具体为：
[0046]将多导联心电图数据切割两次，得到两个多导联子心电图数据，保证切割后的两个多导联子心电图数据有重合心搏；
[0047]利用特征提取模型，从两个多导联子心电图数据中分别得到各自的心搏的分导联绝对特征；
[0048]基于有重合心搏的两个多导联子心电图数据的分导联绝对特征，计算每个重合心搏的均方误差mse，具体为：
[0049][0050]其中，C1
g...

【技术保护点】

【技术特征摘要】
1.一种基于谱空间表示学习的ECG波形聚类方法，其特征在于，包括：获取多导联心电图数据，构建由多导联心电图数据及心搏的分导联谱空间相对特征组成的训练数据集；基于训练数据集，以多导联心电图数据为输入、心搏的分导联绝对特征为输出，对构建的特征提取模型进行训练，学习谱空间表示，其中，基于心搏的分导联谱空间相对特征与分导联绝对特征的欧式距离构建损失函数，对特征提取模型的训练进行约束；将待聚类的多导联心电图数据输入到训练好的特征提取模型中，得到心搏的分导联绝对特征；对得到的心搏的分导联绝对特征进行聚类，得到心搏的聚类结果。2.如权利要求1所述的一种基于谱空间表示学习的ECG波形聚类方法，其特征在于，所述训练数据集的构建步骤为：对多导联心电图数据中的心搏位置进行人工识别，按照时间顺序排列识别出来的心搏位置，构建心搏位置数据集合；基于心搏位置数据集合，对每个导联的心搏两两做IOU计算，构建每个导联的邻接矩阵；根据每个导联的邻接矩阵，计算心搏的分导联谱空间相对特征；多导联心电图数据与最终计算得到的心搏的分导联谱空间相对特征，组成训练数据集。3.如权利要求2所述的一种基于谱空间表示学习的ECG波形聚类方法，其特征在于，所述对每个导联的心搏两两做IOU计算，具体为：对每个导联上的任意两个心搏数据进行补位操作；通过计算补位后的两个心搏数据之间的IOU值，构建每个导联的连接矩阵，IOU值的具体计算方式为：体计算方式为：其中，pre、pos表示两个心搏数据，pre_beat、pos_beat分别表示补位后的两个心搏数据，pre_up、pos_up分别表示pre_beat、pos_beat小于零的部分置为0的波形，pre_down、pos_down分别表示pre_beat、pos_beat大于零的部分置为0的波形。
4.如权利要求2所述的一种基于谱空间表示学习的ECG波形聚类方法，其特征在于，所述计算心搏的分导联谱空间相对特征，具体为：。根据每个导联的邻接矩阵，计算度矩阵，由邻接矩阵矩阵和度矩阵，计算得到拉普拉斯矩阵；对拉普拉斯矩阵进行标准化，并对标准化后的拉普拉斯矩阵特征分解,得到多个特征值及对应的初级特征向量；将多个特征值进行升序排列，并根据特征值的顺序对特征向量进行排序，利用二阶差分方法寻找升序排列的拐点对应的初级特征值，记拐点对应的初级特征值在特征值升序排列中的序号为K；取排序后的前K个特征向量，作为心搏的相对特征，即得到当前导联中心搏的分导联谱空间相对特征。5.如权利要求1所述的一种基于谱空间表示学习的ECG波形聚类方法，其特征在于，所述分导联绝对特征的计算方式为：通过卷积神经网络提取多导联心电图数据中的初级特征；依据导联及心搏位置，得到每个导联上每个心搏对应的初级特征；将每个导联上的心搏特征取平均值，得到每个导联上的心搏绝对特征，构成心搏的分导联绝对特征。6.如权利要求1所述的一种基于谱空间表示学习的ECG波形聚类方法，...

【专利技术属性】
技术研发人员：张伯政，吴军，李彦芳，高希余，王玉杰，房志文，
申请(专利权)人：山东晨泽人工智能研究院有限公司，
类型：发明
国别省市：