CN114376571B - 基于单色光动脉血管图像的血液含氧量状况测量方法及装置 - Google Patents
基于单色光动脉血管图像的血液含氧量状况测量方法及装置 Download PDFInfo
- Publication number
- CN114376571B CN114376571B CN202111491945.4A CN202111491945A CN114376571B CN 114376571 B CN114376571 B CN 114376571B CN 202111491945 A CN202111491945 A CN 202111491945A CN 114376571 B CN114376571 B CN 114376571B
- Authority
- CN
- China
- Prior art keywords
- arterial
- image sequence
- arterial blood
- image
- vessel
- 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.)
- Active
Links
- 239000008280 blood Substances 0.000 title claims abstract description 82
- 210000004369 blood Anatomy 0.000 title claims abstract description 82
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 title claims abstract description 77
- 239000001301 oxygen Substances 0.000 title claims abstract description 77
- 229910052760 oxygen Inorganic materials 0.000 title claims abstract description 77
- 238000000691 measurement method Methods 0.000 title description 4
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 91
- 230000003287 optical effect Effects 0.000 claims abstract description 45
- 238000000034 method Methods 0.000 claims abstract description 38
- 108010064719 Oxyhemoglobins Proteins 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 33
- 230000010349 pulsation Effects 0.000 claims abstract description 28
- 230000008602 contraction Effects 0.000 claims abstract description 8
- 108010054147 Hemoglobins Proteins 0.000 claims description 47
- 102000001554 Hemoglobins Human genes 0.000 claims description 47
- 238000004364 calculation method Methods 0.000 claims description 37
- 238000013139 quantization Methods 0.000 claims description 16
- 230000024883 vasodilation Effects 0.000 claims description 16
- 230000010339 dilation Effects 0.000 claims description 15
- 206010047139 Vasoconstriction Diseases 0.000 claims description 9
- 230000025033 vasoconstriction Effects 0.000 claims description 9
- 238000002835 absorbance Methods 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 5
- 238000011002 quantification Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 6
- 230000005540 biological transmission Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000005498 polishing Methods 0.000 description 2
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010241 blood sampling Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000840 electrochemical analysis Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 230000005541 medical transmission Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000006213 oxygenation reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
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
- A61B5/14551—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 for measuring blood gases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明公开了基于单色光动脉血管图像的血液含氧量状况测量方法及装置,其包括以下步骤:1)将待采集部位放入采集模块中,采集模块的上方设有单色光二极管,下方设有摄像装置,通过采集模块获取多组单色光透射后的光信号,并将多组光信号转化为多组图像数字信号;2)对获得的所有图像数字信号进行图像增强处理和二值化处理;3)根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程;4)根据动脉血管的扩张光程和收缩光程计算血氧饱和度和氧合血红蛋白的浓度,反应人体的血液氧含量状况。本发明为血氧检测设备的小型化提供了可能,为无创氧合血红蛋白浓度的检测提供了新的检测方案。
Description
技术领域
本发明涉及图像识别及处理技术领域,尤其涉及一种基于单色光动脉血管图像的血液含氧量状况测量方法及装置。
背景技术
无创血氧饱和度的测量原理是动脉搏动引起血液对光的吸收的变化,并将光的变化转换为电信号,并对电信号进行处理,进而可计算出血氧饱和度。血氧饱和度一般用血液中氧合血红蛋白占总血红蛋白的百分比来表示,一般通过监测动脉的功能血氧饱和度(SpO2)来反应人体的血液氧含量状况,其表达式为:
其中,SpO2为血氧饱和度;为动脉血管中氧合血红蛋白的浓度;CHb为动脉血管中还原血红蛋白的浓度。
目前医用的透射式无创血氧饱和度的测量方案是采用双波长异步打光方案,如图1所示,手指上壁固定了两个并排放置的发光二极管,波长分别为660nm的红光和904nm的红外光,下壁是一个光电转换器件,将透射光的光强转换为电流信号。通过控制两个LED轮换发光,并对各自发光时输出的电流信号进行I/V变换,放大然后采集,对采集后的信号做处理,得到对应的脉搏波数据,最后计算得到SpO2的值。该方案获得血氧饱和度的方法如中国专利CN112903609A公开的无需校正的双波长静脉血氧饱和度测量方法,包括以下步骤:通过DAC控制红光与红外光二极管发光,使用光频传感器接收透射光,将透光信号转换为对应的频率信号f,然后将得到的频率数据f和f′通过MCU上传到计算机中,再经过推导得到的公式计算即可得出血氧饱和度。
然而,现有的血氧饱和度测量技术需要采用双波长异步打光方案,同时需要复杂的光电转换装置和信号处理单元,计算过程比较繁琐,计算精度不高。采集模块需要复杂的光电转换装置,设备难以做到小型化,设备成本较高。
发明内容
本发明的目的是提供一种基于单色光动脉血管图像的血液含氧量状况测量方法及装置,以解决光信号处理时计算过程繁琐、精度不高、采集模块成本较高等问题。
为了达到目的,本发明提供的技术方案为:
本发明涉及的一种基于单色光动脉血管图像的血液含氧量状况测量方法,其包括以下步骤:
1)将待采集部位放入采集模块中,该采集模块包括光源和摄像装置,所述光源为单波长光源,所述摄像装置采集光透射待采集部位后的光信号,并将所述光信号转化为图像数字信号;
2)对获得的所有图像数字信号进行图像增强处理和二值化处理;
3)根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程;
4)根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度及氧合血红蛋白的浓度,用于反应人体的血液氧含量状况。
优选地,所述的步骤2)中对获得的某一图像数字信号进行图像增强处理和二值化处理的具体步骤包括:
2.1)计算某一图像序列中的每一像素的曲率值,曲率计算公式如下:
其中,
公式中,Cij为某一图像序列的第i行第j列像素的曲率值;Imgij为某一图像序列的第i行第j列像素的灰度值;pij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导;qij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导;rij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导;tij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导;sij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导;
2.2)计算某一图像序列中的全部像素的曲率均值,曲率均值的计算公式如下:
其中,Cavg为某一图像序列中的全部像素的曲率均值;r为某一图像序列的图像行数;c为某一图像序列的图像列数;
2.3)对动脉血管的二值化处理,二值化计算方式为:
其中,fij为某一图像序列的第i行第j列像素的二值化的数值;α为曲率均值二值化系数,取值范围为[0,1]。
优选地,所述的曲率均值二值化系数α的取值为0.5。
优选地,所述的步骤3)中是根据增强处理和二值化处理后图像数字信号的二值化有效面积来计算扩张光程和收缩光程的,具体计算公式为:
Smax=max{S1,S2,S3,...,SA} (10)
Smin=min{S1,S2,S3,...,SA} (11)
Lmax=Smax (13)
Lmin=Smin (16)
公式中,Sm为第m个图像序列的动脉血管二值化有效面积;为第m个图像序列的第i行第j列像素的二值化数值;Smax为动脉血管扩张时的二值化有效面积;Smin为动脉血管收缩时的二值化有效面积;A为获取的图像数字信号的数量;Lmax为动脉血管扩张时的光程;Lmin为动脉血管收缩时的光程。此步骤有利于降低动脉血氧饱和度的计算误差。
优选地,所述的步骤4)根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度的具体步骤包括:
4.1)计算动脉血管扩张时的透射光的强度,计算公式如下:
公式中,IA为动脉血管扩张时的透射光的强度;I0为动脉血管扩张时的入射光的强度;γ为组织内的非脉动成分及静脉血的吸光系数;εHbO2为表示动脉血管中氧合血红蛋白的吸光系数;CHbO2为表示动脉血管中氧合血红蛋白的浓度;εHb为表示动脉血管中还原血红蛋白的吸光系数;CHb为表示动脉血管中还原血红蛋白的浓度;为动脉扩张图像序列的第i行第j列像素的灰度值;/>为动脉扩张图像序列的第i行第j列像素的二值化的数值;r1为动脉扩张图像序列的图像行数;c1为动脉扩张图像序列的图像列数;
4.2)计算动脉血管收缩时的透射光的强度,计算公式如下:
公式中,IB为动脉血管收缩时的透射光的强度;I0’为动脉血管收缩时的入射光的强度;为动脉收缩图像序列的第i行第j列像素的灰度值;/>为动脉收缩图像序列的第i行第j列像素的二值化的数值;r2为动脉收缩图像序列的图像行数;c2为动脉收缩图像序列的图像列数;
4.3)对公式(12)~(14)进行整合和取对数,得到公式(18):
4.4)对公式(15)~(17)进行整合和取对数,得到公式(19):
4.5)对公式(18)~(19)进行整合,得到公式(20)和公式(21),用于计算动脉血管的搏动信息的量化参数:
CHbo2εHbO2+CHbεHb=R (20)
公式中,R为动脉血管的搏动信息的量化参数,通过公式(21)计算具体值;
4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度,计算公式为:
其中,C0为表示动脉血管中血红蛋白的浓度:
CHbO2+CHb=C0 (23)
公式中,CHbO2为表示动脉血管中氧合血红蛋白的浓度;CHb为表示动脉血管中还原血红蛋白的浓度;
4.7)基于动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度,计算公式为:
其中,C0为动脉血管中血红蛋白的浓度与动脉血管的搏动信息的量化参数R存在函数关系,函数关系为:
R=f(C0) (25)
公式中,f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0的函数关系。
优选地,所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0存在线性关系,即:
R=kC0+B (26)
其中,k为线性系数;B为漂移常量;
将公式(26)代入公式(22)和公式(24),可得到动脉血管中氧合血红蛋白的浓度和动脉血氧饱和度的计算公式:
优选地,所述的步骤1)中单色光二极管发出的单色光的波长为800nm~950nm。
优选地,所述的步骤1)中摄像装置的帧率为60帧/s,手指放置时间为至少1s,获得的图像数字信号至少为60张。
本发明还涉及基于单色光动脉血管图像的血液含氧量状况测量装置,其包括:
采集模块,该采集模块包括光源和摄像装置,所述光源为单波长光源,所述摄像装置采集光透射待采集部位后的光信号,并将所述光信号转化为图像数字信号;
图像处理模块,用于对获得的所有图像数字信号进行图像增强处理和二值化处理;
光程计算模块,用于根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程;
动脉血氧饱和度及氧合血红蛋白浓度计算模块,用于根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度和氧合血红蛋白浓度,用于反应人体的血液氧含量状况。
采用本发明提供的技术方案,与现有技术相比,具有如下有益效果:
1、本发明采用一种基于单色光的动脉血管图像的动脉血氧饱和度的计算方法即可完成动脉血氧饱和度的测量,不需要复杂的光电转换装置,能够简化计算过程,提高计算精度,同时为血氧检测设备的小型化和低成本化提供了一种可能。
2、本发明采用的基于单色光的动脉血管图像的动脉血氧饱和度的计算方法不仅可以测得动脉血氧饱和度,还可测出动脉血管中的氧合血红蛋白的浓度的绝对值,为一种无创检测方法,为无创氧合血红蛋白的浓度的检测提供了一种新的检测方案。
附图说明
图1是现有的透射式双波长动脉血氧饱和度测量方法的原理图;
图2是基于单色光动脉血管图像的血液含氧量状况测量方法的流程图;
图3是采集图像数字信号的示意图;
图4是基于单色光动脉血管图像的血液含氧量状况测量装置的框架图。
具体实施方式
为进一步了解本发明的内容,结合实施例对本发明作详细描述,以下实施例用于说明本发明,但不用来限制本发明的范围。
实施例一
参照附图2所示,本实施例涉及的基于单色光动脉血管图像的血液含氧量状况测量方法,其包括以下步骤:
1)参照附图3所示,将手指放入采集模块中,该采集模块的上方设有单色光二极管作为光源,下方设有摄像装置,手指放在单色光二极管和摄像装置之间,获取多组单色光透射后的光信号,并将多组光信号转化为多组图像数字信号;单色光二极管发出的单色光的波长为800nm~950nm,摄像装置的帧率为N帧/s,手指放置时间至少Ts,获得的图像数字信号为NT张;
本实施例中,N为1帧/s,T为1s,得到60张图像数字信号,获得的手指的图像数字信号集合{fvImg}为:
{fvImg}={fvImg1,fvImg2,fvImg3,...,fvImgNT} (30)。
需要说明的是,在该实施例中,单色光二极管放在待采集部位的上方,摄像装置设置于待采集部位的下方;在其他实施方案中,也可将单色光二极管和摄像装置分别置于目标部位的左右两侧,或将单色光二极管和摄像装置放置于目标部位的同一侧。
需要说明的是,在该实施例中,采集部位是手指,在其他实施方案中,采集部位也可以是手掌或面部等其它包含动脉血管的其它部位。
2)对获得的某一图像数字信号进行图像增强处理和二值化处理,其具体步骤包括:
2.1)计算某一图像序列中的每一像素的曲率值,曲率计算公式如下:
其中,
公式中,Cij为某一图像序列的第i行第j列像素的曲率值;Imgij为某一图像序列的第i行第j列像素的灰度值;pij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导;qij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导;rij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导;tij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导;sij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导;
2.2)计算某一图像序列中的全部像素的曲率均值,曲率均值的计算公式如下:
其中,Cavg为某一图像序列中的全部像素的曲率均值;r为某一图像序列的图像行数;c为某一图像序列的图像列数;
2.3)对动脉血管的二值化处理,二值化计算方式为:
其中,fij为某一图像序列的第i行第j列像素的二值化的数值;α为曲率均值二值化系数,取值范围为[0,1],本实施例中,α的取值为0.5;
3)根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程,其计算公式为:
Smax=max{S1,S2,S3,....,SA} (10)
Smin=min{S1,S2,S3,...,SA} (11)
Lmax=Smax (13)
Lmin=Smin (16)
公式中,Sm为第m个图像序列的动脉血管二值化有效面积;为第m个图像序列的第i行第j列像素的二值化数值;Smax为动脉血管扩张时的二值化有效面积;Smin为动脉血管收缩时的二值化有效面积,A为获取的图像数字信号的数量,本实施例A为60张;Lmax为动脉血管扩张时的光程;Lmin为动脉血管收缩时的光程。此步骤有利于降低动脉血氧饱和度的计算误差。
4)根据动脉血管的扩张光程和收缩光程计算出动脉血氧饱和度及氧合血红蛋白的浓度,用于反应人体的血液氧含量状况,其具体步骤包括:
4.1)计算动脉血管扩张时的透射光的强度,计算公式如下:
公式中,IA为动脉血管扩张时的透射光的强度;I0为动脉血管扩张时的入射光的强度;γ为组织内的非脉动成分及静脉血的吸光系数;εHbO2为表示动脉血管中氧合血红蛋白的吸光系数;CHbO2为表示动脉血管中氧合血红蛋白的浓度;εHb为表示动脉血管中还原血红蛋白的吸光系数;CHb为表示动脉血管中还原血红蛋白的浓度;为动脉扩张图像序列的第i行第j列像素的灰度值;/>为动脉扩张图像序列的第i行第j列像素的二值化的数值;r1为动脉扩张图像序列的图像行数;c1为动脉扩张图像序列的图像列数;
4.2)计算动脉血管收缩时的透射光的强度,计算公式如下:
公式中,IB为动脉血管收缩时的透射光的强度;I0’为动脉血管收缩时的入射光的强度;为动脉收缩图像序列的第i行第j列像素的灰度值;/>为动脉收缩图像序列的第i行第j列像素的二值化的数值;r2为动脉收缩图像序列的图像行数;c2为动脉收缩图像序列的图像列数;
4.3)对公式(12)~(14)进行整合和取对数,得到公式(18):
4.4)对公式(15)~(17)进行整合和取对数,得到公式(19):
4.5)对公式(18)~(19)进行整合,得到公式(20)和公式(21),用于计算动脉血管的搏动信息的量化参数:
CHbO2εHbo2+CHbεHb=R (20)
公式中,R为动脉血管的搏动信息的量化参数,通过公式(21)计算具体值;
4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度:
其中,C0为表示动脉血管中血红蛋白的浓度,其计算公式为:
CHbO2+CHb=C0 (23)
公式中,CHbO2为表示动脉血管中氧合血红蛋白的浓度;CHb为表示动脉血管中还原血红蛋白的浓度;
4.7)将公式(23)代入公式(29)中,即获得动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度的计算公式(24),采用该公式计算动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度;
其中,C0为动脉血管中血红蛋白的浓度,其与动脉血管的搏动信息的量化参数R存在函数关系,函数关系为:
R=f(C0) (25)
公式中,f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0的函数关系,本实施例中,所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0存在线性关系,即:
R=kC0+B (26)
其中,k为线性系数;B为漂移常量;
将公式(26)代入公式(22)和公式(24),可得到动脉血管中氧合血红蛋白的浓度和动脉血氧饱和度的计算公式:
用动脉血管中氧合血红蛋白的浓度和动脉血氧饱和度来反应人体的血液氧含量状况。
效果实施例
本实施例随机抽取200名测试者,采用三种测试方法对200名测试者的动脉血氧饱和度或/和氧合血红蛋白浓度进行测量。
方法一:通过传统的血氧饱和度测量方法进行人体动脉采血,利用血气分析仪进行电化学分析,测出血氧分压PO2计算出血氧饱和度。
方法二:借助现有发明专利CN108113682B涉及的测量含氧血红蛋白和脱氧血红蛋白的装置进行氧合血红蛋白浓度,即利用血红蛋白在不同氧合状态下对近红外光具有的不同吸收谱这一特性,对人体组织中氧合血红蛋白进行定量检测。
方法三:采用本发明涉及的装置采集图像并计算出血氧饱和度和氧合血红蛋白浓度。
将方法三中获取的200组血氧饱和度与方法一中获取的对应测试者的血氧饱和度进行比较,将方法三中获取的氧合血红蛋白浓度与方法二中获取的对应测试者的氧合血红蛋白浓度进行比较,当血氧饱和度误差小于1%且氧合血红蛋白浓度误差小于0.5mol/L,则说明方法三中使用的基于单色光动脉血管图像的血液含氧量状况测量装置的测试结果准确,否则说明该装置的测试结果不准确。通过对200位测试的测试,计算出采用本发明涉及的基于单色光动脉血管图像的血液含氧量状况测量装置的测试准确率高达99.5%,因此,证明本发明涉及的基于单色光动脉血管图像的血液含氧量状况测量方法及装置是有效的。
实施例二
参照附图4所述,基于单色光动脉血管图像的血液含氧量状况测量装置包括:
采集模块,,该采集模块包括光源和摄像装置,所述光源为单波长光源,所述摄像装置采集光透射待采集部位后的光信号,并将所述光信号转化为图像数字信号;所述的采集模块用于实现实施例一步骤1)的功能;
图像处理模块,用于对获得的所有图像数字信号进行图像增强处理和二值化处理;图像处理模块用于实现实施例一步骤2)的功能。
光程计算模块,用于根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程,光程计算模块用于实现实施例一步骤3)的功能;
4)动脉血氧饱和度及氧合血红蛋白浓度计算模块,用于根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度和氧合血红蛋白浓度,用于反应人体的血液氧含量状况,动脉血氧饱和度及氧合血红蛋白浓度计算模块用于实现实施例一步骤4)的功能。
以上结合实施例对本发明进行了详细说明,但所述内容仅为本发明的较佳实施例,不能被认为用于限定本发明的实施范围。凡依本发明申请范围所作的均等变化与改进等,均应仍属于本发明的专利涵盖范围之内。
Claims (5)
1.一种基于单色光动脉血管图像的血液含氧量状况测量方法,其特征在于:其包括以下步骤:
1)将待采集部位放入采集模块中,该采集模块包括光源和摄像装置,所述光源为单波长光源,所述摄像装置采集光透射待采集部位后的光信号,并将所述光信号转化为图像数字信号;
2)对获得的所有图像数字信号进行图像增强处理和二值化处理,具体步骤包括:2.1)计算某一图像序列中的每一像素的曲率值,曲率计算公式如下:
其中,
公式中,Cij为某一图像序列的第i行第j列像素的曲率值;Imgij为某一图像序列的第i行第j列像素的灰度值;pij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导;qij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导;rij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导;tij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导;sij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导;
2.2)计算某一图像序列中的全部像素的曲率均值,曲率均值的计算公式如下:
其中,Cavg为某一图像序列中的全部像素的曲率均值;r为某一图像序列的图像行数;c为某一图像序列的图像列数;
2.3)对动脉血管的二值化处理,二值化计算方式为:
其中,fij为某一图像序列的第i行第j列像素的二值化的数值;α为曲率均值二值化系数,取值范围为[0,1];
3)根据增强处理和二值化处理的图像数字信号的二值化有效面积计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程,具体计算公式为:
Smax=max{S1,S2,S3,...,SA} (10)
Smin=min{S1,S2,S3,...,SA} (11)
Lmax=Smax (13)
Lmin=Smin (16)
公式中,Sm为第m个图像序列的动脉血管二值化有效面积;为第m个图像序列的第i行第j列像素的二值化数值;Smax为动脉血管扩张时的二值化有效面积;Smin为动脉血管收缩时的二值化有效面积;A为获取的图像数字信号的数量;Lmax为动脉血管扩张时的光程;Lmin为动脉血管收缩时的光程;
4)根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度及氧合血红蛋白的浓度,用于反应人体的血液氧含量状况,具体步骤包括:
4.1)计算动脉血管扩张时的透射光的强度,计算公式如下:
公式中,IA为动脉血管扩张时的透射光的强度;I0为动脉血管扩张时的入射光的强度;γ为组织内的非脉动成分及静脉血的吸光系数;εHbO2为表示动脉血管中氧合血红蛋白的吸光系数;CHbO2为表示动脉血管中氧合血红蛋白的浓度;εHb为表示动脉血管中还原血红蛋白的吸光系数;CHb为表示动脉血管中还原血红蛋白的浓度;为动脉扩张图像序列的第i行第j列像素的灰度值;/>为动脉扩张图像序列的第i行第j列像素的二值化的数值;r1为动脉扩张图像序列的图像行数;c1为动脉扩张图像序列的图像列数;
4.2)计算动脉血管收缩时的透射光的强度,计算公式如下:
公式中,IB为动脉血管收缩时的透射光的强度;I0’为动脉血管收缩时的入射光的强度;为动脉收缩图像序列的第i行第j列像素的灰度值;/>为动脉收缩图像序列的第i行第j列像素的二值化的数值;r2为动脉收缩图像序列的图像行数;c2为动脉收缩图像序列的图像列数;
4.3)对公式(12)~(14)进行整合和取对数,得到公式(18):
4.4)对公式(15)~(17)进行整合和取对数,得到公式(19):
4.5)对公式(18)~(19)进行整合,得到公式(20)和公式(21),用于计算动脉血管的搏动信息的量化参数:
CHbO2εHbO2+CHbεHb=R (20)
公式中,R为动脉血管的搏动信息的量化参数,通过公式(21)计算具体值;
4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度:
其中,C0为表示动脉血管中血红蛋白的浓度,其计算公式为:
CHbO2+CHb=C0 (23)
公式中,CHbO2为表示动脉血管中氧合血红蛋白的浓度;CHb为表示动脉血管中还原血红蛋白的浓度;
4.7)基于动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度,计算公式为:
其中,C0为动脉血管中血红蛋白的浓度与动脉血管的搏动信息的量化参数R存在函数关系,函数关系为:
R=f(C0) (25)
公式中,f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0的函数关系,所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0存在线性关系,即:
R=kC0+B (26)
其中,k为线性系数;B为漂移常量;
将公式(26)代入公式(22)和公式(24),可得到动脉血管中氧合血红蛋白的浓度和动脉血氧饱和度的计算公式:
2.根据权利要求1所述的基于单色光动脉血管图像的血液含氧量状况测量方法,其特征在于:所述的曲率均值二值化系数α的取值为0.5。
3.根据权利要求1所述的基于单色光动脉血管图像的血液含氧量状况测量方法,其特征在于:所述的步骤1)中单色光二极管发出的单色光的波长为800nm~950nm。
4.根据权利要求1所述的基于单色光动脉血管图像的血液含氧量状况测量方法,其特征在于:所述的步骤1)中摄像装置的帧率为60帧/s,手指放置时间为至少1s,获得的图像数字信号至少为60张。
5.一种基于单色光动脉血管图像的血液含氧量状况测量装置,其特征在于:其包括:
采集模块,该采集模块包括光源和摄像装置,所述光源为单波长光源,所述摄像装置采集光透射待采集部位后的光信号,并将所述光信号转化为图像数字信号;图像处理模块,用于对获得的所有图像数字信号进行图像增强处理和二值化处理,具体步骤包括:
2.1)计算某一图像序列中的每一像素的曲率值,曲率计算公式如下:
其中,
公式中,Cij为某一图像序列的第i行第j列像素的曲率值;Imgij为某一图像序列的第i行第j列像素的灰度值;pij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导;qij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导;rij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导;tij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导;sij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导;
2.2)计算某一图像序列中的全部像素的曲率均值,曲率均值的计算公式如下:
其中,Cavg为某一图像序列中的全部像素的曲率均值;r为某一图像序列的图像行数;c为某一图像序列的图像列数;
2.3)对动脉血管的二值化处理,二值化计算方式为:
其中,fij为某一图像序列的第i行第j列像素的二值化的数值;α为曲率均值二值化系数,取值范围为[0,1];
光程计算模块,用于根据增强处理和二值化处理的图像数字信号的二值化有效面积计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程,具体计算公式为:
Smax=max{S1,S2,S3,...,SA} (10)
Smin=min{S1,S2,S3,...,SA} (11)
Lmax=Smax (13)
Lmin=Smin (16)
公式中,Sm为第m个图像序列的动脉血管二值化有效面积;为第m个图像序列的第i行第j列像素的二值化数值;Smax为动脉血管扩张时的二值化有效面积;Smin为动脉血管收缩时的二值化有效面积;A为获取的图像数字信号的数量;Lmax为动脉血管扩张时的光程;Lmin为动脉血管收缩时的光程;
动脉血氧饱和度及氧合血红蛋白浓度计算模块,用于根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度和氧合血红蛋白浓度,用于反应人体的血液氧含量状况,具体步骤包括:
4.1)计算动脉血管扩张时的透射光的强度,计算公式如下:
公式中,IA为动脉血管扩张时的透射光的强度;I0为动脉血管扩张时的入射光的强度;γ为组织内的非脉动成分及静脉血的吸光系数;εHbO2为表示动脉血管中氧合血红蛋白的吸光系数;CHbO2为表示动脉血管中氧合血红蛋白的浓度;εHb为表示动脉血管中还原血红蛋白的吸光系数;CHb为表示动脉血管中还原血红蛋白的浓度;为动脉扩张图像序列的第i行第j列像素的灰度值;/>为动脉扩张图像序列的第i行第j列像素的二值化的数值;r1为动脉扩张图像序列的图像行数;c1为动脉扩张图像序列的图像列数;
4.2)计算动脉血管收缩时的透射光的强度,计算公式如下:
公式中,IB为动脉血管收缩时的透射光的强度;I0’为动脉血管收缩时的入射光的强度;为动脉收缩图像序列的第i行第j列像素的灰度值;/>为动脉收缩图像序列的第i行第j列像素的二值化的数值;r2为动脉收缩图像序列的图像行数;c2为动脉收缩图像序列的图像列数;
4.3)对公式(12)~(14)进行整合和取对数,得到公式(18):
4.4)对公式(15)~(17)进行整合和取对数,得到公式(19):
4.5)对公式(18)~(19)进行整合,得到公式(20)和公式(21),用于计算动脉血管的搏动信息的量化参数:
CHbO2εHbO2+CHbεHb=R (20)
公式中,R为动脉血管的搏动信息的量化参数,通过公式(21)计算具体值;
4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度:
其中,C0为表示动脉血管中血红蛋白的浓度,其计算公式为:
CHbO2+CHb=C0 (23)
公式中,CHbO2为表示动脉血管中氧合血红蛋白的浓度;CHb为表示动脉血管中还原血红蛋白的浓度;
4.7)基于动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度,计算公式为:
其中,C0为动脉血管中血红蛋白的浓度与动脉血管的搏动信息的量化参数R存在函数关系,函数关系为:
R=f(C0) (25)
公式中,f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0的函数关系,所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C0存在线性关系,即:
R=kC0+B (26)
其中,k为线性系数;B为漂移常量;
将公式(26)代入公式(22)和公式(24),可得到动脉血管中氧合血红蛋白的浓度和动脉血氧饱和度的计算公式:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111491945.4A CN114376571B (zh) | 2021-12-08 | 2021-12-08 | 基于单色光动脉血管图像的血液含氧量状况测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111491945.4A CN114376571B (zh) | 2021-12-08 | 2021-12-08 | 基于单色光动脉血管图像的血液含氧量状况测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114376571A CN114376571A (zh) | 2022-04-22 |
CN114376571B true CN114376571B (zh) | 2023-09-15 |
Family
ID=81195381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111491945.4A Active CN114376571B (zh) | 2021-12-08 | 2021-12-08 | 基于单色光动脉血管图像的血液含氧量状况测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114376571B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004194908A (ja) * | 2002-12-18 | 2004-07-15 | Hamamatsu Photonics Kk | 血液測定装置 |
CN103479363A (zh) * | 2013-09-30 | 2014-01-01 | 深圳市倍轻松科技股份有限公司 | 血液中血氧饱和度的测量方法及系统 |
CN111991004A (zh) * | 2020-08-20 | 2020-11-27 | 广州医软智能科技有限公司 | 血氧饱和度测量设备、测量方法和测量装置 |
CN112858196A (zh) * | 2020-10-20 | 2021-05-28 | 重庆大学 | 一种三波长静脉血氧浓度测量的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180333088A1 (en) * | 2017-05-17 | 2018-11-22 | Microsoft Technology Licensing, Llc | Pulse Oximetry Capturing Technique |
-
2021
- 2021-12-08 CN CN202111491945.4A patent/CN114376571B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004194908A (ja) * | 2002-12-18 | 2004-07-15 | Hamamatsu Photonics Kk | 血液測定装置 |
CN103479363A (zh) * | 2013-09-30 | 2014-01-01 | 深圳市倍轻松科技股份有限公司 | 血液中血氧饱和度的测量方法及系统 |
CN111991004A (zh) * | 2020-08-20 | 2020-11-27 | 广州医软智能科技有限公司 | 血氧饱和度测量设备、测量方法和测量装置 |
CN112858196A (zh) * | 2020-10-20 | 2021-05-28 | 重庆大学 | 一种三波长静脉血氧浓度测量的方法 |
Non-Patent Citations (5)
Title |
---|
"Continuous blood oxygen saturation detection with single‑wavelength photoacoustics";Qiwen Peng等;Progress in Biomedical Optics and Imaging ‑ Proceedings of SPIE;93231R-1-93231R-7 * |
"Single-Wavelength Blood Oxygen Saturation Sensing With Combined Optical Absorption and Scattering";Fei Gao等;IEEE SENSORS JOURNAL;第第16卷卷(第第7期期);第1943-1948页 * |
"single-wavelength oxygen saturation detection fusing optical absorption and scattering properties:a phantom study";MA YAO等;光谱学与光谱分析;第第37卷卷(第第8期期);第2652-2656页 * |
"无创血氧饱和度的测量 及无线监测系统的研制";谢正洋等;中国仪器仪表(第第4期期);第65-67、90页 * |
"非侵入式脉搏血氧检测系统设计";张振等;中国集成电路(第第168期期);第15-20页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114376571A (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101347334B (zh) | 血氧饱和度测量方法和装置 | |
CN105050492B (zh) | 用于确定对象的生命体征的设备和方法 | |
CN105009173B (zh) | 用于确定对象的动脉血氧饱和度的设备和方法 | |
CN100335001C (zh) | 使用光测量技术的诊断装置 | |
JP6194105B2 (ja) | 遠隔フォトプレチスモグラフィ波形を取得するための改良された信号選択 | |
CN100450437C (zh) | 低灌注下测量血氧的方法 | |
US5127406A (en) | Apparatus for measuring concentration of substances in blood | |
US9591999B2 (en) | Determination of tissue oxygenation in vivo | |
CN113080872A (zh) | 用于确定对象的生命体征的设备和方法 | |
JP2011502716A (ja) | 脈動性生物測定信号の処理方法および装置 | |
CN104968259A (zh) | 用于确定对象的生命体征信息的系统和方法 | |
CN107320112B (zh) | 微循环的多参数成像检测方法及装置 | |
US20070299323A1 (en) | Apparatus for measuring one or more physiological functions of a body and a method using the same | |
US20080144004A1 (en) | Optical Spectrophotometer | |
CN112294282A (zh) | 基于rppg的情绪检测装置的自标定方法 | |
WO2022237171A1 (zh) | 一种基于ccd的双波长红外血氧检测系统 | |
CN110604584A (zh) | 用于血氧测量系统的信号检测方法 | |
CN114376571B (zh) | 基于单色光动脉血管图像的血液含氧量状况测量方法及装置 | |
EP3520696A1 (en) | Blood glucose measurement device, blood glucose calculation method, and blood glucose calculation program | |
Vogel et al. | In-ear heart rate monitoring using a micro-optic reflective sensor | |
CN113974617B (zh) | 基于组织血氧宽场成像的血氧检测方法及系统 | |
RU2294141C1 (ru) | Устройство для измерения уровня оксигенации и частоты пульса | |
CN114387212A (zh) | 基于双色光动脉血管图像的动脉血氧饱和度测量方法及装置 | |
CN111631733B (zh) | 一种动脉血液光谱检测方法及装置 | |
CN111759292B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |