CN113171059B - 多模态监测信息的术后end风险预警及相关设备 - Google Patents
多模态监测信息的术后end风险预警及相关设备 Download PDFInfo
- Publication number
- CN113171059B CN113171059B CN202110447993.7A CN202110447993A CN113171059B CN 113171059 B CN113171059 B CN 113171059B CN 202110447993 A CN202110447993 A CN 202110447993A CN 113171059 B CN113171059 B CN 113171059B
- Authority
- CN
- China
- Prior art keywords
- data
- time period
- clinical data
- early warning
- clinical
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 73
- 230000002980 postoperative effect Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 claims abstract description 51
- 238000012545 processing Methods 0.000 claims description 34
- 238000012549 training Methods 0.000 claims description 23
- 238000003070 Statistical process control Methods 0.000 claims description 21
- 238000011282 treatment Methods 0.000 claims description 14
- 238000012795 verification Methods 0.000 claims description 14
- 238000000605 extraction Methods 0.000 claims description 12
- 230000002537 thrombolytic effect Effects 0.000 claims description 9
- 238000005457 optimization Methods 0.000 claims description 8
- 238000001356 surgical procedure Methods 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 7
- 206010061216 Infarction Diseases 0.000 claims description 6
- 239000003146 anticoagulant agent Substances 0.000 claims description 6
- 230000007574 infarction Effects 0.000 claims description 6
- 230000036772 blood pressure Effects 0.000 claims description 5
- 238000004590 computer program Methods 0.000 claims description 5
- 229940079593 drug Drugs 0.000 claims description 5
- 239000003814 drug Substances 0.000 claims description 5
- 229910003798 SPO2 Inorganic materials 0.000 claims description 4
- 101100478210 Schizosaccharomyces pombe (strain 972 / ATCC 24843) spo2 gene Proteins 0.000 claims description 4
- 210000004369 blood Anatomy 0.000 claims description 4
- 239000008280 blood Substances 0.000 claims description 4
- 238000001990 intravenous administration Methods 0.000 claims description 4
- 238000010845 search algorithm Methods 0.000 claims description 4
- 206010002091 Anaesthesia Diseases 0.000 claims description 3
- 206010003658 Atrial Fibrillation Diseases 0.000 claims description 3
- 206010069729 Collateral circulation Diseases 0.000 claims description 3
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 claims description 3
- 208000032843 Hemorrhage Diseases 0.000 claims description 3
- 206010020772 Hypertension Diseases 0.000 claims description 3
- 208000034800 Leukoencephalopathies Diseases 0.000 claims description 3
- 208000031481 Pathologic Constriction Diseases 0.000 claims description 3
- 230000037005 anaesthesia Effects 0.000 claims description 3
- 230000000702 anti-platelet effect Effects 0.000 claims description 3
- 230000010100 anticoagulation Effects 0.000 claims description 3
- 238000004820 blood count Methods 0.000 claims description 3
- 208000029078 coronary artery disease Diseases 0.000 claims description 3
- 239000008103 glucose Substances 0.000 claims description 3
- 230000006872 improvement Effects 0.000 claims description 3
- 210000000265 leukocyte Anatomy 0.000 claims description 3
- 210000000440 neutrophil Anatomy 0.000 claims description 3
- 230000036262 stenosis Effects 0.000 claims description 3
- 208000037804 stenosis Diseases 0.000 claims description 3
- 238000010200 validation analysis Methods 0.000 claims description 3
- 238000005406 washing Methods 0.000 claims description 3
- 101000637031 Homo sapiens Trafficking protein particle complex subunit 9 Proteins 0.000 claims 1
- 108010089417 Sex Hormone-Binding Globulin Proteins 0.000 claims 1
- 208000007536 Thrombosis Diseases 0.000 claims 1
- 102100031926 Trafficking protein particle complex subunit 9 Human genes 0.000 claims 1
- 238000004886 process control Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 18
- 239000000523 sample Substances 0.000 description 17
- 230000008569 process Effects 0.000 description 9
- 238000013528 artificial neural network Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 230000006866 deterioration Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000029058 respiratory gaseous exchange Effects 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 4
- 230000004913 activation Effects 0.000 description 4
- 230000004872 arterial blood pressure Effects 0.000 description 4
- 238000004140 cleaning Methods 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 210000005036 nerve Anatomy 0.000 description 4
- 238000002560 therapeutic procedure Methods 0.000 description 4
- 230000005856 abnormality Effects 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000013135 deep learning Methods 0.000 description 3
- 238000010837 poor prognosis Methods 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 208000028399 Critical Illness Diseases 0.000 description 2
- 238000012276 Endovascular treatment Methods 0.000 description 2
- 208000032382 Ischaemic stroke Diseases 0.000 description 2
- 208000006011 Stroke Diseases 0.000 description 2
- 229940127219 anticoagulant drug Drugs 0.000 description 2
- 238000013176 antiplatelet therapy Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000003727 cerebral blood flow Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 230000000926 neurological effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000000250 revascularization Effects 0.000 description 2
- 238000012552 review Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 206010008111 Cerebral haemorrhage Diseases 0.000 description 1
- 208000002403 Encephalocele Diseases 0.000 description 1
- 206010038563 Reocclusion Diseases 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 210000000709 aorta Anatomy 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013481 data capture Methods 0.000 description 1
- 230000006837 decompression Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001787 epileptiform Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005713 exacerbation Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000036403 neuro physiology Effects 0.000 description 1
- 230000007658 neurological function Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 230000002087 whitening effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7271—Specific aspects of physiological measurement analysis
- A61B5/7275—Determining trends in physiological measurement data; Predicting development of a medical condition based on physiological measurements, e.g. determining a risk factor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Public Health (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明涉及多模态监测信息的术后END风险预警及相关设备。设备包括:存储器和处理器;所述处理器用于执行以下操作:获取待测用户在当前时间周期的临床数据;将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;利用统计过程控制方法获取所述临床数据的统计值,将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确实现下一时间周期的术后END风险预警。本申请实现了实时更新术后END预警。
Description
技术领域
本发明涉及生物信息技术领域,更具体地,涉及多模态监测信息的急性缺血性卒中血管内治疗术后早期神经功能恶化自动化预警、电子设备和计算机可读存储介质。
背景技术
尽管目前血管内治疗已获得全球多项国际组织及国际指南的最高级治疗推荐。但并非使所有患者均能从中获益,约40%接受血管内治疗患者术后最终仍发展为不良预后。很多患者在术后出现非预期病情恶化,最终表现为不良预后。早期神经功能恶化(EarlyNeurological Deterioration,END)是血管内治疗术后最常见,最有代表性的不良事件。目前最广泛使用的END定义为血管再通后24h小时内出现的神经功能恶化,且NIHSS评分较基线增加大于等于4分。END症状性脑出血、再闭塞、无效再通、脑疝等是END发生的常见原因,但具体发生机制仍不清楚。及时预测END的发生,提前采取积极的对应干预措施,能够逆转患者结局。早期神经功能恶化是指患者接受血管再通治疗(包括静脉溶栓与血管内治疗)后短期内发生的神经功能好转后恶化,最终表现为预后不良,甚至死亡。
临床发现一些血管内治疗术后END患者不仅神经功能监测可出现明显异常,也可出现特征性的基础体征信号变化,和SND(神经功能的突发恶化)发生密切相关。这些复合生理信号的改变有时提示END的出现时间。但需要注意的是这些异常生理监测信号的敏感性虽然显著,但预测END的特异性不高,很多重症患者会因其它***疾病问题出现类似改变。基于现有技术目前END预测方面研究结果显示单纯基于术前术中情况、临床影像信息,或单纯依靠电生理监测信息都不能很好的综合预判。基于静态的临床信息、影像信息、同时整合术后连续多模监护信息可以更好指导对于血管内治疗术后END进行实时预测。然而在实际临床这些复杂的信息、多模式监护数据等多形式数据虽被实时监测并有限记录,但其中仅很小部分被筛选用来临床分析。同时这些多维、复杂的数据往往需高年资医师凭借多年临床经验才可部分分析,且时常会出现误判的情况。此外重症患者个体化趋势很明显,基于普通病患的评价标准有时不适用,我们很难使用统一标准来判断这些海量数据,更加剧了分析的难度。
发明内容
鉴于上述问题,本发明利用深基于卷积理论的深度自编码器提取复杂多维信息与电生理信号的时空表征,并利用统计过程控制方法(Statistical Process Control,SPC)对提取的电生理特征进行监测,以达到实时更新END预警。
一种基于多模态监测信息的术后END风险预警设备,所述设备包括:存储器和处理器;
所述存储器用于存储程序指令;
所述处理器用于调用程序指令,当程序指令被执行时,用于执行以下操作:
获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;
利用统计过程控制方法获取所述临床数据的统计值;
将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警;
所述基于卷积理论的深度自编码器的训练方式和各个时间周期的预警阈值的确定方式包括:
获取术后无END患者在各个时间周期的临床数据作为训练样本,获取术后有END患者在各个时间周期的临床数据作为验证样本;
利用训练样本训练基于卷积理论的深度自编码器;利用验证样本验证所述采用基于卷积理论的深度自编码器,并对每个时间周期的验证样本的输入数据和输出数据均利用统计过程控制方法进行处理,得到每个时间周期对应的预警阈值。
所述基线数据包括基线信息、既往史信息、基线实验室指标、基线影像信息、术中情况等方面的特征变量;
进一步,基线信息(B)包括年龄,性别、基线NIHSS评分,基线血压、发病时间;既往史信息(H)包括高血压、冠心病、房颤、既往卒中,既往用药史(抗血小板治疗,抗凝治疗);基线实验室指标(L)包括入院血糖、白细胞计数、中性粒计数、PLT计数、LDL);基线影像信息(I)包括梗死体积、Mismatch体积、ASPECT评分,梗死部位、责任血管、闭塞程度、早期征象、侧枝循环、脑白质病变程度、微出血情况等;术中情况(E)包括是否静脉溶栓治疗,麻醉方式、血管内治疗方式(支架取栓、抽吸取栓、支架成形等),取栓次数,术中抗凝、抗血小板应用,术后TICI分级、残余狭窄,术后NIHSS评分、手术时间等。
进一步,动态监测数据包括ECG(e)、RESP(r)、NIBP(nb)、ABP(ab)、SPO2(s)、PULSE(p)、EEG(eeg)、TCD(tcd)等。
进一步,对临床数据进行特征提取前先进行数据清洗,数据清洗包括对缺失数据处理和对离散及噪音数据的处理。
所述数据清洗包括对缺失数据处理和对离散及噪音数据的处理。所述对缺失数据处理可基于临床分析及模型参数优化完成;所述对离散及噪音数据的处理可采用Binning方法处理;
优选的,对临床数据进行数据清洗后,对所有数据进行归一化处理。
优选的,基于卷积理论的深度自编码器对术后无END患者组和术后有END患者组高维度的基线数据、高数据量的动态监测数据进行特征提取;
优选的,通过基于卷积理论的深度自编码器对高维度的基线数据、高数据量的动态监测数据进行特征提取,进而获得可整合的修正信息编码。
优选的,通过基于卷积理论的深度自编码器对临床数据进行特征提取,进而获得可整合的修正信息编码。
深度自动编码器(Deep Auto-encoder,DAE)主要用于完成转换学习任务,可完成
无监督学习及非线性特征提取。其基本思想就是直接使用一层或者多层的神经网络对输入
数据进行映射,得到输出向量,作为从输入数据提取出的特征(见图)。DAE是一种利用无监
督预训练和***性参数优化的多层非线性网络从无类标数据中提取高维复杂输入数据的
分层特征,并得到原始数据的分布式特征表示的深度学习神经网络结构。DAE由编码器、解
码器和隐含层组成。
常见DAE不能有效解决繁杂数据中的池化与白化问题,并且大量冗余参数被强迫参与计算使得运算效率较低,而基于卷积理论的DAE是可用于处理本申请中多模态数据的,该结构利用重要的局部特征重构原始数据,且输入数据的所有局部特征共享权值矩阵,因而该类DAE的隐含层能完整保存受局部空间限制的边缘特征。基于卷积理论的DAE用基本模块的线性组合来拟合多模态数据信号,显著提升了信号综合识别的速度和准确率。
优选的,采用启发式搜索算法优化自动编码器结构。
优选的,采用启发式搜索算法,选取自动预警模型的结构和超参数,对自动预警模 型进行优化。
优选的,使用统计过程控制中的指数加权移动平均控制图(ExponentiallyWeighted Moving-A verage,EWMA)进行处理,得到每个时间周期对应的预警阈值,侦测多模态监测信号的微小偏移,监测指标变异时自动预警;
进一步,深度自编码器模型输出的值越大,提示具有更大的概率发生END。
优选的,对术后有END患者组临床数据进行特征提取,改进所述深度自编码器模型,获取监测异常参数,提示术后END的发生;
优选的,定义参数RE为监测异常参数,通过训练后的深度编码器模型计算REs值(i=1~N),该值越大时提示此时具有更大概率发生END。
优选的,计算统计过程控制的指数加权移动平均控制图(EWMA)统计量,得到样本均值的加权平均值。
优选的,根据验证样本数据,计算统计过程控制的指数加权移动平均控制图统计量,得到预警阈值。
优选的,获取术后有END患者和术后无END患者在各个时间周期的临床数据作为验证样本改进预警阈值。通过正常和异常样本数据进一步优化预警阈值。
优选的,每个时间周期对应的预警阈值是不尽相同的。
EMWA通过给与观察值不同的权数,按照不同权数求得移动平均值,确定预测值的
方法。采用EMWA更符合人体生物学信号的改变特点,该方法既可以结合患者基线情况及基
础监测指标,又可以满足观察监测指标的近期观察值对于预测值有更大影响,更能反映近
期变化的趋势。运用该技术解决了危重症患者多模监测中突变信号的智能理解问题。
一种基于多模态监测信息的术后END风险预警方法,包括:
获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;
利用统计过程控制方法获取所述临床数据的统计值;
将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
一种基于多模态监测信息的术后END风险预警装置,包括:
获取单元,用于获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
处理单元,用于将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据,利用统计过程控制方法获取所述临床数据的统计值;
预测单元,用于将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现上述的术后END风险预警方法。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获取其他的附图。
图1是本发明实施例提供的一种基于多模态监测信息的术后END风险预警方法的示意流程图;
图2是本发明实施例提供的一种基于多模态监测信息的术后END风险预警装置的示意框图;
图3是本发明实施例提供的深度自编码器内部结构图;
图4是本发明实施例提供的模型训练与改进原理图;
图5是本发明实施例提供的预警模型构建及与预警示意图;
图6是床旁多模监测示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
在本发明的说明书和权利要求书及上述附图中的描述的一些流程中,包含了按照特定顺序出现的多个操作,但是应该清楚了解,这些操作可以不按照其在本文中出现的顺序来执行或并行执行,操作的序号如101、102等,仅仅是用于区分开各个不同的操作,序号本身不代表任何的执行顺序。另外,这些流程可以包括更多或更少的操作,并且这些操作可以按顺序执行或并行执行。需要说明的是,本文中的“第一”、“第二”等描述,是用于区分不同的消息、设备、模块等,不代表先后顺序,也不限定“第一”和“第二”是不同的类型。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域技术人员在没有作出创造性劳动前提下所获取的所有其他实施例,都属于本发明保护的范围。
图1是本发明实施例提供的一种基于多模态监测信息的术后END风险预警方法的示意流程图,具体地,所述方法包括如下步骤:
101:获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
本发明实施例中,所述基线数据包括基线信息、既往史信息、基线实验室指标、基线影像信息、术中情况等方面的数据。
在一个实施例中,患者入选标准:年龄≥18岁;急性缺血性卒中患者;发病至到院时间≤24小时;MRA或CTA或DSA提示大动脉急性闭塞;接受血管内治疗。
在一个实施例中,基线信息、既往史信息、基线实验室指标、基线影像信息、术中情况等临床信息包括基本信息、人口学特征、院前急救(静脉溶栓、血管内治疗)、既往史、家族史、既往用药、入院查体、入院诊断、基线NIHSS评分、住院期间辅助检查、住院期间治疗、最终诊断、出院带药、出院NIHSS评分等、发生END患者的发生时间、程度、原因以及随访信息;影像学信息包括发病24h内完成影像数据。
在一个实施例中,影像学信息包括1)术前影像:基线CT、多模CT/MRI等;2)术中影像:DSA;3)术后24小时复查影像:术后即刻CT、术后24小时复查MRI或CT。
在一个实施例中,所述基线信息(B)包括年龄,性别、基线NIHSS评分,基线血压、发病时间;既往史信息(H)包括高血压、冠心病、房颤、既往卒中,既往用药史(抗血小板治疗,抗凝治疗);基线实验室指标(L)包括入院血糖、白细胞计数、中性粒计数、PLT计数、LDL);基线影像信息(I)包括梗死体积、Mismatch体积、ASPECT评分,梗死部位、责任血管、闭塞程度、早期征象、侧枝循环、脑白质病变程度、微出血情况等;术中情况(E)包括是否静脉溶栓治疗,麻醉方式、血管内治疗方式(支架取栓、抽吸取栓、支架成形等),取栓次数,术中抗凝、抗血小板应用,术后TICI分级、残余狭窄,术后NIHSS评分、手术时间等。
在一个实施例中,基本监护参数及神经电生理监测信号包括入组患者术后24小时的连续心电、呼吸、血压、脉氧等基本生命体征监护信息;术后24h脑电图、脑血流监测信息。某些实施例中,采用图6所示床旁多模监测基本监护参数及神经电生理监测信号。
在一个实施例中,基本监护参数及神经电生理监测信号包括ECG(e)、RESP(r)、NIBP(nb)、ABP(ab)、SPO2(s)、PULSE(p)、EEG(eeg)、TCD(tcd)等。
在一个实施例中,临床信息通过基于标准化设计的电子数据捕获***(EDC)实现在线或离线录入数据库。影像学资料统一保存成DICOM格式,影像学结果判读采取中心化盲法判读的方法。重症监护记录单、操作记录、出院小结、病历首页或门急诊及其他临床信息:上传照片形式。其他检查:脑血流、脑电图需保存为EDF格式或ASCII格式,所有纳入训练集中的信息及数据格式须经后台指控合格。
在一个实施例中,基本监护参数包括:ECG(心电)、RESP(呼吸)、NIBP(无创血压)或ABP(动脉血压)、SPO2(血氧)、PULSE(脉搏)。上述监测指标的参数要求(包括波形、增益、滤波等)按照标准ICU监护模式要求进行,信号不准确情况可进行床边矫正,所有参数的原始数据用于深度学习。
在一个实施例中,采用可移动式视频脑电图对于患者进行床旁持续EEG监测,监测时间要求24h以上。按照国际10-20***安装16导电极(颅骨减压患者例外)。最终脑电图结果分析由两位以上专职判读人员进行并标定异常事件(如癫痫样放电、NCS等),具体脑电图判读应采用美国临床神经生理学会重症监护脑电图标准术语进行。
在一个实施例中,床旁经颅多普勒超声(Transcranial Doppler sonography,TCD)检查采用经颅多普勒超声仪配备2.0MHz脉冲波多普勒超声探头。根据TCD技术规范,制定适宜的输出功率,选取取样容积10-15mm。根据实际情况调整增益、标尺、基线等。TCD检查需至少包括MCA、ICA、V A、BA探测。除非窗位不佳,应进行各窗位探查。所有技术操作应由专职TCD技术员进行,需按预先设计的CRF表单内容进行检测,最终由具有判读资质技术人员进行结果解释。
在一个实施例中,对所述临床数据进行数据清洗;
在一个实施例中,数据清洗主要从两个方面进行,一方面对于缺失数据处理,可基于临床分析及模型参数优化完成;另一方面对于离散及“噪音”采用Binning方法处理。
在一个实施例中,数据清洗步骤,第一步对于缺失数据处理,可基于临床分析及模型参数优化完成,第二步对于离散及“噪音”采用Binning方法处理,第三步对于所有数据进行归一化处理。
102:将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;
利用统计过程控制方法获取所述临床数据的统计值;
在一个实施例中,临床数据中含有高维度的临床及影像信息和高数据量的监测信号,通过基于卷积理论的深度自动编码器对术后无END患者和术后有END患者临床进行特征提取。
在一个实施例中,为了优化基于卷积理论的深度自动编码器结构,利用启发式搜索算法,使获取的编码最优。
解码和编码具体流程见如下公式(2-5)。
x=[B1,B2,B3..,H,H2,H3…,L1,L2,L3…,I,I2,I3…,E1,E2,E3…,e1,e2,e3…,r1,r2,r3…,nb,ab1,ab2,ab3…,s1,s2,s3…,p1,p2,p3…,eeg1,eeg2,eeg3…,tcd1,tcd2,tcd3…]T (1)
y=f(x) (2)
g(.)代表y的解码过程,采用(4)sigmoid函数作为神经网络激活函数。
公式(5)描述了深度自编码器(DAE)训练过程中获取最优参数的方法。其中Wl和bl为l层的权值和偏差,L是DAE的层数,N为数据集的样本量。上述深度编码器的内部结构图见图3。
深度自动编码器(Deep Auto-encoder,DAE)是一种利用无监督预训练和***性参数优化的多层非线性网络从无类标数据中提取高维复杂输入数据的分层特征,并得到原始数据的分布式特征表示的深度学习神经网络结构。DAE由编码器、解码器和隐含层组成。编码器是输入y到隐含表示h的映射,表示为:
h=f(x)=St(W+bb) (6)
其中,St是非线性激活函数,一般为逻辑函数,其表达式为:
解码器函数g(h)将隐含层数据映射回重构y,表示为:
y=g(h)=Sg(W′h+by) (8)
其中,Sg是解码器的激活函数,一般为线性函数或sigmoid函数。训练DAE的过程是在训练样本集D上寻找参数θ={W,by,bh}的最小化重构误差,重构误差的表达式为:
JAE=∑X∈DL(x,g(f(x))) (9)
其中,L为重构误差函数,一般可用平方误差函数或交叉熵损失函数,二者分别表示为:
L(x,y)=||x-y||2 (10)
其中,平方误差用于线性Sg,交叉熵损失函数用于sigmoid。
基于卷积理论的DAE是可用于处理多模态数据的神经网络,该结构利用重要的局部特征重构原始数据,且输入数据的所有局部特征共享权值矩阵,因而该类DAE的隐含层能完整保存受局部空间限制的边缘特征。
对于单通道的输入x其k阶特征映射的隐含表示是:
hk=σ(x*Wk+bk) (12)
其中,σ是激活函数,一般采用sigmoid函数,“*”表示2D卷积运算,Wk表示权值矩阵,bk表示偏置向量,卷积运算表达式为:
该类DAE的重构函数为:
其中,c为每个数据通道的偏置,H是隐含特征映射集,是权值矩阵的批处理,权值更新的规则为随机梯度下降。值得注意的是,在该类DAE中,一个隐含映射对应一个偏置值,偏置向量b对整个映射均有效,并且每个映射负责捕捉数据的一个特征,便于对神经网络进行预训练和精雕,有效缩短了特征提取的时间,简化了特征提取的过程,实现了数据特征的分层提取。基于卷积理论的DAE用基本模块的线性组合来拟合多模态数据信号,显著提升了信号综合识别的速度和准确率。目前,该类DAE可完成目标辨识、动态跟随和模拟等任务,有效解决了原有DAE处理多模态监测数据时出现的识别速度慢,准确率低,需要大量类标数据等问题。
在一个实施例中,选取其中正常组患者(术后未发生END)临床数据为模型训练集,将训练集临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;利用统计过程控制方法获取每个时间周期的样本的输入数据和输出数据进行处理,构建初步预警模型。
在一个实施例中,使用统计过程控制的指数加权移动平均控制图,利用验证集患者(术后有END患者和术后无END患者)临床数据对构建的初步预警模型进行验证及改进,获得预警模型。我们使用术后有END患者相关临床数据并定义参数RE为监测异常的参数,用其提示END发生。通过训练后的深度编码器模型计算REs值(i=1~N),该值越大时提示此时具有更大概率发生END,同时我们使用统计过程控制的指数加权移动平均控制图(Exponentially Weighted Moving-Average,EWMA)估计发生END的控制极限,并得到预警阈值,当模型达到极限时进行报警(END预警),提示未来某时间END将会发生。本节段特征提取与模型训练具体流程见图4。
EMWA通过给与观察值不同的权数,按照不同权数求得移动平均值,确定预测值的方法。采用EMWA更符合人体生物学信号的改变特点,该方法既可以结合患者基线数据及动态监测数据,又可以满足观察监测数据的近期观察值对于预测值有更大影响,更能反映近期变化的趋势。指数加权移动平均控制图定义公式为:
Zi=λxi+(1-λ)Zi-1 (15)
其中常数λ的取值范围为0<λ≤1,Zi是EWMA统计量,即所有之前样本均值的加权平均值。
该公式的初始值Z0(当i=1时)取流程的目标值(即在Z0=μ0),有时也用初始数据的均值作为初始值,即Z0=XBar。由于EWMA是所有之前和当前样本的加权平均,因此对数据的正态性假设很不敏感,因此用于单个观测值是十分理想的。
如果观测值xi是独立的随机变量,方差为σ2,那么Zi的方差为
因此,EWMA控制图的纵轴是Zi,横轴是样本序号或时间,中心线和控制限的计算公式为:
Center line=μ0
注意到上述公式中的(1-λ)2i部分,当i逐渐增大时,(1-λ)2i将很快收敛到0,因此当i增大时,UCL和LCL将稳定到下面两个值,此为EWMA的控制限收敛:
在一个实施例中,每个时间周期可以是1分钟-3小时内任意时间段,示例性的,每个时间周期可以是1分钟、2分钟、5分钟、10分钟、15分钟、20分钟、25分钟、30分钟、1小时、1.5小时、2小时、2.5小时、3小时等。
图5是本发明实施例提供的预警模型构建及与预警示意图。
在一个实施例中,基于卷积理论的深度自编码器的训练方式和各个时间周期的预警阈值的确定方式包括:
获取术后无END患者在各个时间周期的临床数据作为训练样本,获取术后有END患者在各个时间周期的临床数据作为验证样本;
利用训练样本训练基于卷积理论的深度自编码器;利用验证样本验证所述采用基于卷积理论的深度自编码器,并对每个时间周期的验证样本的输入数据和输出数据均利用统计过程控制方法进行处理,得到每个时间周期对应的预警阈值,构建预警模型。
取待测用户在当前时间周期的临床数据,将所述临床数据输入构建预警模型,所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
103:将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
在一个实施例中,术后24小时通过内网服务器连续自动获取动态监测信号,输入预警模型,实现连续实时END预警。
图2是本发明实施例提供的一种基于多模态监测信息的术后END风险预警装置的示意框图。
一种基于多模态监测信息的术后END风险预警装置,包括:
获取单元,用于获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
处理单元,用于将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据,利用统计过程控制方法获取所述临床数据的统计值;
预测单元,用于将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
一种基于多模态监测信息的术后END风险预警方法,包括:
获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;
利用统计过程控制方法获取所述临床数据的统计值;
将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警
一种计算机可读存储介质,所述计算机可读取存储介质存储有计算机程序,所述计算机程序被处理器执行实现上述基于多模态监测信息的术后END风险预警方法。
本验证实施例的验证结果表明,为适应症分配固有权重相对于默认设置来说可以适度改善本方法的性能。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的***,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的***,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个***,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:只读存储器(ROM,Read Only Memory)、随机存取存储器(RAM,RandomAccess Memory)、磁盘或光盘等。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上对本发明所提供的一种计算机设备进行了详细介绍,对于本领域的一般技术人员,依据本发明实施例的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (12)
1.一种基于多模态监测信息的术后END风险预警设备,所述设备包括:存储器和处理器;
所述存储器用于存储程序指令;
所述处理器用于调用程序指令,当程序指令被执行时,用于执行以下操作:
获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;
利用统计过程控制方法获取所述临床数据的统计值;
将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警;
所述基于卷积理论的深度自编码器的训练方式和各个时间周期的预警阈值的确定方式包括:
获取术后无END患者在各个时间周期的临床数据作为训练样本,获取术后有END患者在各个时间周期的临床数据作为验证样本;
利用训练样本训练基于卷积理论的深度自编码器;利用验证样本验证所述基于卷积理论的深度自编码器,并对每个时间周期的验证样本的输入数据和输出数据均利用统计过程控制方法进行处理,得到每个时间周期对应的预警阈值。
2.根据权利要求1所述的设备,其特征在于,采用启发式搜索算法优化深度自动编码器结构。
3.根据权利要求1所述的设备,其特征在于,所述对每个时间周期的验证样本的输入数据和输出数据均利用统计过程控制方法进行处理是对每个时间周期的验证样本的输入数据和输出数据均利用统计过程控制的指数加权移动平均控制图进行处理,得到每个时间周期对应的预警阈值。
4.根据权利要求1所述的设备,其特征在于,所述基线数据包括基线信息、既往史信息、基线实验室指标、基线影像信息、术中情况方面的数据;动态监测数据包括术后基本监护参数及神经电生理监测信号数据。
5.根据权利要求4所述的设备,其特征在于,所述基线信息包括年龄,性别、基线 NIHSS评分,基线血压、发病时间;既往史信息包括高血压、冠心病、房颤、既往卒中,既往用药史;基线实验室指标包括入院血糖、白细胞计数、中性粒计数、PLT 计数、LDL;基线影像信息包括梗死体积、Mismatch 体积、ASPECT 评分,梗死部位、责任血管、闭塞程度、早期征象、侧枝循环、脑白质病变程度、微出血情况;术中情况包括是否静脉溶栓治疗,麻醉方式、血管内治疗方式,取栓次数,术中抗凝、抗血小板应用,术后 TICI 分级、残余狭窄,术后 NIHSS 评分、手术时间。
6.根据权利要求4所述的设备,其特征在于,所述动态监测数据包括ECG、RESP、NIBP、ABP、SPO2、PULSE、EEG、TCD方面的数据。
7.根据权利要求1所述的设备,其特征在于,对临床数据进行特征提取前先进行数据清洗,数据清洗包括对缺失数据处理和对离散及噪音数据的处理。
8.根据权利要求7所述的设备,其特征在于,所述对缺失数据处理基于临床分析及模型参数优化完成;所述对离散及噪音数据的处理采用 Binning 方法处理。
9.根据权利要求7所述的设备,其特征在于,对临床数据进行数据清洗后,对所有数据进行归一化处理。
10.根据权利要求1所述的设备,其特征在于,获取术后有END患者和术后无END患者在各个时间周期的临床数据作为验证样本改进预警阈值。
11.一种基于多模态监测信息的术后END风险预警装置,包括:
获取单元,用于获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
处理单元,用于将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据,利用统计过程控制方法获取所述临床数据的统计值;
预测单元,用于将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
12.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现基于多模态监测信息的术后END风险预警方法,所述方法包括:
获取待测用户在当前时间周期的临床数据,所述临床数据包括基线数据、所述当前时间周期内得到的动态监测数据,和历史时间周期内得到的动态监测数据;
将所述临床数据输入基于卷积理论的深度自编码器,得到所述临床数据对应的特征数据;
利用统计过程控制方法获取所述临床数据的统计值;
将所述临床数据的统计值与当前时间周期对应的预警阈值进行比较,根据比较结果确定下一时间周期的术后END风险预警。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110447993.7A CN113171059B (zh) | 2021-04-25 | 2021-04-25 | 多模态监测信息的术后end风险预警及相关设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110447993.7A CN113171059B (zh) | 2021-04-25 | 2021-04-25 | 多模态监测信息的术后end风险预警及相关设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113171059A CN113171059A (zh) | 2021-07-27 |
CN113171059B true CN113171059B (zh) | 2022-03-01 |
Family
ID=76925245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110447993.7A Active CN113171059B (zh) | 2021-04-25 | 2021-04-25 | 多模态监测信息的术后end风险预警及相关设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113171059B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115641916B (zh) * | 2022-09-29 | 2024-05-03 | 杭州准芯生物技术有限公司 | 一种分子生物学检测方法、计算机设备及可读存储介质 |
CN116719558B (zh) * | 2023-08-09 | 2024-06-21 | 北京比瓴科技有限公司 | gitlab基线检查方法、装置、设备及可读存储介质 |
CN118039070B (zh) * | 2024-04-11 | 2024-06-28 | 四川省肿瘤医院 | 一种用于介入手术的临床护理*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110556178A (zh) * | 2018-05-30 | 2019-12-10 | 西门子医疗有限公司 | 用于医学疗法规划的决策支持*** |
CN110763997A (zh) * | 2019-11-04 | 2020-02-07 | 华北电力大学(保定) | 一种同步电机定子早期故障预警方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10485486B2 (en) * | 2014-11-18 | 2019-11-26 | Baylor College Of Medicine | Clinical metric for predicting onset of cardiorespiratory deterioration in patients |
-
2021
- 2021-04-25 CN CN202110447993.7A patent/CN113171059B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110556178A (zh) * | 2018-05-30 | 2019-12-10 | 西门子医疗有限公司 | 用于医学疗法规划的决策支持*** |
CN110763997A (zh) * | 2019-11-04 | 2020-02-07 | 华北电力大学(保定) | 一种同步电机定子早期故障预警方法 |
Non-Patent Citations (3)
Title |
---|
人工智能技术在神经重症监护中的应用;聂曦明 等;《临床荟萃》;20180405;第33卷(第4期);全文 * |
基于多模态传感数据的诊断预测方法与应用;王丹妮;《中国优秀硕士学位论文全文数据库》;20210115;第2021年卷(第01期);全文 * |
基于异常检测的血透患者死亡风险预测模型研究;娄国锋;《中国优秀硕士学位论文全文数据库》;20210215;第2021年卷(第02期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113171059A (zh) | 2021-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113171059B (zh) | 多模态监测信息的术后end风险预警及相关设备 | |
CN109069081B (zh) | 用于预测、筛查和监测脑病/谵妄的设备、***和方法 | |
JP6857612B2 (ja) | 心血管劣化の警告スコア | |
US11139048B2 (en) | Discovering novel features to use in machine learning techniques, such as machine learning techniques for diagnosing medical conditions | |
Clifford et al. | Robust parameter extraction for decision support using multimodal intensive care data | |
RU2629799C2 (ru) | Оценка и наблюдение острого повреждения легких (ali)/синдрома острой дыхательной недостаточности (ards) | |
EP1133256B1 (en) | A diagnostic tool using a predictive instrument | |
US20220093215A1 (en) | Discovering genomes to use in machine learning techniques | |
RU2657384C2 (ru) | Способ и система неинвазивной скрининговой оценки физиологических параметров и патологий | |
EP1277433B1 (en) | Adjustable coefficients to customize predictive instruments | |
KR102439082B1 (ko) | 딥러닝 기반 심전도 분석데이터를 이용한 일반인의 간질환 보유 가능성 예측 방법 | |
CN115299887B (zh) | 一种动态代谢功能的检测量化方法和*** | |
CN115024725A (zh) | 融合心理状态多参数检测的肿瘤治疗辅助决策*** | |
CN116864104A (zh) | 基于人工智能的慢性血栓栓塞性肺动脉高压风险分级*** | |
Acharya et al. | A systems approach to cardiac health diagnosis | |
Panagiotou et al. | A multi: modal decision making system for an ambient assisted living environment | |
Sahu et al. | Application of fuzzy logic to healthcare industry | |
US20230293103A1 (en) | Analysis device | |
CN113208609A (zh) | 心电信息管理*** | |
Coiera | Automated signal interpretation | |
TWI803893B (zh) | 敗血症之人工智慧輔助醫療診斷方法及其系統 | |
JP7420266B2 (ja) | 分析装置 | |
Bellos et al. | Categorization of COPD patient's health level through the use of the CHRONIOUS wearable platform | |
Cauli | Predicting Myocardial Injury After Noncardiac Surgery | |
US20240206821A1 (en) | Apparatus and a method for predicting a physiological indicator |
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 |