CN102499694B - 用于消除血氧饱和度监测干扰的方法 - Google Patents
用于消除血氧饱和度监测干扰的方法 Download PDFInfo
- Publication number
- CN102499694B CN102499694B CN2011102829311A CN201110282931A CN102499694B CN 102499694 B CN102499694 B CN 102499694B CN 2011102829311 A CN2011102829311 A CN 2011102829311A CN 201110282931 A CN201110282931 A CN 201110282931A CN 102499694 B CN102499694 B CN 102499694B
- Authority
- CN
- China
- Prior art keywords
- signal
- infrared light
- hongguang
- component
- ppg
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
用于消除血氧饱和度监测干扰的方法,它对采集到的红光、红外光PPG信号进行数字滤波,去除直流成分中的低频干扰信号,用以计算准确的直流分量值,再通过计算PPG信号包络提取交流分量,以分离出与干扰相关的信号作为参考信号进行自适应滤波,消除交流成分中频率与目标信号混叠的运动干扰,计算得到准确的PPG交流分量值,最后根据朗伯比尔定律计算出消除干扰后的血氧饱和度。本发明方法能很好地抑制运动干扰,计算得到准确的血氧饱和度值,此外,本发明方法直接从红光和红外光的PPG信号中提取出参考信号进行自适应滤波,构建方法简单,不需要额外装置提供参考信号,也无需进行多次自适应滤波,算法简单,计算量小,利于实现监测设备小型化和佩戴化。
Description
技术领域
本发明涉及一种干扰消除方法,特别是一种血氧饱和度监测的干扰消除方法。
背景技术
与传统的医院环境血氧饱和度监测相比,动态环境人体血氧饱和度监测具有更广泛的应用前景,它在院前急救、战场士兵状态监测、高原病预防、航天员生理状态监测等领域都能发挥积极的作用。但由于受运动干扰的影响,动态环境血氧饱和度监测结果的准确性、稳定性通常较差,限制了它的发展与应用。
针对上述问题,近年来国内外出现了不少抑制运动干扰的方法,有的对探头结构进行改进,有的通过算法进行信号处理,这些方法都起到了消除或降低运动干扰影响的作用。但对于动态环境监测,其运动通常较为剧烈,现存的抗运动干扰方法在处理这种强烈运动干扰时,还存在许多不足和局限性。
例如通过增加设备模块,引入监测部位运动信息进行自适应滤波,或者增设多波长光源进行监测完成干扰抵消,都能实现运动干扰的消除,但这些方法增设了硬件模块,物理结构复杂,不易满足动态环境监测对设备小型化和佩戴化的要求。
再如通过特征识别法,删除包含有运动干扰的信号,以达到排除运动引起的计算误差的目的,但这种方法仅仅在微弱运动或偶发性运动造成的干扰方面有较好的抗干扰性能,而对于较强的连续性运动干扰效果不明显。由Massimo公司提出的离散饱和度变换算法(DST),通过获取的红光与红外光信号计算产生一组带有待定参数的参考信号进行自适应滤波,在0~100%血氧饱和度的范围内进行扫描,寻找出当前实际的血氧饱和度值。这种方法具有极低的假阳性率,是一种较为理想的抗运动干扰血氧饱和度测量方法。但该算法获取一次血氧值需要在0~100%的血氧饱和度范围内进行多次自适应滤波,算法复杂,计算量大,受动态环境监测设备微处理器计算能力及其功耗的限制,实现有很大困难。
发明内容
本发明的目的就是提供一种用于消除血氧饱和度监测干扰的方法,它根据人体PPG信号的特点,从交流分量入手,有效的抑制运动干扰,准确计算出血氧饱和度值。
本发明的目的是通过这样的技术方案实现的,其步骤为:
4)根据朗伯比尔定律计算得到光密度比值R对应相应的血氧饱和度,
进一步,步骤3)中所述的自适应消噪算法为:
进一步,步骤1)所述得到PPG信号后,对PPG信号的高频成分进行低通平滑滤波,去除信号毛刺。
进一步,步骤2)中所述去除直流分量的低频干扰具体方法是:通过低通平滑滤波后的PPG信号以单个波为单位计算信号包络线,通过上、下包络线计算出包络线的均值,用包络线均值对信号进行处理,去除低频干扰,再利用处理后的信号计算两种光的PPG信号的直流分量。
进一步,包络线计算方法为:搜寻PPG信号中的极值点,以极值点为原始数据,以采样点个数为数据量进行插值,得到PPG信号的包络线。
由于采用了上述技术方案,本发明具有如下的优点:
算法能很好的抑制运动干扰,计算得到准确的血氧饱和度值,弥补了传统PPG信号处理中无法有效去除运动干扰的缺陷,算法直接从红光和红外光的PPG信号中提取出参考信号进行自适应滤波,构建方法简单,不需要额外装置提供参考信号。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和权利要求书来实现和获得。
附图说明
本发明的附图说明如下。
图1为本发明自适应算法原理图。
图2为本发明参考信号获取的原理图。
图3为本发明包络线获取示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
用于消除血氧饱和度监测干扰的方法,其步骤为:
根据朗伯比尔定律计算得到光密度比值R对应相应的血氧饱和度,
红光、红外光通过人体反射后进行信号采集,采集到的信号P可分为直流分量和交流分量,在直流分量和交流分量中含有干扰信号和目标信号,根据朗伯比尔定律,完全去除直流分量的干扰信号和交流分量的干扰信号,可以准确的计算出相应的血氧饱和度。加载在信号交流分量上的干扰信号属于高频干扰,其频率带包含目标信号频率到更高频的频带。其中与目标信号频带重叠的部分不易消除。加载在信号直流分量上的干扰信号属于低频干扰,主要表现为信号基线的漂移,可用频率滤波的方式消除干扰。
如图1所示,步骤3)中所述的自适应消噪算法为:
血氧饱和度可通过两光源产生的PPG信号的交直流比值关系确定。对于动态血氧饱和度监测,若能消除PPG信号中交流成分与直流成分中混杂的运动干扰,便可计算得到准确的血氧饱和度值。其中直流成分处于低频段,可通过数字滤波对其进行处理,而交流成分中的运动干扰与信号频带重叠,无法通过简单的滤波实现信号去噪。因此,设计采用最小均方LMS法对PPG交流分量中的运动干扰进行自适应消噪。
根据假设,应为恒定值,表现为一直流信号。但当有运动干扰引入时,红光和红外光的交流信号、都会随着光路的变化而产生波动。因氧合血红蛋白和还原血红蛋白对红光、红外光的吸光系数存在差异,使得两光的交流信号波动程度不同,这时两光交流信号之差为:
差值会随着干扰的的引入而出现波动,表征了波动的程度。中直流分量为目标信号差值,即为一直流信号,而交流信号为噪声引起的波动。波动信号与干扰信号相关而与目标信号无关,使用这个波动信号作为参考信号,如图2所示,分别对红光与红外光的交流成分、进行自适应滤波去除直流信号,即可得到。
步骤1)所述得到PPG信号后,对PPG信号的高频成分进行低通平滑滤波,去除信号毛刺。对高频成分进行低通平滑滤波是对信号做出预处理,消除信号中的毛刺,使信号最后得消噪效果更好。
如图3所示,步骤2)中所述去除直流分量的低频干扰具体方法是:通过低通平滑滤波后的PPG信号以单个波为单位计算信号包络线,通过上、下包络线计算出包络线的均值,用包络线均值对信号进行处理,去除低频干扰,再利用处理后的信号计算两种光的PPG信号的直流分量。
包络线计算方法为:搜寻PPG信号中的极值点,以极值点为原始数据,以采样点个数为数据量进行插值,得到PPG信号的包络线。
要实现对PPG信号交流分量的运动干扰自适应消噪,参考信号的获取尤为重要。包络线是一种有效反映信号整体变化的方法,本文将其用于提取脉搏波中干扰信号,作为自适应对消器的参考信号。采用包络线方法同时带了以下两个优点:
扩充可用于计算的数据量:通常情况下,通过峰峰值计算PPG信号的交流信号至少需要一个完整的波形,也即是一个PPG波形仅能获取一个可用于计算的交流信号值。这样得到的数据量太少,自适应滤波需要很长的时间才能收敛。而先计算出包络线,通过上下包络线之差得到交流分量,这样可使交流信号数据量扩展为采样点个数。
更好体现交流分量的变化趋势:每个脉搏波的交流分量之间加入了大量数据点,能够更为细致的表现信号波动的情况,有助于参考信号的构建。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (4)
1. 用于消除血氧饱和度监测干扰的方法,其特征在于,其步骤为:
4)根据朗伯比尔定律计算得到光密度比值R对应相应的血氧饱和度,
步骤3)中所述的自适应消噪算法为:
2. 如权利要求1所述的用于消除血氧饱和度监测干扰的方法,其特征在于:步骤1)所述得到红光、红外光的PPG信号P后,对红光、红外光的PPG信号P的高频成分进行低通平滑滤波,去除信号毛刺。
3. 如权利要求2所述的用于消除血氧饱和度监测干扰的方法,其特征在于,步骤2)中所述去除直流分量的低频干扰具体方法是:
通过低通平滑滤波后的红光、红外光的PPG信号P以单个波为单位计算信号包络线,通过上、下包络线计算出包络线的均值,用包络线均值对信号进行处理,去除低频干扰,再利用处理后的信号计算红光、红外光的PPG信号P的直流分量。
4. 如权利要求3所述的用于消除血氧饱和度监测干扰的方法,其特征在于,包络线计算方法为:搜寻红光、红外光的PPG信号P中的极值点,以极值点为原始数据,以采样点个数为数据量进行插值,得到红光、红外光的PPG信号P的包络线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102829311A CN102499694B (zh) | 2011-09-22 | 2011-09-22 | 用于消除血氧饱和度监测干扰的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102829311A CN102499694B (zh) | 2011-09-22 | 2011-09-22 | 用于消除血氧饱和度监测干扰的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102499694A CN102499694A (zh) | 2012-06-20 |
CN102499694B true CN102499694B (zh) | 2013-11-27 |
Family
ID=46211876
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102829311A Expired - Fee Related CN102499694B (zh) | 2011-09-22 | 2011-09-22 | 用于消除血氧饱和度监测干扰的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102499694B (zh) |
Families Citing this family (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103027690B (zh) * | 2012-11-28 | 2014-08-27 | 华中科技大学 | 一种基于自相关建模法的低灌注血氧饱和度测量方法 |
CN104095640A (zh) * | 2013-04-03 | 2014-10-15 | 达尔生技股份有限公司 | 血氧饱和度检测方法及装置 |
US20160066863A1 (en) * | 2013-04-05 | 2016-03-10 | Nitto Denko Corporation | METHOD AND APPARATUS FOR DETERMINING SpO2 OF A SUBJECT FROM AN OPTICAL MEASUREMENT |
CN103698410B (zh) * | 2013-12-19 | 2016-01-20 | 江苏大学 | 一种导波信号包络线的绘制方法 |
CN103860152B (zh) * | 2014-04-02 | 2017-01-04 | 辛勤 | 一种对脉搏波信号进行处理的方法 |
WO2015176999A1 (en) * | 2014-05-22 | 2015-11-26 | Koninklijke Philips N.V. | Method and apparatus for optical sensing of tissue variation at increased accuracy |
EP3203912B1 (en) | 2014-10-10 | 2021-03-10 | Medtor LLC | System and method for a non-invasive medical sensor |
CN105099573B (zh) * | 2014-11-21 | 2018-10-26 | 深圳迈瑞生物医疗电子股份有限公司 | 一种干扰信号处理方法、装置及医疗检测设备 |
DE102015103894A1 (de) * | 2015-03-17 | 2016-09-22 | Fritz Stephan Gmbh Medizintechnik | Beatmungsgeräte sowie Steuerverfahren für Beatmungsgeräte |
CN107666860B (zh) * | 2015-06-03 | 2021-06-29 | 皇家飞利浦有限公司 | 光电体积描记装置 |
CN105322974B (zh) * | 2015-10-22 | 2017-11-28 | 深圳市美的连医疗电子股份有限公司 | 一种干扰信号消除方法及用户终端 |
TWI657794B (zh) | 2017-01-09 | 2019-05-01 | 財團法人工業技術研究院 | 生理資訊偵測裝置及使用其之生理資訊偵測方法 |
CN109009159B (zh) * | 2018-07-24 | 2020-06-02 | 华中科技大学 | 一种适用于可穿戴血氧芯片的抗运动干扰加速电路 |
WO2020258022A1 (zh) * | 2019-06-25 | 2020-12-30 | 深圳市蓝瑞格生物医疗科技有限公司 | 一种提升测量系统测量精度的方法及系统 |
CN112587133A (zh) * | 2020-12-09 | 2021-04-02 | 深圳市奋达智能技术有限公司 | 一种血氧饱和度的测量方法 |
CN113367676B (zh) * | 2021-06-07 | 2022-06-21 | 广东电网有限责任公司 | 一种用于可穿戴式设备的ppg运动噪声滤波方法及装置 |
CN113545777A (zh) * | 2021-08-31 | 2021-10-26 | 维沃移动通信有限公司 | 血氧检测方法、血氧检测装置、电子设备和可读存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0283289A2 (en) * | 1987-03-20 | 1988-09-21 | C.R. Bard, Inc. | Excitation and detection apparatus for remote sensor connected by optical fiber |
US5081998A (en) * | 1989-09-01 | 1992-01-21 | Critikon, Inc. | Optically stabilized infrared energy detector |
CN101632588A (zh) * | 2009-08-04 | 2010-01-27 | 中国人民解放军第四军医大学 | 消除运动干扰的血氧饱和度测量方法 |
-
2011
- 2011-09-22 CN CN2011102829311A patent/CN102499694B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0283289A2 (en) * | 1987-03-20 | 1988-09-21 | C.R. Bard, Inc. | Excitation and detection apparatus for remote sensor connected by optical fiber |
US5081998A (en) * | 1989-09-01 | 1992-01-21 | Critikon, Inc. | Optically stabilized infrared energy detector |
CN101632588A (zh) * | 2009-08-04 | 2010-01-27 | 中国人民解放军第四军医大学 | 消除运动干扰的血氧饱和度测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102499694A (zh) | 2012-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102499694B (zh) | 用于消除血氧饱和度监测干扰的方法 | |
CN110840445B (zh) | 一种动态心电信号的自动降噪方法 | |
CN105078505B (zh) | 生理信号处理方法及处理装置 | |
Han et al. | Electrocardiogram signal denoising based on a new improved wavelet thresholding | |
Wei et al. | Study on conditioning and feature extraction algorithm of photoplethysmography signal for physiological parameters detection | |
CN107949321A (zh) | 时域干扰去除和改进的心率测量跟踪机理 | |
CN107464226B (zh) | 一种基于改进二维经验模态分解算法的图像去噪方法 | |
CN106889984B (zh) | 一种心电信号自动降噪方法 | |
CN105997043B (zh) | 一种基于腕式可穿戴设备的脉率提取方法 | |
CN110664390A (zh) | 基于腕带式ppg和深度学习的心率监测系统及方法 | |
CN105662375A (zh) | 非接触式检测生命体征信号的方法和装置 | |
CN104305992B (zh) | 一种交互式胎儿心电快速自动提取方法 | |
US20160029968A1 (en) | Tracking slow varying frequency in a noisy environment and applications in healthcare | |
CN111160090B (zh) | 一种bcg信号降噪方法及系统 | |
CN103860152B (zh) | 一种对脉搏波信号进行处理的方法 | |
CN113349752A (zh) | 一种基于传感融合的可穿戴设备实时心率监测方法 | |
CN103230267A (zh) | 一种抗运动干扰的脉率提取方法 | |
CN113848544A (zh) | 基于多普勒雷达的人体存在感应检测方法、装置、存储介质 | |
CN112370036A (zh) | 基于级联型rls自适应滤波的ppg心率提取装置和方法 | |
CN107184187A (zh) | 基于DTCWT‑Spline的脉搏波信号去噪处理方法 | |
Zhang et al. | ECG signals denoising method based on improved wavelet threshold algorithm | |
CN109758140A (zh) | 一种心率检测方法和装置 | |
CN112587133A (zh) | 一种血氧饱和度的测量方法 | |
Wang et al. | Research on denoising algorithm for ECG signals | |
TWI505816B (zh) | 血氧飽和度檢測方法及裝置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20131127 Termination date: 20170922 |
|
CF01 | Termination of patent right due to non-payment of annual fee |