CN109567788A - 一种去除振铃的心电信号滤波方法 - Google Patents
一种去除振铃的心电信号滤波方法 Download PDFInfo
- Publication number
- CN109567788A CN109567788A CN201811443998.7A CN201811443998A CN109567788A CN 109567788 A CN109567788 A CN 109567788A CN 201811443998 A CN201811443998 A CN 201811443998A CN 109567788 A CN109567788 A CN 109567788A
- Authority
- CN
- China
- Prior art keywords
- signal
- filtering
- mutation
- section
- replacement
- 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.)
- Granted
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
-
- 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]
- A61B5/319—Circuits for simulating ECG signals
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开一种去除振铃的心电信号滤波方法,包括以下步骤:步骤S1、采集原始心电信号,并识别所述原始心电信号中的突变段;步骤S2、截取所述突变段得到突变信号,通过替换段替换所述原始心电信号中的突变段,得到替换信号;步骤S3、分别对所述突变信号以及替换信号进行工频滤波,得到突变滤波信号和替换滤波信号;步骤S4、拼接所述突变滤波信号和替换滤波信号,得到完整的心电滤波信号。本发明可有效去除心电滤波过程中的振铃现象。
Description
技术领域
本发明涉及心电信号滤波技术领域,具体涉及一种去除振铃的心电信号滤波方法。
背景技术
心电信号相对于电网的工频干扰来说是极其微弱的。滤除工频干扰是心电信号处理中需要解决的一个重要问题。目前,数字滤波技术已在工频干扰的处理上得到较多的应用,但现有的数字陷波器在滤除工频干扰时往往会产生“振铃”现象,从而影响数据的分析处理的准确性。心电信号发生了突变,经过工频滤波器,如50Hz或者60Hz的陷波器,在突变信号后会产生来回震荡的信号现象,在心电信号处理中称之为振铃现象。该种突变信号包括QRS波群高尖、起搏器信号、高大T波或由外界干扰引入的异常高尖信号。特别的,在儿童心电图中,由于儿童心电的特有属性,幅值高,信号上升沿短,振铃现象在儿童心电图十分常见,该振铃现象干扰的医师的心电图判读,对临床上的应用极为不利。
发明内容
本发明的目的在于克服上述技术不足,提供一种去除振铃的心电信号滤波方法,解决现有技术中滤除工频干扰时往往会产生振铃现象的技术问题。
为达到上述技术目的,本发明的技术方案提供一种去除振铃的心电信号滤波方法,包括以下步骤:
步骤S1、采集原始心电信号,并识别所述原始心电信号中的突变段;
步骤S2、截取所述突变段得到突变信号,通过替换段替换所述原始心电信号中的突变段,得到替换信号;
步骤S3、分别对所述突变信号以及替换信号进行工频滤波,得到突变滤波信号和替换滤波信号;
步骤S4、拼接所述突变滤波信号和替换滤波信号,得到完整的心电滤波信号。
与现有技术相比,本发明的有益效果包括:本发明首先识别并截取突变信号,通过替换段替换原始信号中的突变段,使得突变段与非突变段分割开,然后分别对突变信号和替换信号进行工频滤波,最后再对滤波后的信号进行拼接得到完整的滤波信号,由于对原始心电信号中的突变段进行了分割与隔离,因此避免了滤波过程中因突变段而产生的振铃现象,去除心电滤波信号中的振铃,提高心电滤波信号的准确性。
附图说明
图1是本发明提供的去除振铃的心电信号滤波方法的流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例1:
如图1所示,本发明的实施例1提供了一种去除振铃的心电信号滤波方法,包括以下步骤:
步骤S1、采集原始心电信号,并识别所述原始心电信号中的突变段;
步骤S2、截取所述突变段得到突变信号,通过替换段替换所述原始心电信号中的突变段,得到替换信号;
步骤S3、分别对所述突变信号以及替换信号进行工频滤波,得到突变滤波信号和替换滤波信号;
步骤S4、拼接所述突变滤波信号和替换滤波信号,得到完整的心电滤波信号。
通过心电图机或者其他心电采集设备采集原始心电信号sig_ori,本实施例中原始心电信号sig_ori的采样率为1000Hz。对原始心电信号sig_ori中突变段进行识别,通过对原始心电信号sig_ori进行检测分析,找出其信号段中的突变段sig_1;:识别出信号的突变段sig_1后,进行信号的分割处理,截取出突变段sig_1,同时,通过替换段替换原始心电信号中的突变段,得到替换信号sig_2。分别对信号sig_1、sig_2进行工频滤波,得到去工频后的信号sig_1’、sig_2’。拼接滤波后信号sig_1’、sig_2’,得到完整的心电滤波信号。
本发明通过分割、分别工频滤波,再拼接的方法,有效消除了心电信号滤波过程中的振铃现象,消除了滤波干扰,提高了心电滤波信号的准确度。本发明可用于QRS复合波、起搏信号、高大T波通过陷波器引起振铃的消除。
优选的,所述步骤S1具体包括:
步骤S11、对所述原始心电信号进行差分计算,得到差分信号;
步骤S12、根据所述差分信号统计所述原始心电信号的上升沿的长度值以及下降沿的长度值;
步骤S13、筛选出上升沿的长度值小于设定长度阈值的信号段以及下降沿的长度值小于所述设定长度阈值的信号段,得到突变段候选区;
步骤S14、筛选出所述突变段候选区中信号幅值大于设定幅值阈值的信号段,得到所述突变段。
在对原始心电信号sig_ori中识别出突变信号段sig_1过程中,对原始心电信号sig_ori中相邻两点进行差分计算,得到差分信号df_sig_ori。通过对差分信号df_sig_ori中每个点进行正负数判断,统计差分信号df_sig_ori中信号上升沿与下降沿的长度值width[i]。计算各个候选区的长度值,若第i个候选区D{i}的长度大于或者等于设定长度阈值,本实施例中为80ms,时间间隔为2ms,即40个采样点,则将该信号段舍去,该信号段不作为突变信号候选段。经过长度筛选后,剩下的突变段候选区域D{m}中,对每个候选区域的point[i]点的幅值进行判断,point[i]幅值小于设定幅度阈值时,该信号段不作为突变段,否则可判定其为突变段。
优选的,所述步骤S12具体为:
步骤S121、遍历所述差分信号,找出上升沿与下降沿的连接点;
步骤S122、以所述连接点为起点分别向前和向后等时间间隔的搜索信号点,直到搜索到的信号点的差分信号符号与前一个信号点的差分信号符号相反,则停止搜索;
步骤S123、根据搜索到的信号点的个数计算上升沿的长度值或下降沿的长度值。
当差分信号df_sig_ori中任意一信号点i的差分信号大小df_sig_ori[i]>0时,表示该点处于信号的上升沿处;当df_sig_ori[i]<0时,表示该点处于信号的下降沿。因此可以通过差分信号符号,即差分信号的正负,判断信号点处于上升沿还是下降沿,从而找到连续的处于上升沿/下降沿的信号点,进而计算出上升沿和下降沿的长度。
具体的,遍历整个信号段的差分信号df_sig_ori,找出上升沿与下降沿连接点point[i],分别从point[i]处往前或往后依次寻找信号点,不妨设定:point[i]前一个点point[i-1]为上升沿点,则继续往前搜索,直至碰到非上升沿点停止;point[i]后一个点point[i+1]为下降沿点,则继续往后搜索,直至碰到非下降沿点停止。按照此方法得到信号sig_ori中一系列的异常信号候选区域D{m},m=1,2,3,...,N。其中,每个候选区域中保存的是信号sig_ori中point[i-j]到point[i+k]的下标。由于信号点是等时间间隔的,因此,可以根据下标的差值与时间间隔的乘积达到上升沿/下降沿的长度值。
优选的,所述设定长度阈值在20ms-40ms之间,所述设定幅值阈值在0.1mv-0.5mv之间。
从大量实验中得知,经过工频滤波后能够产生振铃现象的信号具有这样的特点:该信号段的上升沿/下降沿一般不超过40ms,信号幅值一般不低于0.1mv。因此,将设定长度阈值设置在20ms-40ms之间较为合适,将设定幅值阈值设置在0.1mv-0.5mv之间较为合适。
优选的,所述步骤S122还包括,搜索到的差分信号符号相同的信号点的个数达到设定阈值时,停止搜索,并判定信号段不是突变段。
当差分信号符号相同的信号点的个数达到设定阈值时,说明该信号段的尺度已经超过设定长度阈值,因此排除该信号段为突变段的可能,停止搜索,不用再等到搜索到差分信号符号相反的信号点再停止搜索,节约搜索时间,提高搜索效率。设定阈值根据设定长度阈值设置,例如设定长度阈值为40ms,信号点的时间间隔为1ms,则设定阈值为40。
优选的,所述替换段具体为:与所述突变段长度一致的直线信号。
利用与突变信号长度一致的直线,首尾连接原始心电信号sig_ori中突变段sig_1的起点与终点,得到替换信号。替换段还可以采用其他类型的信号,本实施例采用直线信号,直线信号便于识别和区分。
优选的,所述步骤S4具体为:将所述突变滤波信号拼接至所述替换滤波信号的对应位置,得到所述心电滤波信号。
对步骤S3中分别进行工频滤波后的信号sig_1’、sig_2’进行拼接,具体的,本实施例中进行工频滤波的工频滤波器为FIR滤波器,具有的线性的时延,通过延时的调整,将sig_1’的位置对应回sig_2’,得到工频滤波后完整的心电滤波信号。
实施例2:
本发明的实施例2提供了一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时,实现以上任一实施例所述的去除振铃的心电信号滤波方法。
本发明提供的计算机存储介质,用于实现上述去除振铃的心电信号滤波方法,因此,上述去除振铃的心电信号滤波方法所具备的技术效果,计算机存储介质同样具备,在此不再赘述。
以上所述本发明的具体实施方式,并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形,均应包含在本发明权利要求的保护范围内。
Claims (8)
1.一种去除振铃的心电信号滤波方法,其特征在于,包括以下步骤:
步骤S1、采集原始心电信号,并识别所述原始心电信号中的突变段;
步骤S2、截取所述突变段得到突变信号,通过替换段替换所述原始心电信号中的突变段,得到替换信号;
步骤S3、分别对所述突变信号以及替换信号进行工频滤波,得到突变滤波信号和替换滤波信号;
步骤S4、拼接所述突变滤波信号和替换滤波信号,得到完整的心电滤波信号。
2.根据权利要求1所述的去除振铃的心电信号滤波方法,其特征在于,所述步骤S1具体包括:
步骤S11、对所述原始心电信号进行差分计算,得到差分信号;
步骤S12、根据所述差分信号统计所述原始心电信号的上升沿的长度值以及下降沿的长度值;
步骤S13、筛选出上升沿的长度值小于设定长度阈值的信号段以及下降沿的长度值小于所述设定长度阈值的信号段,得到突变段候选区;
步骤S14、筛选出所述突变段候选区中信号幅值大于设定幅值阈值的信号段,得到所述突变段。
3.根据权利要求2所述的去除振铃的心电信号滤波方法,其特征在于,所述步骤S12具体为:
步骤S121、遍历所述差分信号,找出上升沿与下降沿的连接点;
步骤S122、以所述连接点为起点分别向前和向后等时间间隔的搜索信号点,直到搜索到的信号点的差分信号符号与前一个信号点的差分信号符号相反,则停止搜索;
步骤S123、根据搜索到的信号点的个数计算上升沿的长度值或下降沿的长度值。
4.根据权利要求3所述的去除振铃的心电信号滤波方法,其特征在于,所述步骤S122还包括,搜索到的差分信号符号相同的信号点的个数达到设定阈值时,停止搜索,并判定信号段不是突变段。
5.根据权利要求2所述的去除振铃的心电信号滤波方法,其特征在于,所述设定长度阈值在20ms-40ms之间,所述设定幅值阈值在0.1mv-0.5mv之间。
6.根据权利要求1所述的去除振铃的心电信号滤波方法,其特征在于,所述替换段具体为:与所述突变段长度一致的直线信号。
7.根据权利要求1所述的去除振铃的心电信号滤波方法,其特征在于,所述步骤S4具体为:将所述突变滤波信号拼接至所述替换滤波信号的对应位置,得到所述心电滤波信号。
8.一种计算机存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时,实现如权利要求1-5任一所述的去除振铃的心电信号滤波方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811443998.7A CN109567788B (zh) | 2018-11-29 | 2018-11-29 | 一种去除振铃的心电信号滤波方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811443998.7A CN109567788B (zh) | 2018-11-29 | 2018-11-29 | 一种去除振铃的心电信号滤波方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109567788A true CN109567788A (zh) | 2019-04-05 |
CN109567788B CN109567788B (zh) | 2021-08-20 |
Family
ID=65925663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811443998.7A Active CN109567788B (zh) | 2018-11-29 | 2018-11-29 | 一种去除振铃的心电信号滤波方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109567788B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0751236A (ja) * | 1993-07-16 | 1995-02-28 | Siemens Elema Ab | フィルタ処理されたecg信号のリンギングを除去する装置 |
CN101808054A (zh) * | 2010-03-26 | 2010-08-18 | 北京天碁科技有限公司 | 信道估计的实现方法和装置 |
CN101953684A (zh) * | 2009-07-17 | 2011-01-26 | 深圳市理邦精密仪器股份有限公司 | 一种消除振铃的工频干扰滤波方法 |
US20110037514A1 (en) * | 2009-08-17 | 2011-02-17 | Edan Instruments, Inc. | Method for the Elimination of Ringing of Power Line Interference Filter |
CN104729591A (zh) * | 2015-01-16 | 2015-06-24 | 合肥工业大学 | 一种基于数据替换的涡街流量计抗低频强瞬态冲击振动的信号处理方法 |
CN104796109A (zh) * | 2015-04-09 | 2015-07-22 | 北京楚捷科技有限公司 | 一种滤波方法及装置 |
CN106037655A (zh) * | 2016-06-17 | 2016-10-26 | 美的集团股份有限公司 | 一种睡眠周期信号的分离及周期计算方法及装置 |
US20170172443A1 (en) * | 2014-02-25 | 2017-06-22 | Centre Hospitalier Regional Universitaire De Lille | Method, device, system and computer programme for filtering an rr series obtained from a cardiac signal with automatic checking of the quality of the rr series |
CN106889984A (zh) * | 2017-01-22 | 2017-06-27 | 河北大学 | 一种心电信号自动降噪方法 |
CN107890347A (zh) * | 2017-09-08 | 2018-04-10 | 曼森伯格(深圳)科技发展有限公司 | 滤除工频干扰信号的控制方法及系统 |
CN108113664A (zh) * | 2017-12-07 | 2018-06-05 | 广州视源电子科技股份有限公司 | 心电信号获取过程中抑制陷波滤波器振铃的方法及系统 |
CN108732558A (zh) * | 2018-03-26 | 2018-11-02 | 东南大学 | 一种基于拷贝信号分段的匹配滤波计算方法 |
-
2018
- 2018-11-29 CN CN201811443998.7A patent/CN109567788B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0751236A (ja) * | 1993-07-16 | 1995-02-28 | Siemens Elema Ab | フィルタ処理されたecg信号のリンギングを除去する装置 |
CN101953684A (zh) * | 2009-07-17 | 2011-01-26 | 深圳市理邦精密仪器股份有限公司 | 一种消除振铃的工频干扰滤波方法 |
US20110037514A1 (en) * | 2009-08-17 | 2011-02-17 | Edan Instruments, Inc. | Method for the Elimination of Ringing of Power Line Interference Filter |
CN101808054A (zh) * | 2010-03-26 | 2010-08-18 | 北京天碁科技有限公司 | 信道估计的实现方法和装置 |
US20170172443A1 (en) * | 2014-02-25 | 2017-06-22 | Centre Hospitalier Regional Universitaire De Lille | Method, device, system and computer programme for filtering an rr series obtained from a cardiac signal with automatic checking of the quality of the rr series |
CN104729591A (zh) * | 2015-01-16 | 2015-06-24 | 合肥工业大学 | 一种基于数据替换的涡街流量计抗低频强瞬态冲击振动的信号处理方法 |
CN104796109A (zh) * | 2015-04-09 | 2015-07-22 | 北京楚捷科技有限公司 | 一种滤波方法及装置 |
CN106037655A (zh) * | 2016-06-17 | 2016-10-26 | 美的集团股份有限公司 | 一种睡眠周期信号的分离及周期计算方法及装置 |
CN106889984A (zh) * | 2017-01-22 | 2017-06-27 | 河北大学 | 一种心电信号自动降噪方法 |
CN107890347A (zh) * | 2017-09-08 | 2018-04-10 | 曼森伯格(深圳)科技发展有限公司 | 滤除工频干扰信号的控制方法及系统 |
CN108113664A (zh) * | 2017-12-07 | 2018-06-05 | 广州视源电子科技股份有限公司 | 心电信号获取过程中抑制陷波滤波器振铃的方法及系统 |
CN108732558A (zh) * | 2018-03-26 | 2018-11-02 | 东南大学 | 一种基于拷贝信号分段的匹配滤波计算方法 |
Non-Patent Citations (1)
Title |
---|
刘忠 等: "心室晚电位检测中消除振铃方法的研究", 《北京生物医学工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109567788B (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101449973B (zh) | 用于心电干扰信号识别的判断指标的生成方法及装置 | |
CN108814590B (zh) | 一种心电qrs波群的检测方法及其心电分析方法 | |
WO2012087702A1 (en) | Adaptive ecg wandering correction | |
CN103185837B (zh) | 一种电力系统频率测量的方法 | |
FI107081B (fi) | Menetelmä ja järjestely osittaispurkauslähteiden lukumäärän selvittämiseksi | |
CN104095627B (zh) | 一种心电图数字化信号质量软判决方法和装置 | |
CN108324271A (zh) | 心电信号识别方法、系统和心电监测设备 | |
Sornmo et al. | Adaptive QRS detection: A study of performance | |
CN110596458B (zh) | Demon谱谐波线谱和基频自动估计方法 | |
CN104382589B (zh) | 基于部分按段重采样的胎儿心电图分离提取方法 | |
CN109567788A (zh) | 一种去除振铃的心电信号滤波方法 | |
US8108164B2 (en) | Determination of a common fundamental frequency of harmonic signals | |
CN111948454B (zh) | 一种抗衰减直流分量干扰的同步相量测量方法及系统 | |
Birkemeier et al. | Pattern recognition techniques for the detection of epileptic transients in EEG | |
CN105262729A (zh) | 木马检测方法及系统 | |
CN103795473B (zh) | 消除工频干扰的方法和系统 | |
CN110179456A (zh) | 心电噪声识别模型训练及心电噪声检测方法、装置 | |
CN103746669A (zh) | 一种切比雪夫低通滤波器的设计方法 | |
Zaitsev | Elimination on power line interference from ECG signal using combined bidirectional narrow-band notch filter | |
CN108490311B (zh) | 基于工频采样的弱冲击信号提取与分离方法 | |
WO2014018662A1 (en) | Method of extracting zero crossing data from full spectrum signals | |
CN110507299A (zh) | 一种心率信号检测装置及方法 | |
CN115754832A (zh) | 变压器绕组振荡波高压开关暂态噪声消除方法及装置 | |
CN111208451A (zh) | 心电导联脱落检测电路、方法和医疗监护设备 | |
CN111693685A (zh) | 血栓弹力信号提取方法、存储介质及计算机设备 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |