[发明专利]基于形态学多分辨率分解的心电信号预处理方法

说明书

技术领域

本发明涉及一种心电信号预处理方法，尤其涉及一种面向远程医疗系统中资源受限心电处理模块中的心电信号预处理。

背景技术

心脏病已经成为人类健康的头号杀手，目前全世界每年约有1750万人死于心脏病，占全部死亡人数的30%。作为诊断心脏生理、病理状态的重要依据，心电（ECG，Electrocardiogram）诊断技术自其十九世纪初诞生以来，已经发展成为一门专业学科，而对于ECG信号的采集、处理、存储、传输、远程实时监控及自动诊断等方面更是当今该领域的研究重点。

ECG是将电极置于体表特定位置所记录到的心脏电活动的图形。典型的心电图如图1所示，主要包括以下几部分，且这几部分为临床心电图的诊断所依据的重要特征参量：P波在QRS波群之前，前半部分主要为右心房的电活动，后半部分主要为左心房的电活动；QRS波群代表电激动通过心室肌；T波代表心室的复极过程，波形比较宽且圆钝；ST段代表心室的所有部分均已除极或除极末期与复极初期相互叠加形成的彼此抵消阶段；RR期间为相邻心动周期R波峰值之间的时间间隔。

ECG信号中主要包含三类干扰：人体呼吸或测量电极移位造成的基线漂移、工频干扰以及肌电或其他原因带来的高频噪声。如果对ECG信号中的干扰不进行预处理，会直接影响到QRS复波检测的正确率，增大医生正确诊断的难度。因此，对ECG信号进行预处理是ECG信号处理的首要任务，是后续处理的基础。

ECG信号中包含的三类干扰具有不同特点。其中，基线漂移通常频率很低，工频干扰固定为50/60Hz，而肌电等干扰频率较高。对于工频干扰，通常利用50/60Hz陷波器进行滤除；其余两种干扰是接下来讨论的重点。

早期采用的噪声抑制技术是基于带通滤波的方法。然而，带通的线性滤波器具有固定的截止频率，会使ST段和QRS复波严重失真。同时不能够自适应地跟踪ECG信号的时变特性，即每心动周期作为一周期的伪随机性。此外，一些自适应滤波技术被用来对ECG信号进行噪声抑制。大多数自适应去噪方法都是基于最小二次均方（LMS）或迭代最小二乘（RLS）准则。它们逐渐缩小输入信号与参考信号之间的均方误差。但某些情况下，难于获得合适的参考信号限制了此类方法的应用。

随着小波分析技术的发展及成熟，由于良好的时频特性，使其在ECG信号预处理领域中获得了广泛认可。基于小波分析的ECG信号预处理方法首先将ECG信号进行小波变换，由于基线漂移、高频干扰与信号所在分量不同，可以通过阈值操作将它们分离。研究的重点在于选择适合ECG信号的小波函数、去噪阈值的自适应选择以及优化快速算法等方面。

基于小波分析技术的心电信号预处理方法，由于不能够从频率上精确地分离信号与基线漂移及高频干扰，因此很容易造成不能够充分滤除基线漂移以及将ECG信号中频率相对高的分量与高频噪声一同滤除，导致波形失真的情况出现。

形态滤波是另外一种有效的非线性ECG信号预处理方法。基本思想是，应用直线型结构元素构成的滤波器校正基线漂移，应用三角型结构元素构成的滤波器滤除高频噪声。研究的重点在于形态滤波器的构成、结构元素的选择等方面。

基于形态滤波的预处理方法，应用直线型结构元素构造多分辨率形态滤波器实质上是应用一系列直线代替原始信号中的曲线，因此会造成一些信号特征的缺失；同时，随着结构元素长度的增加，滤波后信号失真也随之增大，这显然不能满足ECG信号预处理的高保真度需求。

为进一步优化性能，近来又提出结合小波变换与形态滤波的预处理方法。基本思想是，应用小波变换方法滤除高频噪声，应用形态滤波方法校正基线漂移。

基于小波变换与形态滤波的预处理方法，这在一定程度上改善了特征波形失真，不能够避免直线型结构元素造成的块状效应。

发明内容

为了解决现有技术中问题，本发明提供了一种基于形态学多分辨率分解的心电信号预处理方法，

步骤1：应用直线型结构元素对(B_o,B_c)的形态滤波器修正基线漂移；

步骤2：应用提升方法构造形态学多分辨率分解心电信号至第二层和第三层；

步骤3：选择第二层与第三层信号分量中QRS复波失真更小的分量作为预处理后的心电信号。

作为本发明的进一步改进，步骤1中，所述结构元素B_o,B_c长度的选择依赖于特征波形的持续时间T_w与心电信号的采样频率F_s。

作为本发明的进一步改进，所述特征波形的持续时间小于0.2s。