CN114376571B - 基于单色光动脉血管图像的血液含氧量状况测量方法及装置 - Google Patents

Abstract

本发明公开了基于单色光动脉血管图像的血液含氧量状况测量方法及装置，其包括以下步骤：1)将待采集部位放入采集模块中，采集模块的上方设有单色光二极管，下方设有摄像装置，通过采集模块获取多组单色光透射后的光信号，并将多组光信号转化为多组图像数字信号；2)对获得的所有图像数字信号进行图像增强处理和二值化处理；3)根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程；4)根据动脉血管的扩张光程和收缩光程计算血氧饱和度和氧合血红蛋白的浓度，反应人体的血液氧含量状况。本发明为血氧检测设备的小型化提供了可能，为无创氧合血红蛋白浓度的检测提供了新的检测方案。

Description

基于单色光动脉血管图像的血液含氧量状况测量方法及装置

技术领域

本发明涉及图像识别及处理技术领域，尤其涉及一种基于单色光动脉血管图像的血液含氧量状况测量方法及装置。

背景技术

无创血氧饱和度的测量原理是动脉搏动引起血液对光的吸收的变化，并将光的变化转换为电信号，并对电信号进行处理，进而可计算出血氧饱和度。血氧饱和度一般用血液中氧合血红蛋白占总血红蛋白的百分比来表示，一般通过监测动脉的功能血氧饱和度(SpO₂)来反应人体的血液氧含量状况，其表达式为：

其中，SpO₂为血氧饱和度；为动脉血管中氧合血红蛋白的浓度；C_Hb为动脉血管中还原血红蛋白的浓度。

目前医用的透射式无创血氧饱和度的测量方案是采用双波长异步打光方案，如图1所示，手指上壁固定了两个并排放置的发光二极管，波长分别为660nm的红光和904nm的红外光，下壁是一个光电转换器件，将透射光的光强转换为电流信号。通过控制两个LED轮换发光，并对各自发光时输出的电流信号进行I/V变换，放大然后采集，对采集后的信号做处理，得到对应的脉搏波数据，最后计算得到SpO₂的值。该方案获得血氧饱和度的方法如中国专利CN112903609A公开的无需校正的双波长静脉血氧饱和度测量方法，包括以下步骤：通过DAC控制红光与红外光二极管发光，使用光频传感器接收透射光，将透光信号转换为对应的频率信号f，然后将得到的频率数据f和f′通过MCU上传到计算机中，再经过推导得到的公式计算即可得出血氧饱和度。

然而，现有的血氧饱和度测量技术需要采用双波长异步打光方案，同时需要复杂的光电转换装置和信号处理单元，计算过程比较繁琐，计算精度不高。采集模块需要复杂的光电转换装置，设备难以做到小型化，设备成本较高。

发明内容

本发明的目的是提供一种基于单色光动脉血管图像的血液含氧量状况测量方法及装置，以解决光信号处理时计算过程繁琐、精度不高、采集模块成本较高等问题。

为了达到目的，本发明提供的技术方案为：

本发明涉及的一种基于单色光动脉血管图像的血液含氧量状况测量方法，其包括以下步骤：

1)将待采集部位放入采集模块中，该采集模块包括光源和摄像装置，所述光源为单波长光源，所述摄像装置采集光透射待采集部位后的光信号，并将所述光信号转化为图像数字信号；

2)对获得的所有图像数字信号进行图像增强处理和二值化处理；

3)根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程；

4)根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度及氧合血红蛋白的浓度，用于反应人体的血液氧含量状况。

优选地，所述的步骤2)中对获得的某一图像数字信号进行图像增强处理和二值化处理的具体步骤包括：

2.1)计算某一图像序列中的每一像素的曲率值，曲率计算公式如下：

其中，

公式中，C_ij为某一图像序列的第i行第j列像素的曲率值；Img_ij为某一图像序列的第i行第j列像素的灰度值；p_ij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导；q_ij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导；r_ij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导；t_ij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导；s_ij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导；

2.2)计算某一图像序列中的全部像素的曲率均值，曲率均值的计算公式如下：

其中，C_avg为某一图像序列中的全部像素的曲率均值；r为某一图像序列的图像行数；c为某一图像序列的图像列数；

2.3)对动脉血管的二值化处理，二值化计算方式为：

其中，f_ij为某一图像序列的第i行第j列像素的二值化的数值；α为曲率均值二值化系数，取值范围为[0，1]。

优选地，所述的曲率均值二值化系数α的取值为0.5。

优选地，所述的步骤3)中是根据增强处理和二值化处理后图像数字信号的二值化有效面积来计算扩张光程和收缩光程的，具体计算公式为：

S_max＝max{S₁，S₂，S₃，...，S_A} (10)

S_min＝min{S₁，S₂，S₃，...，S_A} (11)

L_max＝S_max (13)

L_min＝S_min (16)

公式中，S_m为第m个图像序列的动脉血管二值化有效面积；为第m个图像序列的第i行第j列像素的二值化数值；S_max为动脉血管扩张时的二值化有效面积；S_min为动脉血管收缩时的二值化有效面积；A为获取的图像数字信号的数量；L_max为动脉血管扩张时的光程；L_min为动脉血管收缩时的光程。此步骤有利于降低动脉血氧饱和度的计算误差。

优选地，所述的步骤4)根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度的具体步骤包括：

4.1)计算动脉血管扩张时的透射光的强度，计算公式如下：

公式中，I_A为动脉血管扩张时的透射光的强度；I₀为动脉血管扩张时的入射光的强度；γ为组织内的非脉动成分及静脉血的吸光系数；ε_HbO2为表示动脉血管中氧合血红蛋白的吸光系数；C_HbO2为表示动脉血管中氧合血红蛋白的浓度；ε_Hb为表示动脉血管中还原血红蛋白的吸光系数；C_Hb为表示动脉血管中还原血红蛋白的浓度；为动脉扩张图像序列的第i行第j列像素的灰度值；/>为动脉扩张图像序列的第i行第j列像素的二值化的数值；r₁为动脉扩张图像序列的图像行数；c₁为动脉扩张图像序列的图像列数；

4.2)计算动脉血管收缩时的透射光的强度，计算公式如下：

公式中，I_B为动脉血管收缩时的透射光的强度；I₀’为动脉血管收缩时的入射光的强度；为动脉收缩图像序列的第i行第j列像素的灰度值；/>为动脉收缩图像序列的第i行第j列像素的二值化的数值；r₂为动脉收缩图像序列的图像行数；c₂为动脉收缩图像序列的图像列数；

4.3)对公式(12)～(14)进行整合和取对数，得到公式(18)：

4.4)对公式(15)～(17)进行整合和取对数，得到公式(19)：

4.5)对公式(18)～(19)进行整合，得到公式(20)和公式(21)，用于计算动脉血管的搏动信息的量化参数：

C_Hbo2ε_HbO2+C_Hbε_Hb＝R (20)

公式中，R为动脉血管的搏动信息的量化参数，通过公式(21)计算具体值；

4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度，计算公式为：

其中，C₀为表示动脉血管中血红蛋白的浓度：

C_HbO2+C_Hb＝C₀ (23)

公式中，C_HbO2为表示动脉血管中氧合血红蛋白的浓度；C_Hb为表示动脉血管中还原血红蛋白的浓度；

4.7)基于动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度，计算公式为：

其中，C₀为动脉血管中血红蛋白的浓度与动脉血管的搏动信息的量化参数R存在函数关系，函数关系为：

R＝f(C₀) (25)

公式中，f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀的函数关系。

优选地，所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀存在线性关系，即：

R＝kC₀+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}＝{fvImg₁,fvImg₂，fvImg₃，...，fvImg_NT} (30)。

需要说明的是，在该实施例中，单色光二极管放在待采集部位的上方，摄像装置设置于待采集部位的下方；在其他实施方案中，也可将单色光二极管和摄像装置分别置于目标部位的左右两侧，或将单色光二极管和摄像装置放置于目标部位的同一侧。

需要说明的是，在该实施例中，采集部位是手指，在其他实施方案中，采集部位也可以是手掌或面部等其它包含动脉血管的其它部位。

2)对获得的某一图像数字信号进行图像增强处理和二值化处理，其具体步骤包括：

2.1)计算某一图像序列中的每一像素的曲率值，曲率计算公式如下：

其中，

公式中，C_ij为某一图像序列的第i行第j列像素的曲率值；Img_ij为某一图像序列的第i行第j列像素的灰度值；p_ij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导；q_ij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导；r_ij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导；t_ij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导；s_ij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导；

2.2)计算某一图像序列中的全部像素的曲率均值，曲率均值的计算公式如下：

其中，C_avg为某一图像序列中的全部像素的曲率均值；r为某一图像序列的图像行数；c为某一图像序列的图像列数；

2.3)对动脉血管的二值化处理，二值化计算方式为：

其中，f_ij为某一图像序列的第i行第j列像素的二值化的数值；α为曲率均值二值化系数，取值范围为[0，1]，本实施例中，α的取值为0.5；

3)根据增强处理和二值化处理的图像数字信号计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程，其计算公式为：

S_max＝max{S₁，S₂，S₃，....，S_A} (10)

S_min＝min{S₁，S₂，S₃，...，S_A} (11)

L_max＝S_max (13)

L_min＝S_min (16)

公式中，S_m为第m个图像序列的动脉血管二值化有效面积；为第m个图像序列的第i行第j列像素的二值化数值；S_max为动脉血管扩张时的二值化有效面积；S_min为动脉血管收缩时的二值化有效面积，A为获取的图像数字信号的数量，本实施例A为60张；L_max为动脉血管扩张时的光程；L_min为动脉血管收缩时的光程。此步骤有利于降低动脉血氧饱和度的计算误差。

4)根据动脉血管的扩张光程和收缩光程计算出动脉血氧饱和度及氧合血红蛋白的浓度，用于反应人体的血液氧含量状况，其具体步骤包括：

4.1)计算动脉血管扩张时的透射光的强度，计算公式如下：

公式中，I_A为动脉血管扩张时的透射光的强度；I₀为动脉血管扩张时的入射光的强度；γ为组织内的非脉动成分及静脉血的吸光系数；ε_HbO2为表示动脉血管中氧合血红蛋白的吸光系数；C_HbO2为表示动脉血管中氧合血红蛋白的浓度；ε_Hb为表示动脉血管中还原血红蛋白的吸光系数；C_Hb为表示动脉血管中还原血红蛋白的浓度；为动脉扩张图像序列的第i行第j列像素的灰度值；/>为动脉扩张图像序列的第i行第j列像素的二值化的数值；r₁为动脉扩张图像序列的图像行数；c₁为动脉扩张图像序列的图像列数；

4.2)计算动脉血管收缩时的透射光的强度，计算公式如下：

公式中，I_B为动脉血管收缩时的透射光的强度；I₀’为动脉血管收缩时的入射光的强度；为动脉收缩图像序列的第i行第j列像素的灰度值；/>为动脉收缩图像序列的第i行第j列像素的二值化的数值；r₂为动脉收缩图像序列的图像行数；c₂为动脉收缩图像序列的图像列数；

4.3)对公式(12)～(14)进行整合和取对数，得到公式(18)：

4.4)对公式(15)～(17)进行整合和取对数，得到公式(19)：

4.5)对公式(18)～(19)进行整合，得到公式(20)和公式(21)，用于计算动脉血管的搏动信息的量化参数：

C_HbO2ε_Hbo2+C_Hbε_Hb＝R (20)

公式中，R为动脉血管的搏动信息的量化参数，通过公式(21)计算具体值；

4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度：

其中，C₀为表示动脉血管中血红蛋白的浓度，其计算公式为：

C_HbO2+C_Hb＝C₀ (23)

公式中，C_HbO2为表示动脉血管中氧合血红蛋白的浓度；C_Hb为表示动脉血管中还原血红蛋白的浓度；

4.7)将公式(23)代入公式(29)中，即获得动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度的计算公式(24)，采用该公式计算动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度；

其中，C₀为动脉血管中血红蛋白的浓度，其与动脉血管的搏动信息的量化参数R存在函数关系，函数关系为：

R＝f(C₀) (25)

公式中，f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀的函数关系，本实施例中，所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀存在线性关系，即：

R＝kC₀+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)计算某一图像序列中的每一像素的曲率值，曲率计算公式如下：

其中，

公式中，C_ij为某一图像序列的第i行第j列像素的曲率值；Img_ij为某一图像序列的第i行第j列像素的灰度值；p_ij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导；q_ij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导；r_ij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导；t_ij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导；s_ij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导；

2.2)计算某一图像序列中的全部像素的曲率均值，曲率均值的计算公式如下：

其中，C_avg为某一图像序列中的全部像素的曲率均值；r为某一图像序列的图像行数；c为某一图像序列的图像列数；

2.3)对动脉血管的二值化处理，二值化计算方式为：

其中，f_ij为某一图像序列的第i行第j列像素的二值化的数值；α为曲率均值二值化系数，取值范围为[0，1]；

3)根据增强处理和二值化处理的图像数字信号的二值化有效面积计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程，具体计算公式为：

S_max＝max{S₁，S₂，S₃，...，S_A} (10)

S_min＝min{S₁，S₂，S₃，...，S_A} (11)

L_max＝S_max (13)

L_min＝S_min (16)

公式中，S_m为第m个图像序列的动脉血管二值化有效面积；为第m个图像序列的第i行第j列像素的二值化数值；S_max为动脉血管扩张时的二值化有效面积；S_min为动脉血管收缩时的二值化有效面积；A为获取的图像数字信号的数量；L_max为动脉血管扩张时的光程；L_min为动脉血管收缩时的光程；

4)根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度及氧合血红蛋白的浓度，用于反应人体的血液氧含量状况，具体步骤包括：

4.1)计算动脉血管扩张时的透射光的强度，计算公式如下：

公式中，I_A为动脉血管扩张时的透射光的强度；I₀为动脉血管扩张时的入射光的强度；γ为组织内的非脉动成分及静脉血的吸光系数；ε_HbO2为表示动脉血管中氧合血红蛋白的吸光系数；C_HbO2为表示动脉血管中氧合血红蛋白的浓度；ε_Hb为表示动脉血管中还原血红蛋白的吸光系数；C_Hb为表示动脉血管中还原血红蛋白的浓度；为动脉扩张图像序列的第i行第j列像素的灰度值；/>为动脉扩张图像序列的第i行第j列像素的二值化的数值；r₁为动脉扩张图像序列的图像行数；c₁为动脉扩张图像序列的图像列数；

4.2)计算动脉血管收缩时的透射光的强度，计算公式如下：

公式中，I_B为动脉血管收缩时的透射光的强度；I₀’为动脉血管收缩时的入射光的强度；为动脉收缩图像序列的第i行第j列像素的灰度值；/>为动脉收缩图像序列的第i行第j列像素的二值化的数值；r₂为动脉收缩图像序列的图像行数；c₂为动脉收缩图像序列的图像列数；

4.3)对公式(12)～(14)进行整合和取对数，得到公式(18)：

4.4)对公式(15)～(17)进行整合和取对数，得到公式(19)：

4.5)对公式(18)～(19)进行整合，得到公式(20)和公式(21)，用于计算动脉血管的搏动信息的量化参数：

C_HbO2ε_HbO2+C_Hbε_Hb＝R (20)

公式中，R为动脉血管的搏动信息的量化参数，通过公式(21)计算具体值；

4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度：

其中，C₀为表示动脉血管中血红蛋白的浓度，其计算公式为：

C_HbO2+C_Hb＝C₀ (23)

公式中，C_HbO2为表示动脉血管中氧合血红蛋白的浓度；C_Hb为表示动脉血管中还原血红蛋白的浓度；

4.7)基于动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度，计算公式为：

其中，C₀为动脉血管中血红蛋白的浓度与动脉血管的搏动信息的量化参数R存在函数关系，函数关系为：

R＝f(C₀) (25)

公式中，f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀的函数关系，所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀存在线性关系，即：

R＝kC₀+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)计算某一图像序列中的每一像素的曲率值，曲率计算公式如下：

其中，

公式中，C_ij为某一图像序列的第i行第j列像素的曲率值；Img_ij为某一图像序列的第i行第j列像素的灰度值；p_ij为某一图像序列的第i行第j列像素在x轴方向上的一阶偏导；q_ij为某一图像序列的第i行第j列像素在y轴方向上的一阶偏导；r_ij为某一图像序列的第i行第j列像素在x轴方向上的二阶偏导；t_ij为某一图像序列的第i行第j列像素在y轴方向上的二阶偏导；s_ij为某一图像序列的第i行第j列像素在x轴方向和y轴方向上的二阶混合偏导；

2.2)计算某一图像序列中的全部像素的曲率均值，曲率均值的计算公式如下：

其中，C_avg为某一图像序列中的全部像素的曲率均值；r为某一图像序列的图像行数；c为某一图像序列的图像列数；

2.3)对动脉血管的二值化处理，二值化计算方式为：

其中，f_ij为某一图像序列的第i行第j列像素的二值化的数值；α为曲率均值二值化系数，取值范围为[0，1]；

光程计算模块，用于根据增强处理和二值化处理的图像数字信号的二值化有效面积计算由动脉血管搏动而产生的动脉血管的扩张光程和收缩光程，具体计算公式为：

S_max＝max{S₁，S₂，S₃，...，S_A} (10)

S_min＝min{S₁，S₂，S₃，...，S_A} (11)

L_max＝S_max (13)

L_min＝S_min (16)

公式中，S_m为第m个图像序列的动脉血管二值化有效面积；为第m个图像序列的第i行第j列像素的二值化数值；S_max为动脉血管扩张时的二值化有效面积；S_min为动脉血管收缩时的二值化有效面积；A为获取的图像数字信号的数量；L_max为动脉血管扩张时的光程；L_min为动脉血管收缩时的光程；

动脉血氧饱和度及氧合血红蛋白浓度计算模块，用于根据动脉血管的扩张光程和收缩光程计算动脉血氧饱和度和氧合血红蛋白浓度，用于反应人体的血液氧含量状况，具体步骤包括：

4.1)计算动脉血管扩张时的透射光的强度，计算公式如下：

公式中，I_A为动脉血管扩张时的透射光的强度；I₀为动脉血管扩张时的入射光的强度；γ为组织内的非脉动成分及静脉血的吸光系数；ε_HbO2为表示动脉血管中氧合血红蛋白的吸光系数；C_HbO2为表示动脉血管中氧合血红蛋白的浓度；ε_Hb为表示动脉血管中还原血红蛋白的吸光系数；C_Hb为表示动脉血管中还原血红蛋白的浓度；为动脉扩张图像序列的第i行第j列像素的灰度值；/>为动脉扩张图像序列的第i行第j列像素的二值化的数值；r₁为动脉扩张图像序列的图像行数；c₁为动脉扩张图像序列的图像列数；

4.2)计算动脉血管收缩时的透射光的强度，计算公式如下：

公式中，I_B为动脉血管收缩时的透射光的强度；I₀’为动脉血管收缩时的入射光的强度；为动脉收缩图像序列的第i行第j列像素的灰度值；/>为动脉收缩图像序列的第i行第j列像素的二值化的数值；r₂为动脉收缩图像序列的图像行数；c₂为动脉收缩图像序列的图像列数；

4.3)对公式(12)～(14)进行整合和取对数，得到公式(18)：

4.4)对公式(15)～(17)进行整合和取对数，得到公式(19)：

4.5)对公式(18)～(19)进行整合，得到公式(20)和公式(21)，用于计算动脉血管的搏动信息的量化参数：

C_HbO2ε_HbO2+C_Hbε_Hb＝R (20)

公式中，R为动脉血管的搏动信息的量化参数，通过公式(21)计算具体值；

4.6)根据公式(20)计算动脉血管中氧合血红蛋白的浓度：

其中，C₀为表示动脉血管中血红蛋白的浓度，其计算公式为：

C_HbO2+C_Hb＝C₀ (23)

公式中，C_HbO2为表示动脉血管中氧合血红蛋白的浓度；C_Hb为表示动脉血管中还原血红蛋白的浓度；

4.7)基于动脉血管中氧合血红蛋白的浓度计算动脉血氧饱和度，计算公式为：

其中，C₀为动脉血管中血红蛋白的浓度与动脉血管的搏动信息的量化参数R存在函数关系，函数关系为：

R＝f(C₀) (25)

公式中，f为动脉血管的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀的函数关系，所述的搏动信息的量化参数R与动脉血管中血红蛋白的浓度C₀存在线性关系，即：

R＝kC₀+B (26)

其中，k为线性系数；B为漂移常量；

将公式(26)代入公式(22)和公式(24)，可得到动脉血管中氧合血红蛋白的浓度和动脉血氧饱和度的计算公式：