CN110710982B - 用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法 - Google Patents

用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法 Download PDF

Info

Publication number
CN110710982B
CN110710982B CN201910988379.4A CN201910988379A CN110710982B CN 110710982 B CN110710982 B CN 110710982B CN 201910988379 A CN201910988379 A CN 201910988379A CN 110710982 B CN110710982 B CN 110710982B
Authority
CN
China
Prior art keywords
wavelength
light
blood
absorbance
hemoglobin concentration
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
Application number
CN201910988379.4A
Other languages
English (en)
Other versions
CN110710982A (zh
Inventor
张珣
梁思阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BOE Technology Group Co Ltd
Original Assignee
BOE Technology Group Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BOE Technology Group Co Ltd filed Critical BOE Technology Group Co Ltd
Priority to CN201910988379.4A priority Critical patent/CN110710982B/zh
Publication of CN110710982A publication Critical patent/CN110710982A/zh
Application granted granted Critical
Publication of CN110710982B publication Critical patent/CN110710982B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/145Measuring 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/1455Measuring 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7246Details of waveform analysis using correlation, e.g. template matching or determination of similarity

Abstract

本发明的实施例提供一种用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法,涉及血液检测领域,可以在无创伤的情况下,提高检测血红蛋白浓度的准确度。一种用于检测血红蛋白浓度的模型的获取方法,包括:建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式,
Figure DDA0002237440310000011
将两组测试数据分别代入第一关系式,计算得到x和y的值;其中,每组测试数据均包括血红蛋白浓度值、血液对第一波长的光的吸光度值、血液对第二波长的光的吸光度值;设x=a,y=b,a和b为定值,将a和b代入第一关系式,得到用于检测血红蛋白浓度的模型,
Figure DDA0002237440310000012

Description

用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的 检测方法
技术领域
本发明涉及血液检测领域,尤其涉及一种用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法。
背景技术
临床中,检测贫血以及指导输血时,均需对血液中的血红蛋白浓度进行检测。现有方法是将血液样品与氰化高铁混合,通过形成氰化高铁血红蛋白标准液,在540nm的光源照射下,检测吸光度,进而计算出对应的血红蛋白浓度。
该方法为有创方法,在获取受检者的血液样品时,轻则导致受检者疼痛,重则会带来感染风险,因此,如何通过无创方法,准确检测血液中的血红蛋白浓度是本领域技术人员亟待解决的问题。
发明内容
本发明的实施例提供一种用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法,可以在无创伤的情况下,提高检测血红蛋白浓度的准确度。
为达到上述目的,本发明的实施例采用如下技术方案:
一方面,提供一种用于检测血红蛋白浓度的模型的获取方法,包括:
建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式,
Figure BDA0002237440290000011
其中,所述第一波长λ1和所述第二波长λ2不同;
Figure BDA0002237440290000012
为血液对第一波长的光的吸光度,
Figure BDA0002237440290000013
为血液对第二波长的光的吸光度;CtHb为所述血红蛋白浓度,包括还原血红蛋白浓度和氧合血红蛋白浓度;
将两组测试数据分别代入第一关系式,计算得到x和y的值;其中,每组测试数据均包括血红蛋白浓度值、血液对第一波长的光的吸光度值、血液对第二波长的光的吸光度值;
设x=a,y=b,a和b为定值,将a和b代入所述第一关系式,得到用于检测所述血红蛋白浓度的模型,
Figure BDA0002237440290000021
可选地,在建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式之后,在将两组测试数据分别代入第一关系式,计算得到x和y的值之前,所述用于检测血红蛋白浓度的模型的获取方法还包括:
获取两组测试数据;
其中,每组测试数据的获取方法包括:获取血液在所述第一波长的光照射下的第一光电容积脉搏描记信号和在第二波长的光照射下的第二光电容积脉搏描记信号;
根据所述第一光电容积脉搏描记信号的最大值
Figure BDA0002237440290000022
和最小值
Figure BDA0002237440290000023
得到血液对所述第一波长的光的吸光度为
Figure BDA0002237440290000024
根据所述第二光电容积脉搏描记信号的最大值
Figure BDA0002237440290000025
和最小值
Figure BDA0002237440290000026
得到血液对所述第二波长的光的吸光度为
Figure BDA0002237440290000027
再一方面,提供一种血红蛋白浓度的检测方法,包括:
获取受检血液在第一波长的光照射下的第一受检光电容积脉搏描记信号和在第二波长的光照射下的第二受检光电容积脉搏描记信号;其中,所述第一波长和所述第二波长不同;
分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和所述受检血液对所述第二波长的光的吸光度;
将所述受检血液对所述第一波长的光的吸光度和对所述第二波长的光的吸光度,代入如上所述的用于检测血红蛋白浓度的模型的获取方法所获得的用于检测所述血红蛋白浓度的模型中,计算得到所述受检血液的血红蛋白浓度。
可选地,分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和所述受检血液对所述第二波长的光的吸光度,包括:
根据所述第一受检光电容积脉搏描记信号的最大值
Figure BDA0002237440290000031
和最小值
Figure BDA0002237440290000032
得到所述受检血液对所述第一波长的光的吸光度为
Figure BDA0002237440290000033
根据所述第二受检光电容积脉搏描记信号的最大值
Figure BDA0002237440290000034
和最小值
Figure BDA0002237440290000035
得到所述受检血液对所述第二波长的光的吸光度为
Figure BDA0002237440290000036
可选地,分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和对所述第二波长的光的吸光度,包括:
分别获取所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号的自然对数,得到对应的第一中间光电容积脉搏描记信号和第二中间光电容积脉搏描记信号;
利用块稀疏贝叶斯对所述第一中间光电容积脉搏描记信号、所述第二中间光电容积脉搏描记信号进行稀疏分解,得到所述第一中间光电容积脉搏描记信号对应的第一频谱、所述第二中间光电容积脉搏描记信号对应的第二频谱;
将所述第一频谱与所述第二频谱叠加得到合成频谱;
获取所述合成频谱的最高峰值点;
获取所述最高峰值点对应的频率,作为最优频率;
获取所述第一频谱中所述最优频率对应的幅值,作为所述受检血液对所述第一波长的光的吸光度;
获取所述第二频谱中所述最优频率对应的幅值,作为所述受检血液对所述第二波长的光的吸光度。
可选地,在获取受检血液在第一波长的光照射下的第一受检光电容积脉搏描记信号和在第二波长的光照射下的第二受检光电容积脉搏描记信号之后,在分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和所述受检血液对所述第二波长的光的吸光度之前,所述血红蛋白浓度的检测方法还包括:
分别对第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号进行带通滤波;
其中,频率范围为0.4~5Hz。
可选地,获取所述合成频谱的最高峰值点,包括:获取所述合成频谱在0.4~5Hz的频率范围内的最高峰值点。
另一方面,提供一种计算机设备,包括存储单元和处理单元;所述存储单元中存储可在所述处理单元上运行的计算机程序;所述处理单元执行所述计算机程序时实现如上所述的用于检测血红蛋白浓度的模型的获取方法或如上所述的血红蛋白浓度的检测方法。
又一方面,提供一种计算机可读介质,其存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如上所述的用于检测血红蛋白浓度的模型的获取方法或如上所述的血红蛋白浓度的检测方法。
本发明实施例提供一种用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法,通过简化血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度之间的关系,建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式
Figure BDA0002237440290000041
然后,通过将两组测试数据分别代入第一关系式,计算得到x和y的值,精确度较高;再设x=a,y=b,a和b为定值,将a和b代入第一关系式,即可得到用于检测血红蛋白浓度的模型,
Figure BDA0002237440290000042
由此,所得的用于检测血红蛋白浓度的模型在后续用于检测时,既对受检者不会产生伤害,又可以简单准确的得到血红蛋白浓度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为一种光电容积脉搏波信号的示意图;
图2为本发明实施例提供的一种用于检测血红蛋白浓度的模型的获取方法的流程示意图;
图3为本发明实施例提供的另一种用于检测血红蛋白浓度的模型的获取方法的流程示意图;
图4为本发明实施例提供的一种血红蛋白浓度的检测方法的流程示意图;
图5为本发明实施例提供的另一种血红蛋白浓度的检测方法的流程示意图;
图6为本发明实施例提供的再一种血红蛋白浓度的检测方法的流程示意图;
图7为一种低灌注度情况下的光电容积脉搏波信号的时域与频域示意图;
图8为本发明实施例提供的又一种血红蛋白浓度的检测方法的流程示意图;
图9为本发明实施例提供的又一种血红蛋白浓度的检测方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在心脏搏动周期内,外周血管中的微动脉、毛细血管和微静脉内流过的血液相应的呈搏动性变化。当心脏收缩时血液容积最大,而在心脏舒张时血液容积最小,血液容积这种脉动性变化一般可通过光电容积探测器获得,所得的波形中含有容积脉搏血流信息。由此,可以通过光学技术获取光线在人体组织中的吸收变化的描记图,以得到光电容积脉搏波信号,该方法称为光电容积脉搏波描记法(Photo Plethysmo Graphy,PPG)。
光电容积脉搏波描记法的工作原理是:将特定强度和波长的光束照射到指端皮肤表面时,光束将通过透射或反射方式传送到光电容积探测器,光电容积探测器根据检测到的光强进行描绘,以描绘出如图1所示的交流形式的光电容积脉搏波信号。
在此过程中,由于受到指端皮肤肌肉组织和血液的吸收衰减作用,光电容积探测器检测到的光强度将会减弱。其中,皮肤、肌肉和组织等对光的吸收在整个血液循环中是保持恒定不变的,经光电容积探测器检测到的就是直流分量;而皮肤内的血液容积在心脏收缩舒张作用下呈搏动性变化。当心脏收缩时外周血管血容量最多,光吸收量也最大,检测到的光强度最小;而在心脏舒张时外周血管血容量最少,光吸收量最小,检测到的光强度最大,这样使得光电容积探测器检测到的光强度随之呈脉动性变化,并具有最大值和最小值。
基于此,首先,根据朗伯比尔定律(Lambert-Beer law)可知,吸光度Aλi的修正公式为:
Figure BDA0002237440290000061
其中,
Figure BDA0002237440290000062
表示光源发出的波长为λi的光的光强,
Figure BDA0002237440290000063
表示光电容积探测器接收到波长为λi的光的光强,DPF表示光传播的差分路径因子,ρ表示光源到光电容积探测器的水平距离,1≤m≤M,且m为整数,M代表吸光物质的总数,
Figure BDA0002237440290000064
表示第m吸光物质的吸光系数,Cm表示第m吸光物质的浓度,G表示背景吸收、散射引起的衰减量,可以代表非动脉成分(例如脂肪、肌肉、骨骼)的吸光度。
其次,一般动脉具有搏动部分和静止部分,而其中动脉的搏动部分会对入射的光的光强产生影响,由此,动脉血液在搏动时的吸光度
Figure BDA0002237440290000071
可表示为:
Figure BDA0002237440290000072
其中,
Figure BDA0002237440290000073
表示动脉血液的最大吸光度,
Figure BDA0002237440290000074
表示动脉血液的最小吸光度,
Figure BDA0002237440290000075
表示光电容积探测器检测到的波长为λi的光的最小光强,
Figure BDA0002237440290000076
表示光电容积探测器检测到的波长为λi的光的最大光强。
然后,动脉血液中的还原血红蛋白和氧合血红蛋白对光的吸收影响较大,由此可认为,吸光物质为还原血红蛋白和氧合血红蛋白。
结合公式(1)和公式(2)可得:
Figure BDA0002237440290000077
其中,DPF1表示光电容积描记信号的最小光强对应的差分路径因子,DPF2表示光电容积描记信号的最大光强对应的差分路径因子,
Figure BDA0002237440290000078
表示还原血红蛋白对应波长为λi的光的吸光系数,
Figure BDA0002237440290000079
表示氧合血红蛋白对应波长为λi的光的吸光系数,CRHb表示还原血红蛋白的浓度,CHbO2表示氧合血红蛋白的浓度。
基于上述,本发明的实施例提供一种用于检测血红蛋白浓度的模型的获取方法,如图2所示,包括:
S10、建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式,
Figure BDA00022374402900000710
其中,第一波长λ1和第二波长λ2不同;
Figure BDA00022374402900000711
为血液对第一波长的光的吸光度,
Figure BDA00022374402900000712
为血液对第二波长的光的吸光度;CtHb为血红蛋白浓度,包括还原血红蛋白浓度和氧合血红蛋白浓度。
S20、将两组测试数据分别代入第一关系式,计算得到x和y的值;其中,每组测试数据均包括血红蛋白浓度值、血液对第一波长的光的吸光度值、血液对第二波长的光的吸光度值。
S30、设x=a,y=b,a和b为定值,将a和b代入第一关系式,得到用于检测血红蛋白浓度的模型,
Figure BDA0002237440290000081
利用两路光源检测血红蛋白浓度,第一光源对应第一波长λ1,第二光源对应第二波长λ2,第一光源和第二光源可以分别为发光二极管(Light Emtting Diode,LED)。
其中,第一光源与光电容积探测器之间的水平距离、第二光源与光电容积探测器之间的水平距离相等,并且该水平距离可以是预先设定好的距离,或者也可以是通过仪器测量处理的距离。
需要说明的是,可以使第一光源和第二光源分时发光,从而使光电容积探测器可以分时接收第一光源和第二光源照射指端皮肤表层后出射的光,以描绘出对应第一波长λ1的第一光电容积脉搏描记信号和对应第二波长λ2的第二光电容积脉搏描记信号。
此时,根据第一光电容积脉搏描记信号和第二光电容积脉搏描记信号得到血红蛋白浓度的过程如下:
第一步:第一光电容积脉搏描记信号的最大值,即为光电容积探测器检测到的第一波长λ1的光的最大光强
Figure BDA0002237440290000082
第一光电容积脉搏描记信号的最小值即为光电容积探测器检测到的第一波长λ1的光的最小光强
Figure BDA0002237440290000083
Figure BDA0002237440290000084
Figure BDA0002237440290000085
代入公式(2),可得:
Figure BDA0002237440290000086
代入公式(3)可得:
Figure BDA0002237440290000087
第二步:第二光电容积脉搏描记信号的最大值,即为光电容积探测器检测到的第二波长λ2的光的最大光强
Figure BDA0002237440290000088
第二光电容积脉搏描记信号的最小值,即为光电容积探测器检测到的第二波长λ2的光的最小光强
Figure BDA0002237440290000089
Figure BDA00022374402900000810
Figure BDA00022374402900000811
代入公式(2),可得:
Figure BDA0002237440290000091
代入公式(3)可得:
Figure BDA0002237440290000092
第三步:将公式(5)的等式两边同乘以
Figure BDA0002237440290000093
可得:
Figure BDA0002237440290000094
将公式(7)等式两边同乘以
Figure BDA0002237440290000095
可得:
Figure BDA0002237440290000096
利用公式(8)减去公式(9):
Figure BDA0002237440290000097
由此,可得血红蛋白浓度的公式为:
Figure BDA0002237440290000098
第四步:获取差分路径因子、还原血红蛋白对应第一波长λ1的光的吸光系数、还原血红蛋白对应第二波长的光λ2的吸光系数、氧合血红蛋白对应第一波长λ1的光的吸光系数、氧合血红蛋白对应第二波长的光λ2的吸光系数,光源到光电容积探测器的水平距离,将其全部数据及计算得到的血液分别对第一波长λ1的光的吸光度和对第二波长λ2的光的吸光度,代入公式(10),然后计算得到血红蛋白浓度。
然而,上述第四步中需要获取的数据较多,过程繁琐,并且,若差分路径因子、还原血红蛋白对应第一波长λ1的光的吸光系数、还原血红蛋白对应第二波长的光λ2的吸光系数、氧合血红蛋白对应第一波长λ1的光的吸光系数、氧合血红蛋白对应第二波长的光λ2的吸光系数,光源到光电容积探测器的水平距离中的任一项产生误差,将导致血红蛋白浓度的结果不准确。
又结合公式(10)可知,在用于检测血红蛋白的第一光源和第二光源出射的光的波长不改变的情况下,
Figure BDA0002237440290000101
Figure BDA0002237440290000102
不会改变,此时,若能通过简单的方法直接得到该两项的值,将简化上述第四步的过程以及简化上述公式(10),在后续检测中,只需计算得到血液分别对第一波长λ1的光的吸光度和对第二波长λ2的光的吸光度两个值,即可代入简化后的公式(10),简单又准确的得到血红蛋白浓度。
基于此,结合公式(10)的代数形式,建立血红蛋白浓度与血液对第一波长λ1的光的吸光度、血液对第二波长λ2的光的吸光度的第一关系式为:
Figure BDA0002237440290000103
然后,利用两组测试数据(血红蛋白浓度值、血液对第一波长的光的吸光度值、血液对第二波长的光的吸光度值)分别代入第一关系式,联立并反向推导出x和y。再设x=a,y=b,将a和b代入第一关系式,即可得到用于检测血红蛋白浓度的模型,
Figure BDA0002237440290000104
由于a和b已知,在后续利用该模型时,计算将简单又准确。
本发明实施例提供一种用于检测血红蛋白浓度的模型的获取方法,通过简化血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度之间的关系,建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式
Figure BDA0002237440290000105
然后,通过将两组测试数据分别代入第一关系式,计算得到x和y的值,精确度较高;再设x=a,y=b,a和b为定值,将a和b代入第一关系式,即可得到用于检测血红蛋白浓度的模型,
Figure BDA0002237440290000111
由此,所得的用于检测血红蛋白浓度的模型在后续用于检测时,既对受检者不会产生伤害,又可以简单准确的得到血红蛋白浓度。
可选地,如图3所示,在上述S10之后,在S20之前,用于检测血红蛋白浓度的模型的获取方法还包括:
S11、获取两组测试数据。
需要说明的是,该两组测试数据可以通过对一个人检测两次,或者,针对两个人分别检测,只要得到两组测试数据即可,具体方式本发明对此不作限定。
其中,如图4所示,每组测试数据的获取方法包括:
S111、获取血液在第一波长的光照射下的第一光电容积脉搏描记信号和在第二波长的光照射下的第二光电容积脉搏描记信号。
S112、根据第一光电容积脉搏描记信号的最大值
Figure BDA0002237440290000112
和最小值
Figure BDA0002237440290000113
得到血液对第一波长的光的吸光度为
Figure BDA0002237440290000114
S113、根据第二光电容积脉搏描记信号的最大值
Figure BDA0002237440290000115
和最小值
Figure BDA0002237440290000116
得到血液对第二波长的光的吸光度为
Figure BDA0002237440290000117
此外,每组测试数据的获取方法还包括:
通过有创方法得到血红蛋白浓度值。
示例的,以下提供多个步骤,以清楚说明用于检测血红蛋白浓度的模型的获取方法。
第一步:建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式,
Figure BDA0002237440290000118
第二步:获取第一组测试数据。
针对第一位测试者P,通过有创方法得到该测试者P的血红蛋白浓度值为CtHb_P,通过光电容积探测器测量得到该测试者P的血液在第一波长λ1的光照射下的第一光电容积脉搏描记信号为PPG1_P和在第二波长λ2的光照射下的第二光电容积脉搏描记信号为PPG2_P
第一光电容积脉搏描记信号PPG1_P的最大值为
Figure BDA0002237440290000121
和最小值为
Figure BDA0002237440290000122
得到测试者P的血液对第一波长λ1的光的吸光度值为
Figure BDA0002237440290000123
第二光电容积脉搏描记信号PPG2_P的最大值
Figure BDA0002237440290000124
和最小值
Figure BDA0002237440290000125
得到测试者P的血液对第二波长λ2的光的吸光度值为
Figure BDA0002237440290000126
第三步:获取第二组测试数据。
针对第二位测试者Q,通过有创方法得到该测试者Q的血红蛋白浓度值为CtHb_Q,通过光电容积探测器测量得到该测试者Q的血液在第一波长λ1的光照射下的第一光电容积脉搏描记信号为PPG1_Q和在第二波长λ2的光照射下的第二光电容积脉搏描记信号为PPG2_Q
第一光电容积脉搏描记信号PPG1_Q的最大值为
Figure BDA0002237440290000127
和最小值为
Figure BDA0002237440290000128
得到测试者Q的血液对第一波长λ1的光的吸光度值为
Figure BDA0002237440290000129
第二光电容积脉搏描记信号PPG2_Q的最大值
Figure BDA00022374402900001210
和最小值
Figure BDA00022374402900001211
得到测试者Q的血液对第二波长λ2的光的吸光度值为
Figure BDA00022374402900001212
第四步:将第一组测试数据代入第一关系式,得到
Figure BDA00022374402900001213
将第二组测试数据代入第一关系式,得到
Figure BDA0002237440290000131
Figure BDA0002237440290000132
根据公式(11)和公式(12),分别计算得到
Figure BDA0002237440290000133
Figure BDA0002237440290000134
第五步:设x=a,y=b,则得到用于检测血红蛋白浓度的模型,
Figure BDA0002237440290000135
其中,
Figure BDA0002237440290000136
本发明的实施例还提供一种血红蛋白浓度的检测方法,如图5所示,包括:
S100、获取受检血液在第一波长的光照射下的第一受检光电容积脉搏描记信号和在第二波长的光照射下的第二受检光电容积脉搏描记信号。
其中,第一波长和第二波长不同。
S200、分别根据第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号,得到受检血液对第一波长的光的吸光度和对第二波长的光的吸光度。
S300、将受检血液对第一波长的光的吸光度和对第二波长的光的吸光度,代入如上所述的用于检测血红蛋白浓度的模型的获取方法所获得的用于检测血红蛋白浓度的模型中,计算得到受检血液的血红蛋白浓度。
本发明的实施例提供一种血红蛋白浓度的检测方法,通过获取受检血液在第一波长的光照射下的第一受检光电容积脉搏描记信号和在第二波长的光照射下的第二受检光电容积脉搏描记信号;分别根据第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号,得到受检血液对第一波长的光的吸光度和对第二波长的光的吸光度;然后,受检血液对第一波长的光的吸光度和对第二波长的光的吸光度,代入如上所述的用于检测血红蛋白浓度的模型的获取方法所获得的用于检测血红蛋白浓度的模型中,即可简单又准确的计算得到受检血液的血红蛋白浓度。
可选地,如图6所示,上述S200包括:
S201、根据第一受检光电容积脉搏描记信号的最大值
Figure BDA0002237440290000141
和最小值
Figure BDA0002237440290000142
得到受检血液对第一波长的光的吸光度为
Figure BDA0002237440290000143
S202、根据第二受检光电容积脉搏描记信号的最大值
Figure BDA0002237440290000144
和最小值
Figure BDA0002237440290000145
得到受检血液对第二波长的光的吸光度为
Figure BDA0002237440290000146
由此,将受检血液对第一波长的光的吸光度
Figure BDA0002237440290000147
以及对第二波长的光的吸光度
Figure BDA0002237440290000148
代入用于检测血红蛋白浓度的模型的获取方法所获得的用于检测所述血红蛋白浓度的模型
Figure BDA0002237440290000149
中,可得受检血液的血红蛋白浓度为
Figure BDA00022374402900001410
针对一般情况,利用S201、S202的方法就可以检测出受检血液的血红蛋白浓度。但是,对贫血或者在输血的受检者进行血红蛋白浓度检测时,如图7所示,由于该受检者出现低灌注状态(即,体内重要脏器微循环血液灌注在短时间内急剧降低的低血流状态),导致光电容积探测器采集的交流信号大幅下降,在噪声的干扰下,有效的信息几乎会被噪声信号淹没,信噪比极低。此时,需要对第一光电容积脉搏描记信号第二光电容积脉搏描记信号进行处理。
可选地,如图8所示,上述S200包括:
S211、分别获取第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号的自然对数,得到对应的第一中间光电容积脉搏描记信号和第二中间光电容积脉搏描记信号。
其中,获取自然对数指的是获取以常数e为底数的对数。
示例的,若第一受检光电容积脉搏描记信号为
Figure BDA00022374402900001411
则对应的第一中间光电容积脉搏描记信号表示为
Figure BDA00022374402900001412
需要说明的是,若第一受检光电容积脉搏描记信号的最大值为
Figure BDA0002237440290000151
和最小值为
Figure BDA0002237440290000152
得到受检血液对第一波长的光的吸光度为
Figure BDA0002237440290000153
Figure BDA0002237440290000154
等于第一中间光电容积脉搏描记信号
Figure BDA0002237440290000155
的幅值(即,
Figure BDA0002237440290000156
的最大值减去最小值)。由此反向推理,只要计算得到第一中间光电容积脉搏描记信号的幅值,即等于直接得到了受检血液对第一波长的光的吸光度。
同理,只要计算得到第二中间光电容积脉搏描记信号的幅值,即等于直接得到了受检血液对第二波长的光的吸光度。
S212、利用块稀疏贝叶斯对第一中间光电容积脉搏描记信号、第二中间光电容积脉搏描记信号进行稀疏分解,得到第一中间光电容积脉搏描记信号对应的第一频谱、所述第二中间光电容积脉搏描记信号对应的第二频谱。
需要说明的是,利用块稀疏贝叶斯可以实现第一中间光电容积脉搏描记信号和第二中间光电容积脉搏描记信号的重构,得到的对应的第一频谱和第二频谱具有较高的频谱分辨率,可以区分出相近的周期性噪声,减小噪声对血红蛋白浓度计算所产生的干扰。
块稀疏贝叶斯表示为Y=Φn,mX+V。
其中,Y指的是第一中间光电容积脉搏描记信号和第二中间光电容积脉搏描记信号组成的N×2维信号向量;X指的是M×2维待求解向量;
Figure BDA0002237440290000157
指的是N×M维的离散傅里叶变换基矩阵,0≤n≤N-1,0≤m≤M-1,n和m均为整数;V为未知的噪声向量。
由此,根据Y和Φn,m,求解X,得到第一中间光电容积脉搏描记信号对应的第一频谱和第二中间光电容积脉搏描记信号对应的第二频谱。
S213、将第一频谱与第二频谱叠加,得到合成频谱。
可以理解的是,将第一频谱与第二频谱叠加,可以抑制非周期噪声。
S214、获取合成频谱的最高峰值点。
可选地,上述S214,包括:
获取合成频谱在0.4~5Hz的频率范围内的最高峰值点。
S215、获取最高峰值点对应的频率,作为最优频率。
S216、获取第一频谱中最优频率对应的幅值,作为受检血液对第一波长的光的吸光度。
S217、获取第二频谱中最优频率对应的幅值,作为受检血液对第二波长的光的吸光度。
可选地,在上述S211之前,如图9所示,血红蛋白浓度的检测方法还包括:
S210、分别对第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号进行带通滤波。
其中,频率范围为0.4~5Hz。
去除0.4~5Hz以外的噪声以及运动干扰,以使得第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号具有初步的稀疏性。
本发明的实施例还提供一种计算机设备,包括存储单元和处理单元;存储单元中存储可在处理单元上运行的计算机程序;处理单元执行计算机程序时实现如上所述的用于检测血红蛋白浓度的模型的获取方法或如上所述的血红蛋白浓度的检测方法。
本发明的实施例还提供一种计算机可读介质,其存储有计算机程序,所述计算机程序被处理器执行时实现如上所述的用于检测血红蛋白浓度的模型的获取方法或如上所述的血红蛋白浓度的检测方法。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (8)

1.一种用于检测血红蛋白浓度的模型的获取方法,其特征在于,包括:
建立血红蛋白浓度与血液对第一波长的光的吸光度、血液对第二波长的光的吸光度的第一关系式,
Figure FDF0000019351490000011
其中,所述第一波长λ1和所述第二波长λ2不同;
Figure FDF0000019351490000012
为血液对所述第一波长的光的吸光度,
Figure FDF0000019351490000013
为血液对所述第二波长的光的吸光度;CtHb为所述血红蛋白浓度,包括还原血红蛋白浓度和氧合血红蛋白浓度;
获取两组测试数据;其中,每组测试数据的获取方法包括:获取血液在所述第一波长的光照射下的第一光电容积脉搏描记信号和在所述第二波长的光照射下的第二光电容积脉搏描记信号;根据所述第一光电容积脉搏描记信号的最大值
Figure FDF0000019351490000014
和最小值
Figure FDF0000019351490000015
得到血液对所述第一波长的光的吸光度为
Figure FDF0000019351490000016
根据所述第二光电容积脉搏描记信号的最大值
Figure FDF0000019351490000017
和最小值
Figure FDF0000019351490000018
得到血液对所述第二波长的光的吸光度为
Figure FDF0000019351490000019
将所述两组测试数据分别代入所述第一关系式,计算得到x和y的值;其中,每组测试数据均包括血红蛋白浓度值、血液对所述第一波长的光的吸光度值、血液对所述第二波长的光的吸光度值;
设x=a,y=b,a和b为定值,将a和b代入所述第一关系式,得到用于检测所述血红蛋白浓度的模型,
Figure FDF0000019351490000021
2.一种血红蛋白浓度的检测方法,其特征在于,包括:
获取受检血液在第一波长的光照射下的第一受检光电容积脉搏描记信号和在第二波长的光照射下的第二受检光电容积脉搏描记信号;其中,所述第一波长和所述第二波长不同;
分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和所述受检血液对所述第二波长的光的吸光度;
将所述受检血液对所述第一波长的光的吸光度和对所述第二波长的光的吸光度,代入如权利要求1所述的用于检测血红蛋白浓度的模型的获取方法所获得的用于检测所述血红蛋白浓度的模型中,计算得到所述受检血液的血红蛋白浓度。
3.根据权利要求2所述的血红蛋白浓度的检测方法,其特征在于,分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和所述受检血液对所述第二波长的光的吸光度,包括:
根据所述第一受检光电容积脉搏描记信号的最大值
Figure FDF0000019351490000022
和最小值
Figure FDF0000019351490000023
得到所述受检血液对所述第一波长的光的吸光度为
Figure FDF0000019351490000031
根据所述第二受检光电容积脉搏描记信号的最大值
Figure FDF0000019351490000032
和最小值
Figure FDF0000019351490000033
得到所述受检血液对所述第二波长的光的吸光度为
Figure FDF0000019351490000034
4.根据权利要求2所述的血红蛋白浓度的检测方法,其特征在于,分别根据所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号,得到所述受检血液对所述第一波长的光的吸光度和对所述第二波长的光的吸光度,包括:
分别获取所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号的自然对数,得到对应的第一中间光电容积脉搏描记信号和第二中间光电容积脉搏描记信号;
利用块稀疏贝叶斯对所述第一中间光电容积脉搏描记信号、所述第二中间光电容积脉搏描记信号进行稀疏分解,得到所述第一中间光电容积脉搏描记信号对应的第一频谱、所述第二中间光电容积脉搏描记信号对应的第二频谱;
将所述第一频谱与所述第二频谱叠加,得到合成频谱;
获取所述合成频谱的最高峰值点;
获取所述最高峰值点对应的频率,作为最优频率;
获取所述第一频谱中所述最优频率对应的幅值,作为所述受检血液对所述第一波长的光的吸光度;
获取所述第二频谱中所述最优频率对应的幅值,作为所述受检血液对所述第二波长的光的吸光度。
5.根据权利要求4所述的血红蛋白浓度的检测方法,其特征在于,分别获取所述第一受检光电容积脉搏描记信号和所述第二受检光电容积脉搏描记信号的自然对数,得到对应的第一中间光电容积脉搏描记信号和第二中间光电容积脉搏描记信号之前,所述血红蛋白浓度的检测方法还包括:
分别对第一受检光电容积脉搏描记信号和第二受检光电容积脉搏描记信号进行带通滤波;
其中,频率范围为0.4~5Hz。
6.根据权利要求4或5所述的血红蛋白浓度的检测方法,其特征在于,获取所述合成频谱的最高峰值点,包括:
获取所述合成频谱在0.4~5Hz的频率范围内的最高峰值点。
7.一种计算机设备,其特征在于,包括存储单元和处理单元;所述存储单元中存储可在所述处理单元上运行的计算机程序;所述处理单元执行所述计算机程序时实现如权利要求1所述的用于检测血红蛋白浓度的模型的获取方法,和/或如权利要求2-6任一项所述的血红蛋白浓度的检测方法。
8.一种计算机可读介质,其存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1所述的用于检测血红蛋白浓度的模型的获取方法,和/或如权利要求2-6任一项所述的血红蛋白浓度的检测方法。
CN201910988379.4A 2019-10-17 2019-10-17 用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法 Active CN110710982B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910988379.4A CN110710982B (zh) 2019-10-17 2019-10-17 用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910988379.4A CN110710982B (zh) 2019-10-17 2019-10-17 用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法

Publications (2)

Publication Number Publication Date
CN110710982A CN110710982A (zh) 2020-01-21
CN110710982B true CN110710982B (zh) 2023-01-13

Family

ID=69211803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910988379.4A Active CN110710982B (zh) 2019-10-17 2019-10-17 用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法

Country Status (1)

Country Link
CN (1) CN110710982B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112545461A (zh) * 2020-12-05 2021-03-26 深圳市美的连医疗电子股份有限公司 一种无创血红蛋白浓度值的检测方法、装置、系统及计算机可读存储介质

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH01174940A (ja) * 1987-12-29 1989-07-11 Shimadzu Corp ヘモグロビン測定法
KR100612827B1 (ko) * 2001-04-19 2006-08-14 삼성전자주식회사 비 침습적인 헤모글로빈 농도와 산소 포화도 모니터링방법 및 장치
JP2005253478A (ja) * 2002-03-18 2005-09-22 Citizen Watch Co Ltd ヘモグロビン分析装置
SE530244C2 (sv) * 2006-05-05 2008-04-08 Hemocue Ab Förfarande och system för kvantitativ hemoglobinbestämning
EP2994042B1 (en) * 2013-05-10 2023-09-27 University Of Utah Research Foundation Devices, systems, and methods for measuring blood loss
CN104931441A (zh) * 2015-06-15 2015-09-23 马国鹭 一种快速测量血红蛋白的方法和装置
CN105067550B (zh) * 2015-07-30 2017-12-26 中北大学 一种基于分块稀疏贝叶斯优化的红外光谱波长选择方法
EP3834710A1 (en) * 2015-09-25 2021-06-16 Sanmina Corporation System and method for health monitoring using a non-invasive, multi-band biosensor
CN108288295A (zh) * 2018-01-30 2018-07-17 深圳大学 基于结构信息的红外小目标图像的快速重构方法及系统
CN109222992B (zh) * 2018-08-24 2021-04-30 河南大学 基于光学反射式原理的血红蛋白浓度无创测量方法
CN109106376B (zh) * 2018-09-20 2022-04-19 京东方科技集团股份有限公司 一种血液中总血红蛋白浓度的检测方法及装置

Also Published As

Publication number Publication date
CN110710982A (zh) 2020-01-21

Similar Documents

Publication Publication Date Title
US6002952A (en) Signal processing apparatus and method
US11375906B2 (en) Method and apparatus for detecting fetal blood oxygen saturation, computer-readable storage medium and computer device
JP5096310B2 (ja) 身体の部位における血液の灌流を決定するための方法及び装置
CN103381094B (zh) 胎儿脉搏血氧饱和度监测系统及方法
Schmitt et al. Optical determination of dental pulp vitality
CN107920786B (zh) 脉搏血氧测定
JP2007209782A (ja) 信号対ノイズ比改善のための信号処理方法及び装置
JP2005516642A6 (ja) 信号対ノイズ比改善のための信号処理方法及び装置
CA2827981A1 (en) Regional saturation determination using photoacoustic technique
JP6667456B2 (ja) 対象のヘマトクリット値を非侵襲的に決定するプロセッサ、プログラム、デバイス及び方法
CN110710982B (zh) 用于检测血红蛋白浓度的模型的获取方法、血红蛋白浓度的检测方法
CN109106376B (zh) 一种血液中总血红蛋白浓度的检测方法及装置
Aarthi et al. Fingertip based estimation of heart rate using photoplethysmography
Johnson et al. A review of photoplethysmography-based physiological measurement and estimation, Part 2: multi-input methods
Song et al. Use of bi-level pulsed frequency-division excitation for improving blood oxygen saturation precision
US10653347B2 (en) Non-invasive measurement device and method for estimating local metabolic parameters
RU2221485C2 (ru) Устройство для неинвазивного измерения насыщения крови кислородом
RU2173082C1 (ru) Способ неинвазивного измерения насыщения крови кислородом
RU2805810C1 (ru) Носимое устройство с функцией определения концентрации гемоглобина, способ и система для определения концентрации гемоглобина
CN112932475B (zh) 血氧饱和度的计算方法、装置、电子设备及存储介质
Jayasree et al. Design and development of a simple hardware setup for sensing blood volume pulse and a PIC microcontroller based heart rate meter
Tatiparti et al. Smart non-invasive hemoglobin measurement using portable embedded technology
Prokop et al. Application of a multi-sensor set for comparative evaluation of the photoplethysmographic waveforms
Veerabathraswamy et al. Simulation of a Non-Contact Triple-Wavelength Reflectance Pulse Oximeter to Measure Oxygen and Carbon Monoxide Saturations
Vijeya Kaveri et al. Investigation on Optical Sensors for Heart Rate Monitoring

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