CN103006238B - 一种测量血氧饱和度的方法 - Google Patents
一种测量血氧饱和度的方法 Download PDFInfo
- Publication number
- CN103006238B CN103006238B CN201310004409.6A CN201310004409A CN103006238B CN 103006238 B CN103006238 B CN 103006238B CN 201310004409 A CN201310004409 A CN 201310004409A CN 103006238 B CN103006238 B CN 103006238B
- Authority
- CN
- China
- Prior art keywords
- absorbance
- oxygen saturation
- blood oxygen
- max
- relative value
- 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
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了一种测量血氧饱和度的方法,通过彩色图像传感器获取透过手指或手掌的每帧RGB图像数据,对每帧RGB图像数据进行三原色分离并累加三原色的强度值,连续测量多帧RGB图像数据,分别形成R、G和B光强值时间序列;从R、G和B光强值时间序列中分别获取光电容积脉搏波,计算光电容积脉搏波的最大值和最小值;通过光电容积脉搏波的最大值和最小值分别计算R、G和B的吸光度值;通过R、G和B的吸光度值获取血氧饱和度。本方法可以不增加和改造任何硬件电路与设备,将手机或照相机等转换成血氧饱和度测试仪,不仅实现一机多能,更能方便、快捷、迅速的测量血氧饱和度,降低了测试成本。
Description
技术领域
本发明属于光电技术与信号处理领域,特别涉及一种测量血氧饱和度的方法。
背景技术
传统的基于光电容积脉搏波的血氧饱和度测量方法,通常采用专门的红光和红外光二极管以及光电传感器来获取光电容积脉搏波,进而通过光谱来得到血氧饱和度。
发明人在实现本发明的过程中,发现现有技术中至少存在以下缺点和不足:
上述操作的过程中,不仅需要专门的电路,而且对器件参数和测量条件均有较严格的要求,增加了测量的成本,限制了实际应用中的范围。
发明内容
本发明提供了一种测量血氧饱和度的方法,本方法降低了测量的成本,实现了可以实时、方便和快捷地对血氧饱和度进行测量,详见下文描述:
一种测量血氧饱和度的方法,所述方法包括以下步骤:
(1)通过彩色图像传感器获取透过手指或手掌的每帧RGB图像数据,对每帧RGB图像数据进行三原色分离并累加三原色的强度值,连续测量多帧RGB图像数据,分别形成R、G和B光强值时间序列;
(2)从所述R、G和B光强值时间序列中分别获取光电容积脉搏波,计算所述光电容积脉搏波的最大值和最小值;
(3)通过所述光电容积脉搏波的最大值和最小值分别计算R、G和B的吸光度值;
(4)通过所述R、G和B的吸光度值获取血氧饱和度。
所述通过所述光电容积脉搏波的最大值和最小值分别计算R、G和B的吸光度值具体为:
对Rmax、Rmin、Gmax、Gmin、Bmax和Bmin分别取对数并计算该颜色的吸光度值,即:
AR=lnRmax-lnRmin;
AG=lnGmax-lnGmin;
AB=lnBmax-lnBmin;
其中,Rmax、Rmin、Gmax、Gmin、Bmax和Bmin分别为R、G和B的光电容积脉搏波最大值和最小值;AR、AG和AB分别为R、G和B的吸光度值。
所述通过所述R、G和B的吸光度值获取血氧饱和度具体包括:
1)通过R、G和B的吸光度值获取G和R的吸光度相对值、B和R的吸光度相对值;
2)通过所述G和R的吸光度相对值、所述B和R的吸光度相对值获取吸光度比值Q,通过所述吸光度比值Q获取血氧饱和度。
所述通过所述R、G和B的吸光度值获取血氧饱和度具体包括:
1)通过R、G和B的吸光度值获取G和R的吸光度相对值、B和R的吸光度相对值;
2)通过所述G和R的吸光度相对值AG/AR、所述B和R的吸光度相对值AB/AR构造成一个2维表格,依据定标得到与血氧饱和度值的关系,获取血氧饱和度值。
本发明提供的技术方案的有益效果是:通过采用彩色图像传感器得到包含光电容积脉搏波的RGB图像数据,从RGB图像数据中获取到G和R的吸光度相对值、B和R的吸光度相对值,进而得到血氧饱和度值。本方法可以不增加和改造任何硬件电路与设备,将手机或照相机等转换成血氧饱和度测试仪,不仅实现一机多能,更能方便、快捷、迅速的测量血氧饱和度,降低了测试成本。
附图说明
图1为本发明提供的一种测量血氧饱和度的方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,兹例举以下实施例,并配合附图详细说明如下:
为了降低测量的成本,实现可以实时、方便和快捷地对血氧饱和度进行测量,本发明实施例提出了一种测量血氧饱和度的方法,参见图1,该方法包括以下步骤:
101:通过彩色图像传感器获取透过手指或手掌的每帧RGB图像数据,对每帧RGB图像数据进行三原色分离并累加三原色的强度值,连续测量多帧RGB图像数据,分别形成R、G和B光强值时间序列[Ri]、[Gi]和[Bi];
具体实现时,将手指或手掌覆盖在手机或照相机等(即含有彩色图像传感器的任意装置均可,本发明实施例对此不做限制)装置的镜头前,通过彩色图像传感器获取透过手指或手掌的每帧RGB图像数据。
对每帧RGB图像数据进行R、G和B的分离,累加R、G和B三原色的强度值。
其中,i代表时间,k,j分别代表横纵坐标,Rj,k表示R分量中的第k行第j个点;Gj,k表示G分量中的第k行第j个点;Bj,k表示B分量中的第k行第j个点。
测量的RGB图像数据的个数根据实际应用中的需要进行设定,本发明实施例对此不做限制。
102:从R、G和B光强值时间序列[Ri]、[Gi]和[Bi]中分别获取光电容积脉搏波,计算光电容积脉搏波的最大值和最小值;
由上述公式得到的R、G和B光强值时间序列[Ri]、[Gi]和[Bi]中,必定包含被测手指的光电容积脉搏波,分别计算R、G和B的光电容积脉搏波中的最大值和最小值:
Rmax=MAX[Ri] (4)
Rmin=MIN[Ri] (5)
Gmax=MAX[Gi] (6)
Gmin=MIN[Gi] (7)
Bmax=MAX[Bi] (8)
Bmin=MIN[Bi] (9)
103:通过光电容积脉搏波的最大值和最小值分别计算R、G和B的吸光度值;
该步骤具体为:对Rmax、Rmin、Gmax、Gmin、Bmax和Bmin分别取对数并计算该颜色的吸光度值,即:
AR=lnRmax-lnRmin (10)
AG=lnGmax-lnGmin (11)
AB=lnBmax-lnBmin (12)
104:通过R、G和B的吸光度值获取血氧饱和度。
其中,该步骤具体为:
1)通过R、G和B的吸光度值获取G和R的吸光度相对值、B和R的吸光度相对值;
即,将AG除以AR得到G和R的吸光度相对值AG/AR,将AB除以AR得到B和R的吸光度相对值AB/AR,将R、G和B三个分量转化为G和B两个分量。
2)通过G和R的吸光度相对值、B和R的吸光度相对值获取吸光度比值Q,通过吸光度比值Q获取血氧饱和度。
即,吸光度比值Q等于B和R的吸光度相对值AB/AR与G和R的吸光度相对值AG/AR的比值,
Q=(AB/AR)/(AG/AR) (13)
吸光度比值Q值与血氧饱和度的关系可以通过大样本统计得到,本发明实施例在此不做赘述。
具体实现时,步骤2)还可以为:通过G和R的吸光度相对值AG/AR、B和R的吸光度相对值AB/AR构造成一个2维表格,依据定标得到与血氧饱和度值的关系,可得到血氧饱和度值。
例如:采用一个高精度的血氧仪测得同一人手指的精准血氧饱和度值并测量相同时刻的吸光度相对值AB/AR与AG/AR,得到AB/AR与AG/AR与精准血氧饱和度之间的定标关系,再一次测量时可根据定标关系得到的血氧饱和度值。
再如:直接利用血氧饱和度定标仪得到AB/AR与AG/AR与精准血氧饱和度之间的定标关系,在测量时就可根据定标关系得到的血氧饱和度值。
综上所述,本方法可以不增加和改造任何硬件电路与设备,将手机或照相机等转换成血氧饱和度测试仪,不仅实现一机多能,更能方便、快捷、迅速的测量血氧饱和度,降低了测试成本。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种测量血氧饱和度的方法,其特征在于,所述方法包括以下步骤:
(1)通过彩色图像传感器获取透过手指或手掌的每帧RGB图像数据,对每帧RGB图像数据进行三原色分离并累加三原色的强度值,连续测量多帧RGB图像数据,分别形成R、G和B光强值时间序列;
(2)从所述R、G和B光强值时间序列中分别获取光电容积脉搏波,计算所述光电容积脉搏波的最大值和最小值;
(3)通过所述光电容积脉搏波的最大值和最小值分别计算R、G和B的吸光度值;
(4)通过所述R、G和B的吸光度值获取血氧饱和度。
2.根据权利要求1所述的一种测量血氧饱和度的方法,其特征在于,所述通过所述光电容积脉搏波的最大值和最小值分别计算R、G和B的吸光度值具体为:
对Rmax、Rmin、Gmax、Gmin、Bmax和Bmin分别取对数并计算该颜色的吸光度值,即:
AR=lnRmax-lnRmin;AG=lnGmax-lnGmin;AB=lnBmax-lnBmin;
其中,Rmax、Rmin、Gmax、Gmin、Bmax和Bmin分别为R、G和B的光电容积脉搏波最大值和最小值;AR、AG和AB分别为R、G和B的吸光度值。
3.根据权利要求1所述的一种测量血氧饱和度的方法,其特征在于,所述通过所述R、G和B的吸光度值获取血氧饱和度具体包括:
1)通过R、G和B的吸光度值获取G和R的吸光度相对值、B和R的吸光度相对值;
2)通过所述G和R的吸光度相对值、所述B和R的吸光度相对值获取吸光度比值Q,通过所述吸光度比值Q获取血氧饱和度。
4.根据权利要求1所述的一种测量血氧饱和度的方法,其特征在于,所述通过所述R、G和B的吸光度值获取血氧饱和度具体包括:
1)通过R、G和B的吸光度值获取G和R的吸光度相对值、B和R的吸光度相对值;
2)通过所述G和R的吸光度相对值AG/AR、所述B和R的吸光度相对值AB/AR构造成一个2维表格,依据定标得到与血氧饱和度值的关系,获取血氧饱和度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310004409.6A CN103006238B (zh) | 2013-01-07 | 2013-01-07 | 一种测量血氧饱和度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310004409.6A CN103006238B (zh) | 2013-01-07 | 2013-01-07 | 一种测量血氧饱和度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103006238A CN103006238A (zh) | 2013-04-03 |
CN103006238B true CN103006238B (zh) | 2014-11-26 |
Family
ID=47955729
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310004409.6A Expired - Fee Related CN103006238B (zh) | 2013-01-07 | 2013-01-07 | 一种测量血氧饱和度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103006238B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6455994B2 (ja) * | 2014-06-30 | 2019-01-23 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | フォトプレチスモグラフィセンサ装置及び方法 |
CN105030277B (zh) * | 2015-03-03 | 2018-02-09 | 青岛美诺泰科信息技术有限公司 | 一种基于智能终端的监测装置 |
CN107072488B (zh) * | 2015-09-24 | 2019-04-19 | Hoya株式会社 | 分析装置 |
CN109247944A (zh) * | 2018-08-30 | 2019-01-22 | 合肥工业大学 | 一种基于低端彩色相机的非接触式血氧饱和度检测方法 |
CN109363630B (zh) * | 2018-09-03 | 2021-08-27 | 浙江大华技术股份有限公司 | 一种生命体征信息测量方法及装置 |
US11375909B2 (en) | 2018-12-05 | 2022-07-05 | Boe Technology Group Co., Ltd. | Method and apparatus for determining physiological parameters of a subject, and computer-program product thereof |
CN114246562A (zh) * | 2021-12-06 | 2022-03-29 | 北京极豪科技有限公司 | 指纹模组、生理特征检测方法和装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7711403B2 (en) * | 2001-04-05 | 2010-05-04 | Rhode Island Hospital | Non-invasive determination of blood components |
JP5029150B2 (ja) * | 2007-06-06 | 2012-09-19 | ソニー株式会社 | 生体情報取得装置および生体情報取得方法 |
US8150501B2 (en) * | 2009-03-26 | 2012-04-03 | Johnson & Johnson Consumer Companies, Inc. | Method for measuring skin erythema |
LV14514B (lv) * | 2010-10-06 | 2012-08-20 | Latvijas Universitāte | Iekārta un metode sirdsdarbības parametru optiskai bezkontakta kontrolei |
JP5303012B2 (ja) * | 2010-12-15 | 2013-10-02 | 富士フイルム株式会社 | 内視鏡システム、内視鏡システムのプロセッサ装置及び内視鏡システムの作動方法 |
-
2013
- 2013-01-07 CN CN201310004409.6A patent/CN103006238B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103006238A (zh) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103006238B (zh) | 一种测量血氧饱和度的方法 | |
US11193830B2 (en) | Spectrocolorimeter imaging system | |
CN103824544B (zh) | Led显示屏的校正方法、装置及系统 | |
US9514535B1 (en) | Color calibration method of camera module | |
US9542758B2 (en) | Image processing apparatus and light source identification method | |
US8587783B2 (en) | Spectral estimation parameter generation device, method of generating spectral estimation parameter, and spectral estimation device | |
CN108737815B (zh) | 一种图像传感器的质量检测方法及系统 | |
CN112073703A (zh) | 一种色彩校正矩阵的调整方法、装置、终端设备及介质 | |
US10835135B2 (en) | Non-contact heartbeat rate measurement system, method and apparatus thereof | |
CN109253862A (zh) | 一种基于神经网络的色度测量方法 | |
CN111466605A (zh) | 基于rgb值测定烟叶烘烤变黄变褐程度的方法、系统及装置 | |
CN109459136B (zh) | 一种色度测量的方法与装置 | |
CN108259865A (zh) | 一种基于单像素探测器的彩色成像方法和系统 | |
CN103340604A (zh) | 黄疸自动分析检测方法及计算机程序产品 | |
CN106162161A (zh) | 一种图像检测方法及显示装置 | |
CN106153193B (zh) | 一种利用多光谱相机双光源响应值获取光谱反射比的方法 | |
CN109887044A (zh) | 一种生殖数据测评方法及系统 | |
US20140218571A1 (en) | Image processing method capable of reducing image noise | |
Koren | Measuring camera shannon information capacity with a siemens star image | |
Nixon et al. | The importance of a device specific calibration for smartphone colorimetry | |
WO2016169055A1 (zh) | 光学相机亮度检测值与色度检测值校准方法及装置 | |
TWI512274B (zh) | 影像式彩色分析儀之校正方法及量測裝置 | |
KR20170100717A (ko) | 분광반사율 추정을 이용한 피부 분석 장치 및 방법 | |
KR101541977B1 (ko) | 표준 색상표 없이 색상 자가 진단이 가능한 스마트기기 및 이를 이용한 진단 방법 | |
CN108362650B (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: 20141126 Termination date: 20220107 |
|
CF01 | Termination of patent right due to non-payment of annual fee |