一种基于光电容积脉搏波实现有创血压波形估计的方法

1.本发明涉及医疗检测技术领域，特别是涉及一种基于光电容积脉搏波实现有创血压波形估计的方法。


背景技术：

2.对血压参数的长期监测有助于改善预后效果，从而降低心血管疾病cvd的死亡率。因此，能够长时间无创监测血压的设备对于cvd患者的诊断治疗非常重要。然而，传统的袖带式血压测量设备尺寸较大、不便携带，也不适用于长期监测。许多研究为了实现便携、舒适和长时间的血压监测，提出了无袖带血压测量技术及设备，其大多使用光学原理测量光电容积脉搏波(photoplethysmography，ppg)信号，从中可导出血氧、心率、血压等生理参数。这种方法能够获取连续动态的血压参数，有效弥补传统袖带血压测量不便、间断测量等不足。许多已有方法结合心电图(electrocardiogram，ecg)提取脉搏传导时间(pulse transienttime，ptt)，即脉搏波在血管中传播的时间间隔，进而基于脉搏波速度(pulse wave velocity，pwv)理论估算血压。但是，迄今为止，很少有成熟的基于pwv理论的无袖带血压监测设备。目前仍未充分证明ptt能够准确估计血压，并且需要频繁使用袖带血压计进行校准。最后，ecg测量需要电极或双手接触以形成闭环通路，在易用性上较差，同时也使其难以实现全天候的连续血压监测。
3.上述基于pwv理论估算血压的方法存在一定的问题，与之相比，从ppg波形中提取形态学参数以建立血压模型可能更容易学习到更深层次的生理关联性。很多研究也改进了信号预处理等方法，从而有效抑制生理信号中常见的基线漂移和异常噪声的干扰。采用深度学习方法直接从原始波形中提取深度特征，从而避免了人工定义特征或浅层学习算法在拟合血压方面的不足。此外，考虑到原始波形中存在更多描述心血管功能的有用信息，人工定义的形态学参数有可能遗漏这些特征，序列学习的方法通过对完整信号进行互相关等运算，因此更具准确估算血压的潜力。


技术实现要素：

4.本发明的目的是提供一种基于光电容积脉搏波实现有创血压波形估计的方法，通过同时优化深度特征和整合原始ppg的形态学特征，为血压波形估计制定一个具有知识蒸馏(knowledge distillation，kd)策略的新型informer架构(kd-informer)，提供了连续的血压波形预测，更能反映动脉血管内血压的连续变化，对心血管事件预测、靶器官损伤具有更高的预测价值。
5.为实现上述目的，本发明提供了如下方案：
6.一种基于光电容积脉搏波实现有创血压波形估计的方法，包括如下步骤：
7.s1，采集监测人体的光电容积脉搏波信号，同时收集连续血压波形作为参考血压信号；
8.s2，对步骤s1中获得的光电容积脉搏波信号及参考血压信号片段进行数据预处
理，得到训练数据；
9.s3，通过步骤s2数据预处理获取的训练数据，计算得到大量的ppg形态学参数，然后采用befs算法筛选具有鉴别意义的信息参数作为深度学习模型的先验信息；
10.s4，搭建用于无袖带血压波形连续估计监测的深度学习模型，得到具有kd策略的新型informer架构kd-informer，kd-informer采用informer作为学习共享高级特征的骨干网络，采用带有压缩和激发模块的resnet结构用于融合先验信息，基于响应的kd策略，将教师informer学到的知识转移到轻量级学生informer上；
11.s5，将步骤s2中数据预处理获取的训练数据按照一定比例划分训练集、验证集和测试集，其中，利用训练集和验证集对步骤s4建立的深度学习模型进行训练和参数选择，测试集用于检验经训练和参数选择后得到的深度学习模型的泛化能力；
12.s6，采集用户的待测光电容积脉搏波信号数据，输入步骤s5得到的深度学习模型中，输出用户的预测血压波形。
13.进一步的，所述步骤s2中，所述数据预处理包括重采样、滤波降噪、信号分割、时间对齐以及先验特征提取操作。
14.进一步的，所述步骤s3中，采用befs算法筛选具有鉴别意义的信息参数作为深度学习模型的先验信息，具体包括：
15.采用wrapper和embedded方法分别过滤出具有前k位的特征子集，用于估计收缩压和舒张压，并将合并后的集合作为基线集；
16.然后，在特征参数的后向消除阶段，依次消除重要性最低的特征，将剩余的特征集输入到回归器以拟合血压值，并根据优化后的回归结果考虑最终的特征集，作为深度学习模型的先验信息。
17.进一步的，所述步骤s5中，按照7:1.5:1.5的比例划分训练集、验证集和测试集。
18.根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明提供的基于光电容积脉搏波实现有创血压波形估计的方法，充分利用transformer模型在时间依赖性分析方面的优势，应用transformer作为网络骨干，结合从ppg信号中提取的形态学特征，并使用具有挤压和激发(squeeze-and-excitation)模块的resnet结构(se-resnet)将其与深度特征融合；引入的先验信息可以促进模型在学习过程中利用更多的鉴别性特征，从而达到更高的预测精度；同时，为了减少内存消耗，实现轻量级计算，本发明采用了基于响应的kd策略，将教师模型学到的知识转移到轻量级的学生模型参数中；此外，本发明通过修正chimerge算法、后向消除特征选择(backward elimination feature selection，befs)和特征融合方法，不仅可以获得有用特征的稀疏子集，而且在学习过程中，神经网络捕获的序列特征也可以被有效地结合起来。
附图说明
19.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
20.图1是本发明基于光电容积脉搏波实现有创血压波形估计的方法的流程图；
21.图2是本发明基于光电容积脉搏波实现有创血压波形估计的方法的具体的实验设计流程图；
22.图3是本发明数据处理中设计的形态学特征的直观展示图；
23.图4是本发明所述深度学习模型的整体网络框架；
24.图5是本发明输入信号示意图；
25.图6是本发明输出结果示意图。
具体实施方式
26.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
27.本发明的目的是提供一种基于光电容积脉搏波实现有创血压波形估计的方法，通过同时优化深度特征和整合原始ppg的形态学特征，为血压波形估计制定一个具有知识蒸馏(knowledge distillation，kd)策略的新型informer架构(kd-informer)，提供了连续的血压波形预测，更能反映动脉血管内血压的连续变化，对心血管事件预测、靶器官损伤具有更高的预测价值。
28.为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
29.如图1所示，本发明提供的基于光电容积脉搏波实现有创血压波形估计的方法，包括如下步骤：
30.s1，采集监测人体的光电容积脉搏波信号，同时收集连续血压波形作为参考血压信号；本发明方法采用mimic(重症监护的多参数智能监测)数据集，包含了20000多条从波士顿贝斯以色列医院内科、外科和重症监护室的病人监护仪上记录的生理数据，每条记录通常包含24至48小时的连续数据；
31.s2，对步骤s1中获得的光电容积脉搏波信号及参考血压信号片段进行数据预处理，得到训练数据；所述数据预处理主要包括重采样、滤波降噪、信号分割、时间对齐以及先验特征提取等操作；
32.s3，通过步骤s2数据预处理获取的训练数据，计算得到大量的ppg形态学参数，然后采用befs算法筛选具有鉴别意义的信息参数作为深度学习模型的先验信息；befs算法分为基线集生成和特征参数的后向消除阶段，wrapper和embedded方法分别过滤出具有前k位的特征子集，用于估计收缩压和舒张压，并将合并后的集合作为基线集；然后，在特征参数的后向消除阶段，依次消除重要性最低的特征，将剩余的特征集输入到回归器以拟合血压值，并根据优化后的回归结果考虑最终的特征集数量；此外，该步骤也考虑了先验信息与神经网络提取的深度信息进行特征融合的方法设计；
33.s4，搭建用于无袖带血压波形连续估计监测的深度学习模型，如图2所示，得到具有kd策略的新型informer架构kd-informer，kd-informer采用informer作为学习共享高级特征的骨干网络，采用带有压缩和激发模块的resnet结构用于融合先验信息，基于响应的kd策略，将教师informer学到的知识转移到轻量级学生informer上；其中，informer是主要
框架，血压波形预测主要关注不同模态生理序列之间的映射关系，因此删除了原始网络中掩码的稀疏概率多头自注意力(multi-head self-attention，mhsa)结构，以避免序列信息泄露；
34.s5，将步骤s2中数据预处理获取的训练数据按照一定比例划分训练集、验证集和测试集，其中，利用训练集和验证集对步骤s4建立的深度学习模型进行训练和参数选择，测试集用于检验经训练和参数选择后得到的深度学习模型的泛化能力；通过不断向系统模型批量送入训练数据并且通过网络训练来更新网络参数，重复上述操作n次后，将得到的参数视为最优参数；
35.s6，采集用户的待测光电容积脉搏波信号数据，输入步骤s5得到的深度学习模型中，输出用户的预测血压波形。
36.其中，基于ppg信号预测血压波形可以被概括为一个实时长序列预测问题。因此，该模型需要通过固定窗长的时滞窗口学习并获取两个序列之间的映射关系。在发明方法中，对于时间窗口t的输入数据ppg序列(即ppg波形)，输出是相应的血压序列(即血压波形)。与单变量预测问题不同，该模型的特征维度不限于单变量的情况。可以选择ppg序列的导数，即ppg的速度序列，ppg的加速度序列。建议的模型接受输入和输出，由公式(1)给出。
37.y＝f(x；m；θ)
ꢀꢀꢀ
(1)
38.其中，f(
·
)是模型函数，代表引入先验知识，为可选参数，比如ppg波形的形态学特征；θ则代表深度学习模型的超参数。
39.图2说明了本发明方法具体的实验设计流程图。数据库中收集的ppg和血压信号长度不一，并受到离群值和基线漂移的干扰，因此需要对信号进行预处理，如图3中的方框所示。本发明方法将信号统一重采样到125hz，这有助于序列的时间对齐。此外，由于小波变换适合分析非平稳信号，在信号突变、压缩重建和信号去噪中具有更好的时频定位。因此，选择sym4小波为基础小波的小波变换进行信号去噪，并根据软阈值函数对噪声的ppg信号进行两阶段独立分解。然后，使用固定尺寸为8秒、滑动步长为3秒的时滞窗口对滤波后的信号进行分割。最后，本发明方法采用最大最小归一化，以确保模型能够快速收敛。mimic数据集由80454个8秒片段组成，平均收缩压为140.21
±
28.97mmhg，舒张压为68.62
±
15.60mmhg。
40.最后，本发明方法从归一化的ppg信号中提取形态学特征，以为模型决策提供先验信息。本发明方法分别从ppg信号中提取了75个可解释的特征，包括时域、频域和非线性特征，如时间参数信号偏度(公式2)，无维度指标边缘因子(公式3)，以及面积参数k值(公式4)。大部分的特征展示在图4中。
[0041][0042]
[0043][0044]
其中，p
max
，p
min
和p
mean
分别代表周期内ppg振幅的最大值、最小值和平均值。
[0045]
此外，本发明方法额外使用了脉率变异性(pulse rate variability，prv)来描述序列窗口内的变化模式，如描述序列无序程度的近似熵(approximate entropy，apen)。对于ppg序列x，使用长度为m＝2的窗口截断来获得嵌入的序列w＝{w(i),i＝1,2,...,n m-1}，相似性容忍度为r＝0.1
×
sd，而apen定义为。
[0046][0047]
其中，表示在相似性标准r下，整个序列所有长度为m的子片段的平均相似性率，计算公式如下：
[0048][0049]
其中，w(j)满足max(|w(i)-w(j)|)≤r)/(n-m 1)。
[0050]
此外，本发明方法对计算出的连续特征进行监督下的chimerge分箱离散化，具有一下优点：1)离散特征对异常数据具有较好的稳定性；2)稀疏向量内积乘法操作比较高效；3)离散特征有助于简化模型，降低过拟合的风险。chimerge依赖于卡方检验，将相邻的卡方值较低的区间(表明它们具有相似的类分布)结合起来，直到满足停止准则。计算方法如公式(7)所示，a
ij
表示第i个区间、第j个类别中的实例的数量，相应地，e
ij
表示a
ij
的预期频率。
[0051][0052]
连续特征集可以表示为其中称为时间t的s维特征集集合，代表第t个特征集的平均血压值。本发明方法做进一步改进，首先将血压值等距地分割；然后计算相应类别的值，合并最小区间，从而实现特征离散化。
[0053]
通过步骤s2数据预处理阶段获取的训练数据计算了大量的ppg形态学参数，然而，过多的特征集合可能包含冗余的信息，仍需要特征选择和融合的操作，以使深度神经网络更好地学习形态相关的先验信息。此外，该步骤也考虑了先验信息与神经网络提取的深度信息进行特征融合的方法设计。本发明方法设计了befs算法用于筛选具有鉴别意义的特征作为深度模型的先验信息。
[0054]
在得到优化的特征集后，本发明方法考虑设计了特征融合策略，将kd-informer编码器的深度特征与传统的形态学特征相结合。经典的融合方法是通过设置特征权重来确定不同特征集的比例，其计算方法如下：
[0055]ff
＝ε
·fm
(1-ε)
·
fdꢀꢀꢀ
(8)
[0056]
其中，ff代表融合特征，fm和fd分别表示形态学特征和深度特征。然而，虽然这种线性融合方法部署简单，但超参数很难根据数据分布自动调整，需要进行大量的实验来选择合适的比例因子。为此，本发明设计了一种新的特征融合方法，可以自适应地调整特征比例
因子。这种方法通过训练多层感知器来更新特征权重，融合操作定义如下。
[0057][0058]
其中，b是偏置值。全连接层的非线性空间映射能够更好地权衡先验知识和来自深层特征的信息。
[0059]
考虑到深度学习模型较高的计算复杂性和较大的存储需求，在步骤s4中，本发明方法采用了基于响应的kd策略，将教师informer学到的知识转移到轻量级的学生informer身上。这种方法的主要思想是在教师模型的指导下训练一个高效的学生模型，以获得相当的准确性。从预训练的教师模型中学到的暗知识可以帮助学生模型模仿教师模型的行为，具体的提炼过程如图4所示。教师模型的软几率(logits)包含了特定类别中的有用信息，用于softmax回归任务。温度参数t可以控制分类任务中每个标签的重要性，例如，较高的温度在特定类别zi上产生较弱的概率分布pi(如公式10所示)。然而，血压预测任务的输出是一个时间序列，而不是分类任务中的概率分布或类别单值。在回归任务中，没有这样的软几率(logits)；因此，本发明方法没有使用温度参数，转而使用教师模型的预测输出作为软标签来指导学生的训练。kd-informer损失可以定义为公式(11)。第一项是蒸馏损失，衡量教师的软预测和学生的预测之间的距离；第二项是学生损失，它衡量学生的预测和参考血压序列之间的距离。超参数用于调整两个损失函数之间的贡献程度。
[0060][0061][0062]
步骤s5对步骤s2预处理之后的数据按照7:1.5:1.5的比例划分训练集、验证集和测试集，其中训练集和验证集对步骤4建立的系统模型进行训练和参数选择，测试集用于最终检验步骤s5保存下来的最优系统模型的泛化能力。数据切分采用个体间分割的方式，确保每组数据不包含同一患者的数据，从而避免了信息泄露。在训练网络的过程中，使用adam优化器来进行参数的更新，学习率为0.001，权值衰减为0.999，动量为0.8。
[0063]
步骤s6中采集用户ppg信号(如图5所示)，将信号波形作为输入深度学习模型中，可自动输出用户的血压信号，如图6所示，实现基于光电容积脉搏波实现有创血压波形的无创估计。
[0064]
本发明提供的基于光电容积脉搏波实现有创血压波形估计的方法，(1)输入仅由脉搏波组成，因此采集电路只需要采集ppg信号，与传统的基于脉搏波传播速度的方法相比，省去了采集心电信号这一步骤，也无需过多的求导计算，这样方便集成到手环等设备中，不需要袖带等血压测量装备，摆脱了袖带的束缚，使得设备更加便携；(2)本发明的血压估计算法能够实现连续的血压波形估计，也可以实现血压的长期监测，并且可用于日常生活中血压的测量，测量时也不会给人体带来创伤和不适的影响：(3)本发明应用transformer作为网络骨干提取ppg信号中的大量信息，结合从ppg信号中提取的形态学特征，引入这些先验信息可以促进模型在学习过程中利用更多的鉴别性特征，输入信号包含的信息更丰富，从而使得测量出的血压结果更加稳定，达到更高的预测精度；(4)本发明为
了减少内存消耗，实现轻量级计算，采用知识蒸馏策略，将教师模型学到的知识转移到轻量级的学生模型参数中，嵌入到实际的测量设备中，可以大大降低设备的成本；(5)本发明方法设计了详细的特征选择策略，通过修正chimerge算法、后向消除特征选择(backward elimination feature selection,befs)和特征融合方法，不仅可以获得有用特征的稀疏子集，而且在学习训练的过程中，神经网络捕获的序列特征也可以被有效地结合起来。
[0065]
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。