一种跳跃突变噪声的动态心电信号基线漂移校正方法
阅读说明：本技术 一种跳跃突变噪声的动态心电信号基线漂移校正方法 (Dynamic electrocardiosignal baseline drift correction method for jump mutation noise ) 是由 谢小云 王英龙 刘辉 舒明雷 周书旺 刘照阳 于 2019-10-25 设计创作，主要内容包括：一种跳跃突变噪声的动态心电信号基线漂移校正方法,通过准确检测心电信号中包含的跳跃突变区域,并分段单独处理,可以改善传统滤波方法处理该类区域时基线提取不准确导致的心电信号变形问题。适用于各种含有基线漂移的心电信号的基线漂移校正处理,尤其是可穿戴式的动态心电信号。得到的基线漂移校正信号的波形明显优于其他方法。(A jump mutation noise dynamic electrocardiosignal baseline drift correction method can solve the problem of electrocardiosignal deformation caused by inaccurate baseline extraction when a traditional filtering method is used for processing a jump mutation region contained in an electrocardiosignal through accurately detecting the jump mutation region and carrying out segmentation and independent processing. The method is suitable for baseline drift correction processing of various electrocardiosignals containing baseline drift, in particular to wearable dynamic electrocardiosignals. The waveform of the resulting baseline drift correction signal is significantly better than other methods.)
技术领域
本发明涉及心电信号噪声处理技术领域，具体涉及一种跳跃突变噪声的动态心电信号基线漂移校正方法。
背景技术
基线漂移是心电信号中频率小于1Hz的常见噪声。基线校正对于正确诊断ST段抬高或降低等心电异常十分重要。目前常用的基线校正方法有高通滤波、中值滤波、小波变换、样条插值、数学形态学滤波、自适应滤波等。
随着Holter等动态心电监测设备的广泛使用、尤其是近年来小型化的穿戴式、便携式心电监测终端的快速发展，运动干扰成为心电信号的主要噪声源之一，基线漂移变得更加严重，不仅幅度增加而且频率增大，与正常心电信号的频率有一定重合。更为严重的是，在较大幅度动作下人体皮肤与电极之间的阻抗会随位移变化出现快速的升高或降低，导致采集的信号幅值在不足100ms的较短时间内产生跳跃突变，变化可达0.5mV以上(如附图1的区域A、B、C所示)。
然而，上述基线校正方法均没有针对信号局部特征单独处理，基线提取中对跳跃突变区域不敏感，不能有效校正突变区域内产生的较为严重的基线漂移现象。只有中值滤波方法能够检测到突变发生的位置，并进行校正。但中值滤波方法对突变区域定位不准确或不能检测完整突变区域并处理。因此，对动态心电信号中跳跃突变区域进行正确识别和有效消除是动态心电信号基线漂移校正阶段一个较大的挑战。
发明内容
本发明为了克服以上技术的不足，提供了一种能够准确定位发生跳跃突变的区域，通过对跳跃突变区域和平缓变化区域分段单独处理，得到较优的基线漂移校正结果的动态心电信号基线漂移校正方法。
本发明克服其技术问题所采用的技术方案是：
一种跳跃突变噪声的动态心电信号基线漂移校正方法，包括如下步骤：
a)利用计算机对心电信号中跳跃突变区域进行检测，得到心电信号中的平缓变化区域和跳跃突变区域；
b)对跳跃突变区域进行单独处理，提取跳跃突变区域内的基线bw_u；
c)提取平缓变化区域内的基线bw_c，利用滑动平均法对提取的基线进行平滑处理；
d)将提取的跳跃突变区域内的基线bw_u与平缓变化区域内的基线bw_c进行合并，合并后得到初始基线bw，得到校正后的数据，跳跃突变区域与平缓变化区域连接点间的误差通过滤波组合的方法进行平滑处理。
进一步的，步骤a)通过以下步骤处理：
a-1)利用计算机对心电信号中含基线漂移和跳跃突变噪声的心电信号y使用宽度为200ms和600ms的中值滤波器进行过滤，得到去除R波影响的基线信号y'；
a-2)通过公式
a-3)定义阀值thr_min，thr_min＝1.1*min(Li)，将信号Li中满足Li＞thr_min的区域加入候选集，不满足Li＞thr_min的区域加入平缓变化区域，计算变换后信号Li的香浓熵H、计算第i个候选区域香浓熵Hi及振幅变化范围Ampi，定义幅度变化阀值为thr_a；
a-4)遍历候选集中每个区域，将满足条件Hi＞H且Ampi＞thr_a的区域加入跳跃突变区域，将不满足条件Hi＞H且Ampi＞thr_a的区域加入平缓变化区域。
进一步的，步骤b)中利用滑动平均法平滑中值滤波基线信号y'，得到跳跃突变区域内的基线bw_u。
进一步的，步骤c)中通过遍历跳跃突变区域，对每个突变区域单独处理，将该区域内基线定义为(1-thr_r)*该区域内的含噪信号，处理后的基线定义为bw_c，其中thr_r为缩放率，
进一步的，步骤d)通过以下步骤处理：
d-1)通过公式
d-2)通过公式计算使用宽度为wf个样本点的平滑滤波方法对初始平滑基线bw1进行平滑处理，得到最终的基线bw'，通过公式y-bw'计算得到最终校正后的信号。
进一步的，步骤a-3)中通过公式
本发明的有益效果是：通过准确检测心电信号中包含的跳跃突变区域，并分段单独处理，可以改善传统滤波方法处理该类区域时基线提取不准确导致的心电信号变形问题。适用于各种含有基线漂移的心电信号的基线漂移校正处理，尤其是可穿戴式的动态心电信号。得到的基线漂移校正信号的波形明显优于其他方法。
附图说明
图1为心电信号中跳跃突变区域的示例图；
图2为本发明的跳跃突变区域的检测流程图；
图3为本发明的平缓变化区域和跳跃变化区域基线提取流程图；
图4为本发明的平滑基线流程图；
图5为本发明的方法流程图。
具体实施方式
下面结合附图1、附图2对本发明做进一步说明。
一种跳跃突变噪声的动态心电信号基线漂移校正方法，包括如下步骤：
a)利用计算机对心电信号中跳跃突变区域进行检测，得到心电信号中的平缓变化区域和跳跃突变区域；
b)对跳跃突变区域进行单独处理，提取跳跃突变区域内的基线bw_u；
c)提取平缓变化区域内的基线bw_c，利用滑动平均法对提取的基线进行平滑处理；
d)将提取的跳跃突变区域内的基线bw_u与平缓变化区域内的基线bw_c进行合并，合并后得到初始基线bw，得到校正后的数据，跳跃突变区域与平缓变化区域连接点间的误差通过滤波组合的方法进行平滑处理。
通过准确检测心电信号中包含的跳跃突变区域，并分段单独处理，可以改善传统滤波方法处理该类区域时基线提取不准确导致的心电信号变形问题。适用于各种含有基线漂移的心电信号的基线漂移校正处理，尤其是可穿戴式的动态心电信号。得到的基线漂移校正信号的波形明显优于其他方法。
进一步的，如附图2所示，步骤a)通过以下步骤处理：
a-1)利用计算机对心电信号中含基线漂移和跳跃突变噪声的心电信号y使用宽度为200ms和600ms的中值滤波器进行过滤，得到去除R波影响的基线信号y'；
a-2)通过公式
a-3)定义阀值thr_min，thr_min＝1.1*min(Li)，将信号Li中满足Li＞thr_min的区域加入候选集，不满足Li＞thr_min的区域加入平缓变化区域，计算变换后信号Li的香浓熵H、计算第i个候选区域香浓熵Hi及振幅变化范围Ampi，定义幅度变化阀值为thr_a；
a-4)遍历候选集中每个区域，将满足条件Hi＞H且Ampi＞thr_a的区域加入跳跃突变区域，将不满足条件Hi＞H且Ampi＞thr_a的区域加入平缓变化区域。
进一步的，如附图3所示，步骤b)中利用滑动平均法平滑中值滤波基线信号y'，得到跳跃突变区域内的基线bw_u。
进一步的，如附图3所示，步骤c)中通过遍历跳跃突变区域，对每个突变区域单独处理，将该区域内基线定义为(1-thr_r)*该区域内的含噪信号，处理后的基线定义为bw_c，其中thr_r为缩放率，
进一步的，如附图4所示，步骤d)通过以下步骤处理：
d-1)通过公式
d-2)通过公式
进一步的，步骤a-3)中通过公式计算香浓熵H，式中将信号按照其幅值分为n段，j＝{0,1,2,...K}，Pj为变换后信号属于第j段的概率值；通过公式
