CN108652642A - 一种血氧饱和度测量方法、设备及存储介质 - Google Patents
一种血氧饱和度测量方法、设备及存储介质 Download PDFInfo
- Publication number
- CN108652642A CN108652642A CN201810181589.8A CN201810181589A CN108652642A CN 108652642 A CN108652642 A CN 108652642A CN 201810181589 A CN201810181589 A CN 201810181589A CN 108652642 A CN108652642 A CN 108652642A
- Authority
- CN
- China
- Prior art keywords
- pulse signal
- blood oxygen
- infrared light
- feux rouges
- light pulse
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/145—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue
- A61B5/1455—Measuring characteristics of blood in vivo, e.g. gas concentration, pH value; Measuring characteristics of body fluids or tissues, e.g. interstitial fluid, cerebral tissue using optical sensors, e.g. spectral photometrical oximeters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Power Engineering (AREA)
- Optics & Photonics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种血氧饱和度测量方法、设备及存储介质。所述方法包括:采样脉搏血氧信号,并获取红光脉搏信号以及红外光脉搏信号;对红光脉搏信号以及红外光脉搏信号进行滤波处理;确定红光脉搏信号以及红外光脉搏信号的直流成分和交流成分;根据红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小;根据比率值R的大小确定采样脉搏血氧值;对采样脉搏血氧值进行频次统计,得到频次最大值,并将频次最大值作为当前血氧测量结果;将当前血氧测量结果作为血氧饱和度输出。本发明提供的一种血氧饱和度测量方法、设备及存储介质可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
Description
技术领域
本发明涉及一种血氧测量领域,尤其涉及一种血氧饱和度测量方法、设备及存储介质。
背景技术
脉搏血氧饱和度检测是利用血红蛋白和氧合血红蛋白对光的吸收量随动脉波动而变化的原理。血氧饱和度计算一般先计算比率值R,再根据R来计算血氧饱和度,血氧饱和度测量的关键也就在于R值的计算。目前对R值的计算广泛采用的是峰-峰值法和频谱法。
峰-峰值法是在计算每个脉搏的红光和红外光的波峰和波谷,因此波峰波谷位置的准确性对计算脉搏血氧饱和度有重要影响,但抗干扰性稍差。频谱法是另一种计算脉搏血氧饱和度的方法,频谱法较峰值法,提高了抗干扰性能,但增加算法复杂度和计算量,不宜在一些可穿戴设备上推广应用。
发明内容
有鉴于此,本发明实施例提供的一种血氧饱和度测量方法、设备及存储介质解决了现有技术中抗干扰能力差、算法复杂且不适宜在可穿戴设备上推广应用的问题。
本发明实施例提供的一种血氧饱和度测量方法、设备及存储介质可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
本发明实施例提供的一种血氧饱和度测量方法,包括:
采样脉搏血氧信号,并获取红光脉搏信号以及红外光脉搏信号;
对所述红光脉搏信号以及红外光脉搏信号进行滤波处理;
确定红光脉搏信号以及红外光脉搏信号的直流成分和交流成分;
根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小;
根据所述比率值R的大小确定所述采样脉搏血氧值;
对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果;
将所述当前血氧测量结果作为血氧饱和度输出。
优选的,所述对所述红光脉搏信号以及红外光脉搏信号进行滤波处理,具体包括:
设置滤波带宽对所述红光脉搏信号以及红外光脉搏信号进行滤波;
根据滤波结果设置红光脉搏信号序列以及红外光脉搏信号序列。
优选的,所述确定所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分,具体包括:
根据原始采样红光脉搏信号、原始采样红外光脉搏信号、滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定直流成分;
根据所述滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定交流成分;
其中,所述原始采样红光脉搏信号以及原始采样红外光脉搏信号为所述初始采样的红光脉搏信号以及红外光脉搏信号。
优选的,所述根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小,具体包括:
确定红外光脉搏信号的交流成分与红光脉搏信号的交流成分的比值,记做A;
确定红光脉搏信号的直流成分与红外光脉搏信号的直流成分的比值,记做B;
根据比率值R计算公式确定比率值R,即R=A*B。
优选的,所述根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小,具体还包括:
选取所述红光脉搏信号序列与红外光脉搏信号序列的相同长度的序列;
确定所述红光脉搏信号序列的最大值以及所述红外光脉搏信号序列的最大值;
根据递推算法得出比率值R的序列值。
优选的,所述根据所述比率值R的大小确定所述采样脉搏血氧值,具体包括:
用直方图来统计某段时间内出现相同比率值R值次数最多计算采样脉搏血氧值。
优选的,所述对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果,具体包括:
将所述采样脉搏血氧值在1-100范围内进行频次统计,得到血氧直方图;
在所述血氧直方图中找出频次最大值,将所述频次最大值作为当前血氧测量结果。
本发明实施例提供的一种血氧饱和度测量方法可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
另外,为实现上述目的,本发明实施例还提供一个设备,所述设备包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的血氧饱和度测量程序,所述血氧饱和度测量程序被所述处理器执行时实现所述的血氧饱和度测量方法的步骤。
此外,为实现上述目的,本发明还提供一种存储介质,所述存储介质上存储有血氧饱和度测量程序,所述血氧饱和度测量程序被处理器执行时实现所述的血氧饱和度测量方法的步骤。
本发明实施例提供的一种设备及存储介质可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
附图说明
图1为本发明实施例一种血氧饱和度测量方法第一实施方式的流程示意图;
图2为本发明实施例一种血氧饱和度测量方法第二实施方式的流程示意图;
图3为本发明实施例一种血氧饱和度测量方法第三实施方式的流程示意图;
图4为本发明实施例一种血氧饱和度测量方法第四实施方式的流程示意图;
图5为本发明实施例一种血氧饱和度测量方法第五实施方式的流程示意图;
图6为本发明实施例一种血氧饱和度测量方法第六实施方式的流程示意图;
图7示出了适于用来实现本发明实施方式的式例性设备的框图。
具体实施方式
下面,结合附图以及具体实施方式,对本发明做进一步描述,需要说明的是,在不相冲突的前提下,以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。
参阅图1,图1为本发明实施例一种血氧饱和度测量方法第一实施方式的流程示意图。
在步骤S101,采样脉搏血氧信号,并获取红光脉搏信号以及红外光脉搏信号。
具体的,采用血氧探测器以及模数转换电路采样所述脉搏血氧信号,并通过红光和红外光的照射来获取红光脉搏信号以及红外光脉搏信号。
在步骤S102,对所述红光脉搏信号以及红外光脉搏信号进行滤波处理。
具体的,通过低通滤波器对所述红光脉搏信号以及红外光脉搏信号进行滤波处理。
在步骤S103,确定红光脉搏信号以及红外光脉搏信号的直流成分和交流成分。
具体的,经过滤波后的红光脉搏信号以及红外光脉搏信号中包括直流成分和交流成分。
在步骤S104,根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小。
具体的,根据比率值R的计算公式:在已知直流成分和交流成分的条件下计算得到比率值R的值。
在步骤S105,根据所述比率值R的大小确定所述采样脉搏血氧值。
在步骤S106,对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果。
在步骤S107,将所述当前血氧测量结果作为血氧饱和度输出。
本发明实施例提供的一种血氧饱和度测量方法可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
参阅图2,图2为本发明实施例一种血氧饱和度测量方法第二实施方式的流程示意图。具体的,图2为图1中对所述红光脉搏信号以及红外光脉搏信号进行滤波处理的流程示意图,具体包括:
在步骤S201,设置滤波带宽对所述红光脉搏信号以及红外光脉搏信号进行滤波。
具体的,设置滤波带宽范围为0.5Hz—7.5Hz。
在步骤S202,根据滤波结果设置红光脉搏信号序列以及红外光脉搏信号序列。
具体的,红光脉搏信号序列为{r1,r2,r3,…,rn},其中n为正整数,红外光脉搏信号序列为{i1,i2,i3,…in}.
经过滤波后的红光脉搏信号序列为{fr1,fr2,fr3,…,frn},经过滤波后的红外光脉搏信号序列为{fi1f,i2,fi3,…fin}。
参阅图3,图3为本发明实施例一种血氧饱和度测量方法第三实施方式的流程示意图,具体的,图3为图1或图2中所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分的流程示意图,具体包括:
在步骤S301,根据原始采样红光脉搏信号、原始采样红外光脉搏信号、滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定直流成分。
具体的,直流成分计算公式为:
其中,r(k)代表红光原始采样信号,fr(k)代表经滤波后的红光输出信号;i(k)代表红外光原始采样信号,fi(k)代表经滤波后的红外光输出信号。
在步骤S302,根据所述滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定交流成分。
具体的,交流成分采用滤波后的输出结果。
进一步地,所述原始采样红光脉搏信号以及原始采样红外光脉搏信号为所述初始采样的红光脉搏信号以及红外光脉搏信号。
参阅图4,图4为本发明实施例一种血氧饱和度测量方法第四实施方式的流程示意图。具体的,图4为图1-3任一项中根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小的流程示意图,具体包括:
在步骤S401,确定红外光脉搏信号的交流成分与红光脉搏信号的交流成分的比值,记做A。
在步骤S402,确定红光脉搏信号的直流成分与红外光脉搏信号的直流成分的比值,记做B。
具体的,
在步骤S403,根据比率值R计算公式确定比率值R,即R=A*B。
进一步地,两路光经骨骼、皮肤及静脉吸收后的输出光的比值即为B,一般情况下人体的骨骼,皮肤及静脉的光吸收特性很少发生变化,因此B值一般可以认为是恒定的,所以R值的大小取决于A值的大小,A值反映是动脉血对两路光吸收的输出结果的比值;
由公式可知,当y不等于0,当Δy<<y,ΔI<<I,可认为上式两边近似相等。
参阅图5,图5为本发明实施例一种血氧饱和度测量方法第五实施方式的流程示意图。具体的,图5为图1-4任一项中根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小的流程示意图,具体包括:
在步骤S501,选取所述红光脉搏信号序列与红外光脉搏信号序列的相同长度的序列。
在步骤S502,确定所述红光脉搏信号序列的最大值以及所述红外光脉搏信号序列的最大值。
在步骤S503,根据递推算法得出比率值R的序列值。
具体的,假定红光交流序列为IrSig={I1,I2,I3,…,IN},红外光交流信号序列为RedSig={R1,R2,R3,…RN},两段信号的长度都为N;从红外光交流值中找到最大值记为Imax,位置标为pmax,因两路光同步采集,所以可以根据pmax位置找到红光序列的最大值,最大值并记为Rmax;两者的比值为
在pmax前一个点处,选择两路光的采样点,记为(Ipmar-1,Rpmar-1),两者的比值为
只要采样率足够高,那么相邻两个采样点的幅度变化一般很小,则可以认为
所以可以得到
再在pmax前两个点处,选择两路光的采样点,记为(Ipmax-2,Rpmax-2),同样可以得到
以此类推,对pmax前N个采样点处,可以得到
Rpmar≈Rpmar-N+1≈…≈Rpar-N (1-3)
由公式(1-1),(1-2)和(1-3),可得出如下关系式
Rpmar≈Rpmar-1≈…≈Rpmar-N
将N个采样点的R值序列记为
(Rpmax,Rpmax-1,…,Rpmax-N) (1-4)。
进一步地,根据人体生理学相关研究表明,人体的血红蛋白和氧合血红蛋白在短时间内一般变化很小,即红光信号和红外光信号的每个采样点对应的血氧基本是相等的;而脉搏波波形的周期变化反映的是透射光通过的路径随心脏收缩和舒张改变而改变,因此可用直方图来统计某段时间内出现相同R值次数最多来计算血氧。
参阅图6,图6为本发明实施例一种血氧饱和度测量方法第六实施方式的流程示意图。图6为图1-5任一项中根据所述比率值R的大小确定所述采样脉搏血氧值,并对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果的流程示意图,具体包括:
在步骤S601,用直方图来统计某段时间内出现相同比率值R值次数最多计算采样脉搏血氧值。
在步骤S602,将所述采样脉搏血氧值在1-100范围内进行频次统计,得到血氧直方图。
在步骤S603,在所述血氧直方图中找出频次最大值,将所述频次最大值作为当前血氧测量结果。
具体的,血氧计算公式采用曲线拟合法得到,假定血氧是R的二次函数,如下所示:Spo2=aR2+bR+c
将(Rpmax,Rpmax-1,…,Rpmax-N)带入(1-4),得出血氧序列如下
(SpO2pmax,SpO2pmax-1,…,SpO2pmax-N)
将所有计算出的血氧值,按照血氧值范围在0至100进行频次统计,可得到血氧直方图,具体包括:
(1)假定采样信号的长度为L,计算L个采样点的红光交流值记为RedAcSig,计算L个采样点的红外光交流值记为IrAcSig;计算L个采样点的红光直流值记为RedDcSig,计算L个采样点的红外光直流值记为IrDcSig;
(2)将上一步计算的红光交流值RedAcSig、红外光交流值IrAcSig、红光直流值RedDcSig及红外光直流值IrDcSig带入公式计算R值;
(3)将(2)计算的R值带入血氧计算公式Spo2=aR2+bR+c,计算出血氧结果记为Spo2Result;
(4)定义一个长度为101的整型数组,记为Array,Array[0]记录血氧等于0的次数,Array[1]记录血氧等于1的次数,Array[n]记录血氧等于n的次数(n小于等于100且大于等于0);
(5)将(3)计算出的血氧结果Spo2Result,按血氧值大小(0~100)存入(4)定义的数组Array;假定Spo2Result(0)等于90,则将Array[90]加1;假定Spo2Result(1)等于91,则将Array[91]加1,以此类推Spo2Result(L)等于n,则将Array[n]加1(n小于等于100且大于等于0);
进一步地,查找直方图的频次最大值,将频次最多点对应的血氧作为当前血氧测量结果的方法,具体包括:
(6)搜索计算出的Array数组极大值记为maxResult,搜索方向为A[100]向Array[0]方向;
(7)保存(6)搜索到极大值maxResult的对应下标(下标范围为0~100)记为outSpo2Array;
(8)从极大值maxResult(可能极大值的个数大于1)中找到的最大值maxv,并保存最大值maxv所对应的下标记为maxind;
(9)从极大值maxResult(可能极大值的个数大于1)中找到的最早的极大值值earlymaxv,并保存最大值earlymaxvmaxv所对应的下标记为earlymaxvmaxind;
(10)如果earlymaxvmaxv与L之比大于20%,则将earlymaxvmaxind(范围为0~100)作为血氧测量结果;如果earlymaxvmaxv与L之比小于20%,则将maxind(范围为0~100)作为血氧测量结果。
综合参阅图2-图6,本发明实施例提供的一种血氧饱和度测量方法可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
参阅图7,图7示出了适于用来实现本发明实施方式的式例性设备的框图。
所述设备包括:处理器(processor)71、存储器(memory)72、通信接口(Communications Interface)73和总线74;其中:
所述处理器71、存储器72、通信接口73通过所述总线74完成相互间的通信;
所述通信接口73用于其他设备之间的信息传输。
所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的血氧饱和度测量方法,具体包括:
采样脉搏血氧信号,并获取红光脉搏信号以及红外光脉搏信号;
对所述红光脉搏信号以及红外光脉搏信号进行滤波处理;
确定红光脉搏信号以及红外光脉搏信号的直流成分和交流成分;
根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小;
根据所述比率值R的大小确定所述采样脉搏血氧值;
对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果;
将所述当前血氧测量结果作为血氧饱和度输出。
进一步地,所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的所述对所述红光脉搏信号以及红外光脉搏信号进行滤波处理,具体包括:
设置滤波带宽对所述红光脉搏信号以及红外光脉搏信号进行滤波;
根据滤波结果设置红光脉搏信号序列以及红外光脉搏信号序列。
进一步地,所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的所述确定所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分,具体包括:
根据原始采样红光脉搏信号、原始采样红外光脉搏信号、滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定直流成分;
根据所述滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定交流成分;
其中,所述原始采样红光脉搏信号以及原始采样红外光脉搏信号为所述初始采样的红光脉搏信号以及红外光脉搏信号。
进一步地,所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的所述根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小,具体包括:
确定红外光脉搏信号的交流成分与红光脉搏信号的交流成分的比值,记做A;
确定红光脉搏信号的直流成分与红外光脉搏信号的直流成分的比值,记做B:
根据比率值R计算公式确定比率值R,即R=A*B。
进一步地,所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的所述根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小,具体还包括:
选取所述红光脉搏信号序列与红外光脉搏信号序列的相同长度的序列;
确定所述红光脉搏信号序列的最大值以及所述红外光脉搏信号序列的最大值;
根据递推算法得出比率值R的序列值。
进一步地,所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的所述根据所述比率值R的大小确定所述采样脉搏血氧值,具体包括:
用直方图来统计某段时间内出现相同比率值R值次数最多计算采样脉搏血氧值。
进一步地,所述处理器71用于调用所述存储器72中的计算机程序,以执行上述方法实施例所提供的所述对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果,具体包括:
将所述采样脉搏血氧值在1-100范围内进行频次统计,得到血氧直方图;
在所述血氧直方图中找出频次最大值,将所述频次最大值作为当前血氧测量结果。
本发明实施例提供的一种设备可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
此外,本发明实施例还提出一种存储介质,所述存储介质上存储有血氧饱和度测量程序,所述血氧饱和度测量程序被处理器执行时实现如下操作:
采样脉搏血氧信号,并获取红光脉搏信号以及红外光脉搏信号;
对所述红光脉搏信号以及红外光脉搏信号进行滤波处理;
确定红光脉搏信号以及红外光脉搏信号的直流成分和交流成分;
根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小;
根据所述比率值R的大小确定所述采样脉搏血氧值;
对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果;
将所述当前血氧测量结果作为血氧饱和度输出。
进一步地,所述血氧饱和度测量程序被处理器执行时实现如下操作:
设置滤波带宽对所述红光脉搏信号以及红外光脉搏信号进行滤波;
根据滤波结果设置红光脉搏信号序列以及红外光脉搏信号序列。
进一步地,所述血氧饱和度测量程序被处理器执行时实现如下操作:
根据原始采样红光脉搏信号、原始采样红外光脉搏信号、滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定直流成分;
根据所述滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定交流成分;
其中,所述原始采样红光脉搏信号以及原始采样红外光脉搏信号为所述初始采样的红光脉搏信号以及红外光脉搏信号。
进一步地,所述血氧饱和度测量程序被处理器执行时实现如下操作:
确定红外光脉搏信号的交流成分与红光脉搏信号的交流成分的比值,记做A;
确定红光脉搏信号的直流成分与红外光脉搏信号的直流成分的比值,记做B:
根据比率值R计算公式确定比率值R,即R=A*B。
进一步地,所述血氧饱和度测量程序被处理器执行时实现如下操作:
选取所述红光脉搏信号序列与红外光脉搏信号序列的相同长度的序列;
确定所述红光脉搏信号序列的最大值以及所述红外光脉搏信号序列的最大值;
根据递推算法得出比率值R的序列值。
进一步地,所述血氧饱和度测量程序被处理器执行时实现如下操作:
用直方图来统计某段时间内出现相同比率值R值次数最多计算采样脉搏血氧值。
进一步地,所述血氧饱和度测量程序被处理器执行时实现如下操作:
将所述采样脉搏血氧值在1-100范围内进行频次统计,得到血氧直方图;
在所述血氧直方图中找出频次最大值,将所述频次最大值作为当前血氧测量结果。
本发明实施例提供的一种存储介质可准确得到脉搏血氧值,且同时具有抗噪声、计算量小等优点,且可适用于可穿戴式设备。
随着科技的发展,计算机程序的传播途径不再受限于有形介质,还可以直接从网络下载,或者采用其他方式获取。因此,本实施例中的计算机可读介质不仅可以包括有形的介质,还可以包括无形的介质。
本实施例的计算机存储介质可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括——但不限于——电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于——无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码,所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
在本发明所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明个实施例中的功能单元可以集成在一个处理单元中,也可以是个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元,可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中,包括若干指令用以使得一台计算机售药机(可以是个人计算机,服务器,或者网络售药机等)或处理器(processor)执行本发明个实施例所述方法的部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、磁碟或者光盘等种可以存储程序代码的介质。
上述实施方式仅为本发明的优选实施方式,不能以此来限定本发明保护的范围,本领域的技术人员在本发明的基础上所做的任何非实质性的变化及替换均属于本发明所要求保护的范围。
Claims (9)
1.一种血氧饱和度测量方法,其特征在于,包括:
采样脉搏血氧信号,并获取红光脉搏信号以及红外光脉搏信号;
对所述红光脉搏信号以及红外光脉搏信号进行滤波处理;
确定红光脉搏信号以及红外光脉搏信号的直流成分和交流成分;
根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小;
根据所述比率值R的大小确定所述采样脉搏血氧值;
对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果;
将所述当前血氧测量结果作为血氧饱和度输出。
2.根据权利要求1所述的一种血氧饱和度测量方法,其特征在于,所述对所述红光脉搏信号以及红外光脉搏信号进行滤波处理,具体包括:
设置滤波带宽对所述红光脉搏信号以及红外光脉搏信号进行滤波;
根据滤波结果设置红光脉搏信号序列以及红外光脉搏信号序列。
3.根据权利要求2所述的一种血氧饱和度测量方法,其特征在于,所述确定所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分,具体包括:
根据原始采样红光脉搏信号、原始采样红外光脉搏信号、滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定直流成分;
根据所述滤波后的红光脉搏信号以及滤波后的红外光脉搏信号确定交流成分;
其中,所述原始采样红光脉搏信号以及原始采样红外光脉搏信号为所述初始采样的红光脉搏信号以及红外光脉搏信号。
4.根据权利要求3所述的一种血氧饱和度测量方法,其特征在于,所述根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小,具体包括:
确定红外光脉搏信号的交流成分与红光脉搏信号的交流成分的比值,记做A;
确定红光脉搏信号的直流成分与红外光脉搏信号的直流成分的比值,记做B;
根据比率值R计算公式确定比率值R,即R=A*B。
5.根据权利要求4所述的一种血氧饱和度测量方法,其特征在于,所述根据所述红光脉搏信号以及红外光脉搏信号的直流成分和交流成分确定比率值R的大小,具体还包括:
选取所述红光脉搏信号序列与红外光脉搏信号序列的相同长度的序列;
确定所述红光脉搏信号序列的最大值以及所述红外光脉搏信号序列的最大值;
根据递推算法得出比率值R的序列值。
6.根据权利要求5所述的一种血氧饱和度测量方法,其特征在于,所述根据所述比率值R的大小确定所述采样脉搏血氧值,具体包括:
用直方图来统计某段时间内出现相同比率值R值次数最多计算采样脉搏血氧值。
7.根据权利要求6所述的一种血氧饱和度测量方法,其特征在于,所述对所述采样脉搏血氧值进行频次统计,得到频次最大值,并将所述频次最大值作为当前血氧测量结果,具体包括:
将所述采样脉搏血氧值在1-100范围内进行频次统计,得到血氧直方图;
在所述血氧直方图中找出频次最大值,将所述频次最大值作为当前血氧测量结果。
8.一种设备,其特征在于,所述设备包括:存储器、处理器及存储在所述存储器上并可在所述处理器上运行的血氧饱和度测量程序,所述血氧饱和度测量程序被所述处理器执行时实现如权利要求1至7中任一项所述的血氧饱和度测量方法的步骤。
9.一种存储介质,其特征在于,所述存储介质上存储有血氧饱和度测量程序,所述血氧饱和度测量程序被处理器执行时实现如权利要求1至7中任一项所述的血氧饱和度测量方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810181589.8A CN108652642A (zh) | 2018-03-06 | 2018-03-06 | 一种血氧饱和度测量方法、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810181589.8A CN108652642A (zh) | 2018-03-06 | 2018-03-06 | 一种血氧饱和度测量方法、设备及存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108652642A true CN108652642A (zh) | 2018-10-16 |
Family
ID=63784244
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810181589.8A Pending CN108652642A (zh) | 2018-03-06 | 2018-03-06 | 一种血氧饱和度测量方法、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108652642A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109497973A (zh) * | 2018-12-24 | 2019-03-22 | 兰州理工大学 | 日常无监督状态下脉搏和血氧信息的检测系统及检测方法 |
CN109924987A (zh) * | 2019-03-14 | 2019-06-25 | 浙江荷清柔性电子技术有限公司 | 反射式血氧仪的标定方法、系统和可读存储介质 |
CN111685726A (zh) * | 2019-03-12 | 2020-09-22 | 华为终端有限公司 | 一种血氧检测方法及装置 |
CN112741627A (zh) * | 2021-01-29 | 2021-05-04 | 南开大学 | 一种基于多重滤波的手腕血氧饱和度监测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1985764A (zh) * | 2005-12-23 | 2007-06-27 | 深圳迈瑞生物医疗电子股份有限公司 | 能消除运动干扰的血氧测量方法及其装置 |
US20100145170A1 (en) * | 2006-09-21 | 2010-06-10 | Starr Life Sciences Corp. | Small Animal Pulse Oximeter User Interface |
CN101933811A (zh) * | 2010-09-03 | 2011-01-05 | 深圳市纽泰克电子有限公司 | 一种血氧饱和度检测方法及系统 |
CN102512178A (zh) * | 2011-12-23 | 2012-06-27 | 深圳市理邦精密仪器股份有限公司 | 一种血氧测量装置和方法 |
US8594759B2 (en) * | 2009-07-30 | 2013-11-26 | Nellcor Puritan Bennett Ireland | Systems and methods for resolving the continuous wavelet transform of a signal |
CN103549942A (zh) * | 2013-10-31 | 2014-02-05 | 广州视源电子科技股份有限公司 | 通过光信号识别脉率和血氧饱和度的方法 |
CN105433940A (zh) * | 2016-01-02 | 2016-03-30 | 无锡桑尼安科技有限公司 | 一种智能化多参数检测方法 |
CN107456219A (zh) * | 2017-09-07 | 2017-12-12 | 成都云卫康医疗科技有限公司 | 一种基于皮尔逊相关系数动态心率和血氧测量方法 |
-
2018
- 2018-03-06 CN CN201810181589.8A patent/CN108652642A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1985764A (zh) * | 2005-12-23 | 2007-06-27 | 深圳迈瑞生物医疗电子股份有限公司 | 能消除运动干扰的血氧测量方法及其装置 |
US20100145170A1 (en) * | 2006-09-21 | 2010-06-10 | Starr Life Sciences Corp. | Small Animal Pulse Oximeter User Interface |
US8594759B2 (en) * | 2009-07-30 | 2013-11-26 | Nellcor Puritan Bennett Ireland | Systems and methods for resolving the continuous wavelet transform of a signal |
CN101933811A (zh) * | 2010-09-03 | 2011-01-05 | 深圳市纽泰克电子有限公司 | 一种血氧饱和度检测方法及系统 |
CN102512178A (zh) * | 2011-12-23 | 2012-06-27 | 深圳市理邦精密仪器股份有限公司 | 一种血氧测量装置和方法 |
CN103549942A (zh) * | 2013-10-31 | 2014-02-05 | 广州视源电子科技股份有限公司 | 通过光信号识别脉率和血氧饱和度的方法 |
CN105433940A (zh) * | 2016-01-02 | 2016-03-30 | 无锡桑尼安科技有限公司 | 一种智能化多参数检测方法 |
CN107456219A (zh) * | 2017-09-07 | 2017-12-12 | 成都云卫康医疗科技有限公司 | 一种基于皮尔逊相关系数动态心率和血氧测量方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109497973A (zh) * | 2018-12-24 | 2019-03-22 | 兰州理工大学 | 日常无监督状态下脉搏和血氧信息的检测系统及检测方法 |
CN111685726A (zh) * | 2019-03-12 | 2020-09-22 | 华为终端有限公司 | 一种血氧检测方法及装置 |
CN109924987A (zh) * | 2019-03-14 | 2019-06-25 | 浙江荷清柔性电子技术有限公司 | 反射式血氧仪的标定方法、系统和可读存储介质 |
CN112741627A (zh) * | 2021-01-29 | 2021-05-04 | 南开大学 | 一种基于多重滤波的手腕血氧饱和度监测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108652642A (zh) | 一种血氧饱和度测量方法、设备及存储介质 | |
CN103186406B (zh) | 用于控制流分析的方法和装置 | |
CN107577773A (zh) | 一种音频匹配方法与装置、电子设备 | |
Reddy et al. | ECG data compression using Fourier descriptors | |
Zhang et al. | VLSI friendly ECG QRS complex detector for body sensor networks | |
CN102499694A (zh) | 用于消除血氧饱和度监测干扰的方法 | |
Lazzara et al. | Functional anatomy of the canine left bundle branch | |
CN105147251A (zh) | 基于多通道sEMG的肌肉疲劳动态预测方法 | |
CN107870316B (zh) | 一种基于时差计算的tdoa定位方法、装置及系统 | |
CN105852850A (zh) | 一种从心电信号中提取呼吸信号的方法及相关装置 | |
CN107510454A (zh) | 基于多阶滤波的肌电信号采集装置及系统 | |
CN109189207A (zh) | 一种可穿戴设备穿戴状态的检测方法及装置 | |
CN108294738A (zh) | 一种血氧频域抗干扰测量方法、设备及存储介质 | |
CN106354386A (zh) | 利用生理信号进行交互的电子设备及方法 | |
CN108652598A (zh) | 提高无创血压测量准确度的方法、设备及存储介质 | |
Kunzmann et al. | Parameter extraction of ECG signals in real-time | |
CN108836305B (zh) | 一种融合巴特沃斯滤波和小波变换的ecg特征提取方法 | |
Chowdhury et al. | Hardware/software co-design of an ECG-PPG preprocessor: A qualitative & quantitative analysis | |
Sasi et al. | FPGA implementation of DC bias removal filters–a case study with an ECG signal | |
Kohama et al. | An efficient RR interval detection for ECG monitoring system | |
KR102210468B1 (ko) | 심전도 신호의 기저선 변동 잡음 제거 방법 및 컴퓨터 프로그램 | |
CN208851487U (zh) | 基于多阶滤波的肌电信号采集装置及系统 | |
El-Mora et al. | Application of particle swarm optimization algorithm to multiuser detection in CDMA | |
Hui et al. | Revisiting FPGA implementation of digital filters and exploring approximate computing on biomedical signals | |
CN112054854A (zh) | 基于图着色方法的无线体域网间干扰减轻方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181016 |
|
RJ01 | Rejection of invention patent application after publication |