[发明专利]一种基于窗口方差变换的穿戴式心电信号R波检测方法

说明书

一种基于窗口方差变换的穿戴式心电信号R波检测方法，通过窗口方差变换抑制对R波检测影响最大的幅值的噪声和波形，利用自适应更新的阈值和RR间期实现了高准确率、抗干扰性强的实时R波检测方法，以信号原始波形形态为出发点，提取了能更好的表示原始信号的特征，适用于多类心电信号的R波检测。有效抑制对R波检测有干扰的大幅值噪声的影响，可以改善传统基于阈值的方法在检测R波时因噪声幅值较大而导致的误检测率上升的问题，适用于各种包含混合叠加噪声的心电信号的R检测，尤其是穿戴式动态心电信号。

技术领域

本发明涉及心电信号波形识别技术领域，具体涉及一种基于窗口方差变换的穿戴式心电信号R波检测方法。

背景技术

QRS波群是心电信号中幅值最大，特征最明显的波形，其形态多呈较强的脉冲状，持续时间在0.06s-0.1s之间。心电信号波形检测中，R波的准确检测是其他波形检测的基础。目前常用的R波识别方法有数字滤波法、小波变换法、自适应阈值及卷积神经网络等。然而，上述检测方法大多是基于高质量的ECG信号设计的，其检测结果易受噪声干扰，鲁棒性差，尤其是在处理易受运动伪影等混合噪声强烈干扰的穿戴式动态心电信号。传统的R波检测算法并不适用于噪声等级较高的心电信号，因此，如何在高干扰性的心电信号中准确、高效识别R波是心电信号波形识别的一个较大挑战。

发明内容

本发明为了克服以上技术的不足，提供了一种高鲁棒性，有效抑制R波检测有干扰的大幅值噪声的影响的基于窗口方差变换的穿戴式心电信号R波检测方法。

本发明克服其技术问题所采用的技术方案是：

一种基于窗口方差变换的穿戴式心电信号R波检测方法，包括如下步骤：

a)利用计算机对穿戴式动态心电信号中噪声进行抑制，得到降噪后的心电信号；

b)利用窗口方差平方变换进行波形增强处理，突出降噪后的心电信号中QRS波群，抑制其他波形及噪声；

c)根据QRS波群的幅值与心电信号中其他波形的差异，进行R波识别和检测，得到R波集合R_S。

进一步的，步骤a)中利用滑动平均滤波法和中值滤波法对含噪的穿戴式动态心电信号y中的噪声进行过滤，通过公式计算得到去除高频噪声和基线漂移的信号式中表示滤波窗口为M的滑动平均滤波法，M为滤波窗口的宽度，bw_n为中值滤波提取的第n个采样点的基线。进一步的，步骤b)包括如下步骤：

b-1)提取降噪后的心电信号的局部极值点，根据局部极值点对降噪后的心电信号按照前后时间窗口大小为w进行窗口划分，得到极值点窗口信号seg_p＝[S₁,S₂,...,S_K]，S_K为第K个极值点窗口信号；

b-2)利用公式计算第K个局部极值点所对应的窗口的方差v_k，式中为去除高频噪声和基线漂移的信号中的第t个采样点的信号，得到极值点窗口信号方差wvt_p＝[v₁,v₂,...,v_K]；

b-3)将相邻的两个局部极值点间的n个采样点划分为一个窗口，将该窗口方差设置为0，非极值点对应的窗口信号seg_s＝[ns₁,ns₂,...,ns_K-1]，窗口信号seg_s对应的方差为wvt_s＝[0,0,...,0]；

b-4)依据窗口内采样点索引大小，合并极值点窗口信号方差wvt_p和非极值点方差wvt_s，得到整个信号对应的窗口方差变换wvt，对于wvt进行平方操作，得到QRS波群被增强的信号swvt。

进一步的，步骤c)包括如下步骤：