CN102885616A - 一种去除脉搏波信号基线漂移的方法 - Google Patents
一种去除脉搏波信号基线漂移的方法 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
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 238000007781 pre-processing Methods 0.000 claims abstract description 3
- 230000035485 pulse pressure Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 238000010827 pathological analysis Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
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
本发明公开了一种去除脉搏波基线漂移的方法,它是先采集脉搏波信号;对采集的脉搏波信号进行预处理,获得脉搏波信号中的交流成分;然后对所获得的一段交流信号进行滑动均值滤波,抑制工频干扰,得到脉搏波原始信号;最后对原始信号进行波谷搜索,识别出该段脉搏波信号的所有波谷的位置;在相邻两个波谷之间进行样条插值,可以拟合出一条漂移基线;将原始信号减去漂移基线,即可得到去除基线漂移的脉搏波信号。本发明通过滑动均值滤波抑制了信号中的工频干扰,运用样条插值拟合的方法找到信号中的基线漂移;相比于传统的方法,本发明具有准确、稳定、实时更新、抗干扰、简单易行等优点;具有较好的实际应用价值。
Description
技术领域
本发明涉及一种去除脉搏波信号中基线漂移的方法,特别适合去除血氧饱和度测试仪中的脉搏波信号基线漂移。
背景技术
脉搏病理诊断历来受到中外医学界的重视,几乎世界上所有的民族都有过“摸脉”诊断的历史。大量临床研究表明,脉搏波中蕴含着丰富的心血管疾病信息。目前,通过对脉搏信号的分析可以方便的估算出血液动力学的各项参数,如心输出量、血管阻力、血管壁弹性和血液黏性等。然而,要想准确地判断心血管等方面的情况和预测心血管疾病发生的可能性,以便及时采取措施有效地减少危险因素,首要条件就是要能获取准确清晰的脉搏信号。脉搏信号是一种微弱的非平稳信号,输出幅度通常在几十微伏到几个毫伏,具有随机性强、信噪比低的特点。因此,在采集过程中容易受到呼吸、工频干扰和由于皮肤接触时用力大小不同而产生的基线漂移等因素的影响。其中,基线漂移是测量中影响波形准确性最大的一个因素,对信号的识别和分析造成重要的影响,而基线漂移基本上是趋势分量和低频分量且因人而异。为避免它对脉搏波信号分析和进一步处理时产生错误的影响,必须要对其进行去除。去除基线漂移常用的方法主要有拟合基线法和滤波法。传统的去除基线漂移的方法有:FIR数字滤波、小波变换法、形态滤波法和拟合基线法等,其中拟合基线法具有简洁、快速、实时性强等特点。
发明内容
本发明的目的是要提供一种去除脉搏波信号基线漂移的方法。本方法改善了脉搏波信号的抗干扰性,提高了去除基线漂移的准确性、实时性和稳定性。
实现本发明目的的技术方案是:
一种去除脉搏波信号基线漂移的方法,包括如下步骤:
(1)采集脉搏波信号;
(2)对采集到的脉搏波信号进行预处理,获得脉搏波信号中的交流成分;
(3)对所获得的一段交流信号抑制工频干扰,获得脉搏波原始信号;
(4)对原始信号进行波谷搜索,识别出所有波谷的位置;
(5)在相邻两个波谷之间进行样条插值,可以拟合出一条漂移基线;
(6)将原始信号减去漂移基线,即可得到去除基线漂移的脉搏波信号。
步骤(2)所述对采集到的脉搏波信号的预处理是对信号的数据类型转换并保存。
步骤(3)所述的交流信号抑制工频干扰,是采用5点平均法的滑动滤波得到脉搏波原始信号。
步骤(5)所述基线漂移的寻找方法是:采用样条插值拟合的方法进行去除,具体是:
(1)准确识别波谷的位置:将原始信号一阶微分,得到微分信号,通过峰值检测,找出波峰点,对应到原始信号上就是脉压上升最快的点;
(2)在相邻的脉压上升最快点之间搜索最小值,即是要找的波谷点;
(3)在相邻两个波谷之间进行样条插植,即可拟合出一条漂移基线。
本发明的优点是:通过滑动均值滤波抑制了信号中的工频干扰,运用样条插植拟合的方法去除了信号中的基线漂移;相比于传统的方法,本发明具有准确、稳定、实时更新、抗干扰、简单易行等优点;具有较好的实际应用价值。本发明的方法,简单易行,算法计算量小,硬件需求低,计算准确,结果可靠。
附图说明
图1是本发明的系统流程图;
图2是本发明中波谷搜索流程图;
图3是本发明一种波谷搜索过程示意图;
图4是本发明另一种波谷搜索过程示意图;
图5是本发明基线折线图;
图6是本发明去除基线漂移效果示意图。
具体实施方式
下面结合实施例对本发明作进一步说明。
本发明公开的去除脉搏波基线漂移的方法,也可以用于任何周期性的生理信号进行基线漂移的去除,例如心电信号,下面结合附图阐释该方法的具体实现方式。
本发明主要包括两个环节,抑制工频干扰和去除基线漂移。抑制工频干扰通过采用滑动均值滤波的方法;去除基线漂移采用样条插值拟合的方法进行去除。
抑制工频干扰具体包括如下一些步骤:
1、采集一段脉搏波数据,采样率为250Hz,采样点数为1000个;
2、选取相邻五个点的坐标,求和平均作为该点的新坐标;
3、如此循环,直到第1000个点;
4、结束,得到抑制工频干扰后的脉搏波信号。
去除基线漂移具体包括如下一些步骤:
1、在去除基线漂移之前,一个重要的步骤就是能够准确地把波谷的位置识别出来。如图3所示,将经过平滑滤波后的信号一阶微分,通过峰值检测,找出波峰点A,B的坐标,对应到原信号上就是脉压上升最快的一点。原信号中,在A→B之间寻求最小值点,得到的就是波谷点a。信号中其它波谷点也可以通过这种方法找到。从图4的识别结果中可以看到,系统准确地识别出了波谷点(a,b,c,d,e)的位置。
2、本系统选用拟合基线的方法去除基线漂移。其基本思想是:首先识别出脉搏波每个周期的基准点(波谷点a、b、c、d、e的坐标),信号首末采样点A,B也作为基准点(其中A,B点的纵坐标分别与a,b点的纵坐标相同);然后在相邻的两个基准点之间样条插值拟合,可以拟合出一条漂移曲线。如图5所示。将原始的脉搏信号曲线减去漂移曲线后,就可得到去除基线漂移的脉搏信号。从图6中可以清晰地看出脉搏信号的基线漂移得到了很好地去除。
Claims (4)
1.一种去除脉搏波信号基线漂移的方法,其特征是:包括如下步骤:
(1)采集脉搏波信号;
(2)对采集到的脉搏波信号进行预处理,获得脉搏波信号中的交流成分;
(3)对所获得的一段交流信号抑制工频干扰,获得脉搏波原始信号;
(4)对原始信号进行波谷搜索,识别出所有波谷的位置;
(5)在相邻两个波谷之间进行样条插值,可以拟合出一条漂移基线;
(6)将原始信号减去漂移基线,即可得到去除基线漂移的脉搏波信号。
2.根据权利要求1所述的方法,其特征在于:步骤(2)所述对采集到的脉搏波信号的预处理是对信号的数据类型转换并保存。
3.根据权利要求1所述的方法,其特征在于:步骤(3)所述的交流信号抑制工频干扰,是采用5点平均法的滑动滤波得到脉搏波原始信号。
4.根据权利要求1所述的方法,其特征在于:步骤(5)所述基线漂移的寻找方法是:采用样条插值拟合的方法,具体是:
(1)准确识别波谷的位置:将原始信号一阶微分,得到微分信号,通过峰值检测,找出波峰点,对应到原始信号上就是脉压上升最快的点;
(2)在相邻的脉压上升最快点之间搜索最小值,即是要找的波谷点;
(3)在相邻两个波谷之间进行样条插植,即可拟合出一条漂移基线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102469965A CN102885616A (zh) | 2012-07-17 | 2012-07-17 | 一种去除脉搏波信号基线漂移的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2012102469965A CN102885616A (zh) | 2012-07-17 | 2012-07-17 | 一种去除脉搏波信号基线漂移的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102885616A true CN102885616A (zh) | 2013-01-23 |
Family
ID=47529455
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012102469965A Pending CN102885616A (zh) | 2012-07-17 | 2012-07-17 | 一种去除脉搏波信号基线漂移的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102885616A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104224140A (zh) * | 2014-09-09 | 2014-12-24 | 桂林电子科技大学 | 一种利用提升小波变换滤除基线漂移的方法 |
CN104266981A (zh) * | 2014-08-28 | 2015-01-07 | 北京环境特性研究所 | 一种太赫兹波吸收光谱识别方法 |
CN104730336A (zh) * | 2013-12-19 | 2015-06-24 | 财团法人资讯工业策进会 | 频率检测装置及频率检测方法 |
US9188612B2 (en) | 2013-12-19 | 2015-11-17 | Institute For Information Industry | Frequency measurement device and frequency measurement method |
WO2017118127A1 (zh) * | 2016-01-05 | 2017-07-13 | 深圳和而泰智能控制股份有限公司 | 一种心跳信号处理方法、装置和系统 |
CN107184187A (zh) * | 2017-07-03 | 2017-09-22 | 重庆大学 | 基于DTCWT‑Spline的脉搏波信号去噪处理方法 |
CN107468221A (zh) * | 2017-09-12 | 2017-12-15 | 济南大学 | 基于金氏脉学的脉搏波信号处理方法及肺癌检测系统 |
CN109498022A (zh) * | 2018-12-29 | 2019-03-22 | 西安理工大学 | 一种基于光电容积脉搏波的呼吸频率提取方法 |
CN110801210A (zh) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | 一种脉搏波信号的滤波方法、装置、可读介质及电子设备 |
CN111528825A (zh) * | 2020-05-14 | 2020-08-14 | 浙江大学 | 一种光电容积脉搏波信号优化方法 |
CN112932410A (zh) * | 2020-10-23 | 2021-06-11 | 深圳忆恩医疗科技有限公司 | 基于切比雪夫多项式拟合的生理包络信号幅值计算方法 |
CN113100727A (zh) * | 2021-05-12 | 2021-07-13 | 深圳市通久电子有限公司 | 实时分析识别脉搏波峰的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101371801A (zh) * | 2008-10-10 | 2009-02-25 | 华南理工大学 | 基于ZigBee的病房无线监护系统 |
CN101991410A (zh) * | 2009-08-31 | 2011-03-30 | 深圳市理邦精密仪器股份有限公司 | 一种脉率搜索和计算方法 |
CN102028457A (zh) * | 2010-11-24 | 2011-04-27 | 北京麦邦光电仪器有限公司 | 脉率测量方法及指环式脉率测量仪 |
CN102085095A (zh) * | 2009-12-07 | 2011-06-08 | 深圳市新元素医疗技术开发有限公司 | 一种检测心电图中st段的方法、系统和心电检测仪器 |
-
2012
- 2012-07-17 CN CN2012102469965A patent/CN102885616A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101371801A (zh) * | 2008-10-10 | 2009-02-25 | 华南理工大学 | 基于ZigBee的病房无线监护系统 |
CN101991410A (zh) * | 2009-08-31 | 2011-03-30 | 深圳市理邦精密仪器股份有限公司 | 一种脉率搜索和计算方法 |
CN102085095A (zh) * | 2009-12-07 | 2011-06-08 | 深圳市新元素医疗技术开发有限公司 | 一种检测心电图中st段的方法、系统和心电检测仪器 |
CN102028457A (zh) * | 2010-11-24 | 2011-04-27 | 北京麦邦光电仪器有限公司 | 脉率测量方法及指环式脉率测量仪 |
Non-Patent Citations (2)
Title |
---|
姚建阁,等: "基于三次样条插值技术的去除胸阻抗基线漂移的方法研究", 《中国医学物理学杂志》 * |
宋维军: "脉搏波自动采集分析的研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104730336B (zh) * | 2013-12-19 | 2017-09-15 | 财团法人资讯工业策进会 | 频率检测装置及频率检测方法 |
CN104730336A (zh) * | 2013-12-19 | 2015-06-24 | 财团法人资讯工业策进会 | 频率检测装置及频率检测方法 |
US9188612B2 (en) | 2013-12-19 | 2015-11-17 | Institute For Information Industry | Frequency measurement device and frequency measurement method |
TWI510790B (zh) * | 2013-12-19 | 2015-12-01 | Inst Information Industry | 頻率檢測裝置及頻率檢測方法 |
CN104266981A (zh) * | 2014-08-28 | 2015-01-07 | 北京环境特性研究所 | 一种太赫兹波吸收光谱识别方法 |
CN104224140A (zh) * | 2014-09-09 | 2014-12-24 | 桂林电子科技大学 | 一种利用提升小波变换滤除基线漂移的方法 |
WO2017118127A1 (zh) * | 2016-01-05 | 2017-07-13 | 深圳和而泰智能控制股份有限公司 | 一种心跳信号处理方法、装置和系统 |
CN107184187B (zh) * | 2017-07-03 | 2019-08-02 | 重庆大学 | 基于DTCWT-Spline的脉搏波信号去噪处理方法 |
CN107184187A (zh) * | 2017-07-03 | 2017-09-22 | 重庆大学 | 基于DTCWT‑Spline的脉搏波信号去噪处理方法 |
CN107468221A (zh) * | 2017-09-12 | 2017-12-15 | 济南大学 | 基于金氏脉学的脉搏波信号处理方法及肺癌检测系统 |
CN109498022A (zh) * | 2018-12-29 | 2019-03-22 | 西安理工大学 | 一种基于光电容积脉搏波的呼吸频率提取方法 |
CN110801210A (zh) * | 2019-11-06 | 2020-02-18 | 心核心科技(北京)有限公司 | 一种脉搏波信号的滤波方法、装置、可读介质及电子设备 |
CN111528825A (zh) * | 2020-05-14 | 2020-08-14 | 浙江大学 | 一种光电容积脉搏波信号优化方法 |
CN112932410A (zh) * | 2020-10-23 | 2021-06-11 | 深圳忆恩医疗科技有限公司 | 基于切比雪夫多项式拟合的生理包络信号幅值计算方法 |
CN113100727A (zh) * | 2021-05-12 | 2021-07-13 | 深圳市通久电子有限公司 | 实时分析识别脉搏波峰的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102885616A (zh) | 一种去除脉搏波信号基线漂移的方法 | |
CN109907752B (zh) | 一种去除运动伪影干扰与心电特征检测的心电诊断与监护系统 | |
CN103027667B (zh) | 脉搏波的特征参数提取 | |
CN103549950B (zh) | 移动心电监测差分阈值检测改进算法 | |
CN101732033A (zh) | 一种人体波形中的特征参数的提取方法和装置 | |
CN106037694A (zh) | 一种基于脉搏波的连续血压测量装置 | |
Liu et al. | A multi-step method with signal quality assessment and fine-tuning procedure to locate maternal and fetal QRS complexes from abdominal ECG recordings | |
CN107468221B (zh) | 基于金氏脉学的脉搏波信号处理方法及肺癌检测系统 | |
CN104161510A (zh) | 一种多级导联心电信号qrs波形识别方法 | |
CN106618542A (zh) | 一种去噪心率信号检测装置及方法 | |
CN107184187B (zh) | 基于DTCWT-Spline的脉搏波信号去噪处理方法 | |
CN107361753A (zh) | 基于脉搏波形特征点的人体健康状态监护方法 | |
CN114532993A (zh) | 一种癫痫患者脑电高频振荡信号的自动检测方法 | |
TW201717845A (zh) | 從具雜訊之心電圖資料決定心跳率的設備、系統和方法 | |
CN104644160A (zh) | 心电图伪差信号识别方法及装置 | |
CN108272451A (zh) | 一种基于改进小波变换的qrs波识别方法 | |
CN108403094A (zh) | 一种识别脉搏波波峰的方法 | |
CN104095640A (zh) | 血氧饱和度检测方法及装置 | |
CN104680186A (zh) | 一种心电信号中st段形态的自动分类方法 | |
CN101637394A (zh) | 一种心音信号定位、分段方法 | |
CN102334986B (zh) | 动脉压信号中重搏切迹点识别方法 | |
CN109330582A (zh) | 基于心电信号分析的心率及其特征指标检测方法 | |
CN107361764A (zh) | 一种心电信号特征波形r波的快速提取方法 | |
CN110327032A (zh) | 一种单导心电信号pqrst波联合精准识别算法 | |
Tekeste et al. | Adaptive ECG interval extraction |
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 |