一种血氧饱和度检测方法、装置,电子设备及存储介质转让专利
申请号 : CN202111312528.9
文献号 : CN113749655B
文献日 : 2022-03-08
发明人 : 张广洁 , 汪志伟
申请人 : 芯原微电子(南京)有限公司
摘要 :
权利要求 :
1.一种血氧饱和度检测方法,其特征在于,包括:获取待检测用户的PPG检测信号;
基于形态学滤波算法对所述PPG检测信号进行处理;
基于处理后的PPG检测信号,确定所述待检测用户的血氧饱和度;
其中,所述基于形态学滤波算法对所述PPG检测信号进行处理包括:基于所述形态学滤波算法对所述PPG检测信号进行平滑降噪;基于所述形态学滤波算法对平滑降噪后的PPG检测信号进行基线校准;
基于所述形态学滤波算法进行平滑降噪的公式为:;
;
;
其中, 表示所述PPG检测信号; 为三角型结构元素; 为直线型结构元素; 表示腐蚀操作; 表示膨胀操作; 和 均为中间值; 表示所述平滑降噪后的PPG检测信号;
其中,所述获取待检测用户的PPG检测信号,包括:获取所述待检测用户在进行PPG检测时预设时长的原始信号;对所述预设时长的原始信号进行有效性验证;当所述预设时长的原始信号为有效信号时,将所述预设时长的原始信号添加至缓冲区;当所述预设时长的原始信号为无效信号时,将所述缓冲区的信号清空;其中,所述缓冲区可存储第一预设容量的信号;当所述缓冲区存储的信号达到第二预设容量时,所述缓冲区内的原始信号为所述PPG检测信号;所述第二预设容量小于所述第一预设容量;
所述预设时长的原始信号包括预设时长的第一原始信号和预设时长的第二原始信号,所述预设时长的第一原始信号及所述预设时长的第二原始信号为所述待检测用户通过不同波长的两种光源进行PPG检测生成的;所述有效性验证包括所述预设时长的第一原始信号与所述预设时长的第二原始信号之间的相关性验证;
其中,相关性验证通过如下公式进行验证:;
;
;
表示所述预设时长的第一原始信号和所述预设时长的第二原始信号之间的相关系数; 表示所述预设时长的第一原始信号; 表示所述预设时长的第二原始信号;
表示所述预设时长的第一原始信号的方差; 表示所述预设时长的第二原始信号的方差; 表示所述预设时长的第一原始信号和所述预设时长的第二原始信号之间的协方差; 表示数学期望;
所述有效性验证还包括信号振幅验证、信号斜率验证;其中,所述信号振幅验证为将所述预设时长的原始信号的最大值减去最小值作为振幅,将该振幅与有效振幅区间进行比对实现验证;所述信号斜率验证为将所述预设时长的原始信号中的一个指定窗口的斜率进行验证;其中,该指定窗口的斜率表示该指定窗口的最后一个采样点的值和第一个采样点的值组成的斜率;
在所述基于形态学滤波算法对所述PPG检测信号进行处理之后,所述方法还包括:对处理后的第一原始信号和处理后的第二原始信号进行相关性验证;确定所述处理后的第一原始信号和所述处理后的第二原始信号为有效信号。
2.根据权利要求1所述的方法,其特征在于,基于所述形态学滤波算法进行基线校准的公式为:
;
;
;
其中, 表示所述平滑降噪后的PPG检测信号, 为直线型结构元素; 为开运算;为闭运算; 和 均为中间值; 表示基线校准后的PPG检测信号。
3.根据权利要求1所述的方法,其特征在于,所述有效性验证还包括信号阈值验证;
所述缓冲区包括第一缓冲区和第二缓冲区;所述第一缓冲区用于存储有效的第一原始信号,所述第二缓冲区用于存储有效的第二原始信号;所述第一缓冲区和所述第二缓冲区存储的为同一时间段所检测到的信号;当所述第一缓冲区存储的信号达到所述第二预设容量,且所述第二缓冲区存储的信号达到所述第二预设容量,所述第一缓冲区内的第一原始信号和所述第二缓冲区内的第二原始信号为所述PPG检测信号。
4.一种血氧饱和度检测装置,其特征在于,包括:获取模块,用于获取待检测用户的PPG检测信号;
滤波模块,用于基于形态学滤波算法对所述PPG检测信号进行处理;
确定模块,用于处理后的PPG检测信号,确定所述待检测用户的血氧饱和度;
其中,所述滤波模块具体用于基于所述形态学滤波算法对所述PPG检测信号进行平滑降噪;基于所述形态学滤波算法对平滑降噪后的PPG检测信号进行基线校准;
基于所述形态学滤波算法进行平滑降噪的公式为:;
;
;
其中, 表示所述PPG检测信号; 为三角型结构元素; 为直线型结构元素; 表示腐蚀操作; 表示膨胀操作; 和 均为中间值; 表示所述平滑降噪后的PPG检测信号;
其中,所述获取模块还具体用于获取所述待检测用户在进行PPG检测时预设时长的原始信号;对所述预设时长的原始信号进行有效性验证;当所述预设时长的原始信号为有效信号时,将所述预设时长的原始信号添加至缓冲区;当所述预设时长的原始信号为无效信号时,将所述缓冲区的信号清空;其中,所述缓冲区可存储第一预设容量的信号;当所述缓冲区存储的信号达到第二预设容量时,所述缓冲区内的原始信号为所述PPG检测信号;所述第二预设容量小于所述第一预设容量;
所述预设时长的原始信号包括预设时长的第一原始信号和预设时长的第二原始信号,所述预设时长的第一原始信号及所述预设时长的第二原始信号为所述待检测用户通过不同波长的两种光源进行PPG检测生成的;所述有效性验证包括所述预设时长的第一原始信号与所述预设时长的第二原始信号之间的相关性验证;
其中,相关性验证通过如下公式进行验证:;
;
;
表示所述预设时长的第一原始信号和所述预设时长的第二原始信号之间的相关系数; 表示所述预设时长的第一原始信号; 表示所述预设时长的第二原始信号;
表示所述预设时长的第一原始信号的方差; 表示所述预设时长的第二原始信号的方差; 表示所述预设时长的第一原始信号和所述预设时长的第二原始信号之间的协方差; 表示数学期望;
所述有效性验证还包括信号振幅验证、信号斜率验证;其中,所述信号振幅验证为将所述预设时长的原始信号的最大值减去最小值作为振幅,将该振幅与有效振幅区间进行比对实现验证;所述信号斜率验证为将所述预设时长的原始信号中的一个指定窗口的斜率进行验证;其中,该指定窗口的斜率表示该指定窗口的最后一个采样点的值和第一个采样点的值组成的斜率;
所述确定模块还用于在所述基于形态学滤波算法对所述PPG检测信号进行处理之后,对处理后的第一原始信号和处理后的第二原始信号进行相关性验证;确定所述处理后的第一原始信号和所述处理后的第二原始信号为有效信号。
5.一种电子设备,其特征在于,包括:处理器和存储器,所述处理器和所述存储器连接;
所述存储器用于存储程序;
所述处理器用于运行存储在所述存储器中的程序,执行如权利要求1‑3中任一项所述的方法。
6.一种计算机可读存储介质,其特征在于,其上存储有计算机程序,所述计算机程序在被计算机运行时执行如权利要求1‑3中任一项所述的方法。
说明书 :
一种血氧饱和度检测方法、装置,电子设备及存储介质
技术领域
背景技术
入血液,再通过光敏传感器将血液反射或者透射出的光线转换为电信号,通过分析电信号
的变化来分析血液容积的变化。在获取到PPG检测信号后,需要对该信号进行降噪和/或基
线漂移等处理,在处理后,即可通过朗博‑比尔定律等计算方式得到活体的血氧饱和度。
差,影响血氧饱和度检测的准确性。
发明内容
号,确定所述待检测用户的血氧饱和度;其中,所述基于形态学滤波算法对所述PPG检测信
号进行处理包括:基于所述形态学滤波算法对所述PPG检测信号进行平滑降噪;基于所述形
态学滤波算法对平滑降噪后的PPG检测信号进行基线校准;基于所述形态学滤波算法进行
平滑降噪的公式为: ; ; ;其中, 表示所述PPG检测
信号; 为三角型结构元素; 为直线型结构元素; 表示腐蚀操作; 表示膨胀操作;
和 均为中间值; 表示所述平滑降噪后的PPG检测信号。
相比于均值滤波、高斯滤波、带通滤波等方法计算更简单。采用形态学滤波算法对PPG检测
信号再进行基线校准,相比于采用巴特沃兹带通滤波校正基线漂移,形态学滤波的方法能
够大大降低计算的复杂度。且巴特沃兹带通滤波需要仔细选择滤波的截止频率,才能对基
线进行有效校正,而形态学滤波只需要选择直线型结构性元素,结构元素长度在信号周期
误差±20个采样点以内的范围都可以较好地进行基线校准,相比于巴特沃兹带通滤波,形
态学滤波算法对参数更加不敏感,可以适用于更多的人群和场景。可见,通过形态学滤波算
法可以降低信号处理的复杂度,提高信号处理效果,进而提高血氧饱和度检测的准确性。
中, 表示所述平滑降噪后的PPG检测信号, 为直线型结构元素; 为开运算; 为闭运
算; 和 均为中间值; 表示基线校准后的PPG检测信号。
所述预设时长的原始信号进行有效性验证;当所述预设时长的原始信号为有效信号时,将
所述预设时长的原始信号添加至缓冲区;当所述预设时长的原始信号为无效信号时,将所
述缓冲区的信号清空;其中,所述缓冲区可存储第一预设容量的信号;当所述缓冲区存储的
信号达到第二预设容量时,所述缓冲区内的原始信号为所述PPG检测信号;所述第二预设容
量小于所述第一预设容量。
的原始信号,可以起到更加精细化的验证,进而提高验证的准确性。
原始信号及所述预设时长的第二原始信号为所述待检测用户通过不同波长的两种光源进
行PPG检测生成的;所述有效性验证包括信号振幅验证、信号斜率验证、信号阈值验证及所
述预设时长的第一原始信号与所述预设时长的第二原始信号之间的相关性验证的至少一
种;所述缓冲区包括第一缓冲区和第二缓冲区;所述第一缓冲区用于存储有效的第一原始
信号,所述第二缓冲区用于存储有效的第二原始信号;所述第一缓冲区和所述第二缓冲区
存储的为同一时间段所检测到的信号;当所述第一缓冲区存储的信号达到所述第二预设容
量,且所述第二缓冲区存储的信号达到所述第二预设容量,所述第一缓冲区内的第一原始
信号和所述第二缓冲区内的第二原始信号为所述PPG检测信号。
信号进行有效地检测。
和处理后的第二原始信号进行相关性验证;确定所述处理后的第一原始信号和所述处理后
的第二原始信号为有效信号。
过滤无效数据,进而提高后续血氧饱和度检测地准确度。
行处理;确定模块,用于处理后的PPG检测信号,确定所述待检测用户的血氧饱和度;
形态学滤波算法进行平滑降噪的公式为: ; ; ;其中,
表示所述PPG检测信号; 为三角型结构元素; 为直线型结构元素; 表示腐蚀操作;
表示膨胀操作; 和 均为中间值; 表示所述平滑降噪后的PPG检测信号。
的程序,执行如上述第一方面实施例和/或结合上述第一方面实施例的一些可能的实现方
式提供的方法。
施例的一些可能的实现方式提供的方法。
附图说明
作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以
根据这些附图获得其他相关的附图。
具体实施方式
是,但不限于个人计算机(Personal Computer,PC)、智能手机、平板电脑、个人数字助理
(Personal Digital Assistant,PDA)、移动上网设备(Mobile Internet Device,MID)等。
服务器可以是但不限于网络服务器、数据库服务器、云服务器或由多个子服务器构成的服
务器集成等。当然,上述列举的设备仅用于便于理解本申请实施例,其不应作为对本实施例
的限定。
置包括至少一个可以软件或固件(Firmware)的形式存储在存储器120中或固化在电子设备
100的操作系统(Operating System,OS)中的软件模块。处理器110用于执行存储器120中存
储的可执行模块,例如,血氧饱和度检测装置所包括的软件功能模块及计算机程序等,以实
现血氧饱和度检测方法。处理器110可以在接收到执行指令后,执行计算机程序。
器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated
Circuit ,ASIC)、分立门或晶体管逻辑器件、分立硬件组件,可以实现或者执行本申请实施
例中的公开的各方法、步骤及逻辑框图。此外,通用处理器可以是微处理器或者任何常规处
理器等。
Memory,PROM)、可擦可编程序只读存储器(Erasable Programmable Read‑Only Memory,
EPROM),以及电可擦编程只读存储器(Electric Erasable Programmable Read‑Only
Memory,EEPROM)。存储器120用于存储程序,处理器110在接收到执行指令后,执行该程序。
可以通过软件、硬件或其组合实现。
法不以图2及以下所示的顺序为限制,该方法包括:步骤S101‑步骤S103。
长的原始信号为有效信号时,将预设时长的原始信号添加至缓冲区;当预设时长的原始信
号为无效信号时,将缓冲区的信号清空。
秒的信号。
小。示例性的,预设时长为一秒。当获取待检测用户在进行PPG检测时一秒的原始信号后,对
这一秒的原始信号进行有效性验证。若这一秒的原始信号验证通过,即这一秒的原始信号
为有效信号,则将该信号添加至缓冲区中。此时判断缓冲区存储的信号是否达到第二预设
容量,若是,则进行步骤S102的信号处理,若否,则继续对下一秒的原始信号进行验证。若这
一秒的原始信号未验证通过,则将缓冲区内的所有信号清空,并继续对下一秒的原始信号
进行验证。
的有效信号,则无法使得缓冲区内为连续的有效信号,因此,需要将缓冲区内之前添加的所
有信号进行删除,并从下一次有效信号开始重新计算。此外,通过设定可存储第一预设容量
的信号的缓冲区可以控制内存的占用,且每次验证预设时长的原始信号,可以起到实时验
证的效果。再者,每次验证预设时长的原始信号,可以起到更加精细化的验证,进而提高验
证的准确性。
所述待检测用户通过不同波长的两种光源进行PPG检测生成的。
带来血管中血液容积变化,会影响不同波长的光被脱氧血红蛋白和氧合血红蛋白吸收的程
度,从而计算出血氧饱和度,所以血氧饱和度检测使用了红光和红外光两种光源。因此,上
述的两种光源为红光和红外光,即第一原始信号为通过红光进行PPG检测生成,而第二原始
信号为通过红外光进行PPG检测生成,第一原始信号所对应的红光的波长选择范围为600nm
到800nm,第二原始信号所对应的红外光的波长选择范围为800nm到1000nm。
存储有效的第二原始信号。第一缓冲区和第二缓冲区存储的为同一时间段所检测到的信
号,且两个缓冲区的容量相同。当第一缓冲区存储的信号达到第二预设容量,且第二缓冲区
存储的信号达到第二预设容量,第一缓冲区内的第一原始信号和第二缓冲区内的第二原始
信号为PPG检测信号。
行第一次有效性验证。其中,第一次有效性验证方式包括信号振幅验证、信号斜率验证、信
号阈值验证及所述预设时长的第一原始信号与所述预设时长的第二原始信号之间的相关
性验证的至少一种。
可以同时包括上述四种验证方式。
间,则确定当前预设时长的信号有效。有效振幅区间可以根据实际情况进行设定,本申请在
数值上不作限定。
处于预设的有效区间,则确定当前预设时长的信号无效,若当前预设时长的原始信号的最
大有效值和最小有效值均处于预设的有效区间,则确定当前预设时长的信号有效。预设的
有效区间可以根据实际情况进行设定,本申请在数值上不作限定。
样点的值和第一个采样点的值组成的斜率。若当前预设时长的原始信号的指定窗口的斜率
大于预设斜率值,则确定当前预设时长的信号无效;若当前预设时长的原始信号的指定窗
口的斜率小于预设斜率值,则确定当前预设时长的信号有效。
的第二原始信号均进行信号振幅验证;若采用信号阈值验证,则需对预设时长的第一原始
信号和预设时长的第二原始信号均进行信号阈值验证;若采用信号斜率验证,则需对预设
时长的第一原始信号和预设时长的第二原始信号均进行信号斜率验证。
于预设值,则说明二者变化趋势相同,二者信号有效,而若预设时长的第一原始信号和预设
时长的第二原始信号之间的相关系数小于预设值,则说明二者变化趋势不同,至少其中一
个信号出现异常,此时确定预设时长的第一原始信号和预设时长的第二原始信号均为无效
信号。
原始信号之间的相关系数; 表示预设时长的第一原始信号; 表示预设时长的第二原始
信号; 表示预设时长的第一原始信号的方差; 表示预设时长的第二原始信号
的方差; 表示预设时长的第一原始信号和预设时长的第二原始信号之间的协方
差; 表示数学期望。由于 的绝对值介于0 1之间,通常情况下, 越接近1,
~
表示 和 相关程度越强,反之则越弱。因此,可以上述预设值可以设定为0.5、0.6。由于上
述公式均为本领域所熟知,本申请不作过多说明。
信号进行有效地检测。
个缓冲区的信号均清空,并重新开始新一轮的验证。若两个缓冲区存储的信号均到达第二
预设容量后,则执行步骤S102。
号,对此,本申请不作限定。
一原始信号;在处理第二原始信号时, 表示第二原始信号); 为三角型结构元素,具体
值可以为[0,0,0,1,2,3,2,1,0,0,0]; 为直线型结构元素,具体可以为长度为11的值全
为0的直线型; 表示腐蚀操作; 表示膨胀操作; 和 均为中间值;先通过三角型结构
元素对PPG检测信号 进行腐蚀操作,再使用直线型结构元素对腐蚀结果进行膨胀操作即
可得到信号 ;先通过三角型结构元素对PPG检测信号 进行膨胀操作,再使用直线型结
构元素对膨胀结果进行腐蚀操作即可得到信号 ; 表示平滑降噪后的PPG检测信号。需要
说明的是, 的长度需要大于等于结构元素的长度,即 的采样点的个数需要大于11。
一原始信号,在处理第二原始信号时, 表示第二原始信号), 为直线型结构元素; 为
开运算; 为闭运算; 和 均为中间值;先使用开运算对PPG检测信号 进行第一层形
态学滤波,再使用闭运算进行第二层形态学滤波即可得到 ;先使用闭运算对PPG检测信号
进行第一层形态学滤波,再使用开运算进行第二层形态学滤波即可得到 ; 表示基
线校准后的PPG检测信号。
沃兹带通滤波需要仔细选择滤波的截止频率,才能对基线进行有效校正,而形态学滤波只
需要选择直线型结构性元素,结构元素长度在信号周期误差±20个采样点以内的范围都可
以较好地进行基线校准,相比于巴特沃兹带通滤波,形态学滤波算法对参数更加不敏感,可
以适用于更多的人群和场景。
信号进行基线校准。
第一原始信号,在处理第二原始信号时, 表示第二原始信号); 为三角型结构元素,具
体值可以为[0,0,0,1,2,3,2,1,0,0,0]; 为直线型结构元素,具体可以为长度为11的值
全为0的直线型; 表示腐蚀操作; 表示膨胀操作; 和 均为中间值;先通过三角型结
构元素对PPG检测信号 进行腐蚀操作,再使用直线型结构元素对腐蚀结果进行膨胀操作
即可得到信号 ;先通过三角型结构元素对PPG检测信号 进行膨胀操作,再使用直线型
结构元素对膨胀结果进行腐蚀操作即可得到信号 ; 表示平滑降噪后的PPG检测信号。
需要说明的是, 的长度需要大于等于结构元素的长度,即 的采样点的个数需要大于
11。 也为直线型结构元素; 为开运算; 为闭运算; 和 均为中间值;先使用开运算
对平滑降噪后的PPG检测信号 进行第一层形态学滤波,再使用闭运算进行第二层形态学
滤波即可得到 ;先使用闭运算对平滑降噪后的PPG检测信号 进行第一层形态学滤波,
再使用开运算进行第二层形态学滤波即可得到 ; 表示基线校准后的PPG检测信号。
滤波算法,对此,本申请不作限定。
处的整个缓冲区的信号为经过上述处理后的PPG检查信号。
后的第一原始信号和处理后的第二原始信号为有效信号,则执行步骤S103,若处理后的第
一原始信号和处理后的第二原始信号为无效信号,则将两个缓冲区清空。通过该方式可以
进一步过滤无效数据,进而提高后续血氧饱和度检测地准确度。
次验证可以进一步地提高验证的准确性和可靠性。
交流分量,最后使用朗博‑比尔定律计算当前血氧饱和度的结果。由于计算血氧饱和度的方
式为本领域所熟知,此处不作过多赘述。
号处理效果,进而提高血氧饱和度检测的准确性。
有效信号进行处理,并继续基于处理后的缓冲区中所有的有效信号确定待检测用户的血氧
饱和度。
信号包括第一原始信号和第二原始信号);然后,对预设时长的原始信号进行第一次有效性
验证,当预设时长的原始信号为有效信号时,将预设时长的原始信号添加至缓冲区;当预设
时长的原始信号为无效信号时,将缓冲区的信号清空。在存入缓冲区后,判断缓冲区存储的
信号是否达到第二预设容量,若否,则继续对新获取的预设时长的原始信号进行第一次有
效验证,若是,则缓冲区内的信号为PPG检测信号,此时基于形态学滤波算法对PPG检测信号
进行平滑降噪;再基于形态学滤波算法对平滑降噪后的PPG检测信号进行基线校准。接着,
再对基线校准后的PPG检测信号进行第二次有效性验证,此处的验证为对整个缓冲区进行
的验证,验证通过后,则可计算出血氧饱和度,若验证未通过,则清空缓冲区,重新获取预设
时长的原始信号开始验证。
效信号时,将所述预设时长的原始信号添加至缓冲区;当所述预设时长的原始信号为无效
信号时,将所述缓冲区的信号清空;其中,所述缓冲区可存储第一预设容量的信号;当所述
缓冲区存储的信号达到第二预设容量时,所述缓冲区内的原始信号为所述PPG检测信号;所
述第二预设容量小于所述第一预设容量。
用户通过不同波长的两种光源进行PPG检测生成的;所述有效性验证包括信号振幅验证、信
号斜率验证、信号阈值验证及所述预设时长的第一原始信号与所述预设时长的第二原始信
号之间的相关性验证的至少一种;所述缓冲区包括第一缓冲区和第二缓冲区;所述第一缓
冲区用于存储有效的第一原始信号,所述第二缓冲区用于存储有效的第二原始信号;所述
第一缓冲区和所述第二缓冲区存储的为同一时间段所检测到的信号;当所述第一缓冲区存
储的信号达到所述第二预设容量,且所述第二缓冲区存储的信号达到所述第二预设容量,
所述第一缓冲区内的第一原始信号和所述第二缓冲区内的第二原始信号为所述PPG检测信
号。
所述处理后的第一原始信号和所述处理后的第二原始信号为有效信号。
程,在此不再赘述。
硬盘、磁带)、光介质(例如DVD)、或者半导体介质(例如固态硬盘Solid State Disk (SSD))
等。
辑功能划分,实际实现时可以有另外的划分方式,又例如,多个单元或组件可以结合或者可
以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间
的耦合或直接耦合或通信连接可以是通过一些通信接口,装置或单元的间接耦合或通信连
接,可以是电性,机械或其它的形式。
个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的
目的。
的关系或者顺序。
任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。