CN108618773A - 一种心电信号去噪方法、装置和一种心电信号采集设备 - Google Patents
一种心电信号去噪方法、装置和一种心电信号采集设备 Download PDFInfo
- Publication number
- CN108618773A CN108618773A CN201710152788.1A CN201710152788A CN108618773A CN 108618773 A CN108618773 A CN 108618773A CN 201710152788 A CN201710152788 A CN 201710152788A CN 108618773 A CN108618773 A CN 108618773A
- Authority
- CN
- China
- Prior art keywords
- signal
- electrocardiosignal
- component
- noise
- intrinsic mode
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Abstract
本发明实施例公开了一种心电信号去噪方法,用于解决现有噪声滤除方法难以完全滤除心电信号中的噪声干扰的问题。本发明实施例方法包括:获取采集的心电信号;对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;对所述各个本征模函数分量进行工频干扰的噪声滤除;对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;对所述各个信号独立分量进行肌电干扰的噪声滤除;对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。本发明实施例还提供一种心电信号去噪装置和一种心电信号采集设备。
Description
技术领域
本发明涉及信号噪声处理技术领域,尤其涉及一种心电信号去噪方法、装置和一种心电信号采集设备。
背景技术
心血管病是一类由心脏、血管和调节血液循环的神经体液组织构成的循环系统的疾病。根据世界卫生组织的统计,心血管疾病是世界上致死率最高的疾病,在2012年全世界有1750万人死于心血管疾病,占各种死因的31%。心血管疾病严重威胁人类,特别是50岁以上中老年人健康,具有高患病率、高致残率和高死亡率的特点。近年来,随着生活压力的增加,工作时间的增加,心血管疾病在青年群体中的发病率越来越高。心电图通过体表间接记录心脏电活动,这些活动反映了心肌在心动周期中的去极化和复极活动,表示通常与体表点击测量的整体心肌的动作电位和时间空间的总和。心电图应用于临床已经有一个世纪的历史,作为一门独立的临床检查诊断学科已自成体系,心电图在心脏疾病的临床诊断中具有重要价值,能为心血管疾病的正确分析、诊断、治疗和监护提供客观指标。心电图检查已经在全世界得到普及,是各级医疗机构一项极为重要的诊断工具,尤其因其是一项无创性检查,价廉方便,能被患者所接受。
随着科学技术的不断发展,心电图的监测方式越来越多样话,除了使用传统的静态心电图机采集心电信号外,运动心电图仪、Holter记录仪、可穿戴心电图仪等设备,由于其更能强大,方便舒适,能够长时间采集心电信号,得到了越来越多的重视。但随着便捷性,舒适度的提升,心电信号在采集过程中更加容易受到噪声的干扰,导致心电信号质量下降。因此,获取高保真的心电信号,提升心电信号质量是心电图的基础,没有良好的信号质量,后续医生或智能算法的分析诊断则无从谈起。
心电信号受到的噪声干扰可以分为工频干扰、基线漂移干扰以及肌电干扰三种形式。其中工频干扰为50Hz正弦信号及其谐波,由空间电磁场引起;基线漂移干扰主要由人体呼吸运动或电极片接触不良导致,属于低频信号,频率范围一般为0-1Hz之间;而肌电干扰是人体肌肉所产生的电信号,可以近似看成高斯白噪声。传统的噪声滤除方法是使用经典滤波器滤除各类噪声,经典滤波器结构简单、易于实现,已经在工业应用上广泛实现。但由于人体心电信号频率范围一般在0-100Hz,该频率处于之前提到的三种噪声的频率范围之内,造成心电信号和噪声频率混叠,因此无法使用经典滤波器完全滤除噪声干扰,若滤波器的截止频率设置过低,则会导致心电信号失真的情况出现。
发明内容
本发明实施例提供了一种心电信号去噪方法、装置和一种心电信号采集设备,能够滤除心电信号中的工频干扰、基线漂移干扰以及肌电干扰的噪声,在保留心电信号原本特征的同时,达到去噪的效果,提升心电信号的质量。
本发明实施例提供的一种心电信号去噪方法,包括:
获取采集的心电信号;
对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;
对所述各个本征模函数分量进行工频干扰的噪声滤除;
对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;
对所述各个信号独立分量进行肌电干扰的噪声滤除;
对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。
可选地,所述对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量具体包括:
将所述心电信号确定为待分解心电信号;
向所述待分解心电信号加入预设的白噪声信号;
对所述待分解心电信号进行经验模态分解,得到各个本征模函数分量;
若所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件,则向所述待分解心电信号加入新的预设白噪声信号,并返回执行对所述待分解心电信号进行经验模态分解的步骤;
若所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件,则对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
可选地,所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件具体为:
所述本征模函数分量的极大值和极小值的数目之和大于或等于预设的数目阈值;
所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件具体为:
所述本征模函数分量的极大值和极小值的数目之和小于预设的数目阈值。
可选地,所述对所述各个信号独立分量进行肌电干扰的噪声滤除具体包括:
滤除所述各个信号独立分量中不包含心电信息的信号;
对滤除信号后的所述各个信号独立分量进行肌电干扰的噪声滤除。
可选地,在滤除所述各个信号独立分量中不包含心电信息的信号之前,还包括:
根据第一公式计算所述信号独立分量的微分信号,所述第一公式为:D(t)=IC(t)-IC(t-1),其中IC为信号独立分量,D为IC的微分信号,t为信号序列的时间点;
根据第二公式对所述微分信号进行非线性滤波,得到特征信号,所述第二公式为其中F(t)为特征信号,k为预设正整数;
统计所述特征信号的值大于预设第一阈值的点数;
若所述点数小于预设的第二阈值,则将与所述特征信号对应的所述信号独立分量确定为不包含有心电信息的信号。
本发明实施例提供的一种心电信号去噪装置,包括:
心电信号获取模块,用于获取采集的心电信号;
模态分解模块,用于对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;
工频干扰滤除模块,用于对所述各个本征模函数分量进行工频干扰的噪声滤除;
成分分析模块,用于对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;
肌电干扰滤除模块,用于对所述各个信号独立分量进行肌电干扰的噪声滤除;
信号重构模块,用于对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。
可选地,所述模态分解模块具体包括:
待分解信号确定单元,用于将所述心电信号确定为待分解心电信号;
白噪声单元,用于向所述待分解心电信号加入预设的白噪声信号;
经验模态分解单元,用于对所述待分解心电信号进行经验模态分解,得到各个本征模函数分量;
返回触发单元,用于若所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件,则向所述待分解心电信号加入新的预设白噪声信号,并返回触发所述经验模态分解单元;
分量均值单元,用于若所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件,则对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
可选地,所述肌电干扰滤除模块具体包括:
第一滤除单元,用于滤除所述各个信号独立分量中不包含心电信息的信号;
第二滤除单元,用于对滤除信号后的所述各个信号独立分量进行肌电干扰的噪声滤除。
可选地,所述心电信号去噪装置还包括:
微分计算模块,用于根据第一公式计算所述信号独立分量的微分信号,所述第一公式为:D(t)=IC(t)-IC(t-1),其中IC为信号独立分量,D为IC的微分信号,t为信号序列的时间点;
非线性滤波模块,用于根据第二公式对所述微分信号进行非线性滤波处理,得到特征信号,所述第二公式为其中F(t)为特征信号,k为预设正整数;
点数统计模块,用于统计所述特征信号的值大于预设第一阈值的点数;
信号确定模块,用于若所述点数小于预设的第二阈值,则将与所述特征信号对应的所述信号独立分量确定为不包含有心电信息的信号。
本发明实施例提供的一种心电信号采集设备,包括上述的心电信号去噪装置。
从以上技术方案可以看出,本发明实施例具有以下优点:
本发明实施例中,首先,获取采集的心电信号;然后,对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;接着,对所述各个本征模函数分量进行工频干扰的噪声滤除;次之,对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;对所述各个信号独立分量进行肌电干扰的噪声滤除;最后,对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。这样可以将心电信号分解为多个分量,在不同分量上进行工频干扰、基线漂移干扰以及肌电干扰的噪声滤除,在保留心电信号原本特征的同时,达到去噪的效果,提升了心电信号的质量。
附图说明
图1为本发明实施例中一种心电信号去噪方法一个实施例流程图;
图2为本发明实施例中一种心电信号去噪方法中步骤102的具体流程示意图;
图3为本发明实施例中一种心电信号去噪方法在一个应用场景下得到的各个本征模函数分量信号示意图;
图4为本发明实施例中一种心电信号去噪方法在一个应用场景下得到的各个信号独立分量的示意图;
图5为本发明实施例中一种心电信号去噪方法中步骤105的具体流程示意图;
图6为本发明实施例中一种心电信号去噪方法在一个应用场景下滤波前和滤波后的心电信号对比示意图;
图7为本发明实施例中一种心电信号去噪装置一个实施例结构图。
具体实施方式
本申请提供了一种基于快速独立成分分析和集合经验模态分解的心电信号去噪方法,该方法可以应用于各种电子设备,这些设备具有心电信号采集、分析模块,用于采集生物心电信号并对其分析。电子设备可以包括但不限于电脑、智能手机、心电图机、可穿戴心电采集设备、监护仪等等。
应用上述方法的电子设备可以通过上述方法进行去噪处理,提升心电信号的信号质量。已知的常用心电信号去噪方法通常是使用了经典滤波器进行处理,该方法通常是在频域上将某频率分量进行滤除,达到去噪的效果。但心电信号的噪声往往与心电信号混叠,使用这种方法无法完全滤除心电噪声;另一方面,在设计滤波器的时候需要非常小心,不恰当的滤波器类型和参数选择,容易改变心电信号的本身特征,影响后续对心电图的分析和诊断。
为了解决上述一个或多个问题,根据本申请实施例的基于快速独立成分分析和集合经验模态分解的心电信号去噪方法,将心电信号分解为多个分量,在不同分量上进行去除或滤波处理处理,在保留心电信号原本特征的同时,达到去噪的效果。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1,本发明实施例中一种心电信号去噪方法一个实施例包括:
101、获取采集的心电信号;
本实施例中,可通过静态心电图机、可穿戴心电采集设备(如智能手环、智能臂环等)、监护仪等设备采集用户的心电信号。一般来说,采集得到的心电信号中包含有工频干扰、基线漂移干扰以及肌电干扰三种形式的噪声。
102、对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;
经验模态分解(EMD)依据数据自身的时间尺度特征来进行信号分解,其本质是对信号进行平稳化处理,将具有不同时间尺度的信号逐级分解开来。EMD已经在不同的工程领域得到了迅速有效的应用,但依然存在不足。当信号的时间尺度存在跳跃性变化时,对信号进行EMD分解,就会出现一个本征模函数分量(IMF)包含不同时间尺度特征成分的情况,称之为模态混叠。
本实施例中,采用集合经验模态分解(EEMD)进行心电信号的分解,其是针对EMD的不足,提出的一种噪声辅助信号处理的方法。
如图2所示,进一步地,上述步骤102具体可以包括:
201、将所述心电信号确定为待分解心电信号;
202、向所述待分解心电信号加入预设的白噪声信号;
203、对所述待分解心电信号进行经验模态分解,得到各个本征模函数分量;
204、判断所述本征模函数分量的极大值和极小值的数目是否满足预设的阈值条件,若否,则执行步骤205,若是,则执行步骤206;
205、向所述待分解心电信号加入新的预设白噪声信号,并返回执行步骤203;
206、对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
对于上述步骤202和步骤205中的“白噪声信号”,这些白噪声信号可以预先设置后,通过随机生成的方式得到。
对于上述步骤203,对待分解心电信号进行EMD分解的过程可以如下:
a)找出加入白噪声信号后的待分解心电信号的所有极大值点Imax,并用三次样条插值函数拟合形成原数据的上包络线;同样,找出待分解心电信号的所有极小值点Imin,并用三次样条插值函数拟合形成原数据的下包络线。
b)计算上下包络线的均值,将上下包络线的均值记为Ml,将待分解心电信号的信号序列x(t)减去平均包络Ml,得到一个新的数据序列IMF,公式如下:
IMF(t)=X(t)-Ml
从而,可以得到各个本征模函数分量IMF。
对于上述步骤204~206,可以为EMD的分解设定终止条件,当达到终止条件时,结束EMD分解,反之,则加入新的白噪声信号,再次进行EMD分解,直到达成终止条件。
本实施例中,设定“所述本征模函数分量的极大值和极小值的数目是否满足预设的阈值条件”作为EMD分解的“终止条件”。当不满足预设的阈值条件时,则向所述待分解心电信号加入新的预设白噪声信号,并返回执行步骤203;当满足预设的阈值条件时,则终止EMD分解,并执行步骤206。
具体地,所述本征模函数分量的极大值和极小值的数目是否满足预设的阈值条件具体为:所述本征模函数分量的极大值和极小值的数目之和是否小于预设的数目阈值。如,定义IMF的极大值与极小值的数目之和如下:
T=Imax+Imin
其中,Imax和Imin分别为IMF的极大值和极小值的数目。
当T大于或等于该数目阈值时,则向该待分解心电信号加入新的白噪声信号,返回进行下一循环的EMD分解。需要说明的是,在每次EMD分解中,均分解得到各个IMF,这些不同循环下的IMF可以暂存至缓存中。
当T小于该数目阈值时,则终止EMD分解,不仅能够在有效滤除心电信号的基线漂移,也能够节省运算资源,避免不必要的计算。
在结束EMD分解之后,执行步骤206,对各次循环的EMD中分解得到的各个IMF进行均值求解。举例说明,下述为三次循环EMD得到的各个IMF:
第一次循环:IMF1(1)、IMF2(1)、IMF3(1)、……、IMF9(1);
第二次循环:IMF1(2)、IMF2(2)、IMF3(2)、……、IMF9(2);
第三次循环:IMF1(3)、IMF2(3)、IMF3(3)、……、IMF9(3);
各个IMF的均值序列为:(IMF1(1)+IMF1(2)+IMF1(3))/3、(IMF2(1)+IMF2(2)+IMF2(3))/3、(IMF3(1)+IMF3(2)+IMF3(3))/3、……、(IMF9(1)+IMF9(2)+IMF9(3))/3,从而,在一个应用场景下,如图3所示,得到所述心电信号的各个IMF分别为:IMF1、IMF2、IMF3、……、IMF9。
由上述内容可知,通过步骤102的EEMD分解,可以将心电信号分解为各个本征模函数分量的同时,滤除其中的基线漂移干扰的噪声。
103、对所述各个本征模函数分量进行工频干扰的噪声滤除;
本实施例中,在各个IMF进行工频滤波的方法可以有多种,可以采用现有技术中移植和未来开发的任何技术来实现,本申请在此方面没有限制。由于心电信号使用EEMD分解成了不同阶数的IMF,将原本尖锐的QRS波群(正常心电图中幅度最大的波群)分解到不同的IMF上,从而降低了分布在不同IMF上QRS的尖锐程度。所以若是在不同IMF上进行工频滤波,比如使用经典的IIR陷波器,可以有效避免由尖锐QRS波群引起的震荡。
104、对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析(FastICA),得到各个信号独立分量;
在上述步骤中,将心电信号中的工频干扰和基线漂移干扰的噪声滤除后,此时的心电信号(暂时被分解为各个IMF)可以看成干净的心电信号与肌电干扰和残余干扰的累加。所以可以使用Fast ICA将相互独立的信号分离,将其中的干扰信号去除,并在其余的信号独立分量(IC)上进行滤波去噪处理。
具体的,使用Fast ICA能够从混合数据中提取出原始的信号独立分量(IC)。该方法能够很好地从观测信号中估计相互统计独立的、被未知因素混合的原始信号。Fast ICA算法有基于四阶累计量、基于似然最大、基于负熵最大等形式。在本实施例中,可以使用任何一种形式进行Fast ICA运算。经过Fast ICA后,可以得到信号独立分量IC和分离矩阵W,IC如图4所示,不同阶数的IC相互独立。
105、对所述各个信号独立分量进行肌电干扰的噪声滤除;
在进行Fast ICA得到各个信号独立分量IC之后,可以对各个IC进行肌电干扰的噪声滤除。如图5所示,具体可以分为以下两步执行:
501、滤除所述各个信号独立分量中不包含心电信息的信号;
502、对滤除信号后的所述各个信号独立分量进行肌电干扰的噪声滤除。
对于上述步骤501,可以对各个IC进行分析,找出其中不含心电信息的信号,将之置为0。具体为,首先找出各个IC中包含心电信息的信号和不包含心电信息的信号。
优选地,可以根据第一公式计算所述信号独立分量IC的微分信号,所述第一公式为:
D(t)=IC(t)-IC(t-1);
其中IC为信号独立分量,D为IC的微分信号,t为信号序列的时间点。
然后,根据第二公式对所述微分信号进行非线性滤波,得到特征信号,所述第二公式为:
其中F(t)为特征信号,k为预设正整数。
在得到F(t)之后,统计所述F(t)的值大于预设第一阈值的点数P,若所述点数P小于预设的第二阈值,则判定该特征信号对应的IC包含有心电信息;若所述点数P小于预设的第二阈值,则判定所述特征信号对应的IC不包含有心电信息。从而可以确定出哪些IC为不包含心电信息的信号,将这些信号置0,也即过滤掉。从而,可以在滤除不同噪声源信号的同时,保留心电信号成份,避免滤波过程噪声心电信号的失真。
对于上述步骤502,对含有心电信息的IC进行滤波处理,这里的滤波可以采用现有技术中移植或未来开发的任何技术来实现,本申请在此方面没有限制。比如在本实施例中,使用了小波阈值去噪法对IC进行滤波。
106、对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。
本实施例中,具体地,可以将上述完成肌电干扰噪声滤除后的各个IC重构为IMF,然后将重构后的IMF再次重构得到去噪后的心电信号。
例如,可以定义经过滤波处理后的信号独立分量成为IC',则可以得到重构后的本征模函数IMF':
IMF'=W-1IC';
接着,重构IMF'可以得到滤波后的心电信号X':
其中n是本征模函数的阶数。
在本实施例中,首先使用集合经验模态分解,将心电信号分解成多阶本征模函数分量,从频域上滤除基线漂移和工频干扰;接着,将信号使用快速独立成份分析分解成相互独立的信号独立分量,从而有效去除肌电干扰带来的噪声源,并滤除心电信号源的噪声。
图6示出了在一个应用场景下滤波前和滤波后的心电信号对比示意图,可见,经过本发明的心电信号去噪方法进行的滤波处理,提升了心电信号质量,在有效滤除心电信号噪声的同时,能够有效的保留心电信号的原本特征。
上面主要描述了一种心电信号去噪方法,下面将对一种心电信号去噪装置进行详细描述。
图7示出了本发明实施例中一种心电信号去噪装置一个实施例结构图。
本实施例中,一种心电信号去噪装置包括:
心电信号获取模块701,用于获取采集的心电信号;
模态分解模块702,用于对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;
工频干扰滤除模块703,用于对所述各个本征模函数分量进行工频干扰的噪声滤除;
成分分析模块704,用于对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;
肌电干扰滤除模块705,用于对所述各个信号独立分量进行肌电干扰的噪声滤除;
信号重构模块706,用于对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。
进一步地,所述模态分解模块具体可以包括:
待分解信号确定单元,用于将所述心电信号确定为待分解心电信号;
白噪声单元,用于向所述待分解心电信号加入预设的白噪声信号;
经验模态分解单元,用于对所述待分解心电信号进行经验模态分解,得到各个本征模函数分量;
返回触发单元,用于若所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件,则向所述待分解心电信号加入新的预设白噪声信号,并返回触发所述经验模态分解单元;
分量均值单元,用于若所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件,则对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
进一步地,所述返回触发单元可以包括:第一触发子单元,用于若所述本征模函数分量的极大值和极小值的数目之和大于或等于预设的数目阈值,则向所述待分解心电信号加入新的预设白噪声信号,并返回触发所述经验模态分解单元;
所述分量均值单元可以包括:第一均值子单元,用于若所述本征模函数分量的极大值和极小值的数目之和小于预设的数目阈值,则对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
进一步地,所述肌电干扰滤除模块具体可以包括:
第一滤除单元,用于滤除所述各个信号独立分量中不包含心电信息的信号;
第二滤除单元,用于对滤除信号后的所述各个信号独立分量进行肌电干扰的噪声滤除。
进一步地,所述心电信号去噪装置还可以包括:
微分计算模块,用于根据第一公式计算所述信号独立分量的微分信号,所述第一公式为:D(t)=IC(t)-IC(t-1),其中IC为信号独立分量,D为IC的微分信号,t为信号序列的时间点;
非线性滤波模块,用于根据第二公式对所述微分信号进行非线性滤波处理,得到特征信号,所述第二公式为其中F(t)为特征信号,k为预设正整数;
点数统计模块,用于统计所述特征信号的值大于预设第一阈值的点数;
信号确定模块,用于若所述点数小于预设的第二阈值,则将与所述特征信号对应的所述信号独立分量确定为不包含有心电信息的信号。
本发明实施例还公开了一种心电信号采集设备,其包括图6对应实施例中描述的任意一种心电信号去噪装置。具体地,该心电信号采集设备可以为静态心电图机、运动心电图仪、Holter记录仪、可穿戴心电图仪、监护仪,等等。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种心电信号去噪方法,其特征在于,包括:
获取采集的心电信号;
对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;
对所述各个本征模函数分量进行工频干扰的噪声滤除;
对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;
对所述各个信号独立分量进行肌电干扰的噪声滤除;
对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。
2.根据权利要求1所述的心电信号去噪方法,其特征在于,所述对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量具体包括:
将所述心电信号确定为待分解心电信号;
向所述待分解心电信号加入预设的白噪声信号;
对所述待分解心电信号进行经验模态分解,得到各个本征模函数分量;
若所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件,则向所述待分解心电信号加入新的预设白噪声信号,并返回执行对所述待分解心电信号进行经验模态分解的步骤;
若所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件,则对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
3.根据权利要求2所述的心电信号去噪方法,其特征在于,所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件具体为:
所述本征模函数分量的极大值和极小值的数目之和大于或等于预设的数目阈值;
所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件具体为:
所述本征模函数分量的极大值和极小值的数目之和小于预设的数目阈值。
4.根据权利要求1至3中任一项所述的心电信号去噪方法,其特征在于,所述对所述各个信号独立分量进行肌电干扰的噪声滤除具体包括:
滤除所述各个信号独立分量中不包含心电信息的信号;
对滤除信号后的所述各个信号独立分量进行肌电干扰的噪声滤除。
5.根据权利要求4所述的心电信号去噪方法,其特征在于,在滤除所述各个信号独立分量中不包含心电信息的信号之前,还包括:
根据第一公式计算所述信号独立分量的微分信号,所述第一公式为:D(t)=IC(t)-IC(t-1),其中IC为信号独立分量,D为IC的微分信号,t为信号序列的时间点;
根据第二公式对所述微分信号进行非线性滤波,得到特征信号,所述第二公式为其中F(t)为特征信号,k为预设正整数;
统计所述特征信号的值大于预设第一阈值的点数;
若所述点数小于预设的第二阈值,则将与所述特征信号对应的所述信号独立分量确定为不包含有心电信息的信号。
6.一种心电信号去噪装置,其特征在于,包括:
心电信号获取模块,用于获取采集的心电信号;
模态分解模块,用于对所述心电信号进行集合经验模态分解,得到滤除基线漂移干扰后的各个本征模函数分量;
工频干扰滤除模块,用于对所述各个本征模函数分量进行工频干扰的噪声滤除;
成分分析模块,用于对滤除工频干扰后的所述各个本征模函数分量进行快速独立成分分析,得到各个信号独立分量;
肌电干扰滤除模块,用于对所述各个信号独立分量进行肌电干扰的噪声滤除;
信号重构模块,用于对滤除肌电干扰后的所述各个信号独立分量进行信号重构,得到去噪后的心电信号。
7.根据权利要求6所述的心电信号去噪装置,其特征在于,所述模态分解模块具体包括:
待分解信号确定单元,用于将所述心电信号确定为待分解心电信号;
白噪声单元,用于向所述待分解心电信号加入预设的白噪声信号;
经验模态分解单元,用于对所述待分解心电信号进行经验模态分解,得到各个本征模函数分量;
返回触发单元,用于若所述本征模函数分量的极大值和极小值的数目不满足预设的阈值条件,则向所述待分解心电信号加入新的预设白噪声信号,并返回触发所述经验模态分解单元;
分量均值单元,用于若所述本征模函数分量的极大值和极小值的数目满足预设的阈值条件,则对各次经验模态分解得到的各个本征模函数分量分别进行均值求解,得到滤除基线漂移干扰后的各个本征模函数分量的均值,并将所述各个本征模函数分量的均值确定为所述心电信号的各个本征模函数分量。
8.根据权利要求6或7所述的心电信号去噪装置,其特征在于,所述肌电干扰滤除模块具体包括:
第一滤除单元,用于滤除所述各个信号独立分量中不包含心电信息的信号;
第二滤除单元,用于对滤除信号后的所述各个信号独立分量进行肌电干扰的噪声滤除。
9.根据权利要求8所述的心电信号去噪装置,其特征在于,所述心电信号去噪装置还包括:
微分计算模块,用于根据第一公式计算所述信号独立分量的微分信号,所述第一公式为:D(t)=IC(t)-IC(t-1),其中IC为信号独立分量,D为IC的微分信号,t为信号序列的时间点;
非线性滤波模块,用于根据第二公式对所述微分信号进行非线性滤波处理,得到特征信号,所述第二公式为其中F(t)为特征信号,k为预设正整数;
点数统计模块,用于统计所述特征信号的值大于预设第一阈值的点数;
信号确定模块,用于若所述点数小于预设的第二阈值,则将与所述特征信号对应的所述信号独立分量确定为不包含有心电信息的信号。
10.一种心电信号采集设备,其特征在于,包括如权利要求6至9中任一项所述的心电信号去噪装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710152788.1A CN108618773A (zh) | 2017-03-15 | 2017-03-15 | 一种心电信号去噪方法、装置和一种心电信号采集设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710152788.1A CN108618773A (zh) | 2017-03-15 | 2017-03-15 | 一种心电信号去噪方法、装置和一种心电信号采集设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108618773A true CN108618773A (zh) | 2018-10-09 |
Family
ID=63686385
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710152788.1A Pending CN108618773A (zh) | 2017-03-15 | 2017-03-15 | 一种心电信号去噪方法、装置和一种心电信号采集设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108618773A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109171706A (zh) * | 2018-09-30 | 2019-01-11 | 南京信息工程大学 | 基于分类匹配和分数阶扩散的心电信号去噪方法和系统 |
CN109620213A (zh) * | 2019-02-25 | 2019-04-16 | 山东大学 | 一种基于多尺度差分特征的心电识别方法及装置 |
CN110051325A (zh) * | 2019-03-29 | 2019-07-26 | 重庆邮电大学 | 基于小波变换及改进eemd的心电信号综合滤波方法 |
CN111466911A (zh) * | 2020-04-13 | 2020-07-31 | 徐州工业职业技术学院 | 一种电子电信号处理系统及方法 |
CN111865727A (zh) * | 2020-08-05 | 2020-10-30 | 江苏电子信息职业学院 | 一种基于物联网的通信延迟降低方法及系统 |
CN111956210A (zh) * | 2020-08-27 | 2020-11-20 | 重庆邮电大学 | 一种基于eemd和信号结构分析的心电信号r波识别方法 |
CN112401906A (zh) * | 2020-11-10 | 2021-02-26 | 河北省科学院应用数学研究所 | 一种基于幅值的自动脑电去伪迹方法 |
CN114386456A (zh) * | 2021-12-15 | 2022-04-22 | 深圳先进技术研究院 | 一种太赫兹时域信号降噪方法、图像重建方法和系统 |
CN117235446A (zh) * | 2023-11-09 | 2023-12-15 | 北京豪能汇新能源科技有限公司 | 一种自动换电站数据智能管理系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101017201A (zh) * | 2007-02-14 | 2007-08-15 | 中国科学院安徽光学精密机械研究所 | 基于经验模态分解的激光雷达信号处理方法 |
CN101869477A (zh) * | 2010-05-14 | 2010-10-27 | 北京工业大学 | 一种自适应脑电信号中眼电伪迹的自动去除方法 |
CN101894097A (zh) * | 2010-07-23 | 2010-11-24 | 哈尔滨工业大学 | 卡尔曼滤波与经验模态分解有机结合的机动目标跟踪方法 |
CN102631195A (zh) * | 2012-04-18 | 2012-08-15 | 太原科技大学 | 人体表面肌电信号单通道盲源分离法 |
CN103106903A (zh) * | 2013-01-11 | 2013-05-15 | 太原科技大学 | 一种单通道盲源分离法 |
CN104636314A (zh) * | 2015-01-30 | 2015-05-20 | 太原科技大学 | 一种单通道混沌信号盲源分离法 |
CN104702244A (zh) * | 2013-12-05 | 2015-06-10 | 中国科学院深圳先进技术研究院 | 基于eemd算法滤除肌电信号检测中工频干扰的自适应滤波器 |
CN104720797A (zh) * | 2015-03-20 | 2015-06-24 | 合肥工业大学 | 一种基于单通道的脑电信号中肌电噪声消除方法 |
CN105030232A (zh) * | 2015-06-30 | 2015-11-11 | 广东工业大学 | 一种心电信号的基线漂移校正方法 |
CN106344005A (zh) * | 2016-10-28 | 2017-01-25 | 张珈绮 | 一种可移动心电图监测系统及监测方法 |
-
2017
- 2017-03-15 CN CN201710152788.1A patent/CN108618773A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101017201A (zh) * | 2007-02-14 | 2007-08-15 | 中国科学院安徽光学精密机械研究所 | 基于经验模态分解的激光雷达信号处理方法 |
CN101869477A (zh) * | 2010-05-14 | 2010-10-27 | 北京工业大学 | 一种自适应脑电信号中眼电伪迹的自动去除方法 |
CN101894097A (zh) * | 2010-07-23 | 2010-11-24 | 哈尔滨工业大学 | 卡尔曼滤波与经验模态分解有机结合的机动目标跟踪方法 |
CN102631195A (zh) * | 2012-04-18 | 2012-08-15 | 太原科技大学 | 人体表面肌电信号单通道盲源分离法 |
CN103106903A (zh) * | 2013-01-11 | 2013-05-15 | 太原科技大学 | 一种单通道盲源分离法 |
CN104702244A (zh) * | 2013-12-05 | 2015-06-10 | 中国科学院深圳先进技术研究院 | 基于eemd算法滤除肌电信号检测中工频干扰的自适应滤波器 |
CN104636314A (zh) * | 2015-01-30 | 2015-05-20 | 太原科技大学 | 一种单通道混沌信号盲源分离法 |
CN104720797A (zh) * | 2015-03-20 | 2015-06-24 | 合肥工业大学 | 一种基于单通道的脑电信号中肌电噪声消除方法 |
CN105030232A (zh) * | 2015-06-30 | 2015-11-11 | 广东工业大学 | 一种心电信号的基线漂移校正方法 |
CN106344005A (zh) * | 2016-10-28 | 2017-01-25 | 张珈绮 | 一种可移动心电图监测系统及监测方法 |
Non-Patent Citations (3)
Title |
---|
KWANG JIN LEE等: "A modified algorithm of the combined ensemble empirical mode decomposition and independent component analysis for the removal of cardiac artifacts from neuromuscular electrical signals", 《PHYSIOLOGICAL MEASUREMENT》 * |
任斌斌: "基于白噪声分离的集合经验模态分解心电信号去噪方法研究", 《生物医学工程学杂志》 * |
马玉润: "ECG预处理与QRS波群检测技术研究", 《中国优秀硕士学位论文全文数据库(电子期刊)信息科技辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109171706A (zh) * | 2018-09-30 | 2019-01-11 | 南京信息工程大学 | 基于分类匹配和分数阶扩散的心电信号去噪方法和系统 |
CN109171706B (zh) * | 2018-09-30 | 2020-12-29 | 南京信息工程大学 | 基于分类匹配和分数阶扩散的心电信号去噪方法和系统 |
CN109620213A (zh) * | 2019-02-25 | 2019-04-16 | 山东大学 | 一种基于多尺度差分特征的心电识别方法及装置 |
CN110051325A (zh) * | 2019-03-29 | 2019-07-26 | 重庆邮电大学 | 基于小波变换及改进eemd的心电信号综合滤波方法 |
CN111466911A (zh) * | 2020-04-13 | 2020-07-31 | 徐州工业职业技术学院 | 一种电子电信号处理系统及方法 |
CN111865727A (zh) * | 2020-08-05 | 2020-10-30 | 江苏电子信息职业学院 | 一种基于物联网的通信延迟降低方法及系统 |
CN111956210A (zh) * | 2020-08-27 | 2020-11-20 | 重庆邮电大学 | 一种基于eemd和信号结构分析的心电信号r波识别方法 |
CN111956210B (zh) * | 2020-08-27 | 2022-06-03 | 重庆邮电大学 | 一种基于eemd和信号结构分析的心电信号r波识别方法 |
CN112401906A (zh) * | 2020-11-10 | 2021-02-26 | 河北省科学院应用数学研究所 | 一种基于幅值的自动脑电去伪迹方法 |
CN114386456A (zh) * | 2021-12-15 | 2022-04-22 | 深圳先进技术研究院 | 一种太赫兹时域信号降噪方法、图像重建方法和系统 |
CN117235446A (zh) * | 2023-11-09 | 2023-12-15 | 北京豪能汇新能源科技有限公司 | 一种自动换电站数据智能管理系统 |
CN117235446B (zh) * | 2023-11-09 | 2024-01-26 | 北京豪能汇新能源科技有限公司 | 一种自动换电站数据智能管理系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108618773A (zh) | 一种心电信号去噪方法、装置和一种心电信号采集设备 | |
Chiang et al. | Noise reduction in ECG signals using fully convolutional denoising autoencoders | |
CN106214145B (zh) | 一种基于深度学习算法的心电图分类方法 | |
Castro et al. | ECG feature extraction using optimal mother wavelet | |
García et al. | A novel wavelet-based filtering strategy to remove powerline interference from electrocardiograms with atrial fibrillation | |
Soerensen et al. | A comparison of IIR and wavelet filtering for noise reduction of the ECG | |
CN111265210A (zh) | 一种基于深度学习的房颤预测装置和设备 | |
Kumar et al. | Performance analysis of different filters for power line interface reduction in ECG signal | |
CN110720894B (zh) | 一种心房扑动的检测方法、装置、设备及存储介质 | |
Di Marco et al. | Evaluation of an algorithm based on single-condition decision rules for binary classification of 12-lead ambulatory ECG recording quality | |
Li et al. | ECG denoising method based on an improved VMD algorithm | |
Malghan et al. | A review on ECG filtering techniques for rhythm analysis | |
Bortolan et al. | Dynamic filtration of high-frequency noise in ECG signal | |
CN112568886A (zh) | 心律的检测方法、电子设备以及计算机可读存储介质 | |
Mortezaee et al. | An improved SSA-based technique for EMG removal from ECG | |
CN113384277B (zh) | 心电数据的分类方法及分类系统 | |
Bajaj et al. | A robust approach to denoise ECG signals based on fractional Stockwell transform | |
CN108154112A (zh) | 一种处理心电数据的方法、处理心电数据的装置及电子设备 | |
Chieng et al. | The study and comparison between various digital filters for ECG de-noising | |
Li et al. | A method for automatic removal of EOG artifacts from EEG based on ICA-EMD | |
CN110558974B (zh) | 一种基于极值能量分解法的心电图信号分析方法 | |
Bhoraniya et al. | Motion artifacts extraction using DWT from ambulatory ECG (A-ECG) | |
Gupta et al. | Exploiting tunable Q-factor wavelet transform domain sparsity to denoise wrist PPG signals | |
Kurniawan et al. | Arrhythmia Classification on Electrocardiogram Signal Using Convolution Neural Network Based on Frequency Spectrum | |
Jokić et al. | An efficient ECG modeling for heartbeat classification |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181009 |
|
RJ01 | Rejection of invention patent application after publication |