CN113611414B - 急性高原病动态监测与性能评价方法及系统 - Google Patents
急性高原病动态监测与性能评价方法及系统 Download PDFInfo
- Publication number
- CN113611414B CN113611414B CN202110693559.7A CN202110693559A CN113611414B CN 113611414 B CN113611414 B CN 113611414B CN 202110693559 A CN202110693559 A CN 202110693559A CN 113611414 B CN113611414 B CN 113611414B
- Authority
- CN
- China
- Prior art keywords
- data
- monitoring
- time
- training
- abnormal
- 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 47
- 201000010099 disease Diseases 0.000 title claims abstract description 39
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 39
- 230000001154 acute effect Effects 0.000 title claims abstract description 20
- 238000011156 evaluation Methods 0.000 title claims abstract description 15
- 206010021143 Hypoxia Diseases 0.000 claims abstract description 72
- 230000007954 hypoxia Effects 0.000 claims abstract description 58
- 238000012549 training Methods 0.000 claims abstract description 53
- 238000007781 pre-processing Methods 0.000 claims abstract description 7
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 34
- 239000001301 oxygen Substances 0.000 claims description 34
- 229910052760 oxygen Inorganic materials 0.000 claims description 34
- 238000000034 method Methods 0.000 claims description 25
- 230000002159 abnormal effect Effects 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 18
- 230000004044 response Effects 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000001228 spectrum Methods 0.000 claims description 12
- 230000008569 process Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 239000008280 blood Substances 0.000 claims description 6
- 210000004369 blood Anatomy 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000036387 respiratory rate Effects 0.000 claims description 5
- 238000010183 spectrum analysis Methods 0.000 claims description 5
- 208000008445 altitude sickness Diseases 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000013450 outlier detection Methods 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 230000002087 whitening effect Effects 0.000 claims description 2
- 230000001146 hypoxic effect Effects 0.000 description 13
- 230000007958 sleep Effects 0.000 description 5
- 230000006978 adaptation Effects 0.000 description 4
- 238000010219 correlation analysis Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 239000007789 gas Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 208000024891 symptom Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000002045 lasting effect Effects 0.000 description 2
- 206010025482 malaise Diseases 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 230000029058 respiratory gaseous exchange Effects 0.000 description 2
- 230000004622 sleep time Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 206010019233 Headaches Diseases 0.000 description 1
- 206010028813 Nausea Diseases 0.000 description 1
- 208000013738 Sleep Initiation and Maintenance disease Diseases 0.000 description 1
- 206010047700 Vomiting Diseases 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 208000022531 anorexia Diseases 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 206010061428 decreased appetite Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 231100000869 headache Toxicity 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 206010022437 insomnia Diseases 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008693 nausea Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000003938 response to stress Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000008673 vomiting Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/0205—Simultaneously evaluating both cardiovascular conditions and different types of body conditions, e.g. heart and respiratory condition
-
- 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/14542—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 for measuring blood gases
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Epidemiology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Veterinary Medicine (AREA)
- Cardiology (AREA)
- Optics & Photonics (AREA)
- Pulmonology (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明公开了一种急性高原病动态监测与性能评价方法及系统,通过采集间歇性低氧训练的生理数据,数据预处理,慢特征提取及频域分析,最终获得高原病预测指标,本发明将性能监控技术结合信号处理应用于低氧耐受能力动态性能评估中,利用间歇性低氧训练时间序列数据的慢速特征及其频谱特征表示人体低氧耐受能力适应程度,提出了一种有效的高原病风险预测指标。
Description
技术领域
本发明属于高原病风险评估与预防的技术领域,具体涉及一种急性高原病动态监测与性能评价方法及系统。
背景技术
急性高原病(acute mountain sickness,AMS)是高原特发病症之一,其发病率高、危害大,严重影响快速进入高海拔地区人群健康和生命安全。当人在短时间内进入更高海拔地区时,低压低氧的环境会导致头痛、厌食、恶心、呕吐、失眠和倦怠不适等症状。AMS一般是暂时的,随着对高原环境的逐步适应,其不适症状会得到缓解。
间歇性低氧训练通过低氧舱、低氧面罩等工具,使受试者间断性地吸入低于正常氧分压的气体,在每两次低氧呼吸刺激的间隔时间内吸入正常氧分压的气体,造成人体适度缺氧,在不对人体造成损伤的条件下,提高人体对缺氧环境的适应能力。间歇性低氧训练过程中监测的连续生理数据隐含人体低氧适应过程的关键特征,通过对受试者间歇性低氧训练数据动态特性的挖掘,可获知其急进高原后患高原病风险。
目前根据可监测生理数据进行高原病风险预测的算法有:
(1)统计学方法
部分研究只在较为重要的时间段对相关生理参数进行测量,并从中提炼出较有代表性的指标,计算方法不会因数据的改变而变化。
(2)回归分析方法
针对多种相关生理变量对急性高原病的预测能力,研究大多采用logistic回归进行对比分析。
(3)相关性分析方法
相关性分析多用来判断生理指标与高原病发病概率的关联。
综上所述,现有高原病预测技术大多仅在特定时间点采集数据,并未对人体生理数据进行动态连续监控,难以充分表达人体低氧应激反应的关键信息。相关算法仅关注人体低氧适应系统的时域特性,未全面挖掘人体动力学系统的时频域特征。
发明内容
有鉴于此,本发明提供了一种急性高原病动态监测与性能评价方法及系统,能够克服现有技术中仅基于常氧环境、仅采用单点测量生理指标、仅关注时域特性等缺点,可在平原环境下对个体急进高原后患高原病概率与严重程度进行有效预测,该指标可深度挖掘间歇性低氧训练中的连续监测生理变量,提取人体复杂动力学系统低氧应激过程隐含的关键动态特征。
实现本发明的技术方案如下:
一种急性高原病动态监测与性能评价方法,包括以下步骤:
步骤一、对若干名训练对象进行多个训练周期的间歇性低氧训练,每个训练周期内,训练对象摄入浓度周期性波动的氧气,在训练过程中实时采集氧气浓度、呼吸频率、心率(HR)和血氧饱和度(SpO2)数据;
步骤二、采用孤立森林算法智能识别训练过程中所采集的异常数据,将异常数据剔除,并采用线性插值法进行填充,得到处理后数据;
步骤三、对处理后数据进行滞后排列得到轨迹矩阵,并对轨迹矩阵进行奇异值分解、白化处理、广义特征值求解的矩阵变换,得到数据在时域上的慢速变化特征;
步骤四、对所述慢速变化特征进行傅里叶变换后,得到其幅值频谱,取幅值最高点所对应的幅值频率作为高原病预测指标。
进一步地,所述高原病预测指标τ具体为:
一种急性高原病动态监测与性能评价系统,包括数据采集模块、数据预处理模块、慢特征提取模块和频谱分析模块;
所述数据采集模块实时采集训练对象在间歇性低氧训练过程中的氧气浓度、呼吸频率、心率(HR)和血氧饱和度(SpO2)数据;
所述数据预处理模块对数据采集模块所采集数据进行离群值检测、异常值剔除以及缺失值填充,获得无噪声数据;
所述慢特征提取模块对经过预处理的数据进行轨迹矩阵构建及慢特征分析,获得时域尺度上的慢速变化特征;
所述频谱分析模块对慢速变化特征进行傅里叶变换并分析其频谱特征,以频谱峰值频率作为高原病预测指标τ。
进一步地,所述间歇性低氧训练中,训练对象摄入浓度周期性波动的氧气。
进一步地,所述高原病预测指标τ具体为:
有益效果:
1、本发明将性能监控技术结合信号处理应用于低氧耐受能力动态性能评估中,利用间歇性低氧训练时间序列数据的慢速特征及其频谱特征表示人体低氧耐受能力适应程度,提出了一种有效的高原病风险预测指标。
2、本发明最突出的特征是引入性能监控与信号处理的思想,以峰值频率作为高原病预测指标,全面分析间歇性低氧训练数据时域、频域动态特性。相比于静态指标蕴含了更丰富的信息,可以实现实时监测与评估。
3、本发明的数据预处理方法考虑多个生理变量之间的相互作用,更有利于去除来自人体内部的干扰。
4、本发明的慢特征分析方法可从变化复杂且迅速的数据中提取缓慢变化趋势,更加符合人体低氧适应程度变化较为缓慢的特点,具有医学和工程上的可解释性。
5、本发明所获得的高原病风险预测指标无需个体进入高原环境,即可预测急性高原病的严重程度,增强了实际监测的便利性与灵活性。
附图说明
图1是本发明基于间歇性低氧训练的急性高原病预测的动态监测与性能评价系统流程框图。
图2是本发明实施例中低氧响应数据段慢特征示意图。
图3是本发明实施例中低氧数据段频谱指标示意图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种急性高原病动态监测与性能评价方法,如图1所示,具体包括以下步骤:
第一步、对12名训练对象进行间歇性低氧训练。以天为周期进行为期十天的连续间歇性低氧训练(即共有10个间歇性低氧训练周期)。每个训练周期内,训练对象摄入浓度周期性波动的氧气,包括持续5分钟的低浓度氧气阶段(11%-14%的氧气)和持续3分钟的高浓度氧气阶段(35%-38%)。高原缺氧器(GO2Altitude hypoxicator)(BiomedtechAustralia Pty.Ltd.,Melbourne,Australia)用于间歇性低氧训练。此低氧训练装置为面罩式装置,根据气体分离原理,可以提供低浓度和高浓度的氧气,并实时采集氧气浓度与呼吸频率数据。此装置配有脉搏血氧计单元用以测量心率(HR)和血氧饱和度(SpO2),采样时间为1秒。
第二步、分析采集到的生理监测数据,首先进行数据预处理。间歇性低氧训练过程中,受试者深呼吸、摘下氧气面罩等不规范实验行为会导致生理监测数据存在一定偏差,生理监控仪器设备等外部干扰会导致生理监测数据存在噪声。为排除噪声对特征准确度的影响,需要删除异常数据并对时间序列空值进行合理填充。另外,受试者每次低氧训练时间长短不一、每两次低氧训练的间隔存在常氧恢复时间,由于本发明尚未发现常氧段数据的高原病预测能力,因此需要从原始序列数据中分离出低氧响应数据段。
本发明使用孤立森林算法自动检测生理监测数据中的离群值。孤立森林是一种适用于连续数据的无监督异常检测方法,在特征空间里,分布稀疏的区域表示事件发生在该区域的概率很低,因而可以认为落在这些区域里的数据是异常的。
定义x=[x1(t),x2(t)…xd(t)t=1,2…n]T为时间t∈[1,n]内动态监控序列数据,其中d为监控生理变量的维度。对数据构建num个iTree组成孤立森林,实现步骤如下在如图):
(1)从x中随机选择ρ个时间点数据作为子样本,放入树的根节点。
(2)随机指定一个维度,在当前节点数据中随机产生一个切割点p,切割点产生于当前节点数据中指定维度的最大值和最小值之间。
(3)以此切割点生成一个超平面,然后将当前节点数据空间划分为2个子空间:把指定维度里小于p的数据放在当前节点的左边,把大于等于p的数据放在当前节点的右边。
(4)在子节点中递归步骤(2)和(3),不断构造新的子节点,直到子节点中只有一个数据(无法再继续切割)或子节点已到达限定高度。
获得num个iTree之后,对于每一个时间点数据,令其遍历每一棵iTree,并计算该点在每棵树的高度平均值。高度平均值低于设定阈值的时间点数据即为异常。
采用孤立森林经过多次迭代识别出生理变量异常时间点集合{Tb},定义{xsp(t):t=1,2…n}为时间t∈[1,n]内SpO2监测数据,删除其在异常时间点上的原始监测值:
采用线性插值法对经过异常值剔除的SpO2监测数据进行空值填充。根据相应时间氧气浓度值,将氧气浓度低于14%的时间段定义为低氧响应数据段,从而将经过离群值处理的某完整训练周期内的SpO2监测数据划分为M个低氧响应数据段。定义{xsp(t):t=1,2…nm}为单个低氧响应数据段在时间t∈[1,nm]上的SpO2监测数据。
第三步、对上述SpO2低氧响应数据段进行慢速变化特征提取。针对单个低氧响应数据段{xsp(t):t=1,2…nm},选择合适的窗口长度L将原始时间序列进行滞后排列得到轨迹矩阵:
通常情况取L<nm/2,本发明取L=50,令K=nm-L+1,则轨迹矩阵X为L×K的矩阵
将轨迹矩阵表示为L维时间序列其中t∈[1,K]为序列所处时间范围。慢特征分析算法旨在找到一系列特征序列s(t)=[s1(t),s2(t)…sL(t)|t=1,2…K]中变化最慢的特征序列si(t)满足:
s.t.<si(t)>t=0
从时间序列xsp(t)提取s(t)的过程表示为
s(t)=Wxsp(t)
其中W=[w1,w2…wL]T为权值矩阵。W的计算方法如下:
z=Λ-1/2UTx=Qxsp
W=PQ=PΛ-1/2UT
其中,Λ、Ω为特征值对角矩阵。综上可得变换后的一系列特征序列s(t)=Wxsp=WQ-1z=Pz。当Λ、Ω特征值由小到大排列时,变化最慢的特征序列为s1(t),如图2。
第四步、对低氧响应数据段慢特征进行频谱分析。对上述慢特征{s1(t):t=1,2…K}进行傅里叶变换:
峰值频率表示时间序列经过傅里叶变换后信号成分所集中的频率点,因此指标τ可以表征低氧响应数据在频域范围内信号的集中带宽,如图3。
第五步、在间歇性低氧训练前后,用相同的低氧训练面罩装置进行一次低氧耐受能力测试并获得低氧耐受能力指标(hypoxia tolerance index,HTI)。HTI为在420秒内训练对象的SpO2保持在85%以上的持续时间,用HTI1和HTI2代表间歇性低氧训练前后的低氧耐受能力水平。定义训练前后低氧耐受能力水平差值ΔHTI=HTI2-HTI1。在最后一个低氧训练周期后,训练对象在9pm被立即送入常压低氧舱(空军医学研究所,北京,中国)进行一夜的睡眠实验。此常压低氧舱模拟海拔4300米地区的大气环境(氧气浓度约为12%)。睡眠指标包括睡眠时间,深睡时间和比例,夜间觉醒时间、生理指标如血氧饱和度(SpO2)和心率(HR)。微动敏感床垫睡眠检测系统RS-611型(新兴阳升科技有限公司,北京,中国)用于睡眠监测。此系统可以在不影响正常睡眠的情况下进行精确的测量。第二日7am,对训练对象进行路易斯湖(LLS)测试。此实施例对上述预测指标实际预测能力进行了性能验证分析与评价。低氧耐受能力指标LLS、HTI2与低氧舱中测得的生理指标具有显著的Spearman相关性,如表1:
表1低氧耐受能力指标与低氧舱生理指标相关性分析
表中r表示Spearman相关性,p表示显著性水平,p<0.05表示有统计学差异,p<0.01表示有显著统计学差异,p<0.001表示有极其显著的统计学差异。
表1中HTI2比HTI1与低氧舱生理指标具有更为显著的相关性,此结果说明人体的低氧耐受能力通过间歇性低氧训练有所提升,因此间歇性低氧训练过程中的监测数据包含了人体低氧耐受能力的关键信息,适合用户高原病预测指标的挖掘。
所有训练对象的间歇性低氧训练数据得出的急性高原病预测指标τ与在低氧舱中测得的生理指标具有显著的Spearman相关性,如表2:
表2急性高原病预测指标τ与低氧舱生理指标相关性分析
表中r表示Spearman相关性,p表示显著性水平,p<0.05表示有统计学差异,p<0.01表示有显著统计学差异,p<0.001表示有极其显著的统计学差异。
表2说明本发明所提出的高原病预测指标τ与判断人体低氧耐受能力的医学指标具有显著相关性,因此指标τ具有可靠的预测作用。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种急性高原病动态监测与性能评价方法,其特征在于,包括以下步骤:
步骤一、对若干名训练对象进行多个训练周期的间歇性低氧训练,每个训练周期内,训练对象摄入浓度周期性波动的氧气,在训练过程中实时采集氧气浓度、呼吸频率、心率和血氧饱和度数据;
步骤二、采用孤立森林算法智能识别训练过程中所采集的异常数据,将异常数据剔除,并采用线性插值法进行填充,得到处理后数据;
孤立森林是一种适用于连续数据的无监督异常检测方法,在特征空间里,分布稀疏的区域表示事件发生在该区域的概率很低,因而可以认为落在这些区域里的数据是异常的;
定义x=[x1(t),x2(t)…xd(t)t=1,2...n]T为时间t∈[1,n]内动态监控序列数据,其中d为监控生理变量的维度;对数据构建num个iTree组成孤立森林,实现步骤如下:
(1)从x中随机选择ρ个时间点数据作为子样本,放入树的根节点;
(2)随机指定一个维度,在当前节点数据中随机产生一个切割点p,切割点产生于当前节点数据中指定维度的最大值和最小值之间;
(3)以此切割点生成一个超平面,然后将当前节点数据空间划分为2个子空间:把指定维度里小于p的数据放在当前节点的左边,把大于等于p的数据放在当前节点的右边;
(4)在子节点中递归步骤(2)和(3),不断构造新的子节点,直到子节点中只有一个数据或子节点已到达限定高度;
获得num个iTree之后,对于每一个时间点数据,令其遍历每一棵iTree,并计算该点在每棵树的高度平均值,高度平均值低于设定阈值的时间点数据即为异常;
采用孤立森林经过多次迭代识别出生理变量异常时间点集合{Tb},定义{xsp(t):t=1,2…n}为时间t∈[1,n]内SpO2监测数据,删除其在异常时间点上的原始监测值:
采用线性插值法对经过异常值剔除的SpO2监测数据进行空值填充;根据相应时间氧气浓度值,将氧气浓度低于14%的时间段定义为低氧响应数据段,从而将经过离群值处理的某完整训练周期内的SpO2监测数据划分为M个低氧响应数据段;定义{xsp(t):t=1,2…nm}为单个低氧响应数据段在时间t∈[1,nm]上的SpO2监测数据;
步骤三、对处理后数据进行滞后排列得到轨迹矩阵,并对轨迹矩阵进行奇异值分解、白化处理、广义特征值求解的矩阵变换,得到数据在时域上的慢速变化特征;
步骤四、对所述慢速变化特征进行傅里叶变换后,得到其幅值频谱,取幅值最高点所对应的幅值频率作为高原病预测指标。
3.一种急性高原病动态监测与性能评价系统,其特征在于,包括数据采集模块、数据预处理模块、慢特征提取模块和频谱分析模块;
所述数据采集模块实时采集训练对象在间歇性低氧训练过程中的氧气浓度、呼吸频率、心率和血氧饱和度数据;
所述数据预处理模块对数据采集模块所采集数据进行离群值检测、异常值剔除以及缺失值填充,获得无噪声数据;
采用孤立森林算法智能识别所采集的异常数据,将异常数据剔除,并采用线性插值法进行填充,得到处理后数据;
孤立森林是一种适用于连续数据的无监督异常检测方法,在特征空间里,分布稀疏的区域表示事件发生在该区域的概率很低,因而可以认为落在这些区域里的数据是异常的;
定义x=[x1(t),x2(t)…xd(t)t=1,2…n]T为时间t∈[1,n]内动态监控序列数据,其中d为监控生理变量的维度;对数据构建num个iTree组成孤立森林,实现步骤如下:
(5)从x中随机选择ρ个时间点数据作为子样本,放入树的根节点;
(6)随机指定一个维度,在当前节点数据中随机产生一个切割点p,切割点产生于当前节点数据中指定维度的最大值和最小值之间;
(7)以此切割点生成一个超平面,然后将当前节点数据空间划分为2个子空间:把指定维度里小于p的数据放在当前节点的左边,把大于等于p的数据放在当前节点的右边;
(8)在子节点中递归步骤(2)和(3),不断构造新的子节点,直到子节点中只有一个数据或子节点已到达限定高度;
获得num个iTree之后,对于每一个时间点数据,令其遍历每一棵iTree,并计算该点在每棵树的高度平均值,高度平均值低于设定阈值的时间点数据即为异常;
采用孤立森林经过多次迭代识别出生理变量异常时间点集合{Tb},定义{xsp(t):t=1,2…n}为时间t∈[1,n]内SpO2监测数据,删除其在异常时间点上的原始监测值:
采用线性插值法对经过异常值剔除的SpO2监测数据进行空值填充;根据相应时间氧气浓度值,将氧气浓度低于14%的时间段定义为低氧响应数据段,从而将经过离群值处理的某完整训练周期内的SpO2监测数据划分为M个低氧响应数据段;定义{xsp(t):t=1,2…nm}为单个低氧响应数据段在时间t∈[1,nm]上的SpO2监测数据;
所述慢特征提取模块对经过预处理的数据进行轨迹矩阵构建及慢特征分析,获得时域尺度上的慢速变化特征;
所述频谱分析模块对慢速变化特征进行傅里叶变换并分析其频谱特征,以频谱峰值频率作为高原病预测指标τ。
4.如权利要求3所述的一种急性高原病动态监测与性能评价系统,其特征在于,所述间歇性低氧训练中,训练对象摄入浓度周期性波动的氧气。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110693559.7A CN113611414B (zh) | 2021-06-22 | 2021-06-22 | 急性高原病动态监测与性能评价方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110693559.7A CN113611414B (zh) | 2021-06-22 | 2021-06-22 | 急性高原病动态监测与性能评价方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113611414A CN113611414A (zh) | 2021-11-05 |
CN113611414B true CN113611414B (zh) | 2023-03-28 |
Family
ID=78303615
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110693559.7A Active CN113611414B (zh) | 2021-06-22 | 2021-06-22 | 急性高原病动态监测与性能评价方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113611414B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114464319B (zh) * | 2022-01-13 | 2022-11-01 | 北京理工大学 | 一种基于慢特征分析和深度神经网络的ams易感性评估系统 |
CN114748862B (zh) * | 2022-03-28 | 2023-04-11 | 北京理工大学 | 一种低氧耐受能力动态性能主动测试系统 |
CN117059270A (zh) * | 2023-08-14 | 2023-11-14 | 北京理工大学 | 一种结合医学先验知识伪标签的急性高原病风险评估系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111812332A (zh) * | 2020-06-23 | 2020-10-23 | 中国人民解放军军事科学院军事医学研究院 | 一种用于检测高原缺氧的生物标志物及其应用 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2937089A1 (en) * | 2014-04-23 | 2015-10-28 | Linde AG | Device for the treatment of, treatment of complications arising from, and/or prevention of respiratory disorders |
CN110364259B (zh) * | 2019-05-30 | 2022-05-31 | 中国人民解放军总医院 | 一种高原疾病预测方法、系统、介质及电子设备 |
CN111513726B (zh) * | 2020-04-07 | 2021-04-13 | 北京理工大学 | 一种基于iht动态性能评估ams风险的系统 |
-
2021
- 2021-06-22 CN CN202110693559.7A patent/CN113611414B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111812332A (zh) * | 2020-06-23 | 2020-10-23 | 中国人民解放军军事科学院军事医学研究院 | 一种用于检测高原缺氧的生物标志物及其应用 |
Non-Patent Citations (1)
Title |
---|
新兴可穿戴生理感知技术助力急性高原病的监测和辅助诊疗;史大威等;《生命科学仪器》;20200625;第18卷(第3期);第3-22页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113611414A (zh) | 2021-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113611414B (zh) | 急性高原病动态监测与性能评价方法及系统 | |
Montesinos et al. | Multi-modal acute stress recognition using off-the-shelf wearable devices | |
US8204583B2 (en) | System for seizure monitoring and detection | |
EP2018825A1 (en) | Methods and device to quantify human physical activity pattern | |
Belgacem et al. | ECG based human authentication using wavelets and random forests | |
US10368792B2 (en) | Method for detecting deception and predicting interviewer accuracy in investigative interviewing using interviewer, interviewee and dyadic physiological and behavioral measurements | |
CN111466876A (zh) | 一种基于fNIRS和图神经网络的阿尔兹海默症辅助诊断系统 | |
CN111407262A (zh) | 基于blstm的睡眠分期方法以及基于blstm进行睡眠分期的装置 | |
CN116895382B (zh) | 基于数据分析的人体健康监测管理系统 | |
Jaafar et al. | Estimation of breathing rate and heart rate from photoplethysmogram | |
CN110623678A (zh) | 一种血糖测量装置及其数据处理方法、存储介质 | |
Karandikar et al. | Detection of sleep apnea events via tracking nonlinear dynamic cardio-respiratory coupling from electrocardiogram signals | |
Bruun et al. | Automatic atrial fibrillation detection: A novel approach using discrete wavelet transform and heart rate variability | |
CN114464319B (zh) | 一种基于慢特征分析和深度神经网络的ams易感性评估系统 | |
CN113842124A (zh) | 一种基于生理健康指标的精神状态预测方法、系统及设备 | |
Bashar et al. | Extraction of heart rate from PPG Signal: a machine learning approach using decision tree regression algorithm | |
Pathan et al. | A machine learning based human activity recognition during physical exercise using wavelet packet transform of ppg and inertial sensors data | |
CN116035536A (zh) | 一种睡眠行为活跃水平检测量化的方法、系统和装置 | |
Haritha et al. | Automating anxiety detection using respiratory signal analysis | |
Akhter et al. | Feature selection for heart rate variability based biometric recognition using genetic algorithm | |
CN118098581A (zh) | 情绪状态的监测方法及系统 | |
CN110956406B (zh) | 一种基于心率变异性的团队协同能力的评估方法 | |
CN105266764A (zh) | 一种中医宗气测评装置 | |
Wang et al. | Heart rate monitoring from wrist-type PPG based on singular spectrum analysis with motion decision | |
Assadi et al. | QRS complex based human identification |
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 |