CN111513679B - 一种基于形态成分分析及稀疏表示的心电噪声去除方法 - Google Patents
一种基于形态成分分析及稀疏表示的心电噪声去除方法 Download PDFInfo
- Publication number
- CN111513679B CN111513679B CN202010360012.0A CN202010360012A CN111513679B CN 111513679 B CN111513679 B CN 111513679B CN 202010360012 A CN202010360012 A CN 202010360012A CN 111513679 B CN111513679 B CN 111513679B
- Authority
- CN
- China
- Prior art keywords
- data
- signal
- formula
- resonance
- signals
- 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.)
- Active
Links
Images
Classifications
-
- 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/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- 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/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/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/361—Detecting fibrillation
-
- 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
-
- 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/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/023—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the heart
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Power Engineering (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
一种基于形态成分分析及稀疏表示的心电噪声去除方法,通过基于不同成分分析及稀疏表示原理的ECG信号修复技术进行说明,此方法不基于傅里叶变换产生的振幅及频谱。它将信号分解成具有不同特性成分,其中一成分是由多个同时连续振动的信号组成,称作为高共振分量。另一成分是由没有特定状态和持续时间的瞬态冲击信号组成,称作为低共振分量。通过基于形态成分分析及稀疏表示的心电噪声去除方法,包括信号形态成分分析,稀疏表示和群稀疏阈值处理。有效去同频带噪声的干扰。
Description
技术领域
本发明属于心电信号处理技术领域,具体涉及一种基于形态成分分析及稀疏表示的心电噪声去除方法。
背景技术
心电信号是人体重要的生理活动信号,信号的波形是评估患者心脏健康状况的重要依据。ECG信号在收集过程中,由于外部对电子环境的干扰和身体表面电位的变化,收集到的原始信号将与多种形态的噪音混合在一起。由于这些噪音的影响,原ECG信号的波形会被模糊甚至会掩盖波形中的关键特征,无法使用。ECG信号是非线性及非平稳信号,由于噪声对信号有用频率的幅值有调幅作用,并且噪声干扰成分通常存在于整个频率范围内,信号的噪音常常被表现出来。一般以频率为基础的方法,并不能有效的去除噪声瞬态成分在同频带的干扰。
发明内容
本发明为了克服以上技术的不足,提供了一种根据信号政党特性的不同结合形态成分分析,有效去除同频带噪声干扰的基于形态成分分析及稀疏表示的心电噪声去除方法。
本发明克服其技术问题所采用的技术方案是:
一种基于形态成分分析及稀疏表示的心电噪声去除方法,包括如下步骤:
a)加载含有噪声心电数据S;
d)将信号xA *与xB *进行叠加,得到降噪后的心电数据S'。
进一步的,步骤a)中计算机加载心电数据的步骤为:将含有高斯白噪声的心电数据以矩阵S的形式输入,矩阵S为n行1列的数据,通过公式N=length(S)提取矩阵S的全部数据,N为载入的数据的长度。
进一步的,步骤b)中共振稀疏分解步骤如下:
b-1)S=SA+SB,SA为心电信号中由多个同时连续振动的信号组成的高共振成分,SB为心电信号中由没有特定状态和持续时间的多重瞬态冲击信号组成的的低共振成分;
b-2)从心电数据S分离信号SA和SB,通过TWQT可调Q因子小波变换算法利用具有不同品质因子的小波基函数构建高共振信号分量的过完备字典P1以及低共振信号分量的过完备字典P2,通过公式S=P1w1+P2w2表示心电数据,w1为SA在过完备字典P1下的变换系数,w2为SB在过完备字典P2下的变换系数;
b-4)通过公式以及计算得到分解后的高共振信号及低共振信号式中为通过SALSA迭代算法对目标函数J(w1,w2)进行迭代时求得的最小时高共振分量的变换矩阵,为通过SALSA迭代算法对目标函数J(w1,w2)进行迭代时求得的最小时低共振分量的变换矩阵。
进一步的,步骤c)的步骤如下:
c-1)高共振信号由高共振信号中有用的心电数据xA以及高共振信号中的噪声数据zA构成,低共振信号由低共振信号中有用的心电数据xB以及低共振信号中的噪声数据zB构成, N为载入的数据的长度,i为数据标号;
c-3)当K=1时,对函数F(xA)使用MM最小优化算法迭代计算,直至函数收敛时返回数据xA *;c-4)通过公式计算信号xB *,式中R(xB)为惩罚函数,I∈{0,...,N-1},J∈{0,...,K-1},j为第i组的系数指标,λB为正则化参数,λB=0.1;
F(xB)使用MM最小优化算法迭代计算,直至函数收敛时返回数据xB *。
进一步的,步骤d)中通过公式S'=xA *+xB *计算叠加后的心电数据S'。
本发明的有益效果是:通过基于不同成分分析及稀疏表示原理的ECG信号修复技术进行说明,此方法不基于傅里叶变换产生的振幅及频谱。它将信号分解成具有不同特性成分,其中一成分是由多个同时连续振动的信号组成,称作为高共振分量。另一成分是由没有特定状态和持续时间的瞬态冲击信号组成,称作为低共振分量。通过基于形态成分分析及稀疏表示的心电噪声去除方法,包括信号形态成分分析,稀疏表示和群稀疏阈值处理。打破了以传统的以频率划分信号的分解方法。根据信号震荡特性的不同,结合形态成分分析的方法,将具有相同震荡特性且形态不同的成分进行区分,分被设定不同的群稀疏阈值(OGS,Overlapping Group Shrinkage)的处理,有效去同频带噪声的干扰。
附图说明
图1为本发明的方法流程图。
具体实施方式
下面结合附图1对本发明做进一步说明。
一种基于形态成分分析及稀疏表示的心电噪声去除方法,包括如下步骤:
a)加载含有噪声心电数据S;
d)将信号xA *与xB *进行叠加,得到降噪后的心电数据S'。
通过基于不同成分分析及稀疏表示原理的ECG信号修复技术进行说明,此方法不基于傅里叶变换产生的振幅及频谱。它将信号分解成具有不同特性成分,其中一成分是由多个同时连续振动的信号组成,称作为高共振分量。另一成分是由没有特定状态和持续时间的瞬态冲击信号组成,称作为低共振分量。通过基于形态成分分析及稀疏表示的心电噪声去除方法,包括信号形态成分分析,稀疏表示和群稀疏阈值处理。打破了以传统的以频率划分信号的分解方法。根据信号震荡特性的不同,结合形态成分分析的方法,将具有相同震荡特性且形态不同的成分进行区分,分被设定不同的群稀疏阈值(OGS,Overlapping GroupShrinkage)的处理,有效去同频带噪声的干扰。
优选的,步骤a)中计算机加载心电数据的步骤为:将含有高斯白噪声的心电数据以矩阵S的形式输入,矩阵S为n行1列的数据,通过公式N=length(S)提取矩阵S的全部数据,N为载入的数据的长度。
优选的,步骤b)中共振稀疏分解步骤如下:
b-1)S=SA+SB,SA为心电信号中由多个同时连续振动的信号组成的高共振成分,SB为心电信号中由没有特定状态和持续时间的多重瞬态冲击信号组成的的低共振成分;
b-2)从心电数据S分离信号SA和SB,通过TWQT可调Q因子小波变换算法利用具有不同品质因子的小波基函数构建高共振信号分量的过完备字典P1以及低共振信号分量的过完备字典P2,通过公式S=P1w1+P2w2表示心电数据,w1为SA在过完备字典P1下的变换系数,w2为SB在过完备字典P2下的变换系数;
b-4)通过公式以及计算得到分解后的高共振信号及低共振信号式中为通过SALSA迭代算法对目标函数J(w1,w2)进行迭代时求得的最小时高共振分量的变换矩阵,为通过SALSA迭代算法对目标函数J(w1,w2)进行迭代时求得的最小时低共振分量的变换矩阵。
优选的,步骤c)的步骤如下:
c-1)高共振信号由高共振信号中有用的心电数据xA以及高共振信号中的噪声数据zA构成,低共振信号由低共振信号中有用的心电数据xB以及低共振信号中的噪声数据zB构成, N为载入的数据的长度,i为数据标号;
c-3)当K=1时,对函数F(xA)使用MM最小优化算法迭代计算,直至函数收敛时返回数据xA *;c-4)通过公式计算信号xB *,式中R(xB)为惩罚函数,I∈{0,...,N-1},J∈{0,...,K-1},j为第i组的系数指标,λB为正则化参数,λB=0.1;
F(xB)使用MM最小优化算法迭代计算,直至函数收敛时返回数据xB *。
优选的,步骤d)中通过公式S'=xA *+xB *计算叠加后的心电数据S'。
Claims (4)
1.一种基于形态成分分析及稀疏表示的心电噪声去除方法,其特征在于,包括如下步骤:
a)加载含有噪声心电数据S;
d)将信号xA *与xB *进行叠加,得到降噪后的心电数据S';
步骤c)包括如下步骤:
c-1)高共振信号由高共振信号中有用的心电数据xA以及高共振信号中的噪声数据zA构成,低共振信号由低共振信号中有用的心电数据xB以及低共振信号中的噪声数据zB构成, N为载入的数据的长度,i为数据标号;
2.根据权利要求1所述的基于形态成分分析及稀疏表示的心电噪声去除方法,其特征在于,步骤a)中计算机加载心电数据的步骤为:将含有高斯白噪声的心电数据以矩阵S的形式输入,矩阵S为n行1列的数据,通过公式N=length(S)提取矩阵S的全部数据,N为载入的数据的长度。
3.根据权利要求1所述的基于形态成分分析及稀疏表示的心电噪声去除方法,其特征在于,步骤b)中共振稀疏分解步骤如下:
b-1)S=SA+SB,SA为心电信号中由多个同时连续振动的信号组成的高共振成分,SB为心电信号中由没有特定状态和持续时间的多重瞬态冲击信号组成的的低共振成分;
b-2)从心电数据S分离信号SA和SB,通过TWQT可调Q因子小波变换算法利用具有不同品质因子的小波基函数构建高共振信号分量的过完备字典P1以及低共振信号分量的过完备字典P2,通过公式S=P1w1+P2w2表示心电数据,w1为SA在过完备字典P1下的变换系数,w2为SB在过完备字典P2下的变换系数;
4.根据权利要求1所述的基于形态成分分析及稀疏表示的心电噪声去除方法,其特征在于:步骤d)中通过公式S'=xA *+xB *计算叠加后的心电数据S'。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010360012.0A CN111513679B (zh) | 2020-04-29 | 2020-04-29 | 一种基于形态成分分析及稀疏表示的心电噪声去除方法 |
NL2026029A NL2026029B1 (en) | 2020-04-29 | 2020-07-09 | Electrocardio noise removing method based on morphological component analysis and sparse representation |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010360012.0A CN111513679B (zh) | 2020-04-29 | 2020-04-29 | 一种基于形态成分分析及稀疏表示的心电噪声去除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111513679A CN111513679A (zh) | 2020-08-11 |
CN111513679B true CN111513679B (zh) | 2021-06-01 |
Family
ID=71903850
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010360012.0A Active CN111513679B (zh) | 2020-04-29 | 2020-04-29 | 一种基于形态成分分析及稀疏表示的心电噪声去除方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111513679B (zh) |
NL (1) | NL2026029B1 (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0850489A (ja) * | 1994-08-05 | 1996-02-20 | Nissan Motor Co Ltd | 吸音構造体 |
CN105700020A (zh) * | 2016-03-23 | 2016-06-22 | 中国石油天然气集团公司 | 一种地震数据随机噪声压制方法及装置 |
CN108805059A (zh) * | 2018-05-29 | 2018-11-13 | 东华大学 | 稀疏正则化滤波与自适应稀疏分解的齿轮箱故障诊断方法 |
CN110598593A (zh) * | 2019-08-29 | 2019-12-20 | 东南大学 | 基于共振稀疏分解和FastICA算法的行星齿轮箱故障诊断方法 |
-
2020
- 2020-04-29 CN CN202010360012.0A patent/CN111513679B/zh active Active
- 2020-07-09 NL NL2026029A patent/NL2026029B1/en not_active IP Right Cessation
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0850489A (ja) * | 1994-08-05 | 1996-02-20 | Nissan Motor Co Ltd | 吸音構造体 |
CN105700020A (zh) * | 2016-03-23 | 2016-06-22 | 中国石油天然气集团公司 | 一种地震数据随机噪声压制方法及装置 |
CN108805059A (zh) * | 2018-05-29 | 2018-11-13 | 东华大学 | 稀疏正则化滤波与自适应稀疏分解的齿轮箱故障诊断方法 |
CN110598593A (zh) * | 2019-08-29 | 2019-12-20 | 东南大学 | 基于共振稀疏分解和FastICA算法的行星齿轮箱故障诊断方法 |
Non-Patent Citations (4)
Title |
---|
Detection of Myocardial Infarction from Multi-lead ECG using Dual-Q Tunable Q-Factor Wavelet Transform;Liu, Jia et al;《2019 41ST ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY (EMBC) 》;20190727;第1496-1499页 * |
Intelligent bearing fault signature extraction via iterative oscillatory behavior based signal decomposition (IOBSD);Shi, Juanjuan et al;《EXPERT SYSTEMS WITH APPLICATIONS》;20160301;第45卷;第40-45页 * |
基于调Q小波变换的心电信号特征量提取方法;李楠等;《计算机科学》;20141130;第41卷(第11A期);第61-64页 * |
李楠等.基于调Q小波变换的心电信号特征量提取方法.《计算机科学》.2014,第41卷(第11A期), * |
Also Published As
Publication number | Publication date |
---|---|
CN111513679A (zh) | 2020-08-11 |
NL2026029B1 (en) | 2021-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109784242B (zh) | 基于一维残差卷积神经网络的脑电信号去噪方法 | |
CN111616697B (zh) | 一种基于新阈值函数小波变换的心电信号去噪算法 | |
Salman et al. | Novel logarithmic reference free adaptive signal enhancers for ECG analysis of wireless cardiac care monitoring systems | |
Tang et al. | ECG de-noising based on empirical mode decomposition | |
CN113887362A (zh) | 一种局部放电信号的特征提取方法 | |
Goel et al. | An improved wavelet-based signal-denoising architecture with less hardware consumption | |
Kasambe et al. | VLSI wavelet based denoising of PPG signal | |
Butt et al. | Denoising practices for electrocardiographic (ECG) signals: a survey | |
Sharma | EEG signal denoising based on wavelet transform | |
CN112998704A (zh) | 一种可穿戴设备血氧饱和度计算方法 | |
CN111657936A (zh) | 基于小波变换和全变差正则化的信号去噪方法 | |
Yadav et al. | Denoising and SNR improvement of ECG signals using wavelet based techniques | |
Chatterjee et al. | Sparsity-based modified wavelet de-noising autoencoder for ECG signals | |
Lu et al. | Model-based ECG denoising using empirical mode decomposition | |
Hou et al. | Deep neural network denoising model based on sparse representation algorithm for ECG signal | |
CN111513679B (zh) | 一种基于形态成分分析及稀疏表示的心电噪声去除方法 | |
Elshekhidris et al. | WAVELET TRANSFORMS FOR EEG SIGNAL DENOISING AND DECOMPOSITION | |
CN116720056A (zh) | 一种基于增强解码的ae-gan的ecg信号重建方法 | |
Zhou et al. | Denoising and baseline correction of ECG signals using sparse representation | |
CN113598785B (zh) | 基于小波域稀疏特性的心电信号去噪方法 | |
Suchetha et al. | Empirical mode decomposition-based subtraction techniques for 50 Hz interference reduction from electrocardiogram | |
CN114757236A (zh) | 基于tqwt与svmd的脑电信号去噪优化方法及系统 | |
Dai et al. | Denoising ECG by a new wavelet threshold function | |
Malleswari et al. | An improved denoising of electrocardiogram signals based on wavelet thresholding | |
Gon et al. | Removal of Noises from an ECG Signal Using an Adaptive S-Median Thresholding Technique |
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 |