变异性心电信号时空特征提取方法技术

本发明专利技术涉及一种变异性心电信号时空特征提取方法，包括以下步骤：S1、获取原始心电信号；S2、构建时空特征提取模型；S3、利用时空特征提取模型提取心电信号的时空特征；S4、构建异常信号检测模型；S5、使用异常信号检测模型对心电信号异常进行检测。本发明专利技术是一种基于深度学习的心电信号异常检测算法，它融合了基于Transformer的模型和基于CNN的模型，可以有效地提取心电信号的时空特征并检测异常。本发明专利技术不仅考虑了导联之间的相互关系，还利用一维卷积和二维卷积提取相关的特征信号进行异常检测，从而提高了检测的精度，并具有更高的准确性、敏感性、特异性和F1评分。特异性和F1评分。特异性和F1评分。

【技术实现步骤摘要】
变异性心电信号时空特征提取方法


[0001]本专利技术涉及一种心电信号采集方法，具体地说是一种变异性心电信号时空特征提取方法。

技术介绍

[0002]时间序列分析是统计学的分支之一，近年来主流的分析方法主要基于机器学习和深度学习。然而机器学习受特征工程的制约，需要人工对大量提取的特征中进行有效的特征选择，从而影响了该方法的效率。因此，基于深度学习的时间序列分析方法应运而生，该方法克服了人为挑选特征的困难，可以自动提取时间序列中所蕴含的特征。如今，该方法已被广泛应用，并已经证明可以很好的提取时间序列中的异常特征。然而，对于时间序列中的空间特征提取这一问题还有待解决。
[0003]近年来，随着人们物质生活水平的提高以及人口老龄化进程的加快，全球心血管疾病的发病率呈逐年上升趋势。目前，心血管疾病已成为世界范围内死亡的首因，全球每年有1300万余人死于心血管疾病。因此，精准地诊断出心脏疾病信号就显得尤为重要。研究表明，心脏发出的电信号不仅具有时间序列特征，还可以反映心脏电活动的空间特征。心电信号(ECG)是记录心电图机采集的心脏在每一个心动周期产生电活动的变化图，其完整地描述了心脏在三维立体空间中的时空活动。不同生理和病理条件下的人群，心脏各部位形态不同，ECG信号及其时空特征自然也存在很大的差别。因此，提取完整的时空特征、捕捉心电信号的特异性是很有必要的。

技术实现思路

[0004]本专利技术的目的就是提供一种变异性心电信号时空特征提取方法，以解决现有心电信号采集方法对空间特征不能有效提取的问题。<br/>[0005]本专利技术的目的是这样实现的：
[0006]一种变异性心电信号时空特征提取方法，包括以下步骤：
[0007]S1、在获取原始心电信号后，先用基于Daubechies 6小波基的离散小波变换去除原始心电信号中的噪声；用Pan
‑
Tompkins算法检测原始心电信号中的R波；再根据一个完整的心跳周期以及信号的采样频率，对每个心跳周期的心电信号采集R波前250个采样点和R波后400个采样点，共651个采样点，构成一个完整的心拍。
[0008]S2、构建时空特征提取模型，所述时空特征提取模型包括Transformer编码器和全连接层；所述Transformer编码器用于提取心电导联之间的相关性信息，所述全连接层为模型的输出层，用于对提取的相关性信息进行融合，并输出表示时空特征的矩阵。
[0009]S3、利用时空特征提取模型提取心电信号的时空特征，所述时空特征提取模型的输入数据为：
[0010][0011]其中，a
i,j
为输入数据中的一个采样点，i∈(Ⅰ,II,Ⅲ,...,V6)为第i导联，j∈(1,2,3,...,651)为第j个心电心跳采样点。
[0012]输出时空特征矩阵为：M∈R3×
651
，具体表示为：
[0013][0014]其中，b
i,j
为时空特征中的一个采样点，是指时空特征矩阵M的三个正交方向的时变坐标，j∈(1,2,3,...,651)是指坐标的长度。
[0015]S4、构建异常信号检测模型，所述异常信号检测模型包括一维卷积层、二维卷积层和输出层；所述一维卷积层用于从时空特征矩阵M的三个时变坐标中提取相关的定时信息；所述二维卷积层是基于残差结构的卷积层，用于对二维特征图进行进一步的提取和融合；所述输出层包括全局平均池化层、全连接层和SoftMax分类器，用于输出最终的异常检测结果。
[0016]S5、使用异常信号检测模型对异常的心电信号进行检测。
[0017]进一步地，本专利技术提取方法还包括以下步骤：
[0018]S6、对时空特征提取模型和异常信号检测模型分别进行训练，PTB数据库中的所有数据按照训练集﹕验证集﹕测试集＝3﹕1﹕1的比例分开，时空特征提取模型和异常信号检测模型在相同的数据集上进行训练。
[0019]进一步地，在训练中引入均方误差损失函数来计算损失，以有效地反映估计值与实际值之间的差异：
[0020][0021]其中，y表示输入数据，表示模型的输出，x
i
和表示数据中的相应元素，n表示数据的长度。
[0022]进一步地，对于异常信号检测模型是以交叉点损失函数来表示预测值与实际值之间的差异：
[0023]L
CE
(P,Q)＝
‑
plog(q)
‑
(1
‑
p)log(1
‑
q)
[0024]其中，p为实际结果，q为预测结果，L
CE
(P,Q)为损失。
[0025]进一步地，所述Transformer编码器包括多头注意块和MLP块；所述MLP块用于融合多头注意提取的线索之间的相关性特征；所述多头注意块的多头注意力是建立在自注意的基础上；所述自注意是对于每个Query向量，计算它与所有Key向量的相似度，然后，将相似度乘以对应的Value向量，得到Query向量的加权和，也就是输出中与Query向量对应的值，其中，Query向量、Key向量和Value向量是通过映射输入获得的矩阵；自我注意的计算过程为：
[0026][0027]其中，QK
T
表示Query向量与Key向量之间的相似度，d
K
表示Query向量、Key向量和Value向量的维度；用相似度QK
T
除以输入维度来归一化相似度QK
T
。
[0028]进一步地，多头注意的计算过程如下：
[0029]MultiHead(X
in
)＝Concat(head1,
…
,head
h
)W
O
[0030]head
i
＝Attention(X
in
W
iQ
,X
in
W
iK
,X
in
W
iV
)
[0031]其中，W
iQ
，W
iK
和W
iV
是可训练投影参数矩阵，Concat()表示连接多个heat的操作，W
O
是可训练参数的权重矩阵。
[0032]进一步地，异常信号检测模型中的所述一维卷积层为三个加权的一维卷积块，先各从时空特征矩阵M的三个时变坐标中的一个提取一维特征，再对提取的这三个一维特征进行扩展，形成二维特征图；所述二维卷积层对二维特征图进行进一步的提取和融合；全局平均池化层获取每个特征地图的全局信息，全连接层实现数据降维，SoftMax分类器将全连接层的输出映射到心电信号异常的检测结果。
[0033]本专利技术是一种基于深度学习的心电信号异常检测算法，它融合了基于Transformer的模型和基于CNN的模型，可以有效地提取心电信号的时空特征并检测异常。与本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种变异性心电信号时空特征提取方法，其特征是，包括以下步骤：S1、在获取原始心电信号后，先用基于Daubechies 6小波基的离散小波变换去除原始心电信号中的噪声；用Pan
‑
Tompkins算法检测原始心电信号中的R波；再根据一个完整的心跳周期以及信号的采样频率，对每个心跳周期的心电信号采集R波前250个采样点和R波后400个采样点，共651个采样点，构成一个完整的心拍；S2、构建时空特征提取模型，所述时空特征提取模型包括Transformer编码器和全连接层；所述Transformer编码器用于提取心电导联之间的相关性信息，所述全连接层为模型的输出层，用于对提取的相关性信息进行融合，并输出表示时空特征的矩阵；S3、利用时空特征提取模型提取心电信号的时空特征，所述时空特征提取模型的输入数据为：其中，a
i,j
为输入数据中的一个采样点，i∈(Ⅰ,II,Ⅲ,...,V6)为第i导联，j∈(1,2,3,...,651)为第j个心电心跳采样点；输出时空特征矩阵为：M∈R3×
651
，具体表示为：其中，b
i,j
为输出时空特征中的一个采样点，是指时空特征矩阵M的三个正交方向的时变坐标，j∈(1,2,3,...,651)是指坐标的长度；S4、构建异常信号检测模型，所述异常信号检测模型包括一维卷积层、二维卷积层和输出层；所述一维卷积层用于从时空特征矩阵M的三个时变坐标中提取相关的定时信息；所述二维卷积层是基于残差结构的卷积层，用于对二维特征图进行进一步的提取和融合；所述输出层包括全局平均池化层、全连接层和SoftMax分类器，用于输出最终的异常检测结果；S5、使用异常信号检测模型对异常的心电信号进行检测。2.根据权利要求1所述的基于云边协同的可撤销的访问控制方法，其特征是，还包括以下步骤：S6、对时空特征提取模型和异常信号检测模型分别进行训练，PTB数据库中的所有数据按照训练集﹕验证集﹕测试集＝3﹕1﹕1的比例分开，时空特征提取模型和异常信号检测模型在相同的数据集上进行训练。3.根据权利要求2所述的变异性心电信号时空特征提取方法，其特征是，在训练中引入均方误差损失函数来计算损失，以有效地反映估计值与实际值之间的差异：其中，y表示输入数据，表示模型的输出，x
i
和表示数据中的相应...

【专利技术属性】
技术研发人员：刘秀玲，熊鹏，张杰烁，杜海曼，李昆霖，何聪，
申请(专利权)人：河北大学，
类型：发明
国别省市：