CN111065323A - 用于确定对象的血压的装置和方法 - Google Patents

用于确定对象的血压的装置和方法 Download PDF

Info

Publication number
CN111065323A
CN111065323A CN201880048936.8A CN201880048936A CN111065323A CN 111065323 A CN111065323 A CN 111065323A CN 201880048936 A CN201880048936 A CN 201880048936A CN 111065323 A CN111065323 A CN 111065323A
Authority
CN
China
Prior art keywords
blood pressure
calibration
subject
features
unit
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
CN201880048936.8A
Other languages
English (en)
Other versions
CN111065323B (zh
Inventor
L·施米特
J·米尔施泰夫
E·布雷施
X·于
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of CN111065323A publication Critical patent/CN111065323A/zh
Application granted granted Critical
Publication of CN111065323B publication Critical patent/CN111065323B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/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/021—Measuring pressure in heart or blood vessels
    • A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
    • A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
    • 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/7221—Determining signal validity, reliability or quality
    • 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/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
    • 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
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/04—Measuring blood pressure
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/58—Testing, adjusting or calibrating the diagnostic device
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2560/00—Constructional details of operational features of apparatus; Accessories for medical measuring apparatus
    • A61B2560/02—Operational features
    • A61B2560/0223—Operational features of calibration, e.g. protocols for calibrating sensors
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
    • A61B2562/02—Details of sensors specially adapted for in-vivo measurements
    • A61B2562/0219—Inertial sensors, e.g. accelerometers, gyroscopes, tilt switches
    • 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/024—Detecting, measuring or recording pulse rate or heart rate
    • A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
    • 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/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Cardiology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Psychiatry (AREA)
  • Artificial Intelligence (AREA)
  • Signal Processing (AREA)
  • Vascular Medicine (AREA)
  • Hematology (AREA)
  • Acoustics & Sound (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明涉及用于确定对象的血压的装置和方法。为了自动触发校准，所述装置包括：传感器信号输入端(31)，其被配置为获得所述对象的动脉脉搏波传感器信号(11)；特征提取单元(32)，其被配置为从所获得的动脉脉搏波传感器信号中提取多个特征(42)；估计单元(33)，其被配置为针对个体提取的特征和/或提取的特征的组确定多个血压估计值(43)，并且根据所述多个血压估计值来确定所述对象的血压(44)；校准单元(34)，其被配置为基于血压参考测量结果(21)来校准所述估计单元(33)；以及校准触发单元(35)，其被配置为在所述多个血压估计值分散超过分散度极限时触发通过所述校准单元(34)进行的校准。

Description

用于确定对象的血压的装置和方法
技术领域
本发明涉及用于确定对象的血压的装置和方法。另外，本发明涉及用于监测对象的血压的装置。
背景技术
血压(BP)测量在所有医院环境中都是常见的，以评价并监测患者的血液动力学状况。对于不稳定的患者，需要目前主要使用有创动脉内压力导管的连续监测。该方法需要经训练的医务人员来应用，并且与患者的感染风险相关联。然而，如果患者的血液动力学稳定性允许，那么无创测量方法(听诊或示波)被应用。这些方法涉及通常被施加在上臂处的袖带。该方法不允许连续地监测血压，并且对于患者来说是相当不舒适的，特别是当在夜间进行利用例如15分钟间隔的自动化抽查测量时。
另外的方法允许连续的但无干扰的血压测量，优选地不将任何外力或压力施加于患者。这些方法通常基于能够被连续地且无创地感测的生理参数。根据周期性波形数据，计算BP替代参数，其允许在连续的或心跳到心跳的基础上导出BP信息。用于BP替代参数的范例是脉搏到达时间(PAT)，其能够在逐个心跳的基础上根据连续的心电图(ECG)和光电体积描记(PPG)波形数据来计算。PAT表示测量位置处的动脉脉搏的到达时间。随着BP增加/降低，脉搏波更快地/更慢地行进通过动脉树，因此导致降低的/增加的PAT。原理上，从电、光学和加速度测量传感器获得的光电体积描记、心电图和心震图波形信号包含丰富的生理信息，能够从其中提取不同的特征用于确定BP的替代。
US2012/0136261A1公开了用于校准生理参数的计算的系统和方法。两种或更多种校准技术可以用来确定生理测量结果与期望的生理参数之间的关系，诸如差分脉搏传播时间(DPTT)与血压之间的关系。不同的校准技术可以以串联方式(一个接着另一个)或以并联方式被使用，具有根据每种校准技术的不同的权重。当生理或其他变化发生时，校准数据可以被存储用于稍后使用，并且新的校准数据可以被生成。
US2017/0042433A1公开了一种血压估计方法，包括测量包括用户的脉搏波信息的生物信号，确定用于血压估计模型的校准方法，使用所确定的校准方法来校准血压估计模型，并且使用经校准的血压估计模型根据生物信号来估计用户的血压。
仍然需要对患者间以及患者内变化性鲁棒的特征选择和/或自动检测何时需要(重新)校准。
发明内容
本发明的目的是提供用于确定对象的血压的装置和方法，其自动检测何时需要(重新)校准。
在本发明的第一方面中，呈现了一种用于确定对象的血压的装置，包括：
传感器信号输入端，其被配置为获得所述对象的动脉脉搏波传感器信号，
特征提取单元，其被配置为从所获得的动脉脉搏波传感器信号中提取多个特征，
估计单元，其被配置为针对个体提取的特征和/或提取的特征的组确定多个血压估计值，并且根据所述多个血压估计值来确定所述对象的血压，
校准单元，其被配置为基于血压参考测量结果来校准所述估计单元，以及
校准触发单元，其被配置为在所述多个血压估计值分散超过分散度极限时触发通过所述校准单元进行的校准。
在本发明的又一方面中，呈现了一种用于监测对象的血压的装置，包括：
动脉脉搏波传感器，其被配置为采集所述对象的动脉脉搏波传感器信号，
血压参考测量单元，其被配置为采集所述对象的血压参考测量结果，以及
本文中公开的装置，其用于基于所采集的动脉脉搏波传感器信号和所采集的血压参考测量结果来确定对象的血压。
在本发明的另外的其他方面中，提供了对应的方法、计算机程序以及非瞬态计算机可读记录介质，所述计算机程序包括程序代码单元，当所述计算机程序在计算机上被执行时，所述程序代码单元用于使所述计算机执行本文中公开的方法的步骤，所述非瞬态计算机可读记录介质在其中存储了计算机程序产品，当被处理器执行时，所述计算机程序产品使本文中公开的方法被执行。
在从属权利要求中限定了本发明的优选实施例。应当理解，要求保护的方法、装置、计算机程序和介质具有特别是如在从属权利要求中限定的和如在本文中公开的类似和/或相同的优选实施例。
本发明基于在用于确定血压的选定特征失去其与血压的关系时利用触发器用于触发(重新)校准(在本文中一般被称为“校准”)的构思。这种情况被自动检测。确定对于针对特征的集合的校准的需要的问题因此通过使用针对个体提取的特征和/或提取的特征的组的多个血压估计值来解决，根据其对象的血压被确定，并且在多个血压估计值分散超过分散度极限时触发校准。
因此，本发明不确定通过将最近的血压测量结果与一个或多个先前的血压测量结果(即在不同时刻获得的血压测量结果)进行比较确定的血压变化是否超过特定阈值(如根据US2012/0136261A1完成的)。另外，本发明不确定从由不同位置处的多个传感器测量的多个生物信号选择的参考生物信号(例如PPG信号)与被预先存储在模板中的生物信号之间的相似度的程度是否小于或等于预定的阈值(如根据US2017/0042433A1完成的)。替代地，本发明确定两个或更多个血压估计值，其中，每个血压估计值针对从(相同的)所获得的动脉脉搏波传感器信号(例如PPG信号)中提取的不同特征或特征的不同组来确定。
例如，如在实施例中提取的，引入特征分散度度量，其测量跨选定特征的集合的BP估计之间的差，由此提供检测何时需要校准的更直接且可靠的方式。所述校准触发单元因此可以被配置为确定特征分散度度量，所述特征分散度度量表示针对所述多个血压估计值之间的差的度量。分散度度量然后可以与分散度阈值(或极限)进行比较以检测是否需要校准。
在实际的实施例中，所述校准触发单元被配置为在预定数量的或所有所述多个血压估计值的差的平均差、最大差或标准偏差超过分散度阈值时触发校准。
由此，所述校准触发单元优选地被配置为使用预定的或时间变化的分散度阈值，特别地其中，所述时间变化的分散度阈值自上一次校准以后随着时间被减小。这确保在一些时间点处校准将会被进行，并且对于校准的需要不被完全抑制。
在一实施方式中，所述校准单元被配置为通过以下操作来确定用于通过所述估计单元针对个体提取的特征和/或提取的特征的组确定所述多个血压估计值的校准值：将提取的特征或提取的特征的组与血压之间的建模的关系拟合到实际提取的特征值或提取的特征值的组和所述血压参考测量结果，所述拟合特别是经由线性或非线性回归进行的。
在一简单实施例中，所述估计单元被配置为通过获取所述多个血压估计值中的一些或全部的平均值或中值来确定所述对象的血压。
优选地，所述估计单元被配置为通过获取所述多个血压估计值中的一些或全部的加权平均值来确定所述对象的血压，其中，用于所述加权平均值的权重是基于对象的血压与其对应特征之间的相关性来确定的，特别地其中，对象的血压与其对应特征之间的所述相关性是基于所获得的血压参考测量结果来确定的。由于真实血压不是已知的，所以血压参考测量结果可以用来评价相关性。示出与血压的良好相关性的特征可以被选择，并且对于每个选定的特征，血压估计值能够被计算(例如经由根据数学模型的回归，再次使用血压参考测量结果)。如果有多个特征被选择，那么多个血压估计值可以被组合成单个血压估计值，这能够经由加权平均或其他手段来完成。
在一备选实施例中，所述估计单元被配置为通过获取所述多个血压估计值中的一些或全部的加权平均值来确定所述对象的血压，其中，用于所述加权平均值的权重通过选择其对应特征具有与所述对象的血压的最好相关性的血压估计值、或依据所述血压估计值的对应特征的信噪比、或基于回归误差来确定。
现有的系统和方法的另一问题是缺少用于对患者间以及患者内变化性鲁棒的特征选择的方法。该问题通过以下实施例来解决，根据该实施例，所述校准单元被配置为执行所述估计单元针对所有提取的特征的初始校准，并且基于对象的血压与对应特征之间的相关性来选择提取的特征的子集，并且所述估计单元被配置为针对所述子集的个体提取的特征和/或提取的特征的组确定所述多个血压估计值。以此方式，提供了自动特征选择，其中对于个体患者，源于大量可能特征的特征的子集被自动地且单独地选择。优选地，选定特征的集合在每次校准期间被自动更新。
由此，所述校准单元可以被配置为在针对一个或多个特征与所述对象的血压的所述相关性的一个或多个相关性值低于相关性阈值时重复所述估计单元针对一些或所有提取的特征的所述初始校准。备选地，所述校准单元可以被配置为在针对一个或多个特征与所述对象的血压的所述相关性的一个或多个相关性值低于相关性阈值时重复所述估计单元针对一些或所有提取的特征的所述校准。
本发明还提供了一种用于监测对象的血压的装置，包括动脉脉搏波传感器、血压参考测量单元、如本文中公开的用于确定对象的血压的装置。所述动脉脉搏波传感器可以包括加速度计、光电体积描记传感器、超声传感器、雷达传感器和生命体征相机中的一个或多个。
附图说明
本发明的这些和其他方面将参考下文描述的(一个或多个)实施例变得显而易见并将参考下文描述的(一个或多个)实施例得以阐述。在以下附图中，
图1示出了根据本发明的用于监测血压的装置和用于确定血压的装置的实施例，
图2示出了指示特性点的PPG波形信号的示意图，
图3示出了根据本发明的方法的第一实施例的流程图，
图4示出了根据本发明的方法的第二实施例的流程图，
图5示出了根据本发明的方法的第三实施例的流程图，
图6示出了时间变化的分散度阈值的示意图，
图7示出了根据本发明的方法的第四实施例的流程图，以及
图8示出了根据本发明的方法的第五实施例的流程图。
具体实施方式
图1示出了根据本发明的用于监测血压的装置1和用于确定对象(例如患者)的血压的装置30的实施例。用于监测血压的装置1包括动脉脉搏波传感器10，其被配置为采集对象的动脉脉搏波传感器信号11(或多个动脉脉搏波传感器信号)。动脉脉搏波传感器10包括加速度计、光电体积描记(PPG)传感器、超声传感器、雷达传感器和生命体征相机中的一个或多个。装置1还包括血压参考测量单元20，其被配置为采集对象的血压参考测量结果21。血压参考测量单元20可以例如是常规的无创BP(例如NIBP)测量设备。装置1还包括用于基于所采集的动脉脉搏波传感器信号11和所采集的血压参考测量结果21来确定对象的血压的装置30。装置30可以以软件和/或硬件方式被实施，例如被实施为适当编码的处理器、CPU或计算机。
如果需要，可以提供用于附接或定位传感器10和/或单元20的装置，诸如带、粘附物等。
用于确定对象的血压的装置30包括传感器信号输入端31，其被配置为获得对象的动脉脉搏波传感器信号11。传感器信号输入端31可以例如是用于检索或接收直接来自动脉脉搏波传感器10或来自存储设备或缓冲器的动脉脉搏波传感器信号11的无线或有线接口(例如USB接口、WLAN接口、网络接口等)。
装置30还包括特征提取单元32，其被配置为从所获得的动脉脉搏波传感器信号11中提取多个特征42。一般存在可以被使用的大量特征，例如所获得的动脉脉搏波传感器信号11的特性点之间的时间、特性点的幅度等。
装置30还包括估计单元33，其被配置为针对个体提取的特征42和/或提取的特征42的组确定多个血压估计值43，并且根据所述多个血压估计值43来确定对象的血压44。校准单元34基于提取的特征42和血压参考测量结果21来校准估计单元33，提取的特征和血压参考测量结果可以被存储在装置30中或例如直接或经由输入单元31从血压参考测量单元20或从存储设备或缓冲器获得(检索或接收)。如果多个血压估计值43分散超过分散度极限，那么校准触发单元35(例如通过触发信号45)触发通过校准单元34进行的校准，所述分散度极限可以是预定的或时间变化的阈值，并且所述分散度极限可以被存储在装置30中或从外部实体获得。
图2示出了PPG波形信号的示意图，表示范例动脉脉搏波传感器信号11，指示特性点(特征)。在该范例中，这些点包括PPG脉搏的最小值(min)、PPG脉搏的最大值(max)、一阶导数的最大值(deri)、重搏切迹的最小值和最大值(a和b)。在下表中示出了能够从PPG波形信号中提取的特征的范例，这些能够在文献中发现：
Figure BDA0002376982760000071
Figure BDA0002376982760000081
(a)时域特征
Figure BDA0002376982760000082
(b)频域特征
Figure BDA0002376982760000083
(c)小波域特征
由于患者间变化性，具体的BP替代参数具有与一些患者中的BP的更强的相关性和与其他患者中的BP更少的相关性。而且，替代参数需要被校准到针对个体患者的血压以便允许定量估计。为了给出范例：患者的身高确定动脉树的长度，并且因此是当使用PAT作为替代参数时贡献于患者间变化性的许多因素之中的一个因素。
用于校准的标准方法如下：定义用于将替代参数与BP相关的数学模型，包含一个或多个未知的校准参数。多对替代参数值和参考BP值(通常利用袖带设备)被测量。校准常数通过将模型拟合到BP替代参数对(回归)来确定。此后，血压能够根据连续测量的替代参数经由定义的数学模型通过利用所确定的校准常数来估计。
校准过程需要被定期地重复(重新校准)以考虑患者内变化性。作为一范例，血管顺应性能够由于某些生理事件或由于药物而改变，使得校准参数需要被更新。此外，患者的血管状态的改变甚至能够影响具体替代参数与BP的相关性，使得重新校准仅部分地考虑患者内互操作性。
因此，对于使用BP替代，存在两个主要挑战，即总是识别与BP很好地相关的良好替代参数，并且对于所有患者，检测经校准的替代参数何时需要重新校准。关于后一挑战的挑战是当BP估计性能下降时需要重新校准。然而，对BP估计性能的直接评价是不可能的，因为参考BP测量结果不可获得。
在一实施例中，根据本发明使用的校准方法包括，获得并存储来自NIBP设备的一个或多个参考血压测量结果，评价每个特征与血压之间的关系(例如利用p值的相关性系数)，基于存储的参考血压测量结果，选择最强特征的子集(与血压具有最好关系的那些)，以及(例如经由线性或非线性回归)基于存储的血压参考测量结果，确定针对每个特征的校准参数。
根据依据本发明的血压估计方法的实施例，对于每个选定的特征，经由函数关系(例如线性模型)通过使用校准参数来估计血压，并且每个选定的特征的血压估计被组合成单个血压估计(替代)，例如通过获取中值或借助于平均。
根据用于自动触发重新校准(即确定是否以及何时重新校准)的方法的实施例，在选定的特征失去其与血压的关系时重新校准被触发。这种情况可以通过将特征分散度度量(FDM)与阈值(其可选地是时间变化的)进行比较来自动检测。FDM表示针对来自选定的个体特征的BP估计之间的差的度量(例如最大差或标准偏差)。
用于重新校准的方法的实施例使用上面描述的校准方法的步骤：选定的特征的集合和其校准参数被动态地更新。用于重新校准的更老的存储的参考血压测量结果的影响可以被降低(例如通过从存储设备移除它们或通过分配适当的加权因子)。
图3示出了根据本发明的方法100的第一实施例的流程图。该实施例涉及用于连续的BP估计的多个特征的自动重新校准。
在第一步骤101中，选择N个PPG信号特征。特征的选择是基于能够被应用于PPG波形信号并且提供特征值的函数。作为一范例，来自上面示出的表的特征time_deri2a(从第一导数的最大值到重搏切迹的最小值的时间)可以被使用。对于每一个脉搏搏动，通过提取PPG波形信号中的与第一导数的最大值到重搏切迹的最小值相对应的时间差来获得特征值。
作为备选选项，除了PPG信号之外的不同生理信号(例如ECG信号)可以被使用。而且，多个生理信号的组合能够被使用。例如，特征PAT从两个不同的生理信号(心电图和光电体积描记图)中提取。此外，对于每个特征，定义根据若干未知的校准参数将特征与BP相关的数学模型。这些能够是线性或非线性模型。例如，线性模型涉及两个校准参数，一个用于斜率，另一个用于偏移。N个特征的有意义的预先选择可以基于已知一般与BP很好相关的特征的公知常识。备选地，特征的预先选择可以基于关于对象的专门数据收集研究。在该具体实施例中，对于连续的BP估计，预先选择的特征的集合不针对个体患者进行调整，而是对于所有患者都是相同的，并且在整个监测程序期间保持固定。
在第二步骤102中，进行N个特征的初始校准。出于该目的，利用血压测量设备(通常是无创设备)来获得多个BP参考测量结果。该设备可以例如是袖带设备或另一经校准的NIBP设备。BP参考测量值和在参考测量时的对应特征值被存储。随后，对于每个特征，通过将特征与BP之间的建模的关系拟合到特征值和BP参考测量值的实际测量对来确定校准参数。通常，这经由回归来完成，取决于建模的关系，回归是线性的或非线性的。
在第三步骤103中，执行连续的BP估计。对于每个特征，从所接收的生理信号中提取特征值。随后，根据每个特征值，通过使用与BP的建模的数学关系和所确定的校准参数来估计对应的BP值。最后，根据所有特征估计的BP值被组合成单个估计的BP值。该组合能够通过各种方法来完成。在优选实施例中，根据所有特征估计的BP值的中值被选择。为了给出备选选项的一些范例，组合可以通过平均或通过加权平均来完成。加权因子可以根据相应的特征与BP之间的相关性来选择，其中，所述相关性根据校准程序期间的测量结果来确定。另一备选方案将会是选择该估计的BP值，其对应于具有最好相关性的特征。
在第四步骤104中，检测对于重新校准的需要(“重新校准条件满足？”)。该步骤可以立即、在单个新的估计的BP值已经被确定的先前步骤103之后或期间、或备选地在几个新的估计的BP值已经被确定之后进行。在图4中图示了用于确定对于重新校准的需要的自动程序200的实施例(即步骤104的实施例)。
根据第一检查201(“自上一次重新校准以后逝去的时间<T_max？”)，在最近的重新校准/校准之后逝去的时间已经超过某一阈值T_max时，触发重新校准。根据第二检查202(“与上一次重新校准相比BP的变化<ΔBP_max？”)，在估计的BP与最近的参考测量结果的BP之间的差已经超过某一阈值ΔBP_max时，触发重新校准。这头两个校准是可选的。它们表示用于确保重新校准在某一时间间隔(例如2小时)内进行至少一次或在显著的BP变化发生时进行重新校准的条件。
该程序的必要元素是要检查针对不同特征的个体估计的BP值是否仍然良好吻合。如果不是这种情况，那么它是不是所有特征都正在准确地估计BP的强指示。因此，某一特征分散度度量(FDM)被定义(步骤203)，并且与某一阈值FDM_max进行比较(步骤204)。假如FDM值超过阈值，那么重新校准条件(“FDM<FDM_max(t)？”)满足(步骤205)，并且重新校准被触发(参见图3)；否则(步骤206)，重新校准不满足。作为备选选项，FDM_max能够是时间变化的。自上一次重新校准/校准以后逝去的时间越多，阈值越低。在图5中图示了时间变化的阈值FDM_max(t)的范例，其中在重新校准之后的某一时间t1之后，FDM_max被线性地降低直至时间t2，在此之后它再次被保持恒定。阈值轨迹的各种其他方式是可能的。
在优选实施例中，FDM被定义为来自所有特征的估计的BP值的范围，即来自所有特征的估计的BP值的最大值与来自所有特征的估计的BP值的最小值之间的差。定义FDM的备选选项能够是所有特征的估计的BP值的标准偏差、四分位距、平均绝对差、中值绝对偏差或平均绝对偏差。
在图3中示出的方法的第五步骤105中，执行重新校准。重新校准与步骤2中的初始校准类似地被执行。通常，在重新校准期间进行的BP参考测量结果的数量与初始校准相比是稍微更少的。
最后，该循环是闭合的，并且循环再次在第三步骤103中开始。
在图6中示出的流程图中图示了根据本发明的方法300的另一实施例，其涉及用于连续的BP估计的多个特征的自动重新校准。然而，与图3中示出的第一实施例相比的实质差异是在预选选择大量M个特征(步骤301，类似于图3中示出的第一实施例的步骤101)之后并且在初始重新校准步骤102之后，在步骤302中选择最好的N个特征的子集，并且步骤303中检查(“特征选择成功”)N个特征的特征选择是否已经成功。该元素的目的是针对患者间变化性改善鲁棒性。通过基于初始校准结果选择最好的N个特征的更小子集，特征的集合针对个体患者进行定制。对于第一实施例，不是这样的情况，其中特征的固定集合用于所有患者。
在图7中示出了用于从大量M个特征中选择最好的N个特征的程序400(即在图6中示出的方法300的步骤302和303的实施例)。图6中的其余步骤与在图3中是完全相同的，并且在第一实施例中进行了描述。
作为用于对M个预先选择的特征进行排序的准则，选择与BP的相关性。基于在校准程序期间获得的特征值和BP参考值，在步骤401中确定针对每个特征的关于BP参考测量结果的相关性系数。此后，在步骤402中选择具有最高相关性系数的N个特征。在步骤403中，检查(“最小相关性系数>C_min？”)最弱特征的相关性系数是否不超过某一最小要求(C_min)，例如C_min＝0.7，特征选择的输出被表示不成功(步骤404)。否则(步骤405)，结果被表示成功。如果特征选择的结果是不成功的，那么初始校准能够被重复，使得更多的BP参考测量结果被获得。
通过选择与BP的高相关性作为用于选择特征的准则，呈现特征与BP之间的线性关系。然而，如果具体特征与BP之间的关系是强的但是非线性的，那么相关性系数能够是非常低的。因此，如果对于一个或多个特征来说用于将特征与BP相关的数学模型是非线性的，那么通过将特征值拟合到在校准期间获得的参考BP值获得的回归误差被选择作为用于对M个预先选择的特征进行排序的备选准则。
作为另一备选方案，特征选择过程能够通过利用测量的特征的信噪比来改善。例如，如果测量的特征的信噪比是低的，那么测量的相关性系数或测量的回归误差能够被适当的加权因子惩罚。这防止意外地具有良好相关性因子或低回归误差的噪声特征被选择。特征的信噪比能够例如借助于标准偏差或样本方差来测量。
图8示出了根据本发明的方法500的另一实施例的流程图。该实施例也涉及用于连续的BP估计的多个特征的自动重新校准。然而，与在图6中示出的实施例相比的实质差异是，从大量M个特征中选择最好的N个特征的子集在每个重新校准步骤之后进行。这被称为动态特征选择。该元素的目的是与在图6中示出的实施例相比针对患者内变化性(即在时间上发生的血管系统的生理变化)改善鲁棒性。通过基于每个重新校准结果选择最好的N个特征的更小子集，特征的集合不仅针对个体患者进行定制，而且针对该患者随着时间的生理变化进行定制。对于在图6中示出的实施例，不是这样的情况，其中针对个体患者的最好的N个特征的选择在初始校准之后仅进行一次。在图8中示出的个体步骤与在图6中示出的那些完全相同，并且在上面图示的实施例中进行了描述。
本发明可以在所有临床环境(例如普通病房)中、在媒介监护中以及在不需要动脉线路的用于患者的手术室和重症监护室中使用。
尽管已经在附图和前面的描述中详细说明和描述了本发明，但是这样的说明和描述应当被认为是说明性或示范性的，而非限制性的；本发明不限于所公开的实施例。本领域技术人员通过研究附图、说明书以及权利要求书，在实践请求保护的发明时能够理解并实现对所公开的实施例的其他变型。
在权利要求中，“包括”一词不排除其他元件或步骤，并且词语“一”或“一个”不排除多个。单个元件或其他单元可以实现在权利要求中记载的若干项的功能。尽管某些措施被记载在互不相同的从属权利要求中，但是这并不指示不能有利地使用这些措施的组合。
计算机程序可以被存储/分布在合适的非瞬态介质上，例如与其他硬件一起或作为其他硬件的部分供应的光学存储介质或固态介质，但是也可以被以其他形式分布，例如经由互联网或其他有线或无线的电信系统。
权利要求中的任何附图标记都不应被解释为对范围的限制。

Claims (15)

1.一种用于确定对象的血压的装置，包括：
传感器信号输入端(31)，其被配置为获得所述对象的动脉脉搏波传感器信号(11)，
特征提取单元(32)，其被配置为从所获得的动脉脉搏波传感器信号中提取多个特征(42)，
估计单元(33)，其被配置为针对个体提取的特征和/或提取的特征的组确定多个血压估计值(43)，并且根据所述多个血压估计值来确定所述对象的血压(44)，
校准单元(34)，其被配置为基于血压参考测量结果(21)来校准所述估计单元(33)，以及
校准触发单元(35)，其被配置为在所述多个血压估计值分散超过分散度极限时触发通过所述校准单元(34)进行的校准。
2.根据权利要求1所述的装置，
其中，所述校准触发单元(35)被配置为确定特征分散度度量，所述特征分散度度量表示针对所述多个血压估计值之间的差的度量。
3.根据权利要求1所述的装置，
其中，所述校准触发单元(35)被配置为在预定数量的或所有所述多个血压估计值的差的平均差、最大差或标准偏差超过分散度阈值时触发校准。
4.根据权利要求3所述的装置，
其中，所述校准触发单元(35)被配置为使用预定的或时间变化的分散度阈值，特别地其中，所述时间变化的分散度阈值自上一次校准以后随着时间被减小。
5.根据权利要求1所述的装置，
其中，所述校准单元(34)被配置为通过以下操作来确定用于通过所述估计单元(33)针对个体提取的特征和/或提取的特征的组确定所述多个血压估计值的校准值：将提取的特征或提取的特征的组与血压之间的建模的关系拟合到实际提取的特征值或提取的特征值的组和所述血压参考测量结果，所述拟合特别是经由线性或非线性回归进行的。
6.根据权利要求1所述的装置，
其中，所述估计单元(33)被配置为通过获取所述多个血压估计值中的一些或全部的平均值或中值来确定所述对象的血压。
7.根据权利要求1所述的装置，
其中，所述估计单元(33)被配置为通过获取所述多个血压估计值中的一些或全部的加权平均值来确定所述对象的血压，其中，用于所述加权平均值的权重基于对象的血压与对应特征之间的相关性来确定，特别地其中，对象的血压与对应特征之间的所述相关性是基于所获得的血压参考测量结果来确定的。
8.根据权利要求1所述的装置，
其中，所述估计单元(33)被配置为通过获取所述多个血压估计值中的一些或全部的加权平均值来确定所述对象的血压，其中，用于所述加权平均值的权重通过选择其对应特征具有与所述对象的血压的最好相关性的血压估计值、或依据所述血压估计值的对应特征的信噪比、或基于回归误差来确定。
9.根据权利要求1所述的装置，
其中，所述校准单元(34)被配置为执行所述估计单元(33)针对所有提取的特征的初始校准，并且基于对象的血压与对应特征之间的相关性来选择提取的特征的子集，并且
其中，所述估计单元(33)被配置为针对所述子集的个体提取的特征和/或提取的特征的组确定所述多个血压估计值。
10.根据权利要求9所述的装置，
其中，所述校准单元(34)被配置为在针对一个或多个特征与所述对象的血压的所述相关性的一个或多个相关性值低于相关性阈值时重复所述估计单元(33)针对一些或所有提取的特征的所述初始校准。
11.根据权利要求9所述的装置，
其中，所述校准单元(34)被配置为在针对一个或多个特征与所述对象的血压的所述相关性的一个或多个相关性值低于相关性阈值时重复所述估计单元(33)针对一些或所有提取的特征的所述校准。
12.一种用于确定对象的血压的方法，包括：
获得所述对象的动脉脉搏波传感器信号(11)，
从所获得的动脉脉搏波传感器信号中提取多个特征(42)，
针对个体提取的特征和/或提取的特征的组确定多个血压估计值(43)，
根据所述多个血压估计值来确定所述对象的血压(44)，
基于血压参考测量结果(21)来校准对所述多个血压估计值的所述确定，并且
在所述多个血压估计值分散超过分散度极限的情况下触发校准。
13.一种用于监测对象的血压的装置，包括：
动脉脉搏波传感器(10)，其被配置为采集所述对象的动脉脉搏波传感器信号(11)，
血压参考测量单元(20)，其被配置为采集所述对象的血压参考测量结果(21)，以及
根据权利要求1所述的装置(30)，其用于基于所采集的动脉脉搏波传感器信号(11)和所采集的血压参考测量结果(21)来确定对象的血压(31)。
14.根据权利要求13所述的装置，
其中，所述动脉脉搏波传感器(10)包括加速度计、光电体积描记传感器、超声传感器、雷达传感器和生命体征相机中的一个或多个。
15.一种包括程序代码单元的计算机程序，当所述计算机程序在计算机上被执行时，所述程序代码单元用于使所述计算机执行根据权利要求12所述的方法的步骤。
CN201880048936.8A 2017-07-21 2018-07-23 用于确定对象的血压的装置和方法 Active CN111065323B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP17182479.0 2017-07-21
EP17182479.0A EP3430991A1 (en) 2017-07-21 2017-07-21 Apparatus and method for determining blood pressure of a subject
PCT/EP2018/069894 WO2019016405A1 (en) 2017-07-21 2018-07-23 APPARATUS AND METHOD FOR DETERMINING THE BLOOD PRESSURE OF A SUBJECT

Publications (2)

Publication Number Publication Date
CN111065323A true CN111065323A (zh) 2020-04-24
CN111065323B CN111065323B (zh) 2023-04-28

Family

ID=59384018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880048936.8A Active CN111065323B (zh) 2017-07-21 2018-07-23 用于确定对象的血压的装置和方法

Country Status (5)

Country Link
US (1) US11583191B2 (zh)
EP (2) EP3430991A1 (zh)
JP (1) JP7269916B2 (zh)
CN (1) CN111065323B (zh)
WO (1) WO2019016405A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112168155A (zh) * 2020-10-28 2021-01-05 广东小天才科技有限公司 一种血压检测方法、可穿戴设备以及计算机可读存储介质
CN112274127A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于单路脉搏波的无创连续血压检测方法、装置
CN112274120A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于单路脉搏波的无创动脉硬化检测方法、装置
CN112274126A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于多路脉搏波的无创连续血压检测方法、装置
CN112274121A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于多路脉搏波的无创动脉硬化检测方法、装置

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20200144688A (ko) * 2019-06-19 2020-12-30 삼성전자주식회사 생체정보 추정 모델의 캘리브레이션 장치 및 방법과, 생체정보 추정 장치
KR20210004376A (ko) 2019-07-04 2021-01-13 삼성전자주식회사 생체정보 추정 모델의 캘리브레이션 장치 및 방법과, 생체정보 추정 장치
EP3785617A1 (en) * 2019-08-26 2021-03-03 Koninklijke Philips N.V. Method and device for calibrating blood pressure measurements based on measured pulse arrival time
US11367525B2 (en) * 2019-12-20 2022-06-21 Covidien Lp Calibration for continuous non-invasive blood pressure monitoring using artificial intelligence
FR3111064B1 (fr) * 2020-06-08 2022-09-09 Panoramic Digital Health procédé de traitement de mesures effectuées par un capteur porté par une personne
US20220175258A1 (en) * 2020-12-07 2022-06-09 Qualcomm Incorporated Non-invasive blood pressure estimation and blood vessel monitoring based on photoacoustic plethysmography
WO2022146878A1 (en) * 2020-12-30 2022-07-07 Suntech Medical, Inc. Systems, methods and apparatus for generating blood pressure estimations using real-time photoplethysmography data
KR102527271B1 (ko) * 2021-05-24 2023-05-02 주식회사 플렉스시스템 레이더 센서를 이용한 비접촉식 혈압 측정 시스템
CN114305364B (zh) * 2022-01-05 2023-08-11 北京科技大学 基于毫米波雷达的血压检测方法、系统及设备
CN115590489B (zh) * 2022-09-28 2024-01-23 中国人民解放军总医院 一种基于调频连续波雷达的无接触血压监测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5533511A (en) * 1994-01-05 1996-07-09 Vital Insite, Incorporated Apparatus and method for noninvasive blood pressure measurement
JP2000033078A (ja) * 1998-05-12 2000-02-02 Nippon Colin Co Ltd 非観血連続血圧推定装置
CN101088455A (zh) * 2006-06-16 2007-12-19 香港中文大学 对无袖带式动脉血压测量装置进行校准的方法
CN101264011A (zh) * 2007-03-12 2008-09-17 三星电子株式会社 无箍带无侵入地测量腕部血压的方法和装置
US20090076398A1 (en) * 2003-07-07 2009-03-19 Nellcor Puritan Bennett Ireland Continuous Non-Invasive Blood Pressure Measurement Apparatus and Methods Providing Automatic Recalibration

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000107146A (ja) 1998-10-07 2000-04-18 Nippon Colin Co Ltd 血圧監視装置
US9687161B2 (en) * 2008-09-30 2017-06-27 Nellcor Puritan Bennett Ireland Systems and methods for maintaining blood pressure monitor calibration
US9066660B2 (en) * 2009-09-29 2015-06-30 Nellcor Puritan Bennett Ireland Systems and methods for high-pass filtering a photoplethysmograph signal
US9451887B2 (en) * 2010-03-31 2016-09-27 Nellcor Puritan Bennett Ireland Systems and methods for measuring electromechanical delay of the heart
US8825428B2 (en) 2010-11-30 2014-09-02 Neilcor Puritan Bennett Ireland Methods and systems for recalibrating a blood pressure monitor with memory
US20120136261A1 (en) * 2010-11-30 2012-05-31 Nellcor Puritan Bennett Llc Systems and methods for calibrating physiological signals with multiple techniques
JP6613555B2 (ja) 2014-02-13 2019-12-04 日本電気株式会社 血圧推定装置、血圧推定方法、血圧推定プログラム、及び、血圧測定装置
EP3033991B1 (en) * 2014-12-15 2018-02-28 Stichting IMEC Nederland System and method for blood pressure estimation
KR102486700B1 (ko) * 2015-08-11 2023-01-11 삼성전자주식회사 혈압 추정 방법 및 장치

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5533511A (en) * 1994-01-05 1996-07-09 Vital Insite, Incorporated Apparatus and method for noninvasive blood pressure measurement
JP2000033078A (ja) * 1998-05-12 2000-02-02 Nippon Colin Co Ltd 非観血連続血圧推定装置
US20090076398A1 (en) * 2003-07-07 2009-03-19 Nellcor Puritan Bennett Ireland Continuous Non-Invasive Blood Pressure Measurement Apparatus and Methods Providing Automatic Recalibration
CN101088455A (zh) * 2006-06-16 2007-12-19 香港中文大学 对无袖带式动脉血压测量装置进行校准的方法
CN101264011A (zh) * 2007-03-12 2008-09-17 三星电子株式会社 无箍带无侵入地测量腕部血压的方法和装置

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112168155A (zh) * 2020-10-28 2021-01-05 广东小天才科技有限公司 一种血压检测方法、可穿戴设备以及计算机可读存储介质
CN112274127A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于单路脉搏波的无创连续血压检测方法、装置
CN112274120A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于单路脉搏波的无创动脉硬化检测方法、装置
CN112274126A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于多路脉搏波的无创连续血压检测方法、装置
CN112274121A (zh) * 2020-10-28 2021-01-29 河北工业大学 一种基于多路脉搏波的无创动脉硬化检测方法、装置

Also Published As

Publication number Publication date
EP3654837A1 (en) 2020-05-27
WO2019016405A1 (en) 2019-01-24
EP3430991A1 (en) 2019-01-23
CN111065323B (zh) 2023-04-28
JP2020527073A (ja) 2020-09-03
US11583191B2 (en) 2023-02-21
US20200229716A1 (en) 2020-07-23
EP3654837B1 (en) 2024-01-24
JP7269916B2 (ja) 2023-05-09

Similar Documents

Publication Publication Date Title
CN111065323B (zh) 用于确定对象的血压的装置和方法
JP6669409B2 (ja) 血圧値を測定するための方法、機器及びコンピュータプログラム
EP3474749B1 (en) Dynamic calibration of a blood pressure measurement device
JP6310401B2 (ja) 生理的リズムを表す信号を処理する方法、システム及びコンピュータプログラム
US10765374B2 (en) Methods and apparatus for adaptable presentation of sensor data
JP6573988B2 (ja) 血流速度を決定する装置及び方法
US20230233152A1 (en) Methods, apparatus and systems for adaptable presentation of sensor data
KR102206785B1 (ko) 심장탄도의 기준점 사이에서 계측된 시간 간격으로부터 대동맥 맥파 전달 시간을 추정하기 위한 방법 및 기구
US11576585B2 (en) Artifact-tolerant pulse rate variability measurement
EP4173555A1 (en) Device, system and method for calibrating a blood pressure surrogate for use in monitoring a subject&#39;s blood pressure
WO2019006631A1 (zh) 质量评估方法及装置、模型建立方法及模块、穿戴设备
Heimark et al. Accuracy of non-invasive cuffless blood pressure in the intensive care unit: Promises and challenges
EP4173556A1 (en) Device, system and method for calibrating a blood pressure surrogate for use in monitoring a subject&#39;s blood pressure
Chen et al. The effect of signal quality on six cardiac output estimators
WO2017171802A1 (en) Aortic stenosis screening
WO2024047651A1 (en) System and method for sensors integration for non-static continuous blood pressure monitoring
WO2021141572A1 (en) Methods and systems for adaptable presentation of sensor data
Foo et al. Variability in heart rate estimates from different pulse oximeters

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