CN106419898A - Method removing electrocardiosignal baseline drift - Google Patents

Method removing electrocardiosignal baseline drift Download PDF

Info

Publication number
CN106419898A
CN106419898A CN201610675959.4A CN201610675959A CN106419898A CN 106419898 A CN106419898 A CN 106419898A CN 201610675959 A CN201610675959 A CN 201610675959A CN 106419898 A CN106419898 A CN 106419898A
Authority
CN
China
Prior art keywords
wavelet
baseline drift
signal
electrocardiosignal
frequency coefficient
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
Application number
CN201610675959.4A
Other languages
Chinese (zh)
Inventor
罗伟
朱国永
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
WUHAN ZONCARE BIO-MEDICAL ELECTRONICS Co Ltd
Original Assignee
WUHAN ZONCARE BIO-MEDICAL ELECTRONICS Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by WUHAN ZONCARE BIO-MEDICAL ELECTRONICS Co Ltd filed Critical WUHAN ZONCARE BIO-MEDICAL ELECTRONICS Co Ltd
Publication of CN106419898A publication Critical patent/CN106419898A/en
Pending legal-status Critical Current

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms

Abstract

The invention discloses a method removing electrocardiosignal baseline drift; the method comprises the following steps: 1, determining a wavelet decomposition layer: using a signal sampling rate and an electrocardiosignal effective frequency range to determine the wavelet decomposition layer; 2, decomposing wavelet: using wavelet analysis to carry out J layer decomposition for the electrocardiosignal according to the wavelet decomposition layer determined in the step1, wherein the wavelet high frequency coefficients of each layer are respectively CD1, CD2,... CDJ, and the lower frequency coefficient is CAJ; 3, wavelet coefficient processing and signal reconstruction: respectively processing high frequency coefficients and low frequency coefficient in the step2, and carrying out wavelet reconstruction according to the processed wavelet coefficients; 4, moving average filtering step: carrying out moving average filtering on reconstructed signal in the step 3, thus obtaining an electrocardiosignal baseline drift value; 5, removing baseline drift: using an original electrocardiosignal signal to subtract the baseline drift value in the step 4, thus obtaining the baseline drift removed signal.

Description

A kind of method removing ECG baseline drift
Technical field
The present invention relates to biomedicine signals noise management technique field, particularly relate to wavelet analysis combined window and move flat The method that all ECG baseline drift is removed in filtering.
Background technology
Electrocardiosignal (Electrocardiogram, ECG) is that the mankind study the earliest and are applied to clinical medical biological electricity One of signal, it reflects heart and produces in excitement, conduct and the Electrical change in recovery process, is a kind of visitor of cardiac electrical activity Seeing and representing, reflecting the duty of heart in different aspects, clinical diagnosis and treatment for heart disease have very Important reference value.Electrocardiosignal is a kind of faint bioelectrical signals, and its amplitude is millivolt (mV) level, it is easy to by outward The impact of boundary's environment.During to electrocardiographic recording, the activity of the person's of being recorded health, breathing, various communication apparatus etc. all can Producing bigger interference to recording process, these interference have a great impact for the correctness of signal detection.
Electrocardiosignal is mainly disturbed by following three types noise:
1 industrial frequency noise:Industrial frequency noise is a kind of interference being produced by electrified system, and the form presenting is similar at pure electrocardio Sinusoidal wave and similar sinusoidal wave waveform occurs on signal.Such noise frequency is generally 50Hz or 60Hz, and its Amplitude is not high, is about the 50% of this electro-cardiologic signal waveforms high-amplitude.
2 myoelectricity noises:Myoelectricity noise is owing to the movement of muscular tone and the person causes, the frequency of this noise like For 0-10kHz, similar to average Gaussian noise, so not being typically it is obvious that this noise like typically shows fainter, can Myoelectricity interference is regarded as the zero-mean band-limited noise of moment generation.
3 baseline drift noises:Baseline drift interference in ECG is common situation in ECG's data compression, by being remembered Causing when record person breathes, having no idea to avoid, the noise amplitude ratio of baseline drift is relatively low, and baseline drift frequency arrives at 0.05Hz Between 2Hz, this curve ratio is shallower simultaneously, with sinusoidal waveform is similar, is difficult to remove.
The noise main method removing three of the above type has:The method removing industrial frequency noise is relatively more, mainly includes:Flat Sliding filtering, IIR filtering, FIR filtering, adaptive-filtering etc., the way of comparative maturity is to utilize notch filter.Every kind of method all has Its respective good and bad point, selects suitable method according to actual application demand in engineer applied;Remove myoelectricity noise, for flesh Electrical noise typically can take the low pass filter taking upper limiting frequency to be 100Hz on high-pass filtering combined with hardware to eliminate; Remove baseline noise, typically have the methods such as traditional high-pass filtering and wavelet transformation, due to needle position misalignment noise frequency ratio relatively Low (generally 0-0.05Hz), filter order required during design high-pass filter is big especially, computationally intensive, takies calculating money Source is many;Equally, when the method utilizing wavelet analysis is carried out, owing to it is based on convolution algorithm, will be by when removing baseline drift After signal decomposition to ten layers (in the case that particularly signal sampling rate is big), computationally intensive, high to calculating resource requirement, unfavorable Process in real-time.
Content of the invention
The technical problem to be solved is to provide a kind of method removing ECG baseline drift, is used for overcoming The technical problem that prior art exists.
It is as follows that the present invention solves the technical scheme that above-mentioned technical problem taked:
A kind of method removing ECG baseline drift, including:
Step 1) determine wavelet decomposition number of plies step, including:By the effective band of the sample rate of signal and electrocardiosignal, Determine the number of plies of wavelet decomposition;
Step 2) carry out wavelet decomposition step, including:According to the wavelet decomposition number of plies determining in step 1, sampling small echo Analyzing and carrying out J layer decomposition to electrocardiosignal, the small echo high frequency coefficient of every layer is respectively:CD1, CD2 ..., CDJ, low frequency coefficient For:CAJ;
Step 3) wavelet coefficient process and signal reconstruction step, including:Respectively by the high frequency coefficient in step 2 and low frequency system Number is processed, and carries out wavelet reconstruction according to the wavelet coefficient after processing;
Step 4) moving average filtering step, including:Average filter is moved to the signal of reconstruct in step 3, obtains The baseline drift value of electrocardiosignal;
Step 5) remove baseline drift step, including:Deducted the needle position misalignment value in step 4 by original electro-cardiologic signals, I.e. obtain the signal after removing baseline drift.
Preferably, step 1) in, when determining the wavelet decomposition number of plies, carry out by following formula:
J=log2(Fs/f)
Wherein:J is the wavelet decomposition number of plies, and Fs is the sample rate of electrocardiosignal, and f is that the effective band of electrocardiosignal is maximum Value.
Preferably, in step 2 and step 3, the wavelet function carrying out wavelet function feedback employing is Haar function.
Preferably, in step 3, process to by the high frequency coefficient in step 2 and low frequency coefficient, including:
Each for small echo floor height frequency filter factor is set to 0 entirely, i.e.:CA1 '=CA2 '=... ...=CAJ '=0;
Low frequency coefficient keeps constant, CAJ '=CAJ.
Preferably, step 3) in, in signal reconstruction step, take Mallat algorithm.
Preferably, in moving average filtering step, specifically include:
Wherein x [n] is input signal, and y [n] is output signal, and n is positive integer, and the M value taking rolling average algorithm is 500.
After this invention takes such scheme, it is determined that the number of plies of wavelet transformation, need not by wavelet decomposition to ten which floor, Greatly reduce computation complexity, process the signal of band of interest targetedly, reduce the shadow of other band noises Ring, improve the quality of signal.
Other features and advantages of the present invention will illustrate in the following description, and, partly become from specification Obtain it is clear that or understood by implementing the present invention.The purpose of the present invention and other advantages can be by the explanations write Structure specifically noted in book, claims and accompanying drawing realizes and obtains.
Brief description
Below in conjunction with the accompanying drawings the present invention is described in detail, so that the above-mentioned advantage of the present invention is definitely.Its In,
Fig. 1 is that the present invention removes in the method for ECG baseline drift s0016lrem number in PTB database The schematic diagram of two lead signals;
Fig. 2 is that the present invention removes the schematic diagram that wavelet decomposition in the method for ECG baseline drift calculates each layer coefficients;
Fig. 3 is the signal that the present invention removes wavelet reconstruction and moving average filtering in the method for ECG baseline drift Figure;
Fig. 4 is the schematic diagram of the signal that the present invention removes in the method for ECG baseline drift after reconstruct;
Fig. 5 present invention removes the schematic diagram of baseline drift signal in the method for ECG baseline drift;
Fig. 6 present invention removes the schematic diagram removing signal after baseline drift in the method for ECG baseline drift;
Fig. 7 is letter after the present invention removes in the method for ECG baseline drift patient's I electrocardiosignal and goes baseline drift Number schematic diagram;
Fig. 8 is letter after the present invention removes in the method for ECG baseline drift patient's II electrocardiosignal and goes baseline drift Number schematic diagram;
Fig. 9 is after the present invention removes in the method for ECG baseline drift patient's III electrocardiosignal and goes baseline drift The schematic diagram of signal.
Detailed description of the invention
Describe embodiments of the present invention in detail below with reference to drawings and Examples, whereby how the present invention is applied Technological means solves technical problem, and reach technique effect realize that process can fully understand and implement according to this.Need explanation As long as not constituting conflict, each embodiment in the present invention and each feature in each embodiment can be combined with each other, The technical scheme being formed is all within protection scope of the present invention.
In addition, can be in the department of computer science of such as one group of computer executable instructions in the step shown in the flow chart of accompanying drawing System performs, and, although show logical order in flow charts, but in some cases, can be to be different from herein Order perform shown or described step.
Specifically, the present invention is directed to problem present in ECG baseline drift denoising, provide a kind of new Solution:A kind of method being combined with moving average filtering by wavelet transformation, removes the baseline drift noise of electrocardiosignal. Comprise the following steps that:
A kind of method removing ECG baseline drift, including:
Step 1) determine wavelet decomposition number of plies step, including:By the effective band of the sample rate of signal and electrocardiosignal, Determine the number of plies of wavelet decomposition;
Step 2) carry out wavelet decomposition step, including:According to the wavelet decomposition number of plies determining in step 1, sampling small echo Analyzing and carrying out J layer decomposition to electrocardiosignal, the small echo high frequency coefficient of every layer is respectively:CD1, CD2 ..., CDJ, low frequency coefficient For:CAJ;
Step 3) wavelet coefficient process and signal reconstruction step, including:Respectively by the high frequency coefficient in step 2 and low frequency system Number is processed, and carries out wavelet reconstruction according to the wavelet coefficient after processing;
Step 4) moving average filtering step, including:Average filter is moved to the signal of reconstruct in step 3, obtains The baseline drift value of electrocardiosignal;
Step 5) remove baseline drift step, including:Deducted the needle position misalignment value in step 4 by original electro-cardiologic signals, I.e. obtain the signal after removing baseline drift.
After this invention takes such scheme, it is determined that the number of plies of wavelet transformation, need not by wavelet decomposition to ten which floor, Greatly reduce computation complexity, process the signal of band of interest targetedly, reduce the shadow of other band noises Ring, improve the quality of signal.
Wherein, in further embodiment, it is preferred that step 1) in, when determining the wavelet decomposition number of plies, by such as following formula Son is carried out:
J=log2(Fs/f)
Wherein:J is the wavelet decomposition number of plies, and Fs is the sample rate of electrocardiosignal, and f is that the effective band of electrocardiosignal is maximum Value.
Preferably, in step 2 and step 3, the wavelet function carrying out wavelet function feedback employing is Haar function.
Preferably, in step 3, process to by the high frequency coefficient in step 2 and low frequency coefficient, including:
Each for small echo floor height frequency filter factor is set to 0 entirely, i.e.:CA1 '=CA2 '=... ...=CAJ '=0;
Low frequency coefficient keeps constant, CAJ '=CAJ.
Preferably, step 3) in, in signal reconstruction step, take Mallat algorithm.
Preferably, in moving average filtering step, specifically include:
Wherein x [n] is input signal, and y [n] is output signal, and n is positive integer, and the M value taking rolling average algorithm is 500.
Wherein, below in conjunction with the accompanying drawings to the present invention be embodied as be described in further detail.
The present invention have chosen German National Weights and Measures Service ecg database (PTB) data (data number is s0016lrem, Two lead) as experimental data, signal sampling rate is 1000Hz, sampling time 10s, as shown in Figure 1.The experimental procedure of this method As follows:
Step 1:Determine the wavelet decomposition number of plies.By the effective band of the sample rate of signal and electrocardiosignal, determine little The number of plies of Wave Decomposition:
J=log2(Fs/f)
Wherein:J is the wavelet decomposition number of plies, and Fs is the sample rate of electrocardiosignal, and f is that the effective band of electrocardiosignal is maximum Value.As:In the present invention, the electrocardiosignal sample rate of research is 1000Hz, and the effective band of electrocardiosignal is 0--40Hz, according to formula (1) can determine that wavelet decomposition layer should be 4 layers.The frequency band comprising in 4th layer of wavelet low frequency coefficient is 0--31.125Hz, contains The useful electrocardiosignal of major part;
Step 2:Carry out wavelet decomposition.Secondly sample rate 1000Hz according to signal studied in the present invention, uses small echo Analyzing and carrying out 4 layers of decomposition to electrocardiosignal, in order to improve computational efficiency, the base chosen in the wavelet function feedback of signal is little Ripple is Haar small echo.The small echo high frequency coefficient of every layer is respectively:CD1, CD2, CD3, CD4, low frequency coefficient is:CA4, decomposable process See Fig. 2;
Step 3:Wavelet coefficient is processed and signal reconstruction.Small echo High frequency filter coefficient is all set to 0, i.e.:CA1 '=CA2 '= CA3 '=CA4 '=0.Retain low frequency coefficient, CA4 '=CA4.Calculated by Mallat and carry out small echo according to process wavelet coefficient Reconstruct.Wherein, wavelet reconstruction is the inverse process of wavelet decomposition, and signal reconstruction flow process is as it is shown on figure 3, signal such as Fig. 4 institute after Chong Gou Showing, in Fig. 3, target signal filter high-frequency noise, retains most electrocardiosignal information and low-frequency noise, this low-frequency noise master The baseline drift of electrocardiosignal to be shown as.
Step 4:Moving average filtering.Filtering method the most frequently used in the method Digital Signal Processing of moving average filtering it One, there is realization simple, the advantages such as amount of calculation is little, excellent result.Setting input signal as x [n], output signal is y [n], moves Dynamic average filter algorithm is represented by:
Sample rate according to signal in the present invention, the M value taking rolling average algorithm is 500.Before ensureing operation efficiency Put, be found through experiments, when M value takes 400 600, go baseline drift best results.Reconstruction signal in step 3 is carried out Moving average filtering, obtains the baseline drift part of electrocardiosignal.
Step 5:Remove baseline drift.Deducted the needle position misalignment in step 4 by original electro-cardiologic signals, i.e. gone Except the signal after baseline drift, as shown in Figure 5.
The real data containing needle position misalignment is finally used to verify this method.Fig. 7--the initial data in 9 is The electrocardio of patient I, II, the III with obvious baseline drift that the IMAC1800 electrocardiograph that Zoncare company produces is gathered Signal, sample rate is 1000Hz, sampling time 9s, passes through-Fig. 7-9 and can be seen that:This method goes baseline drift respond well, Signal distortion is low.
The frequency band range that first present invention combines electrocardiosignal according to electrocardiosignal sample rate determines number of plies wavelet decomposition layer Number, carries out J time and decomposes to electrocardiosignal, and the small echo high frequency coefficient of every layer is respectively:CD1, CD2 ..., CDJ, low frequency coefficient is: CAJ;By small echo high frequency coefficient CD1, the CD2 of every layer ..., CDJ sets to 0, and retains low frequency coefficient CAJ, then carries out wavelet reconstruction, will The signal obtaining after carrying out wavelet reconstruction, through moving average filter, obtains the baseline drift part of electrocardiosignal, finally leads to Cross original electro-cardiologic signals and deduct needle position misalignment, i.e. obtain the signal after removing baseline drift.
The remarkable result of this algorithm is:Determine the number of plies of wavelet transformation, need not by wavelet decomposition to ten which floor, greatly Reduce computation complexity, process the signal of band of interest targetedly, reduce the impact of other band noises, improve The quality of signal;Additionally, the wavelet basis function chosen during carrying out the wavelet function feedback of signal is Haar small echo Function, improves operation efficiency.Process to layer coefficients each after wavelet decomposition, it will be apparent that improve efficiency of algorithm.Carry out baseline During the extraction of drift line, have employed moving average filter, utilize it to realize simple, the advantages such as amount of calculation is little, excellent result.Make Obtain this method to be applicable to process in real time.Germany's electrocardiogram (ECG) data (PTB) and actual electrocardiogram (ECG) data is finally used to remove base to this method Line drift performance is verified, result shows that this method goes baseline drift effect good, and distorted signals is low.
It should be noted that for said method embodiment, in order to be briefly described, therefore it is all expressed as a series of Combination of actions, but those skilled in the art should know, the application is not limited by described sequence of movement because According to the application, some step can use other orders or carry out simultaneously.Secondly, those skilled in the art also should know Knowing, embodiment described in this description belongs to preferred embodiment, involved action and module not necessarily the application Necessary.
Those skilled in the art are it should be appreciated that embodiments herein can be provided as method, system or computer program Product.Therefore, the application can use complete hardware embodiment, complete software implementation or the reality in terms of combining software and hardware Execute the form of example.
And, the application can use and can use at one or more computers wherein including computer usable program code The upper computer program implemented of storage medium (including but not limited to magnetic disc store, CD-ROM, optical memory etc.) Form.
Finally it should be noted that:The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, Although being described in detail the present invention with reference to previous embodiment, for a person skilled in the art, it still may be used Modify with the technical scheme described in foregoing embodiments, or equivalent is carried out to wherein portion of techniques feature. All within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made, should be included in the present invention's Within protection domain.

Claims (6)

1. the method removing ECG baseline drift, it is characterised in that include:
Step 1) determine wavelet decomposition number of plies step, including:By the effective band of the sample rate of signal and electrocardiosignal, determine Go out the number of plies of wavelet decomposition;
Step 2) carry out wavelet decomposition step, including:According to the wavelet decomposition number of plies determining in step 1, sampling wavelet analysis Carry out J layer decomposition to electrocardiosignal, the small echo high frequency coefficient of every layer is respectively:CD1, CD2 ..., CDJ, low frequency coefficient is: CAJ;
Step 3) wavelet coefficient process and signal reconstruction step, including:Respectively the high frequency coefficient in step 2 and low frequency coefficient are entered Row is processed, and carries out wavelet reconstruction according to the wavelet coefficient after processing;
Step 4) moving average filtering step, including:Average filter is moved to the signal of reconstruct in step 3, obtains electrocardio The baseline drift value of signal;
Step 5) remove baseline drift step, including:Deducted the needle position misalignment value in step 4 by original electro-cardiologic signals, to obtain final product Arrive the signal after removing baseline drift.
2. the method for removal ECG baseline drift according to claim 1, it is characterised in that step 1) in, determine During the wavelet decomposition number of plies, carry out by following formula:
J=log2(Fs/f)
Wherein:J is the wavelet decomposition number of plies, and Fs is the sample rate of electrocardiosignal, and f is the effective band maximum of electrocardiosignal.
3. the method for removal ECG baseline drift according to claim 2, it is characterised in that step 2 and step 3 In, the wavelet function carrying out wavelet function feedback employing is Haar function.
4. the method for removal ECG baseline drift according to claim 1, it is characterised in that in step 3, to by step High frequency coefficient in rapid 2 and low frequency coefficient are processed, including:
Each for small echo floor height frequency filter factor is set to 0 entirely, i.e.:CA1 '=CA2 '=...=CAJ '=0;
Low frequency coefficient keeps constant, CAJ '=CAJ.
5. the method for removal ECG baseline drift according to claim 1, it is characterised in that step 3) in, signal In reconstruction step, take Mallat algorithm.
6. the method for removal ECG baseline drift according to claim 1, it is characterised in that moving average filtering walks In Zhou, specifically include:
y [ i ] = 1 M Σ j = 0 M - 1 x ( i + j )
Wherein x [n] is input signal, and y [n] is output signal, and n is positive integer, and the M value taking rolling average algorithm is 500.
CN201610675959.4A 2016-08-12 2016-08-16 Method removing electrocardiosignal baseline drift Pending CN106419898A (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610664421 2016-08-12
CN2016106644213 2016-08-12

Publications (1)

Publication Number Publication Date
CN106419898A true CN106419898A (en) 2017-02-22

Family

ID=58182453

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610675959.4A Pending CN106419898A (en) 2016-08-12 2016-08-16 Method removing electrocardiosignal baseline drift

Country Status (1)

Country Link
CN (1) CN106419898A (en)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107260163A (en) * 2017-07-31 2017-10-20 广东南方电信规划咨询设计院有限公司 Wireless remote cardioelectric monitor system and method
CN107951482A (en) * 2017-11-13 2018-04-24 浙江好络维医疗技术有限公司 A kind of electrocardiosignal baseline based on wavelet decomposition and spline interpolation evens up method
CN109044335A (en) * 2018-07-17 2018-12-21 西安交通大学 A kind of cardiac function evaluation method based on distance speech stimulation
CN110693486A (en) * 2019-09-27 2020-01-17 武汉中旗生物医疗电子有限公司 Electrocardiogram abnormity labeling method and device
CN111166322A (en) * 2020-01-17 2020-05-19 哈尔滨工业大学 Electrocardiosignal characteristic wave extraction method
CN111772617A (en) * 2019-04-04 2020-10-16 成都心吉康科技有限公司 ECG display filtering method, device and storage medium
CN113288158A (en) * 2021-05-27 2021-08-24 河北省科学院应用数学研究所 Method, device and equipment for removing baseline drift and high-frequency noise
CN113303809A (en) * 2021-05-27 2021-08-27 河北省科学院应用数学研究所 Method, device, equipment and storage medium for removing baseline drift and high-frequency noise
CN113827253A (en) * 2020-06-04 2021-12-24 阿里巴巴集团控股有限公司 Computing device and method for removing noise from electroencephalogram signals

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102247143A (en) * 2011-06-03 2011-11-23 吉林大学珠海学院 Integratable fast algorithm for denoising electrocardiosignal and identifying QRS waves
CN102274020A (en) * 2011-06-30 2011-12-14 东北大学 Low-power consumption portable electrocardiograph monitor and control method thereof
CN104298879A (en) * 2014-10-15 2015-01-21 四川东方电气自动控制工程有限公司 Wavelet packet decomposition based gas turbine engine combustion signal analysis method
CN104367316A (en) * 2014-11-13 2015-02-25 重庆邮电大学 Electrocardiosignal denoising method based on morphological filtering and lifting wavelet transformation
CN104473635A (en) * 2014-12-16 2015-04-01 重庆邮电大学 Left-right hand motor imagery electroencephalogram characteristic extraction method mixing wavelet and common spatial pattern

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102247143A (en) * 2011-06-03 2011-11-23 吉林大学珠海学院 Integratable fast algorithm for denoising electrocardiosignal and identifying QRS waves
CN102274020A (en) * 2011-06-30 2011-12-14 东北大学 Low-power consumption portable electrocardiograph monitor and control method thereof
CN104298879A (en) * 2014-10-15 2015-01-21 四川东方电气自动控制工程有限公司 Wavelet packet decomposition based gas turbine engine combustion signal analysis method
CN104367316A (en) * 2014-11-13 2015-02-25 重庆邮电大学 Electrocardiosignal denoising method based on morphological filtering and lifting wavelet transformation
CN104473635A (en) * 2014-12-16 2015-04-01 重庆邮电大学 Left-right hand motor imagery electroencephalogram characteristic extraction method mixing wavelet and common spatial pattern

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107260163A (en) * 2017-07-31 2017-10-20 广东南方电信规划咨询设计院有限公司 Wireless remote cardioelectric monitor system and method
CN107951482A (en) * 2017-11-13 2018-04-24 浙江好络维医疗技术有限公司 A kind of electrocardiosignal baseline based on wavelet decomposition and spline interpolation evens up method
CN109044335A (en) * 2018-07-17 2018-12-21 西安交通大学 A kind of cardiac function evaluation method based on distance speech stimulation
CN109044335B (en) * 2018-07-17 2020-11-10 西安交通大学 Heart function evaluation method based on instantaneous sound stimulation
CN111772617A (en) * 2019-04-04 2020-10-16 成都心吉康科技有限公司 ECG display filtering method, device and storage medium
CN110693486A (en) * 2019-09-27 2020-01-17 武汉中旗生物医疗电子有限公司 Electrocardiogram abnormity labeling method and device
CN111166322A (en) * 2020-01-17 2020-05-19 哈尔滨工业大学 Electrocardiosignal characteristic wave extraction method
CN113827253A (en) * 2020-06-04 2021-12-24 阿里巴巴集团控股有限公司 Computing device and method for removing noise from electroencephalogram signals
CN113288158A (en) * 2021-05-27 2021-08-24 河北省科学院应用数学研究所 Method, device and equipment for removing baseline drift and high-frequency noise
CN113303809A (en) * 2021-05-27 2021-08-27 河北省科学院应用数学研究所 Method, device, equipment and storage medium for removing baseline drift and high-frequency noise
CN113288158B (en) * 2021-05-27 2022-12-20 河北省科学院应用数学研究所 Method, device and equipment for removing baseline drift and high-frequency noise

Similar Documents

Publication Publication Date Title
CN106419898A (en) Method removing electrocardiosignal baseline drift
CN109907752B (en) Electrocardiogram diagnosis and monitoring system for removing motion artifact interference and electrocardio characteristic detection
Sayadi et al. Multiadaptive bionic wavelet transform: Application to ECG denoising and baseline wandering reduction
Alfaouri et al. ECG signal denoising by wavelet transform thresholding
CN103405227B (en) Double-layer morphological filter based electrocardiosignal preprocessing method
CN110680308B (en) Electrocardiosignal denoising method based on fusion of improved EMD and threshold method
CN102973264B (en) Electrocardiosignal preprocessing method based on morphological multiresolution decomposition
Patro et al. De-noising of ECG raw signal by cascaded window based digital filters configuration
CN105852850A (en) Method and related device for extracting respiratory signals from ECG (electrocardiogram) signals
CN113384277B (en) Electrocardiogram data classification method and classification system
Tabakov et al. Online digital filter and QRS detector applicable in low resource ECG monitoring systems
Abdelmounim et al. Electrocardiogram signal denoising using discrete wavelet transform
CN113180685B (en) Electrocardio abnormity discrimination system and method based on morphological filtering and wavelet threshold
Boucheham et al. Piecewise linear correction of ECG baseline wander: a curve simplification approach
CN116584960A (en) Diaphragmatic electromyographic signal noise reduction method
Ivanov et al. Convolutional autoencoder for filtering of power-line interference with variable amplitude and frequency: Study of 12-lead PTB-XL ECG database
CN113598785A (en) Electrocardiosignal denoising method based on wavelet domain sparse characteristic
Suchetha et al. Denoising and arrhythmia classification using EMD based features and neural network
Saha et al. Noise suppressing cascaded IIR elliptic filter design for ecg signals
Gupta et al. Hardware implementation and reduction of artifacts from ECG signal
Moein et al. Intelligent ECG signal noise removal using PSONN
CN114027847B (en) Electrocardiosignal analysis method based on time-frequency analysis
Srinivasulu et al. Novel method to find the parameter for noise removal from multi-channel ecg waveforms
Narla et al. Noise Removal in ECG Signals Utilizing Fully Convolutional DAE
CN116269429A (en) Self-supervision electrocardio-noise ratio evaluation method and system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20170222

RJ01 Rejection of invention patent application after publication