CN102885616A - Method for removing baseline drift from pulse wave signal - Google Patents
Method for removing baseline drift from pulse wave signal Download PDFInfo
- Publication number
- CN102885616A CN102885616A CN2012102469965A CN201210246996A CN102885616A CN 102885616 A CN102885616 A CN 102885616A CN 2012102469965 A CN2012102469965 A CN 2012102469965A CN 201210246996 A CN201210246996 A CN 201210246996A CN 102885616 A CN102885616 A CN 102885616A
- Authority
- CN
- China
- Prior art keywords
- signal
- pulse wave
- baseline
- drift
- wave signal
- 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
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
The invention discloses a method for removing baseline drift from a pulse wave. The method comprises the following steps of: acquiring a pulse wave signal; pre-treating the acquired pulse wave signal to obtain an alternating current component in the pulse wave signal; performing sliding mean filtering on a section of obtained alternating current signal, and inhibiting the power frequency interference to obtain an original signal of the pulse wave; finally, searching wave troughs of the original signal to recognize the positions of all wave troughs of the pulse wave signal; performing spline interpolation between two adjacent wave troughs to fit a drifting baseline; and substracting the drifting baseline from the original signal to obtain a pulse wave signal in which the baseline drift is removed. According to the method, the power frequency interference in the signal is inhibited through sliding mean filtering, and the baseline drift in the signal is found by a spline interpolation fitting method; compared with the traditional method, the invention has the advantages that the method is accurate, stable, anti-jamming, simple, easy, and updated in real time; and the method has good practical application value.
Description
Technical field
The present invention relates to a kind of method of removing baseline drift in the pulse wave signal, be particularly suitable for removing the pulse wave signal baseline drift in the blood oxygen saturation tester.
Background technology
The pulse pathological diagnosis always is subject to the attention of China and foreign countries' medical circle, and almost there was the history of " feeling pulse " diagnosis in all nationality in the world.A large amount of clinical researches show, are containing abundant cardiovascular disease information in the pulse wave.At present, can estimate easily hemodynamic parameters by the analysis to pulse signal, such as cardiac output, vascular resistance, vessel wall elasticity and blood stickiness etc.Yet, want to judge exactly situation and the pathogenetic probability of predicting cardiovascular disease of the aspects such as cardiovascular, in order in time take measures effectively to reduce risk factor, the most important condition wants can obtain accurately clearly pulse signal exactly.Pulse signal is a kind of faint non-stationary signal, and output amplitude arrives several millivolts at tens microvolts usually, has the advantages that randomness is strong, signal to noise ratio is low.Therefore, in gatherer process, breathed easily, power frequency is disturbed and the impact of the factors such as baseline drift that produce of firmly varying in size during owing to contact skin.Wherein, baseline drift is a factor that affects waveform accuracy maximum in measuring, and identification and the analysis of signal caused important impact, and baseline drift is trend component and low frequency component basically and varies with each individual.For avoiding it to the pulse wave signal analysis and producing wrong impact when further processing, must remove it.Remove baseline drift method commonly used and mainly contain match base-line method and filter method.The method of traditional removal baseline drift has: FIR digital filtering, Wavelet Transform, shape filtering method and match base-line method etc., wherein match base-line method has the characteristics such as succinct, quick, real-time.
Summary of the invention
The objective of the invention is to provide a kind of method of removing the pulse wave signal baseline drift.This method has been improved the anti-interference of pulse wave signal, has improved accuracy, real-time and the stability of removing baseline drift.
The technical scheme that realizes the object of the invention is:
A kind of method of removing the pulse wave signal baseline drift comprises the steps:
(1) gathers pulse wave signal;
(2) pulse wave signal that collects is carried out pretreatment, obtain the alternating component in the pulse wave signal;
(3) the one section AC signal that obtains is suppressed power frequency and disturb, obtain the pulse wave primary signal;
(4) primary signal is carried out the trough search, identify the position of all troughs;
(5) between adjacent two troughs, carry out spline interpolation, can simulate a drift baseline;
(6) primary signal is deducted the drift baseline, can obtain removing the pulse wave signal of baseline drift.
The described pretreatment to the pulse wave signal that collects of step (2) is to the data type conversion of signal and preserves.
The described AC signal of step (3) suppresses power frequency to be disturbed, and is to adopt the glide filter of 5 averaging method to obtain the pulse wave primary signal.
The finding method of the described baseline drift of step (5) is: adopt the method for spline interpolation match to remove, specifically:
(1) accurately identify the position of trough: with primary signal single order differential, obtain differential signal, detect by peak value, find out wave crest point, corresponding on the primary signal is exactly the pulse pressure the fastest point that rises;
(2) searching for minima between adjacent pulse pressure rises the most as soon as possible, namely is the trough point that will look for;
(3) between adjacent two troughs, carry out the batten transplanting, can simulate a drift baseline.
Advantage of the present invention is: suppressed the power frequency in the signal by the slip mean filter and disturbed, used the method for batten transplanting match to remove the baseline drift in the signal; Than traditional method, the present invention has accurately, stablizes, real-time update, the advantage such as anti-interference, simple; Has preferably actual application value.Method of the present invention, simple, the algorithm amount of calculation is little, and hsrdware requirements are low, calculate accurately reliable results.
Description of drawings
Fig. 1 is system flow chart of the present invention;
Fig. 2 is trough search routine figure among the present invention;
Fig. 3 is a kind of trough search procedure of the present invention sketch map;
Fig. 4 is the another kind of trough search procedure of the present invention sketch map;
Fig. 5 is baseline broken line graph of the present invention;
Fig. 6 is that the present invention removes baseline drift effect sketch map.
The specific embodiment
The invention will be further described below in conjunction with embodiment.
The method of removal pulse wave disclosed by the invention baseline drift also can be used for the removal that any periodic physiological signal carry out baseline drift, and electrocardiosignal for example is below in conjunction with the specific implementation of accompanying drawing explaination the method.
The present invention mainly comprises two links, suppresses power frequency and disturbs and remove baseline drift.Suppressing power frequency disturbs by adopting the method for slip mean filter; Removing baseline drift adopts the method for spline interpolation match to remove.
Suppress the power frequency interference and specifically comprise the following step:
1, gather one section pulse wave data, sample rate is 250Hz, and sampling number is 1000;
2, choose the coordinate of adjacent five points, sum-average arithmetic is as the new coordinate of this point;
3, so circulation is until the 1000th point;
4, finish the pulse wave signal after the power frequency that is inhibited is disturbed.
Remove baseline drift and specifically comprise the following step:
1, before removing baseline drift, an important step is exactly can be exactly the location recognition of trough out.As shown in Figure 3, will detect by peak value through the signal single order differential behind the smothing filtering, find out wave crest point A, the coordinate of B, corresponding on the original signal is exactly a bit that pulse pressure rises the fastest.In the original signal, seek minimum point between A → B, that obtain is exactly trough point a.Other trough point also can find by this method in the signal.Can see that from the recognition result of Fig. 4 system has identified the position of trough point (a, b, c, d, e) exactly.
2, native system selects the method for match baseline to remove baseline drift.Its basic thought is: at first identify the datum mark (coordinate of trough point a, b, c, d, e) in each cycle of pulse wave, signal first and last sampled point A, B also as datum mark (A wherein, the vertical coordinate that B is ordered respectively with a, the vertical coordinate that b is ordered is identical); Then spline interpolation match between two adjacent datum marks can simulate a drift curve.As shown in Figure 5.With original pulse signal curve deduct the drift curve after, just can obtain removing the pulse signal of baseline drift.The baseline drift that can clearly find out pulse signal from Fig. 6 has obtained removing well.
Claims (4)
1. a method of removing the pulse wave signal baseline drift is characterized in that: comprise the steps:
(1) gathers pulse wave signal;
(2) pulse wave signal that collects is carried out pretreatment, obtain the alternating component in the pulse wave signal;
(3) the one section AC signal that obtains is suppressed power frequency and disturb, obtain the pulse wave primary signal;
(4) primary signal is carried out the trough search, identify the position of all troughs;
(5) between adjacent two troughs, carry out spline interpolation, can simulate a drift baseline;
(6) primary signal is deducted the drift baseline, can obtain removing the pulse wave signal of baseline drift.
2. method according to claim 1 is characterized in that: the described pretreatment to the pulse wave signal that collects of step (2) is to the data type conversion of signal and preserves.
3. method according to claim 1 is characterized in that: the described AC signal of step (3) suppresses power frequency to be disturbed, and is to adopt the glide filter of 5 averaging method to obtain the pulse wave primary signal.
4. method according to claim 1, it is characterized in that: the finding method of the described baseline drift of step (5) is: adopt the method for spline interpolation match, specifically:
(1) accurately identify the position of trough: with primary signal single order differential, obtain differential signal, detect by peak value, find out wave crest point, corresponding on the primary signal is exactly the pulse pressure the fastest point that rises;
(2) searching for minima between adjacent pulse pressure rises the most as soon as possible, namely is the trough point that will look for;
(3) between adjacent two troughs, carry out the batten transplanting, can simulate a drift baseline.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2012102469965A CN102885616A (en) | 2012-07-17 | 2012-07-17 | Method for removing baseline drift from pulse wave signal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2012102469965A CN102885616A (en) | 2012-07-17 | 2012-07-17 | Method for removing baseline drift from pulse wave signal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN102885616A true CN102885616A (en) | 2013-01-23 |
Family
ID=47529455
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2012102469965A Pending CN102885616A (en) | 2012-07-17 | 2012-07-17 | Method for removing baseline drift from pulse wave signal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102885616A (en) |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104224140A (en) * | 2014-09-09 | 2014-12-24 | 桂林电子科技大学 | Method for filtering baseline drift by using lifting wavelet transformation |
| CN104266981A (en) * | 2014-08-28 | 2015-01-07 | 北京环境特性研究所 | Terahertz wave absorption spectrum recognition method |
| CN104730336A (en) * | 2013-12-19 | 2015-06-24 | 财团法人资讯工业策进会 | Frequency detection device and frequency detection method |
| US9188612B2 (en) | 2013-12-19 | 2015-11-17 | Institute For Information Industry | Frequency measurement device and frequency measurement method |
| WO2017118127A1 (en) * | 2016-01-05 | 2017-07-13 | 深圳和而泰智能控制股份有限公司 | Heartbeat signal processing method, device and system |
| CN107184187A (en) * | 2017-07-03 | 2017-09-22 | 重庆大学 | Pulse Wave Signal Denoising processing method based on DTCWT Spline |
| CN107468221A (en) * | 2017-09-12 | 2017-12-15 | 济南大学 | Pulse wave signal processing method and lung cancer detection system based on Jin Shi pulse studies |
| CN109498022A (en) * | 2018-12-29 | 2019-03-22 | 西安理工大学 | A kind of respiratory rate extracting method based on photoplethysmographic |
| CN110801210A (en) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | Pulse wave signal filtering method and device, readable medium and electronic equipment |
| CN111528825A (en) * | 2020-05-14 | 2020-08-14 | 浙江大学 | Photoelectric volume pulse wave signal optimization method |
| CN112932410A (en) * | 2020-10-23 | 2021-06-11 | 深圳忆恩医疗科技有限公司 | Physiological envelope signal amplitude calculation method based on Chebyshev polynomial fitting |
| CN113100727A (en) * | 2021-05-12 | 2021-07-13 | 深圳市通久电子有限公司 | Method for analyzing and identifying pulse wave crest in real time |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101371801A (en) * | 2008-10-10 | 2009-02-25 | 华南理工大学 | Sickroom wireless monitoring system based on ZigBee |
| CN101991410A (en) * | 2009-08-31 | 2011-03-30 | 深圳市理邦精密仪器股份有限公司 | Pulse rate searching and calculating method |
| CN102028457A (en) * | 2010-11-24 | 2011-04-27 | 北京麦邦光电仪器有限公司 | Pulse rate measuring method and ring type pulse rate measuring meter |
| CN102085095A (en) * | 2009-12-07 | 2011-06-08 | 深圳市新元素医疗技术开发有限公司 | Method, system and electrocardioscanner for detecting ST segment in electrocardiogram |
-
2012
- 2012-07-17 CN CN2012102469965A patent/CN102885616A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101371801A (en) * | 2008-10-10 | 2009-02-25 | 华南理工大学 | Sickroom wireless monitoring system based on ZigBee |
| CN101991410A (en) * | 2009-08-31 | 2011-03-30 | 深圳市理邦精密仪器股份有限公司 | Pulse rate searching and calculating method |
| CN102085095A (en) * | 2009-12-07 | 2011-06-08 | 深圳市新元素医疗技术开发有限公司 | Method, system and electrocardioscanner for detecting ST segment in electrocardiogram |
| CN102028457A (en) * | 2010-11-24 | 2011-04-27 | 北京麦邦光电仪器有限公司 | Pulse rate measuring method and ring type pulse rate measuring meter |
Non-Patent Citations (2)
| Title |
|---|
| 姚建阁，等: "基于三次样条插值技术的去除胸阻抗基线漂移的方法研究", 《中国医学物理学杂志》 * |
| 宋维军: "脉搏波自动采集分析的研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104730336B (en) * | 2013-12-19 | 2017-09-15 | 财团法人资讯工业策进会 | Frequency detection device and frequency detection method |
| CN104730336A (en) * | 2013-12-19 | 2015-06-24 | 财团法人资讯工业策进会 | Frequency detection device and frequency detection method |
| US9188612B2 (en) | 2013-12-19 | 2015-11-17 | Institute For Information Industry | Frequency measurement device and frequency measurement method |
| TWI510790B (en) * | 2013-12-19 | 2015-12-01 | Inst Information Industry | Frequency measurement device and frequency measurement method |
| CN104266981A (en) * | 2014-08-28 | 2015-01-07 | 北京环境特性研究所 | Terahertz wave absorption spectrum recognition method |
| CN104224140A (en) * | 2014-09-09 | 2014-12-24 | 桂林电子科技大学 | Method for filtering baseline drift by using lifting wavelet transformation |
| WO2017118127A1 (en) * | 2016-01-05 | 2017-07-13 | 深圳和而泰智能控制股份有限公司 | Heartbeat signal processing method, device and system |
| CN107184187B (en) * | 2017-07-03 | 2019-08-02 | 重庆大学 | Pulse Wave Signal Denoising processing method based on DTCWT-Spline |
| CN107184187A (en) * | 2017-07-03 | 2017-09-22 | 重庆大学 | Pulse Wave Signal Denoising processing method based on DTCWT Spline |
| CN107468221A (en) * | 2017-09-12 | 2017-12-15 | 济南大学 | Pulse wave signal processing method and lung cancer detection system based on Jin Shi pulse studies |
| CN109498022A (en) * | 2018-12-29 | 2019-03-22 | 西安理工大学 | A kind of respiratory rate extracting method based on photoplethysmographic |
| CN110801210A (en) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | Pulse wave signal filtering method and device, readable medium and electronic equipment |
| CN111528825A (en) * | 2020-05-14 | 2020-08-14 | 浙江大学 | Photoelectric volume pulse wave signal optimization method |
| CN112932410A (en) * | 2020-10-23 | 2021-06-11 | 深圳忆恩医疗科技有限公司 | Physiological envelope signal amplitude calculation method based on Chebyshev polynomial fitting |
| CN113100727A (en) * | 2021-05-12 | 2021-07-13 | 深圳市通久电子有限公司 | Method for analyzing and identifying pulse wave crest in real time |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102885616A (en) | Method for removing baseline drift from pulse wave signal | |
| CN103027667B (en) | Characteristic parameter extraction of pulse wave | |
| CN109907752B (en) | Electrocardiogram diagnosis and monitoring system for removing motion artifact interference and electrocardio characteristic detection | |
| CN103584854B (en) | Extraction method of electrocardiosignal R waves | |
| CN101828918B (en) | Electrocardiosignal R peak detection method based on waveform characteristic matching | |
| CN109124610B (en) | Anti-interference method and device for non-invasive blood pressure measurement | |
| Chanwimalueang et al. | Enabling R-peak detection in wearable ECG: Combining matched filtering and Hilbert transform | |
| CN101732033A (en) | Method and device for extracting characteristic parameter in human body waveform | |
| CN103549950A (en) | Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring | |
| CN107622259B (en) | T wave detection method, electrocardiogram data analysis method and device | |
| CN109381181A (en) | The end-point detecting method of electrocardiosignal signature waveform | |
| CN105147269A (en) | Noninvasive continuous blood pressure measuring method | |
| CN110313903B (en) | Pulse wave frequency domain characteristic parameter extraction method and device | |
| CN104095627B (en) | A kind of electrocardiogram digitized signal quality soft decision method and device | |
| CN109745035A (en) | Electro-cardiologic signal waveforms detection method | |
| CN104068841B (en) | A kind of measuring method and device measuring Indices of Systolic Time parameter | |
| CN109009086A (en) | A kind of R wave of electrocardiosignal detection method and system | |
| CN107361764A (en) | A kind of rapid extracting method of electrocardiosignal signature waveform R ripples | |
| CN111528821A (en) | Method for identifying characteristic points of counterpulsation waves in pulse waves | |
| CN112656392A (en) | Low-power PPG heart rate calculation method and electronic equipment | |
| CN108042107A (en) | A kind of PPG signals puppet difference correcting method | |
| CN103083011B (en) | Electrocardio R wave crest real-time locating method assisted by chest impedance second order difference image | |
| CN102217931A (en) | Method and device for acquiring heart rate variation characteristic parameter | |
| CN101897578B (en) | Method for segmenting arterial pressure signal by beats | |
| CN102334986B (en) | Method for identifying dicrotic notch point in arterial tension signal |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20130123 |