本发明专利技术提供了一种基于形态学多分辨率分解的心电信号预处理方法，步骤1：应用直线型结构元素对的形态滤波器修正基线漂移；步骤2：应用提升方法构造形态学多分辨率分解心电信号至第二层和第三层；步骤3：选择第二层与第三层信号分量中QRS复波失真更小的分量作为预处理后的心电信号。本发明专利技术应用直线型结构元素对修正基线漂移，能够在较大程度上减小波形失真。应用形态学多分辨率分解方法滤除肌电干扰等高频噪声，性能优于基于小波变换及形态滤波的预处理方法，且方法简单、易于实现。采用提升方法构造多层形态学分解，由于选择了适合心电信号的预测、更新算子，不但进一步优化了去噪性能，而且减小了心电信号特征波的失真。
【技术实现步骤摘要】
本专利技术涉及一种心电信号预处理方法，尤其涉及一种面向远程医疗系统中资源受限心电处理模块中的心电信号预处理。
技术介绍
心脏病已经成为人类健康的头号杀手，目前全世界每年约有1750万人死于心脏病，占全部死亡人数的30%。作为诊断心脏生理、病理状态的重要依据，心电(ECG，Electrocardiogram)诊断技术自其十九世纪初诞生以来，已经发展成为一门专业学科，而对于ECG信号的采集、处理、存储、传输、远程实时监控及自动诊断等方面更是当今该领域的研究重点。ECG是将电极置于体表特定位置所记录到的心脏电活动的图形。典型的心电图如图I所示，主要包括以下几部分，且这几部分为临床心电图的诊断所依据的重要特征参量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信号预处理的高保真度需求。为进一步优化性能，近来又提出结合小波变换与形态滤波的预处理方法。基本思想是，应用小波变换方法滤除高频噪声，应用形态滤波方法校正基线漂移。基于小波变换与形态滤波的预处理方法，这在一定程度上改善了特征波形失真，不能够避免直线型结构元素造成的块状效应。
技术实现思路
为了解决现有技术中问题，本专利技术提供了一种，步骤I :应用直线型结构元素对(B。，B。)的形态滤波器修正基线漂移；步骤2 :应用提升方法构造形态学多分辨率分解心电信号至第二层和第三层；步骤3 :选择第二层与第三层信号分量中QRS复波失真更小的分量作为预处理后的心电信号。作为本专利技术的进一步改进，步骤I中，所述结构元素B。，B。长度的选择依赖于特征波形的持续时间Tw与心电信号的采样频率Fs。作为本专利技术的进一步改进，所述特征波形的持续时间小于O. 2s。作为本专利技术的进一步改进，步骤I中形态滤波器的数学表达如式I所示fb=f0 ο B0 · Bc (I)f = f0-fb (2)其中，ο代表形态学中的开运算， 代表形态学中的闭运算；f。是原始的ECG信号，fb是检测到的基线漂移，f是f。抑制噪声和去除基线漂移后的ECG信号；B。用来移除信号的波峰，Bc用来移除余下的波谷。作为本专利技术的进一步改进，Bc的长度大于B。。作为本专利技术的进一步改进，步骤2中，首先应用最小长度的直线型结构元素对信号进行分解，得到信号分量与细节分量；然后应用预测算子、更新算子得到新的信号分量与细节分量，完成第一层分解；然后应用更大长度的直线型结构元素对第一层信号分量进行分解，得到第二层上的信号分量与细节分量，随后应用预测算子、更新算子得到新的信号分量与细节分量，完成第二层分解；然后应用更大长度的直线型结构元素对第二层信号分量进行分解，得到第三层上的信号分量与细节分量，随后应用预测算子、更新算子得到新的信号分量与细节分量，完成第三层分解。作为本专利技术的进一步改进，令预测算子为该样本与前一个样本的差，更新算子为左右相邻样本的三次样条差值。3/5页作为本专利技术的进一步改进，根据提升小波的构造原理，依据分解、预测及更新三步骤，在尺度j上应用尺度分量h与细节分量Yj构造新的尺度分量X, J与新的细节分量I' j，如式3与4所示j' J=Yj- Ji (Xj)(3)X' J=Xj- λ (y1 j)(4)其中λ Iffj — Vj, IVj — Wj ;则尺度j及尺度j_l上重建公式如式5与6所示Xj = χ) +My'j)￠5) xJA =+ ωι (+ ^(Xj-))KL)其中与 ^为尺度j上的尺度合成算子与细节合成算子。作为本专利技术的进一步改进，令预测算子为该样本与前一个样本的差，更新算子为左右相邻样本的三次样条差值，如式7与8所示Ji (X) (η) =χ(η)-χ(η+1)(7)λ (y) (η) =spline (χ (η_1)，χ (η+1))(8)在此基础上，应用构造尺度分析算子V与细节分析算子 ;如式9和10所示；…,)=x, =MF(Xi)ω) (yj)=y ,-MF1 (λ' )(10)其中，定义形态滤波器MFdf)如式11所示，其中h为长度为j+Ι的直线型结构元Μ/'χ/> 丄(/.。。Bi)MD素；本专利技术应用直线型结构元素对修正基线漂移，相较于基于小波变换的方法，能够在较大程度上减小波形失真。应用形态学多分辨率分解方法滤除肌电干扰等高频噪声，性能优于基于小波变换及形态滤波的预处理方法，且方法简单、易于实现。采用提升方法构造多层形态学分解，由于选择了适合心电信号的预测、更新算子，不但进一步优化了去噪性能，而且减小了心电信号特征波的失真。附图说明图图图号;图图图图图图具体实施例方式下面结合附图对本专利技术做进一步说明。，该方法主要工作过程描述如下1、应用直线型结构元素对(B。，B。)的形态滤波器修正基线漂移，结构元素B。，Bc长度的选择依赖于本文档来自技高网...
【技术保护点】
一种基于形态学多分辨率分解的心电信号预处理方法，其特征在于：步骤1：应用直线型结构元素对(Bo,Bc)的形态滤波器修正基线漂移；步骤2：应用提升方法构造形态学多分辨率分解心电信号至第二层和第三层；步骤3：选择第二层与第三层信号分量中QRS复波失真更小的分量作为预处理后的心电信号。
【技术特征摘要】
【专利技术属性】
技术研发人员：张钦宇，张璞，郑石，
申请(专利权)人：哈尔滨工业大学深圳研究生院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。