CN116421457A - 一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法 - Google Patents
一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法 Download PDFInfo
- Publication number
- CN116421457A CN116421457A CN202310206600.2A CN202310206600A CN116421457A CN 116421457 A CN116421457 A CN 116421457A CN 202310206600 A CN202310206600 A CN 202310206600A CN 116421457 A CN116421457 A CN 116421457A
- Authority
- CN
- China
- Prior art keywords
- perfusion pressure
- coronary perfusion
- prediction
- cardiopulmonary resuscitation
- coronary
- 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.)
- Pending
Links
- 230000010412 perfusion Effects 0.000 title claims abstract description 93
- 238000002680 cardiopulmonary resuscitation Methods 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000006835 compression Effects 0.000 claims abstract description 30
- 238000007906 compression Methods 0.000 claims abstract description 30
- 230000002068 genetic effect Effects 0.000 claims abstract description 13
- 238000000605 extraction Methods 0.000 claims abstract description 11
- 230000008569 process Effects 0.000 claims description 21
- 238000012549 training Methods 0.000 claims description 10
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 5
- 230000004872 arterial blood pressure Effects 0.000 claims description 4
- 230000001746 atrial effect Effects 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims 1
- 210000004351 coronary vessel Anatomy 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 10
- 238000010801 machine learning Methods 0.000 description 7
- 238000013186 photoplethysmography Methods 0.000 description 5
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 206010049418 Sudden Cardiac Death Diseases 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 230000033764 rhythmic process Effects 0.000 description 3
- 208000014221 sudden cardiac arrest Diseases 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 210000000349 chromosome Anatomy 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012806 monitoring device Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000010187 selection method Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000003066 decision tree Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000009396 hybridization Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000007477 logistic regression Methods 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 230000013011 mating Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000013642 negative control Substances 0.000 description 1
- 230000004962 physiological condition Effects 0.000 description 1
- 230000035790 physiological processes and functions Effects 0.000 description 1
- 239000013641 positive control Substances 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000029058 respiratory gaseous exchange Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000000087 stabilizing effect Effects 0.000 description 1
- 230000001839 systemic circulation Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61H—PHYSICAL THERAPY APPARATUS, e.g. DEVICES FOR LOCATING OR STIMULATING REFLEX POINTS IN THE BODY; ARTIFICIAL RESPIRATION; MASSAGE; BATHING DEVICES FOR SPECIAL THERAPEUTIC OR HYGIENIC PURPOSES OR SPECIFIC PARTS OF THE BODY
- A61H31/00—Artificial respiration or heart stimulation, e.g. heart massage
- A61H31/004—Heart stimulation
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/25—Bioelectric electrodes therefor
- A61B5/279—Bioelectric electrodes therefor specially adapted for particular uses
- A61B5/28—Bioelectric electrodes therefor specially adapted for particular uses for electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- 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/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
-
- 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/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- 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/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2131—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on a transform domain processing, e.g. wavelet transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- 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)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Artificial Intelligence (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Pulmonology (AREA)
- General Physics & Mathematics (AREA)
- Fuzzy Systems (AREA)
- Physical Education & Sports Medicine (AREA)
- Epidemiology (AREA)
- Pain & Pain Management (AREA)
- Emergency Medicine (AREA)
- Optics & Photonics (AREA)
- Rehabilitation Therapy (AREA)
- Genetics & Genomics (AREA)
- Computational Linguistics (AREA)
Abstract
本公开属于心肺复苏技术领域,具体涉及一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法,包括:获取模块,其被配置为获取心电信号和光电容积脉搏波信号;提取模块,其被配置为基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;预测模块,其被配置为根据所提取的信号特征和实时冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。本公开基于心电信号和脉搏信号预测冠状动脉灌注压,根据所得到的冠状动脉灌注压预测值实时调节心肺复苏的按压深度,提高闭环自动胸外按压控制的鲁棒性和自适应性。
Description
技术领域
本公开属于心肺复苏技术领域,具体涉及一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息,不必然构成在先技术。
心脏骤停是指心脏突然停止泵血,同时引起全身循环中断,呼吸和意识丧失的一种恶性疾病。心肺复苏(Cardiopulmonary resuscitation,简称CPR)可以在患者心脏骤停时维持一定的血液循环,是心脏骤停的关键抢救方法。但是徒手胸外按压对施救人员的体力需求较大,且较难达到按压的标准要求。市面上的一些自动胸外按压仪虽然可以按照标准要求进行按压,但无法根据患者生理状况进行自适应的调整。
发明内容
为了解决上述问题,本公开提出了一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法,基于心电信号和脉搏信号预测冠状动脉灌注压(Coronary PerfusionPressure,简称CPP),根据所得到的冠状动脉灌注压预测值实时调节心肺复苏的按压深度,提高闭环自动胸外按压控制的鲁棒性和自适应性。
根据一些实施例,本公开的第一方案提供了一种冠状动脉灌注压的无创预测系统,采用如下技术方案:
一种冠状动脉灌注压的无创预测系统,包括:
获取模块,其被配置为获取心电信号和光电容积脉搏波信号;
提取模块,其被配置为基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;
预测模块,其被配置为根据所提取的信号特征和冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。
作为进一步的技术限定,在所述提取模块中,通过离散小波变换,对所述获取模块所得到的信号进行不同频带的分解,分别计算不同频带的细节系数和信号特征,根据所得到的不同的细节系数得到不同频带的信号特征,完成所获取信号的特征提取。
作为进一步的技术限定,在所述预测模块中,实时监测有创动脉压和有创右房压,根据所监测的有创动脉压的波谷点的幅值与对应时刻的有创右房压,得到该时刻的冠状动脉灌注压,所得到的该时刻的冠状动脉灌注压为实时冠状动脉灌注压,根据所得到的实时冠状动脉灌注压进行无创预测模型的训练。
作为进一步的技术限定,在所述预测模块中,在每一次迭代过程中,计算种群的信息复杂度分数,结束标准为迭代次数达到设定值;迭代结束后,寻找迭代过程中信息复杂度分数最小的种群,作为特征选择的结果。
进一步的,计算模型在测试集上的性能,使用平均绝对误差、均方误差、均方根误差、拟合优度以及信息复杂度分数;均绝对误差、均方误差和均方根误差越小,模型性能越优;拟合优度越接近1,模型性能越优;信息复杂度分数越小,模型性能越优。
进一步的,在模型训练过程中,每一次迭代都会有一个信息复杂度分数,在有限的迭代次数里,当模型的复杂度和预测精度达到相对平衡时,会出现信息复杂度的最小值;即得到信息复杂度的极小值。
根据一些实施例,本公开的第二方案提供了一种冠状动脉灌注压的无创预测方法,采用如下技术方案:
一种冠状动脉灌注压的无创预测方法,包括:
获取心电信号和光电容积脉搏波信号;
基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;
根据所提取的信号特征和冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。
根据一些实施例,本公开的第三方案提供了一种基于无创预测冠状动脉灌注压的心肺复苏系统,采用如下技术方案:
一种基于无创预测冠状动脉灌注压的心肺复苏系统,包括:
冠状动脉灌注压的无创预测模块,其被配置为采用了如第一方案中所述的冠状动脉灌注压的无创预测系统,获取冠状动脉灌注压预测值;
控制模块,其被配置为根据所获取的冠状动脉灌注压预测值和冠状动脉灌注压,通过PID控制,调节优化按压参数;
心肺复苏模块,其被配置为根据所得到的按压参数,调节心肺复苏按压深度,完成心肺复苏的自适应调节。
作为进一步的技术限定,在所述控制模块中,将所得到的当前冠状动脉灌注压值与设定的冠状动脉灌注压标准值相减,得到当前误差;若PID控制器的输出为正,则控制心肺复苏仪的按压深度增大;若PID控制器的输出为负则控制按压深度减小,输出值越大,增大和减小速度越快,实现按压参数的调节。
作为进一步的技术限定,在所述心肺复苏模块中,根据实时预测的冠状动脉灌注压值与设定冠状动脉灌注压参考值之间的误差,结合PID控制器,自适应调节心肺复苏仪的按压深度,实现预测冠状动脉灌注压值与设定冠状动脉灌注压值的逐步趋近,完成自适应心肺复苏的调节。
与现有技术相比,本公开的有益效果为:
本公开使用基于机器学习的CPP预测模型,使用无创的心电信号(Electrocardiogram,简称ECG)和光电容积脉搏波信号(Photoplethysmography,简称PPG)信号预测CPP,进而实现按压深度和按压频率的深入指导,实现个体化按压。
本公开在CPP预测模型训练的过程中使用遗传算法(Genetic Algorithm,简称GA)和信息复杂度(Information complexity,简称ICOMP)的混合机器学习模型选择特征,在降低模型复杂度的同时,保持较高的模型训练效果。
本公开在自适应调节心肺复苏的过程中,利用PID控制器,根据冠状动脉灌注压预测值自动调整自动胸外按压器的PID控制参数,使整个胸外按压闭环控制系统达到最佳的复苏效果。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解,本公开的示意性实施例及其说明用于解释本公开,并不构成对本公开的不当限定。
图1是本公开实施例一中的冠状动脉灌注压的无创预测系统的流程图;
图2是本公开实施例一中的遗传算法的流程图;
图3是本公开实施例一中的特征子集和模型选择过程的流程图;
图4是本公开实施例三中的基于无创预测冠状动脉灌注压的心肺复苏系统示意图;
图5是本公开实施例三中的心肺复苏的流程图;
图6是本公开实施例三中的PID控制系统的结构图;
其中,1、血氧饱和度探头;2、第一心电电极;3、第二心电电极;4、第三心电电极;5、第四心电电极;6、第五心电电极;7、自动胸外按压心肺复苏仪;8、生理信号监测装置;9、生理信号分析装置。
具体实施方式
下面结合附图与实施例对本公开作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本公开提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本公开的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本公开中,术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系,只是为了便于叙述本公开各部件或元件结构关系而确定的关系词,并非特指本公开中任一部件或元件,不能理解为对本公开的限制。
本公开中,术语如“固接”、“相连”、“连接”等应做广义理解,表示可以是固定连接,也可以是一体地连接或可拆卸连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员,可以根据具体情况确定上述术语在本实公开中的具体含义,不能理解为对本公开的限制。
在不冲突的情况下,本公开中的实施例及实施例中的特征可以相互组合。
实施例一
本公开实施例一介绍了一种冠状动脉灌注压的无创预测系统。
一种冠状动脉灌注压的无创预测系统,包括:
获取模块,其被配置为获取心电信号和光电容积脉搏波信号;
提取模块,其被配置为基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;
预测模块,其被配置为根据所提取的信号特征和冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。
如图1所示,使用心电信号和脉搏波信号预测CPP:在获得ECG和PPG信号之后,对其进行预处理操作,提取特征;在进行特征选择之后进入已训练好的机器学习模型,最后输出预测的CPP结果,得到CPP预测值。
在本实施例中,采用离散小波变换(Discrete Wavelet Transform,DWT)进行信号的特征提取。下面以ECG信号为例展开详细介绍:
DWT将一个平方可积的x(t)∈V0分解为详细的Wj和近似的Vj子级数。在这个框架中,投影到Vj上的x(t)在尺度j上的近似信息记为:
其中,φj,p是比例函数,aj,p=<x(t),φj,p(t)>被称为近似系数,其中<·>是内积算子。同样地,投影到Wj上的x(t)在尺度j上的详细信息被定义为:
因此,DWT通过g(n)和h(n)将信号x(t)分解为近似系数aj,p(低频分量)和细节系数dj,p(高频分量)系数:
在所提出的特征提取过程中,初始ECG信号被分为3秒段。经过划分过程后,通过DWT分解为高频带(细节系数,dj)和低频带(近似系数,aj),直到第五级分解(j=1,2,...,5)。因此,可以从每个层次上的这些子带中分别计算出许多显著的特征。
在本实施例中,为减少数据冗余,所有的细节系数(dj,j=1,2,...,5)和最后一个分解水平(第5级)的近似系数都被用于特征提取过程。每个子带可以分别计算出各种特征,如曲线长度,能量,最大值,最小值,中位数,平均值,熵,范围,峭度,偏度,趋势,过零点数等12个特征;可从第一个细节系数(d1)计算出来。这样,就可以从其他子带中计算出相同的特征。
同样地,PPG可以通过上述过程类似的获取特征。
对于预测模型的标签,使用通过ABP信号和RAP信号计算出的CPP值。找到每个片段的ABP信号的波谷点,将该点的幅度值减去对应时间点的RAP的值;每个片段的CPP值为计算出的CPP平均值;将该平均值作为预测模型的标签。通过分割数据进行交叉验证,之后进行特征子集和模型选择过程,最后根据计算出的ICOMP评分和其他性能指标选择最佳模型。
信息复杂度标准ICOMP旨在惩罚处理评估估计参数之间相关性的复杂性增加,而不是仅仅惩罚不必要的参数。对于单变量和多变量模型,ICOMP的结构为:
其中,n是样本的数目,误差方差是由实际组标签和预测组标签之间的均方根误差来估计的(实际组标签为经过实验获得的CPP值,预测组标签为经过模型预测得到的CPP值,通过计算两者的均方根误差得到误差方差),/>为模型参数向量的稳定和平滑凸和协方差矩阵的估计量:
通过稳定和平滑凸和协方差矩阵的估计量,它能够克服病态协方差矩阵,除了选择最优核函数外,还能使误差最小化。
为降低系统输入的维数和模型的复杂性,本实施例使用了一种新的特征子集选择方法,并将此过程集成到机器学习模型的训练过程中。基本上,所使用的特征子集选择方法是基于遗传算法GA,其中ICOMP作为适应度函数。
GA根据每一代的种群大小产生一定数量的特征子集,这些子集的创建是通过使用一些特殊的操作符来执行的,如交叉、突变、迁移和选择;在与这些特征子集对应的数据集上估计机器学习模型。对具有不同特征子集的估计模型的性能进行了评分。这些分数显示了特征子集对拟合优度的贡献。此外,它们还可以用来在再生过程中为下一代创建新的特征子集。
GA循环中的算子如图2所示,详细步骤为:
a.随机创建一个初始种群,包括被称为染色体的特征子集,它只包含二进制数{0,1}。例如,如果特征Fj包含在特征子集,则二进制数为1,否则为0。
b.为所有的子集(染色体)分别训练任何机器学习模型,并根据ICOMP标准进行评分。
c.将一些得分最高的有利子集(近1%)发送给下一代。
d.利用启发式或随机算子,如随机普遍抽样,将所青睐的特征子集(父本)收集到交配池中。(此过程称为选择)
e.将突变和交叉过程应用于选定的父本,并为下一代创建新的特征子集(子集)。(c,d和e阶段被称为繁殖过程)。
f.如果满足其中一个停止标准,则终止该程序;否则,返回第三阶段。
g.保持最佳功能。
为保证获得最佳的特征子集和回归模型,本实施例中的过程可以重复多次。从本质上讲,具有特征选择的训练过程也可能估计一些备选模型。为在这些替代方案中选择最佳模型,可以考虑某些独特的性能统计,如平均绝对误差,均方误差,均方根误差和拟合优度,以及ICOMP标准和其他标准。特征子集和模型选择的流程图如图3所示。
特征子集选择结束的标准是ICOMP分数出现极小值。模型选择标准是该模型在测试集上的性能最优;可以选择的模型包括支持向量回归(核函数有线性核、二次核、三次核、径向核以及高斯核),集成学习(通过自适应增强(Adaboost)与决策树和K近邻等弱学习器进行杂交),以及逻辑回归。
在每一次迭代过程中,计算种群的ICOMP分数,结束标准为迭代次数达到设定值。结束之后,寻找迭代过程中ICOMP分数最小的种群,作为特征选择的结果。
计算模型在测试集上的性能,使用平均绝对误差,均方误差,均方根误差,拟合优度以及ICOMP分数。理论上讲,均绝对误差,均方误差和均方根误差越小,模型性能越优;拟合优度越接近1,模型性能越优;ICOMP分数越小,模型性能越优。由上述结果确定最优模型,作为模型选择的结果。
本实施例在CPP预测模型训练的过程中使用遗传算法(Genetic Algorithm,简称GA)和信息复杂度(Information complexity,简称ICOMP)的混合机器学习模型选择特征,在降低模型复杂度的同时,保持较高的模型训练效果。
实施例二
本公开实施例二介绍了一种冠状动脉灌注压的无创预测方法,采用了实施例一中所介绍的冠状动脉灌注压的无创预测系统。
一种冠状动脉灌注压的无创预测方法,包括:
获取心电信号和光电容积脉搏波信号;
基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;
根据所提取的信号特征和实时冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。
详细步骤与实施例一提供的冠状动脉灌注压的无创预测系统相同,在此不再赘述。
实施例三
本公开实施例三介绍了一种基于无创预测冠状动脉灌注压的心肺复苏系统。
一种基于无创预测冠状动脉灌注压的心肺复苏系统,其特征在于,包括:
冠状动脉灌注压的无创预测模块,其被配置为采用了实施例一中的冠状动脉灌注压的无创预测系统,获取冠状动脉灌注压预测值;
控制模块,其被配置为根据所获取的冠状动脉灌注压预测值和实时冠状动脉灌注压,通过PID控制,调节优化按压参数;
心肺复苏模块,其被配置为根据所得到的按压参数,调节心肺复苏按压深度,完成心肺复苏的自适应调节。
如图4所示,本实施例通过血氧饱和度探头1实时监测患者的血氧饱和度,在第一心电电极2、第二心电电极3、第三心电电极4、第四心电电极5和第五心电电极6的作用下获取ABP和RAP,通过生理信号监测装置8获取光电容积脉搏波信号和心电信号;在生理信号分析装置9的作用下基于实时ECG和PPG信号,通过预测模型获取的CPP值来判断按压质量并反馈至心肺复苏仪。
如图5所示,在发现患者之后,对患者进行心律监测,若患者心律正常,即无需对患者进行操作;若患者心律异常,将心肺复苏仪放置在患者身上,并启动。在整个心肺复苏过程中,实时监测患者的生理信号,通过预测模型获取CPP值,通过PID控制系统调节按压参数,直至患者成功复苏。
本实施例采用PID控制系统对按压深度进行自动闭环控制调节,如图6所示,通过PID系统的输出u(t)调节心肺复苏仪的按压深度d(t),调整按压深度,从而调整生理状态,形成自适应心肺复苏系统,可最大限度地提高闭环自动胸外按压控制的鲁棒性和自适应性;在获得当前CPP值c(t)后,与设定的CPP标准值r(t)相减,得到当前误差e(t),PID控制器的输出u(t)为正控制心肺复苏仪的按压深度d(t)增大;为负则控制按压深度减小。
以上所述仅为本公开的优选实施例而已,并不用于限制本公开,对于本领域的技术人员来说,本公开可以有各种更改和变化。凡在本公开的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。
上述虽然结合附图对本公开的具体实施方式进行了描述,但并非对本公开保护范围的限制,所属领域技术人员应该明白,在本公开的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。
Claims (10)
1.一种冠状动脉灌注压的无创预测系统,其特征在于,包括:
获取模块,其被配置为获取心电信号和光电容积脉搏波信号;
提取模块,其被配置为基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;
预测模块,其被配置为根据所提取的信号特征和冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。
2.如权利要求1中所述的冠状动脉灌注压的无创预测系统,其特征在于,在所述提取模块中,通过离散小波变换,对所述获取模块所得到的信号进行不同频带的分解,分别计算不同频带的细节系数和信号特征,根据所得到的不同的细节系数得到不同频带的信号特征,完成所获取信号的特征提取。
3.如权利要求1中所述的冠状动脉灌注压的无创预测系统,其特征在于,在所述预测模块中,实时监测有创动脉压和有创右房压,根据所监测的有创动脉压的波谷点的幅值与对应时刻的有创右房压,得到该时刻的冠状动脉灌注压,所得到的该时刻的冠状动脉灌注压为实时冠状动脉灌注压,根据所得到的实时冠状动脉灌注压进行无创预测模型的训练。
4.如权利要求1中所述的冠状动脉灌注压的无创预测系统,其特征在于,在所述预测模块中,在每一次迭代过程中,计算种群的信息复杂度分数,结束标准为迭代次数达到设定值;迭代结束后,寻找迭代过程中信息复杂度分数最小的种群,作为特征选择的结果。
5.如权利要求4中所述的冠状动脉灌注压的无创预测系统,其特征在于,计算模型在测试集上的性能,使用平均绝对误差、均方误差、均方根误差、拟合优度以及信息复杂度分数;均绝对误差、均方误差和均方根误差越小,模型性能越优;拟合优度越接近1,模型性能越优;信息复杂度分数越小,模型性能越优。
6.如权利要求4中所述的冠状动脉灌注压的无创预测系统,其特征在于,在模型训练过程中,每一次迭代都会有一个信息复杂度分数,在有限的迭代次数里,当模型的复杂度和预测精度达到相对平衡时,会出现信息复杂度的最小值;即得到信息复杂度的极小值。
7.一种冠状动脉灌注压的无创预测方法,其特征在于,包括:
获取心电信号和光电容积脉搏波信号;
基于离散小波变换对所获取的信号进行不同频带的分解,提取所获取信号的特征;
根据所提取的信号特征和冠状动脉灌注压,以信息复杂度为适应度函数,基于遗传算法构建无创预测模型,得到冠状动脉灌注压预测值,完成冠状动脉灌注压的无创预测。
8.一种基于无创预测冠状动脉灌注压的心肺复苏系统,其特征在于,包括:
冠状动脉灌注压的无创预测模块,其被配置为采用了如权利要求1-6中任一项所述的冠状动脉灌注压的无创预测系统,获取冠状动脉灌注压预测值;
控制模块,其被配置为根据所获取的冠状动脉灌注压预测值和冠状动脉灌注压,通过PID控制,调节优化按压参数;
心肺复苏模块,其被配置为根据所得到的按压参数,调节心肺复苏按压深度,完成心肺复苏的自适应调节。
9.如权利要求8中所述的基于无创预测冠状动脉灌注压的心肺复苏系统,在所述控制模块中,将所得到的当前冠状动脉灌注压值与设定的冠状动脉灌注压标准值相减,得到当前误差;若PID控制器的输出为正,则控制心肺复苏仪的按压深度增大;若PID控制器的输出为负则控制按压深度减小,输出值越大,增大和减小速度越快,实现按压参数的调节。
10.如权利要求8中所述的基于无创预测冠状动脉灌注压的心肺复苏系统,其特征在于,在所述心肺复苏模块中,根据实时预测的冠状动脉灌注压值与设定冠状动脉灌注压参考值之间的误差,结合PID控制器,自适应调节心肺复苏仪的按压深度,实现预测冠状动脉灌注压值与设定冠状动脉灌注压值的逐步趋近,完成自适应心肺复苏的调节。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310206600.2A CN116421457A (zh) | 2023-02-28 | 2023-02-28 | 一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310206600.2A CN116421457A (zh) | 2023-02-28 | 2023-02-28 | 一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116421457A true CN116421457A (zh) | 2023-07-14 |
Family
ID=87091550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310206600.2A Pending CN116421457A (zh) | 2023-02-28 | 2023-02-28 | 一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116421457A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117137492A (zh) * | 2023-11-01 | 2023-12-01 | 山东大学齐鲁医院 | 冠状动脉血流异常检测系统、存储介质及终端 |
-
2023
- 2023-02-28 CN CN202310206600.2A patent/CN116421457A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117137492A (zh) * | 2023-11-01 | 2023-12-01 | 山东大学齐鲁医院 | 冠状动脉血流异常检测系统、存储介质及终端 |
CN117137492B (zh) * | 2023-11-01 | 2024-02-09 | 山东大学齐鲁医院 | 冠状动脉血流异常检测系统、存储介质及终端 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20180146929A1 (en) | Device for predicting ventricular arrhythmia and method therefor | |
EP2544124A1 (en) | Methods and systems for non-invasive measurement of glucose levels | |
CN105473065B (zh) | 确定在cpr期间的自主循环恢复 | |
US20160256063A1 (en) | Analyte assessment and arrhythmia risk prediction using physiological electrical data | |
EP2544125A1 (en) | Methods and systems for non-invasive measurement of blood pressure | |
CN116421457A (zh) | 一种基于无创预测冠状动脉灌注压的心肺复苏系统及方法 | |
EP2512325A2 (en) | Processing physiological sensor data using a physiological model combined with a probabilistic processor | |
Amann et al. | Reduction of CPR artifacts in the ventricular fibrillation ECG by coherent line removal | |
CN110812699A (zh) | 远程除颤防护系统及其控制方法 | |
Rheinberger et al. | Removal of CPR artifacts from the ventricular fibrillation ECG by adaptive regression on lagged reference signals | |
Mateo et al. | Radial basis function neural networks applied to efficient QRST cancellation in atrial fibrillation | |
Mathews et al. | Leveraging a discriminative dictionary learning algorithm for single-lead ECG classification | |
CN111832537B (zh) | 异常心电信号识别方法及异常心电信号识别装置 | |
CN113967022A (zh) | 一种基于个体自适应的运动想象脑电特征表征方法 | |
Aamir et al. | Automatic Heart Disease Detection by Classification of Ventricular Arrhythmias on ECG Using Machine Learning. | |
Alonso et al. | A machine learning framework for pulse detection during out-of-hospital cardiac arrest | |
CN110974213A (zh) | 基于深度堆栈网络的心电信号识别方法 | |
CN116236162A (zh) | 一种睡眠分期方法、装置、设备及存储介质 | |
Ming et al. | Detection of shockable rhythm during chest compression based on machine learning | |
CN115670902B (zh) | 一种心肺复苏的智能决策系统及方法 | |
US11589824B2 (en) | System and method for estimating the brain blood volume and/or brain blood flow and/or depth of anesthesia of a patient | |
CN114082169A (zh) | 基于脑电信号的伤残手软体康复机器人运动想象识别方法 | |
EP3014501B1 (en) | Measuring respiration | |
Saxena et al. | Extraction of various features of ECG signal | |
Valencia et al. | Multiscale regularity analysis of the Heart Rate Variability: stratification of cardiac death risk |
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 |