CN116978570A - 一种在线实时患者危重程度评估及生命体征参数预测系统 - Google Patents

一种在线实时患者危重程度评估及生命体征参数预测系统 Download PDF

Info

Publication number
CN116978570A
CN116978570A CN202311238621.9A CN202311238621A CN116978570A CN 116978570 A CN116978570 A CN 116978570A CN 202311238621 A CN202311238621 A CN 202311238621A CN 116978570 A CN116978570 A CN 116978570A
Authority
CN
China
Prior art keywords
patient
vital sign
sign parameter
sampling
death
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.)
Granted
Application number
CN202311238621.9A
Other languages
English (en)
Other versions
CN116978570B (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.)
Zhejiang Lab
Original Assignee
Zhejiang Lab
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 Zhejiang Lab filed Critical Zhejiang Lab
Priority to CN202311238621.9A priority Critical patent/CN116978570B/zh
Publication of CN116978570A publication Critical patent/CN116978570A/zh
Application granted granted Critical
Publication of CN116978570B publication Critical patent/CN116978570B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/60ICT specially adapted for the handling or processing of patient-related medical or healthcare data for patient-specific data, e.g. for electronic patient records
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Abstract

本发明公开了一种在线实时患者危重程度评估及生命体征参数预测系统,包括:患者生命体征参数获取模块;患者生理状态特征获取模块;采样实例获取模块,用于根据患者生理状态特征获取人体状态统计分布参数,配置累积测量权重,根据人体状态统计分布参数和累积测量权重获取当前的人体累积状态统计分布,对当前的人体累积状态统计分布进行采样,得到采样实例;死亡风险预测模块,用于获取患者死亡风险分布,患者死亡风险分布通过死亡率在每一区间的发生概率以及每一区间的期望死亡率表示;将死亡率在每一区间的发生概率与每一区间的期望死亡率的乘积之和作为当前患者的期望死亡率;将死亡率大于高死亡概率阈值的发生概率之和作为高死亡发生概率。

Description

一种在线实时患者危重程度评估及生命体征参数预测系统
技术领域
本发明属于医疗健康信息技术领域,尤其涉及一种在线实时患者危重程度评估及生命体征参数预测系统。
背景技术
死亡率是目前死亡风险评估最常用的指标,但是单纯的死亡率作为死亡风险评价指标并不全面,例如潜在猝死风险人群的死亡风险应包括引发猝死的疾病发生概率以及该疾病的死亡率两个方面。
其次,当前现有的死亡率预测手段多是基于患者进入ICU后一段时间内(一般为24小时或48小时)的监测数据进行一次死亡率预测,然而患者的生命体征是处于实时变化的,该方法并不能实时动态的反应患者的实际死亡风险。
再者,当前实时死亡率预测的方法在训练中将患者的任意时间的生存状况都以最终存活情况作为标签,但是在临床实际中患者的生理状态是动态实时变化的,最终的存活与否并不意味着患者之前的死亡风险都是固定不变的。
此外临床中由于采样周期等因素往往会出现临床数据缺失,现有的方法往往采用插值的方法对缺失数据做填充,再通过神经网络模型预测死亡率。但是在临床实际中患者的生理状态并不是平稳的,插值产生的数据不可避免的会和实际真实情况存在偏差,导致预测结果不准确。
发明内容
有鉴于此,本发明提供了一种在线实时患者危重程度评估及生命体征参数预测系统。
第一方面,本发明实施例提供了一种在线实时患者危重程度评估及生命体征参数预测系统,所述系统包括:
患者生命体征参数获取模块,用于实时采集患者生命体征参数;
患者生理状态特征获取模块,用于根据采集的患者生命体征参数设置第一患者生命体征参数,对第一患者生命体征参数特征进行特征提取得到患者生理状态特征;
采样实例获取模块,用于根据患者生理状态特征获取人体状态统计分布参数,配置累积测量权重,根据人体状态统计分布参数和累积测量权重获取当前的人体累积状态统计分布,对当前的人体累积状态统计分布进行采样,得到采样实例;
死亡风险预测模块,通过对各采样实例评估死亡率并进行分段统计得到患者死亡风险分布,患者死亡风险分布通过死亡率在每一区间的发生概率以及每一区间的期望死亡率表示;将死亡率在每一区间的发生概率与每一区间的期望死亡率的乘积之和作为当前患者的期望死亡率;设置高死亡概率阈值,将死亡率大于高死亡概率阈值的发生概率之和作为高死亡发生概率。
第二方面,本发明实施例提供了一种电子设备,包括存储器和处理器,所述存储器与所述处理器耦接;其中,所述存储器用于存储程序数据,所述处理器用于执行所述程序数据以实现上述的一种在线实时患者危重程度评估及生命体征参数预测系统。
第三方面,本发明实施例提供了一种计算机可读存储介质,其上存储有计算机程序,所述程序被处理器执行时实现上述的一种在线实时患者危重程度评估及生命体征参数预测系统。
与现有技术相比,本发明的有益效果为:
1、本发明通过实时采集患者生命体征参数,可以在线实时评估患者死亡风险。
2、本发明通过对各采样实例评估死亡率并进行分段统计得到患者死亡风险分布,根据患者死亡风险分布得到当前患者的期望死亡率以及高死亡发生概率。本发明考虑到在临床实际中患者的生理状态是动态实时变化的,最终的存活与否并不意味着患者之前的死亡风险都是固定不变的。通过患者死亡风险分布得到的期望死亡率及高死亡发生概率可以评估患者实时的死亡风险。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种患者危重程度评估及生命体征参数预测系统的示意图;
图2为本发明实施例提供的患者生理状态特征获取模块的示意图;
图3为本发明实施例提供的采样实例获取模块的示意图;
图4为本发明实施例提供的患者死亡风险分布的示意图;
图5为本发明实施例提供的一种电子设备的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
需要说明的是,在不冲突的情况下,下述的实施例及实施方式中的特征可以相互组合。
如图1所示,本发明实施例提供了一种在线实时患者危重程度评估及生命体征参数预测系统,所述系统包括:
患者生命体征参数获取模块,用于实时采集患者生命体征参数,并获取患者生命体征参数测量标识和患者生命体征参数统计向量。
所述患者生命体征参数包括但不限于心率HR、血氧饱和度SpO2、呼吸率RespR、体温Temp、收缩压Syst、舒张压Dias、平均压Mean等。
其中,如图2所示,获取患者生命体征参数测量标识和患者生命体征参数统计向量包括:获取前m时间段内各生命体征参数的测量标识,在本实例中m取1小时。在测量时间段内判断患者生命体征参数是否被测量,当患者生命体征参数没有被测量,则将患者生命体征参数对应的测量标识记为0;当患者生命体征参数被测量,则将患者生命体征参数对应的测量标识记为1。
记生命体征参数的测量标识M=[MHR,MSpO2,MRespR,MTemp,MSyst,MDias,MMean]。
式中,MHR为心率HR对应的测量标识,MSpO2为血氧饱和度SpO2对应的测量标识,MRespR为呼吸率RespR对应的测量标识,MTemp为体温Temp对应的测量标识,MSyst为收缩压Syst对应的测量标识,MDias为舒张压Dias对应的测量标识,MMean为平均压Mean对应的测量标识。
其中,所述统计向量包括平均值、标准差、最大值和最小值。
具体地,记患者生命体征参数统计向量V=[VHR,VSpO2,VRespR,VTemp,VSyst,VDias,VMean]。
式中,VHR为心率HR对应的统计向量,VSpO2为血氧饱和度SpO2对应的统计向量,VRespR为呼吸率RespR对应的统计向量,VTemp为体温Temp对应的统计向量,VSyst为收缩压Syst对应的统计向量,VDias为舒张压Dias对应的统计向量,VMean为平均压Mean对应的统计向量。
进一步地,当患者生命体征参数没有被测量,患者生命体征参数对应的测量标识为0时,该患者生命体征参数对应的统计向量记为默认值[-1,-1,-1,-1]。
进一步地,当患者生命体征参数被测量,获取患者生命体征参数对应的统计向量包括:当患者生命体征参数被测量,患者生命体征参数对应的测量标识为1时,先对患者生命体征参数进行截断归一化处理,再计算归一化后的患者生命体征参数对应的统计向量[平均值,标准差,最大值,最小值]。
进一步地,对患者生命体征参数进行截断归一化处理包括:
根据各患者生命体征参数的极限取值范围对各患者生命体征参数的测量值进行饱和截断,若患者生命体征参数的测量值大于该患者生命体征参数的最大极限值,则将该测量值设置为最大极限值;若患者生命体征参数的测量值小于该患者生命体征参数的最小极限值,则将该测量值设置为最小极限值。
需要说明的是,心率HR的极限取值范围可为0—300次/分钟,血氧饱和度SpO2的极限取值范围可为0—100%,呼吸率RespR的极限取值范围可为0—50次/分钟,体温Temp的极限取值范围可为0—50℃,收缩压Syst的极限取值范围可为0—200mmHg,舒张压Dias的极限取值范围为0—200 mmHg,平均压Mean的极限取值范围为0—200 mmHg。
然后再进行归一化,公式如下:
式中,为归一化后的患者生命体征参数,X为患者生命体征参数的测量值, LowerLimit为患者生命体征参数的最小极限值, UpperLimit为患者生命体征参数的最大极限值。
患者生理状态特征获取模块,用于根据采集的患者生命体征参数设置第一患者生命体征参数,对第一患者生命体征参数特征进行特征提取得到患者生理状态特征;
其中,根据采集的患者生命体征参数设置第一患者生命体征参数具体包括:对每一患者生命体征参数测量标识与对应的第二患者生命体征参数特征的乘积之和并进行归一化缩放得到第一患者生命体征参数特征,所述第二患者生命体征参数特征通过对患者生命体征参数统计向量进行特征提取得到。
具体地,将第二患者生命体征参数特征记为E=[EHR,ESpO2,ERespR,ETemp,ESyst,EDias,EMean]。
式中,EHR为心率HR对应的高维嵌入特征,ESpO2为血氧饱和度SpO2对应的高维嵌入特征,ERespR为呼吸率RespR对应的高维嵌入特征,ETemp为体温Temp对应的高维嵌入特征,ESyst为收缩压Syst对应的高维嵌入特征,EDias为舒张压Dias对应的高维嵌入特征,EMean为平均压Mean对应的高维嵌入特征。
在本实例中,通过特征提取网络(在本实例中,采用卷积神经网络)对患者生命体征参数统计向量V进行特征提取,得到第二患者生命体征参数特征(即患者生命体征参数高维嵌入特征)。表达式如下:
EHR=sigmoid(WHR×VHR+bHR)
ESpO2=sigmoid(WSpO2×VSpO2+bSpO2)
ERespR=sigmoid(WRespR×VRespR+bRespR)
ETemp=sigmoid(WTemp×VTemp+bTemp)
ESyst=sigmoid(WSyst×VSyst+bSyst)
EDias=sigmoid(WDias×VDias+bDias)
EMean=sigmoid(WMean×VMean+bMean)
式中,WHR、WSpO2、WRespR、WTemp、WSyst、WDias、WMean为连接权重,bHR、bSpO2、bRespR、bTemp、bSyst、bDias、bMean为偏置参数。
对每一患者生命体征参数测量标识与对应的第二患者生命体征参数特征的乘积之和并进行归一化缩放得到第一患者生命体征参数特征,表达式如下:
式中,G为归一化增益,归一化增益G为每一患者生命体征参数测量标识之和的倒数,表达式如下:
G=1/(MHR+MSpO2+MRespR+MTemp+MSyst+MDias)
需要说明的是,若每一患者生命体征参数测量标识均为0,则归一化增益G为0。
对第一患者生命体征参数特征进行特征提取得到患者生理状态特征S,表达式如下:
式中,WE为连接权重,bE为偏置参数。
采样实例获取模块,用于根据患者生理状态特征获取人体状态统计分布参数,配置累积测量权重,根据人体状态统计分布参数和累积测量权重获取当前的人体累积状态统计分布,对当前的人体累积状态统计分布进行采样,得到采样实例。
具体地,如图3所示,首先根据患者生理状态特征S获取人体状态统计分布参数,其中,所述人体状态统计分布参数采用高斯分布,通过高斯分布参数估计网络拟合得到μinstant和/>,该高斯分布参数估计网络可采用卷积神经网络结构,表达式如下:
μinstant=Wμ_instant×S+bμ_instant
σinstant=Wσ_instant×S+bσ_instant
式中,μinstant为人体状态统计分布的数学期望,为人体状态统计分布的方差,Wμ_instant、Wσ_instant为连接权重,bμ_instant、bσ_instant为偏置参数。
配置累积测量权重α,包括:当人体状态统计分布参数为第一次测量时,将累积测量权重α配置为1;否则,通过神经网络(在本实例中采用卷积神经网络)配置累积测量权重,表达式如下:
式中,Watten为连接权重,batten为偏置参数。
估计当前的人体累积状态统计分布参数,其中当前的人体累积状态统计分布可采用高斯分布,表达式如下:
其中,为前一次估计的人体累积状态统计分布参数。
对当前的人体累积状态统计分布采样,表达式如下:
ISi=u+σ×εi,i∈[0,K)
式中,ISi为采样实例,K为采样实例的个数,εi为符合标准正态分布的随机数。
死亡风险预测模块,通过对各采样实例评估死亡率并进行分段统计得到患者死亡风险分布,患者死亡风险分布通过死亡率在每一区间的发生概率以及每一区间的期望死亡率表示;将死亡率在每一区间的发生概率与每一区间的期望死亡率的乘积之和作为当前患者的期望死亡率;设置高死亡概率阈值,将死亡率大于高死亡概率阈值的发生概率之和作为高死亡发生概率。
通过对各采样实例评估死亡率并进行分段统计得到患者死亡风险分布包括:利用卷积神经网络对各采样实例ISi进行死亡率评估,表达式如下:
yi=sigmoid(WIS×ISi+bIS),i∈[0,K)
式中,WIS为连接权重,bIS为偏置参数。
对各采样实例死亡率进行分段区间统计,计算各死亡率区间段的实例个数,进而计算各死亡率段发生的概率即患者死亡风险分布。
示例性地,在本实例中,按10%为间隔进行统计,得到患者死亡风险分布,如图4所示,表达式为:
P = [f0:p0,f1:p1,f2:p2,f3:p3,f4:p4,f5:p5,f6:p6,f7:p7,f8:p8,f9:p9];
式中,p0表示死亡率在0%到10%之间的发生概率,f0表示该区段的期望死亡率为5%;p1表示死亡率在10%到20%之间的发生概率,f1表示该区段的期望死亡率为15%;p2表示死亡率在20%到30%之间的发生概率,f2表示该区段的期望死亡率为25%;p3表示死亡率在30%到40%之间的发生概率,f3表示该区段的期望死亡率为35%;p4表示死亡率在40%到50%之间的发生概率,f4表示该区段的期望死亡率为45%;p5表示死亡率在50%到60%之间的发生概率,f5表示该区段的期望死亡率为55%;p6表示死亡率在60%到70%之间的发生概率,f6表示该区段的期望死亡率为65%;p7表示死亡率在70%到80%之间的发生概率,f7表示该区段的期望死亡率为75%;p8表示死亡率在80%到90%之间的发生概率,f8表示该区段的期望死亡率为85%;p9表示死亡率在90%到100%之间的发生概率,f9表示该区段的期望死亡率为95%。
计算当前患者的期望死亡率,所述当前患者的期望死亡率为死亡率在每一区间的发生概率与每一区间的期望死亡率的乘积之和,表达式为:
式中,j为间隔区间的个数。
计算高死亡发生概率,包括:设置高死亡概率阈值,将死亡率大于高死亡概率阈值的发生概率之和作为高死亡发生概率。
进一步地,所述死亡风险预测模块还包括:
若当前时间段内存在未测量的患者生命体征参数,则利用卷积神经网络对各采样实例对应的未知的患者生命体征参数统计向量进行预测;表达式如下:
VHR_recons_i=sigmoid(WHR_recons×ISi+bHR_recons)
VSpO2_recons_i=sigmoid(WSpO2_recons×ISi+bSpO2_recons)
VRespR_recons_i=sigmoid(WRespR _recons×ISi+bRespR _recons)
VTemp_recons_i=sigmoid(WTemp_recons×ISi+bTemp_recons)
VSyst_recons_i=sigmoid(WSyst_recons×ISi+bSyst_recons)
VDias_recons_i=sigmoid(WDias_recons×ISi+bDias_recons)
VMean_recons_i=sigmoid(WMean_recons×ISi+bMean_recons)
式中,WHR_recons、WSpO2_recons、WRespR_recons、WTemp_recons、WSyst_reconst、WDias_recons、WMean_recons为连接权重,bHR_recons、bSpO2_recons、bRespR_recons、bTemp_recons、bSyst_recons、bDias_recons、bMean_recons为偏置参数。
将各实例的未知生命体征统计向量预测值进行平均作为当前时间段未知的患者生命体征统计向量期望值;表达式如下:
Vrecons=[VHR_recons,VSpO2_recons,VRespR_recons,VTemp_recons,VSyst_recons,VDias_recons,VMean_recons]。
根据患者生命体征参数的极限取值范围、未知的患者生命体征统计向量期望值得到当前时间段内未知的患者生命体征参数期望值。
式中,为未知的患者生命体征参数期望值,X为根据未知的患者生命体征统计向量期望值得到的未知的患者生命体征参数对应的初值,XLowerLimit为患者生命体征参数的最小极限值,XUpperLimit为患者生命体征参数的最大极限值。
所述死亡风险预测模块还包括:对高死亡率风险事件进行预警;具体地,包括:
设置期望死亡率阈值,当期望死亡率大于期望死亡率阈值时,触发第一预警信号。
设置高死亡发生概率阈值,当高死亡发生概率大于高死亡发生概率阈值时,触发第二预警信号。示例性地,设置死亡率大于50%的发生概率阈值为30%,计算死亡率大于50%的发生概率为,当死亡率大于50%的发生概率大于30%时,触发第二预警信号。
预测到的未测量的患者生命体征参数统计向量中的平均值超过患者生命体征参数的报警取值范围,触发第三预警信号。在本实例中,所述报警取值范围设置为:心率HR的预设范围可为50—120次/分钟,血氧饱和度SpO2的预设范围可为90—100%,呼吸率RespR的预设范围可为12—20次/分钟,体温Temp的预设范围可为36—37℃,收缩压Syst的预设范围可为90—140mmHg,舒张压Dias的预设范围为60—90 mmHg,平均压Mean的预设范围为70—110 mmHg。
需要说明的是,本发明能够进行更全面的死亡风险预警,包含期望死亡率、高死亡风险发生概率和未知生命体征参数异常。
训练模块,用于构建患者生命体征参数序列样本,设置损失函数,通过损失函数训练患者生理状态特征获取模块、采样实例获取模块和死亡风险预测模块;所述损失函数为死亡率预测损失、生命体征参数重构损失、实例约束损失的加权之和。
具体地,首先建立包含N个ICU患者生命体征参数序列样本,生命体征包括但不限于心率XHR=[XHR_1,XHR_2,…,XHR_T],血氧饱和度XSpO2=[XSpO2_1,XSpO2_2,…,XSpO2_T],呼吸率XRespR=[XRespR_1,XRespR_2,…,XRespR_T],体温XTemp=[XTemp_1,XTemp_2,…,XTemp_T],收缩压XSyst=[XSyst_1,XSyst_2,…,XSyst_T],舒张压XDias=[XDias_1,XDias_2,…,XDias_T],平均压XMean=[XMean_1,XMean_2,…,XMean_T]。
对任一患者,以m时间(1小时)为时间间隔将该患者的各生命体征序列分为T段,计 算各时间段内各生命体征参数的测量标识Mn_t及统计值向量Vn_t,K个采样实例ISn_t_i,各实 例ISn_t_i的死亡率yn_t_i,该时间段的患者的期望死亡率及生命体征参数统计向量期望 Vn_recons_t
其中,死亡率预测损失根据死亡风险预测模块得到的各采样实例的死亡率以及对应的死亡率标签得到;表达式如下:
式中,为该实例的死亡率标签,该死亡率标签定义为:
需要说明的是,本发明考虑到在临床实际中患者的生理状态是动态实时变化的,最终的存活与否并不意味着患者之前的死亡风险都是固定不变的。因此对各时间段内采用非固定的死亡与否标记。
生命体征参数重构损失根据患者生命体征参数获取模块得到的患者生命体征参数测量标识和患者生命体征参数统计向量、死亡风险预测模块得到的未知的患者生命体征统计向量期望值得到;表达式如下:
实例约束损失根据采样实例获取模块得到的采样实例以及对应的采样实例的期望分布得到;表达式如下:
式中,g(ISn_t_i)为采样实例ISn_t_i的期望分布。
所述采样实例的期望分布为:当采样实例处于采样分布约束范围内,采样实例的期望等于采样实例;当采样实例大于采样分布约束范围上限,采样实例的期望为第一系数与采样实例的乘积,且第一系数小于1;当采样实例小于采样分布约束范围下限,采样实例的期望为第二系数与采样实例的乘积,且第二系数大于1。
在本实例中,将采样实例ISn_t_i的期望分布g(ISn_t_i)定义为:
式中,Uplimit(t)为采样分布约束范围上限,Lowlimit(t)为采样分布约束范围下限,t表示表示距离死亡的时间。
其中,对于存活的患者样本,采样分布约束范围上限设置为第一预设值,采样分布约束范围下限设置为第二预设值,在本实例中采样分布约束范围上限设置为1,采样分布约束范围下限设置为0;对于死亡的患者样本,采样分布约束范围上限和采样分布约束范围下限根据距离死亡的时间进行设置。
在本实例中,采样分布约束范围上限和采样分布约束范围下限根据距离死亡的时间定义为:
所述损失函数为死亡率预测损失、生命体征参数重构损失、实例约束损失的加权之和,表达式如下:
L(x1,x2,…,xn│θ)=ω0×L0(x1,x2,…,xn│θ)+ω1×L1(x1,x2,…,xn│θ)+ω2×L2(x1,x2,…,xn│θ)
式中,ω0、ω1、ω2为权重系数,该权重系数可定义为1/3,也可通过网格搜索获取最优值。
与前述一种在线实时患者危重程度评估及生命体征参数预测系统的实施例相对应,本发明还提供了一种电子设备的实施例。
参见图5,本发明实施例提供的一种电子设备,包括一个或多个处理器,用于实现上述实施例中的一种在线实时患者危重程度评估及生命体征参数预测系统。
本发明电子设备的实施例可以应用在任意具备数据处理能力的设备上,该任意具备数据处理能力的设备可以为诸如计算机等设备或装置。装置实施例可以通过软件实现,也可以通过硬件或者软硬件结合的方式实现。以软件实现为例,作为一个逻辑意义上的装置,是通过其所在任意具备数据处理能力的设备的处理器将非易失性存储器中对应的计算机程序指令读取到内存中运行形成的。从硬件层面而言,如图5所示,为本发明电子设备所在任意具备数据处理能力的设备的一种硬件结构图,除了图5所示的处理器、内存、网络接口、以及非易失性存储器之外,实施例中装置所在的任意具备数据处理能力的设备通常根据该任意具备数据处理能力的设备的实际功能,还可以包括其他硬件,对此不再赘述。
上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程,在此不再赘述。
对于装置实施例而言,由于其基本对应于方法实施例,所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本发明方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
本发明实施例还提供一种计算机可读存储介质,其上存储有程序,该程序被处理器执行时,实现上述实施例中的一种在线实时患者危重程度评估及生命体征参数预测系统。
所述计算机可读存储介质可以是前述任一实施例所述的任意具备数据处理能力的设备的内部存储单元,例如硬盘或内存。所述计算机可读存储介质也可以是任意具备数据处理能力的设备,例如所述设备上配备的插接式硬盘、智能存储卡(Smart Media Card,SMC)、SD卡、闪存卡(Flash Card)等。进一步的,所述计算机可读存储介质还可以既包括任意具备数据处理能力的设备的内部存储单元也包括外部存储设备。所述计算机可读存储介质用于存储所述计算机程序以及所述任意具备数据处理能力的设备所需的其他程序和数据,还可以用于暂时地存储已经输出或者将要输出的数据。
本领域技术人员在考虑说明书及实践这里公开的内容后,将容易想到本申请的其它实施方案。本申请旨在涵盖本申请的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本申请的一般性原理并包括本申请未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的。
应当理解的是,本申请并不局限于上面已经描述并在附图中示出的精确结构,并且可以在不脱离其范围进行各种修改和改变。

Claims (10)

1.一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,所述系统包括:
患者生命体征参数获取模块,用于实时采集患者生命体征参数;
患者生理状态特征获取模块,用于根据采集的患者生命体征参数设置第一患者生命体征参数,对第一患者生命体征参数特征进行特征提取得到患者生理状态特征;
采样实例获取模块,用于根据患者生理状态特征获取人体状态统计分布参数,配置累积测量权重,根据人体状态统计分布参数和累积测量权重获取当前的人体累积状态统计分布,对当前的人体累积状态统计分布进行采样,得到采样实例;
死亡风险预测模块,通过对各采样实例评估死亡率并进行分段统计得到患者死亡风险分布,患者死亡风险分布通过死亡率在每一区间的发生概率以及每一区间的期望死亡率表示;将死亡率在每一区间的发生概率与每一区间的期望死亡率的乘积之和作为当前患者的期望死亡率;设置高死亡概率阈值,将死亡率大于高死亡概率阈值的发生概率之和作为高死亡发生概率。
2.根据权利要求1所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,所述患者生命体征参数包括:心率、血氧饱和度、呼吸率、体温、收缩压、舒张压、平均压。
3.根据权利要求1所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,所述患者生命体征参数获取模块还包括:获取患者生命体征参数测量标识和患者生命体征参数统计向量;
所述患者生理状态特征获取模块还包括:对每一患者生命体征参数测量标识与对应的第二患者生命体征参数特征的乘积之和并进行归一化缩放得到第一患者生命体征参数特征;所述第二患者生命体征参数特征通过对患者生命体征参数统计向量进行特征提取得到。
4.根据权利要求3所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,获取患者生命体征参数测量标识和患者生命体征参数统计向量包括:
当患者生命体征参数没有被测量时,则将患者生命体征参数对应的测量标识记为0,统计向量设置为默认值;
当患者生命体征参数被测量时,则将患者生命体征参数对应的测量标识记为1;先对患者生命体征参数进行截断归一化处理,再计算截断归一化后的患者生命体征参数对应的统计向量,所述统计向量为[平均值, 标准差, 最大值, 最小值];
其中,对患者生命体征参数进行截断归一化处理包括:根据各患者生命体征参数的极限取值范围对各患者生命体征参数的测量值进行饱和截断,若患者生命体征参数的测量值大于患者生命体征参数的最大极限值,则将患者生命体征参数的测量值设置为最大极限值;若患者生命体征参数的测量值小于患者生命体征参数的最小极限值,则将患者生命体征参数的测量值设置为最小极限值。
5.根据权利要求1所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,根据患者生理状态特征获取人体状态统计分布参数,配置累积测量权重,根据人体状态统计分布参数和累积测量权重获取当前的人体累积状态统计分布,对当前的人体累积状态统计分布进行采样,得到采样实例包括:
通过高斯分布参数估计网络对生理状态特征进行拟合,得到人体状态统计分布参数;
当人体状态统计分布参数为第一次测量时,将累积测量权重配置为1;否则,通过神经网络配置累积测量权重;
根据人体状态统计分布参数、累积测量权重、前一次估计的人体累积状态统计分布获取当前的人体累积状态统计分布;
预设采样实例的个数,对当前的人体累积状态统计分布进行采样,得到采样实例。
6.根据权利要求1所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,所述死亡风险预测模块还包括:
若当前时间段内存在未测量的患者生命体征参数,则对各采样实例对应的未知的患者生命体征参数统计向量进行预测;
将各实例的未知生命体征统计向量预测值进行平均作为当前时间段未知的患者生命体征统计向量期望值;
根据患者生命体征参数的极限取值范围、未知的患者生命体征统计向量期望值得到当前时间段内未知的患者生命体征参数期望值。
7.根据权利要求1所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,所述系统还包括:
训练模块,用于设置损失函数,通过损失函数训练患者生理状态特征获取模块、采样实例获取模块和死亡风险预测模块;所述损失函数为死亡率预测损失、生命体征参数重构损失、实例约束损失的加权之和;
其中,死亡率预测损失根据死亡风险预测模块得到的各采样实例的死亡率以及对应的死亡率标签得到;
生命体征参数重构损失根据患者生命体征参数获取模块得到的患者生命体征参数测量标识和患者生命体征参数统计向量、死亡风险预测模块得到的未知的患者生命体征统计向量期望值得到;
实例约束损失根据采样实例获取模块得到的采样实例以及对应的采样实例的期望分布得到。
8.根据权利要求7所述的一种在线实时患者危重程度评估及生命体征参数预测系统,其特征在于,所述采样实例的期望分布为:
当采样实例处于采样分布约束范围内,采样实例的期望等于采样实例;
当采样实例大于采样分布约束范围上限,采样实例的期望为第一系数与采样实例的乘积,且第一系数小于1;
当采样实例小于采样分布约束范围下限,采样实例的期望为第二系数与采样实例的乘积,且第二系数大于1;
其中,对于存活的患者样本,采样分布约束范围上限设置为第一预设值,采样分布约束范围下限设置为第二预设值;对于死亡的患者样本,采样分布约束范围上限和采样分布约束范围下限根据距离死亡的时间进行设置。
9.一种电子设备,包括存储器和处理器,其特征在于,所述存储器与所述处理器耦接;其中,所述存储器用于存储程序数据,所述处理器用于执行所述程序数据以实现上述权利要求1-8任一项所述的一种在线实时患者危重程度评估及生命体征参数预测系统。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述程序被处理器执行时实现如权利要求1-8中任一所述的一种在线实时患者危重程度评估及生命体征参数预测系统。
CN202311238621.9A 2023-09-25 2023-09-25 一种在线实时患者危重程度评估及生命体征参数预测系统 Active CN116978570B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202311238621.9A CN116978570B (zh) 2023-09-25 2023-09-25 一种在线实时患者危重程度评估及生命体征参数预测系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202311238621.9A CN116978570B (zh) 2023-09-25 2023-09-25 一种在线实时患者危重程度评估及生命体征参数预测系统

Publications (2)

Publication Number Publication Date
CN116978570A true CN116978570A (zh) 2023-10-31
CN116978570B CN116978570B (zh) 2024-02-06

Family

ID=88477108

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202311238621.9A Active CN116978570B (zh) 2023-09-25 2023-09-25 一种在线实时患者危重程度评估及生命体征参数预测系统

Country Status (1)

Country Link
CN (1) CN116978570B (zh)

Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140371610A1 (en) * 2013-06-14 2014-12-18 Massachusetts Institute Of Technology Method and apparatus for beat-space frequency domain prediction of cardiovascular death after acute coronary event
CN107742534A (zh) * 2010-03-15 2018-02-27 新加坡保健服务集团有限公司 患者存活性预测系统
CN109195507A (zh) * 2016-07-08 2019-01-11 爱德华兹生命科学公司 低血压分析参数的预测加权
WO2019215742A1 (en) * 2018-05-09 2019-11-14 B.G. Negev Technologies And Applications Ltd., At Ben-Gurion University System and method for hypotensive episode prediction
CN110944577A (zh) * 2017-07-27 2020-03-31 长桑医疗(海南)有限公司 一种血氧饱和度的检测方法与系统
CN111009327A (zh) * 2019-12-19 2020-04-14 京东方科技集团股份有限公司 一种风险预测方法、装置及系统、介质
US20200161000A1 (en) * 2017-06-02 2020-05-21 University Of Florida Research Foundation, Incorporated Method and apparatus for prediction of complications after surgery
CN111329459A (zh) * 2020-03-06 2020-06-26 北京化工大学 一种改良的患者生命体征监测预警评分方法
CN111542892A (zh) * 2017-11-09 2020-08-14 罗伯特·博世有限公司 生命体征监视系统和方法
CN112466469A (zh) * 2020-12-08 2021-03-09 杭州脉兴医疗科技有限公司 一种重大危机及死亡风险预测方法
CN112967803A (zh) * 2021-01-29 2021-06-15 成都一尧科技有限公司 基于集成模型的急诊患者早期死亡率预测方法及系统
CN114420231A (zh) * 2022-01-14 2022-04-29 东南大学 一种可解释的急性肾损伤持续预警方法、系统、存储介质及电子设备
CN114758786A (zh) * 2022-04-14 2022-07-15 军事科学院系统工程研究院卫勤保障技术研究所 基于无创参数的创伤后失血性休克动态早期预警系统
KR102526577B1 (ko) * 2022-09-27 2023-04-27 주식회사 비투엔 환자 관리를 위한 인공지능 기반의 위험 예측 시스템 및 방법
CN116525117A (zh) * 2023-07-04 2023-08-01 之江实验室 一种面向数据分布漂移检测与自适应的临床风险预测系统

Patent Citations (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107742534A (zh) * 2010-03-15 2018-02-27 新加坡保健服务集团有限公司 患者存活性预测系统
US20140371610A1 (en) * 2013-06-14 2014-12-18 Massachusetts Institute Of Technology Method and apparatus for beat-space frequency domain prediction of cardiovascular death after acute coronary event
CN109195507A (zh) * 2016-07-08 2019-01-11 爱德华兹生命科学公司 低血压分析参数的预测加权
US20200161000A1 (en) * 2017-06-02 2020-05-21 University Of Florida Research Foundation, Incorporated Method and apparatus for prediction of complications after surgery
CN110944577A (zh) * 2017-07-27 2020-03-31 长桑医疗(海南)有限公司 一种血氧饱和度的检测方法与系统
CN111542892A (zh) * 2017-11-09 2020-08-14 罗伯特·博世有限公司 生命体征监视系统和方法
WO2019215742A1 (en) * 2018-05-09 2019-11-14 B.G. Negev Technologies And Applications Ltd., At Ben-Gurion University System and method for hypotensive episode prediction
CN111009327A (zh) * 2019-12-19 2020-04-14 京东方科技集团股份有限公司 一种风险预测方法、装置及系统、介质
CN111329459A (zh) * 2020-03-06 2020-06-26 北京化工大学 一种改良的患者生命体征监测预警评分方法
CN112466469A (zh) * 2020-12-08 2021-03-09 杭州脉兴医疗科技有限公司 一种重大危机及死亡风险预测方法
CN112967803A (zh) * 2021-01-29 2021-06-15 成都一尧科技有限公司 基于集成模型的急诊患者早期死亡率预测方法及系统
CN114420231A (zh) * 2022-01-14 2022-04-29 东南大学 一种可解释的急性肾损伤持续预警方法、系统、存储介质及电子设备
CN114758786A (zh) * 2022-04-14 2022-07-15 军事科学院系统工程研究院卫勤保障技术研究所 基于无创参数的创伤后失血性休克动态早期预警系统
KR102526577B1 (ko) * 2022-09-27 2023-04-27 주식회사 비투엔 환자 관리를 위한 인공지능 기반의 위험 예측 시스템 및 방법
CN116525117A (zh) * 2023-07-04 2023-08-01 之江实验室 一种面向数据分布漂移检测与自适应的临床风险预测系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RUEY-KAI SHEU: "Multi-Modal Data Analysis for Pneumonia Status Prediction Using Deep Learning (MDA-PSP)", 《DIAGNOSTICS》 *
朱志勇;陈一昕;李建功;: "大数据技术在医疗急重症领域的应用", 邮电设计技术, no. 08, pages 28 - 31 *

Also Published As

Publication number Publication date
CN116978570B (zh) 2024-02-06

Similar Documents

Publication Publication Date Title
US10722179B2 (en) Residual-based monitoring of human health
Clifton et al. A large-scale clinical validation of an integrated monitoring system in the emergency department
US11580432B2 (en) System monitor and method of system monitoring to predict a future state of a system
CN111657888A (zh) 重度急性呼吸窘迫综合症预警方法及系统
US20220165417A1 (en) Population-level gaussian processes for clinical time series forecasting
WO2016057806A1 (en) Weaning readiness indicator, sleeping status recording device, and air providing system applying nonlinear time-frequency analysis
CN116978570B (zh) 一种在线实时患者危重程度评估及生命体征参数预测系统
CN116942149B (zh) 基于毫米波雷达的腰椎监控方法、装置、设备及存储介质
AU2021236462A1 (en) Residual-based monitoring of human health
CN112998717B (zh) 一种量化人体生理状态的分析方法和设备
CN113223713A (zh) 基于风险评估的动态监护方法、系统以及存储介质
CN111466877A (zh) 一种基于lstm网络的氧减状态预测方法
Dey et al. Patient Health Observation and Analysis with Machine Learning and IoT Based in Realtime Environment
US20140031658A1 (en) System and method for monitoring glucose or other compositions in an individual
CN112971791B (zh) 一种个体化的生理状态监测分析方法和设备
CN112971792B (zh) 一种基于连续生理数据的个体化状态监测分析方法和设备
CN112957017B (zh) 一种昏迷患者的实时体征监测方法及系统
Al Alwani et al. Neonatal events recognition using LBP descriptor and wavelet thresholding technique
Islam S et al. An efficient compression technique for Foetal phonocardiogram signals in remote healthcare monitoring systems
CN112971742A (zh) 一种个体化的生理状态监测分析方法和设备
CN114271804A (zh) 基于消费级设备的心率状态辅助检测系统、心脏状态监控系统、方法、存储介质及终端
CN112992379A (zh) 一种个体化的生理状态监测分析方法和设备
CN116246793A (zh) 基于可穿戴设备的用户健康画像生成方法及装置
CN116965785A (zh) 一种生命体征监测分析方法、系统、终端及存储介质
CN116823752A (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