CN102271574B - 推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法 - Google Patents

推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法 Download PDF

Info

Publication number
CN102271574B
CN102271574B CN200880132612.9A CN200880132612A CN102271574B CN 102271574 B CN102271574 B CN 102271574B CN 200880132612 A CN200880132612 A CN 200880132612A CN 102271574 B CN102271574 B CN 102271574B
Authority
CN
China
Prior art keywords
blood pressure
reverse
value
casp
waveform
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
Application number
CN200880132612.9A
Other languages
English (en)
Other versions
CN102271574A (zh
Inventor
陈崇铭
蔡岳辉
白伟民
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.)
Kin capital technology (Suzhou Industrial Park) Limited
Original Assignee
Kin Capital Technology (suzhou Industrial Park) Ltd
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 Kin Capital Technology (suzhou Industrial Park) Ltd filed Critical Kin Capital Technology (suzhou Industrial Park) Ltd
Publication of CN102271574A publication Critical patent/CN102271574A/zh
Application granted granted Critical
Publication of CN102271574B publication Critical patent/CN102271574B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/024—Detecting, measuring or recording pulse rate or heart rate
    • 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
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Cardiology (AREA)
  • Engineering & Computer Science (AREA)
  • Heart & Thoracic Surgery (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)
  • Public Health (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

用于通过以下步骤而推导中心主动脉收缩压力的方法、系统以及计算机可读介质：将预定数量的血压测量值的集合的顺序逆转，以获得逆转血压集；对逆转血压集求平均，从而平均集代表移动平均波形；使逆转血压集和移动平均波形重叠；识别逆转动脉波形和移动平均波形上的交点；以及将中心主动脉收缩压力设定为逆转血压集中最接近交点的逆转血压值。

Description

推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法
技术领域
本发明涉及推导中心主动脉收缩压力值的方法和分析动脉波形数据以推导中心主动脉收缩压力值的方法。
背景技术
以下对本发明的背景的讨论意在促进对本发明的理解。然而，应当意识到，如在本申请的优先权日的任何管辖之下，该讨论并非将任何所提及的材料承认或接纳为公布、已知或作为常见的一般知识的一部分。
在发达国家，心脏疾病是一种严重的健康问题。潜在的心脏疾病的一项重要指标是从“正常”值的指定范围的血压变化。
一种常见的确定血压值的方法是使用压力袖带在肱动脉处得到收缩血压读数和舒张血压读数这两者。这些值通常被认为指示包括主动脉压力值的其他动脉中的一般血压。虽然该假设在历史上有益于指示潜在的心脏疾病，但最近的研究表明，以该方式测量的“正常的”肱血压值可能掩盖异常的中心主动脉收缩压力值。
Atcor Medical Pty Ltd采用了一种针对该问题的解决方案。Atcor的解决方案使用一般传递函数公式来将径向压力波形转换成中心主动脉血压波形。然而，该一般传递公式假设所有主动脉都相同，该公式基于通过对患者的交叉代表进行测试而确定的相关值，并且，因此，一般假设可能在计算来自不属于交叉代表性样本的范围的患者的中心主动脉压力的计算中将错误引入。该一般传递公式还涉及获得复杂的导数，在将被处理的数据在非常短的时间帧内显著增加时，该导数可能转化为显而易见的时滞。
发明内容
遍及该文献，除非另有指明与此相反，否则，术语“包括”、“包含”等被解释为非穷尽的，或者，换句话说，被解释为意指“包括但不限于”。
根据本发明的第一方面，存在着一种推导中心主动脉收缩压力的方法，包括步骤：
a.获得预定数量的血压测量值的集合；
b.将预定数量的血压测量值的集合的顺序逆转，以获得逆转血压集，所述逆转血压集代表根据下面的公式的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n代表所述预定数量；
c.确定整数区间值；
d.从所述逆转血压集中的第q个逆转血压测量值开始对逆转血压集中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
e.将平均值存储在平均集中；
f.以q值为1开始并且每次增加1来重复步骤d和e，直到q值加上所述整数区间值等于所述集合中的血压测量值的所述预定数量为止，从而值的平均集代表移动平均波形；
g.识别所述逆转动脉波形和所述移动平均波形的交点；以及
h.将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA。
由于逆转动脉波形和移动平均波形之间的交点典型地在逆转动脉波形上血压测量值密度可能更高的区域处，因而这可能具有减少任何时滞影响的优点。而且，从不作关于交叉代表性样本的一般血压波形的假设的血压测量值中推导中心主动脉收缩压力。这可以减少错误的机会。
根据本发明的另一个方面，存在着一种推导中心主动脉收缩压力的方法，包括步骤：
a)接收动脉波形数据集；
b)将所述动脉波形数据集中的所述动脉波形分成具有预定数量的血压测量值的代表性集合；
c)将具有所述预定数量的血压测量值的集合的顺序逆转，以获得逆转血压测量值集，所述逆转血压集代表根据下面的公式的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n是所述预定数量；
d)针对正在被处理的集合而确定整数区间值；
e)从所述集合中的第q个逆转血压测量值开始对所述集合中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
f)将所述平均值存储在平均数据集中；所述平均数据集代表移动平均波形；
g)识别所述逆转动脉波形和所述移动平均波形上的交点；以及
h)将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA。
其中，针对所述动脉波形数据集中的每个动脉波形而重复步骤b至h，并且，对于每个这样的重复，以q值从1开始并且每次增加1来进一步重复步骤e和f，直到q值加上所述整数区间值等于正在被处理的集合中的逆转血压测量值的所述预定数量为止。
根据本发明的另一个方面，存在着一种分析动脉波形数据以推导中心主动脉收缩压力值的方法，包括步骤：
a)接收动脉波形数据集，其中，所述数据集中的每个动脉波形包括具有预定数量的血压测量值的代表性集合；
b)将血压测量值的所述代表性集合的顺序逆转，以获得逆转血压测量值集，所述逆转血压集代表根据下面的公式的逆转动脉波形：
rbp[1...n]＝bp[n...1]
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n是所述预定数量；
c)针对正在被处理的集合而确定整数区间值；
d)从所述逆转血压集中的第q个逆转血压测量值开始对所述集中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
e)将所述平均值存储在平均数据集中，所述平均数据集代表移动平均波形；
f)识别所述逆转血压集和所述移动平均波形上的交点；以及
g)将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA，
其中，针对所述动脉波形数据集中的每个动脉波形而重复步骤b至g，并且，对于每个这样的重复，以q值从1开始并且每次增加1来进一步重复步骤d和e，直到q值加上所述整数区间值等于正在被处理的所述集中的血压测量值的所述预定数量为止。
优选地，设定所述中心主动脉收缩压力(CASP最终)还包括在所述逆转动脉波形上定位与所述点CASPA相邻的两个点(CASPA1和CASPA2)的步骤；根据所述公式而确定所述CASP最终：
CASP最终＝(CASPA1+CASPA2+CASPA)/3。
优选地，压力测量值的集合基本上等于来自所述动脉波形的值的均匀分布。
优选地，确定所述动脉波形的持续时间以及从中的血压测量值的所述预定数量的步骤是根据下面的公式而确定的：
预定数量＝sr×t
其中，
sr＝用于记录所述集合中的血压测量值的测量设备的以Hz为单位的采样率；并且
t＝以秒为单位的所述动脉波形的持续时间。
优选地，所述整数区间值是所述采样率的除法运算。
在优选的实施例中，整数区间值是所述采样率除以5。
在优选的实施例中，整数区间值是所述预定数量的除法运算。
在优选的实施例中，整数区间值在如下确定的边界的范围内：
i范围＝n/(t×v)±(n/(t×30))
其中，
n是所述集合中的血压测量值的预定数量；
t是所述波形的持续时间(以秒为单位)；并且
v是预定的除数。
优选地，所述预定除数是5。
在优选的实施例中，整数区间值等于60除以预定除数。
在优选的实施例中，整数区间值(i)是12。
在优选的实施例中，所述集合中的血压测量值的预定数量等于或大于12。
在优选的实施例中，所述集合中的血压测量值的预定数量为至少15。
根据本发明的另一个方面，存在着一种用于推导中心主动脉收缩压力的系统，包括：
动脉波形测量设备；以及
处理单元，
其中，所述动脉波形测量设备以预定区间得到血压测量值，直到所得到的血压测量值的集合代表至少一个动脉波形为止，代表一个动脉波形的血压测量值的集合然后被通信至所述处理单元，所述处理单元：
a)将所述血压集的顺序逆转，以获得逆转血压测量集，所述逆转血压集代表根据下面的公式而确定的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n是所述预定数量；
b)确定整数区间值；
c)从所述逆转血压集中的第q个逆转血压测量值开始对逆转血压集中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
d)将平均值存储在平均集中；以及
e)识别所述逆转血压集和所述平均集上的交点；以及
f)将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA，
其中，以q值从1开始并且每次增加1来重复步骤c至d，直到q值加上所述整数区间值等于所述集合中的血压测量值的数量为止。
优选地，设定中心主动脉收缩压力(CASP最终)还包括在所述逆转动脉波形上定位与所述点CASPA相邻的两个点(CASPA1和CASPA2)的步骤；根据所述公式而确定所述CASP最终：
CASP最终＝(CASPA1+CASPA2+CASPA)/3。
根据本发明的另一个方面，存在着一种计算机可读介质，在其上记录有软件，以推导中心主动脉收缩压力，所述软件包括：
用于将具有预定数量的血压测量值的集合的顺序逆转以获得逆转血压测量值的集合的模块，逆转血压测量值代表根据下面的公式而确定的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n是所述预定数量；
第一存储模块，其用于将所述预定数量的逆转血压测量值存储在逆转血压集中；
用于确定整数区间值的模块；
求平均模块，其从所述逆转血压集中的第q个逆转血压测量值开始对所述集中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
第二存储模块，其用于将平均值存储在平均集中，所述平均集代表平均波形；
识别模块，其用于识别所述逆转动脉波形和所述平均波形上的交点；以及
结果模块，其用于将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA。
其中，以q值从1开始并且每次增加1来重复所述求平均模块和所述第二存储模块的功能，直到q值加上所述整数区间值等于所述集合中的血压测量值的所述预定数量为止。
根据本发明的另一个方面，存在着一种计算机可读介质，在其上记录有软件以分析动脉波形数据从而推导中心主动脉收缩压力值，所述软件包括：
通信模块，其用于接收动脉波形数据集；
设定模块，其用于将所述动脉波形数据集中的所述动脉波形分成具有预定数量的血压测量值的代表性集合；
用于将所述预定数量的血压测量值的顺序逆转以获得逆转血压测量值的集合的模块，所述逆转血压测量值代表根据下面的公式而确定的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n是所述预定数量；
区间模块，其用于针对正在被处理的所述逆转血压集而确定整数区间值；
求平均模块，其从所述逆转血压集中的第f个逆转血压测量值开始对所述集中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
第一存储模块，其用于将平均值存储在平均集中；
识别模块，其用于识别所述逆转动脉波形和所述平均波形上的交点；
结果模块，其将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点处的所述逆转血压值的逆转血压值CASPA；以及
第二存储模块，其用于将CASP最终值存储在与正在被处理的所述动脉波形在所述动脉波形数据集中所占据的位置相对应的位置处，其中，针对所述动脉波形数据集中的每个动脉波形而重复所述设定模块、区间模块、求平均模块、第一存储模块以及第二存储模块的功能，并且，针对每个这样的重复，以所述f值从1开始并且每次增加1来进一步重复所述求平均模块和第一存储模块的功能，直到所述f值加上所述整数区间值等于正在被处理的所述集合中的血压测量值的预定数量为止。
根据本发明的另一方面，存在着一种计算机可读介质，在其上记录有软件以分析动脉波形数据从而推导中心主动脉收缩压力值，所述软件包括：
通信模块，其用于接收动脉波形数据集，其中，所述数据集中的每个动脉波形包括具有预定数量的血压测量值的代表性集合；
用于将所述预定数量的血压测量值的顺序逆转以获得逆转血压测量值的集合的模块，所述逆转血压测量值代表根据下面的公式而确定的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表血压测量值的集合，
rbp代表逆转血压测量值的集合，
并且，n是所述预定数量；
区间模块，其用于针对正在被处理的所述逆转血压集而确定整数区间值；
求平均模块，其用于从所述集合中的第q个逆转血压测量值开始对所述集合中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
第一存储模块，其用于将平均值存储在平均集中；
识别模块，其用于识别所述逆转动脉波形和所述平均波形上的交点；
结果模块，其将所述中心主动脉收缩压力(CASP最终)设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA；以及
第二存储模块，其用于存储CASP最终值，
其中，针对所述动脉波形数据集中的每个动脉波形而重复所述区间模块、求平均模块、第一存储模块以及第二存储模块的功能，并且，针对每个这样的重复，以q值从1开始并且每次增加1来进一步重复所述求平均模块和第一存储模块的功能，直到q值加上所述整数区间值等于正在被处理的所述集合中的血压测量值的预定数量为止。
优选地，设定中心主动脉收缩压力(CASP最终)还包括在所述逆转动脉波形上定位与所述点CASPA相邻的两个点(CASPA1和CASPA2)的步骤；根据所述公式而确定所述CASP最终：
CASP最终＝(CASPA1+CASPA2+CASPA)/3。
附图说明
现在，仅经由示例并参考附图而描述本发明，其中：
图1是根据本发明的第一实施例的确定动脉的中心主动脉压力值的方法的流程图；
图2是根据本发明的第二实施例的分析动脉波形数据集以产生相应的中心主动脉压力数据集的方法的流程图；
图3a是根据本发明的第三实施例的包括测试值以产生相应的中心主动脉收缩压力数据集的表；
图3b是描绘逆转动脉波形和移动平均波形以及它们的交点的图，能够通过当前发明的任何实施例对其处理以得到中心主动脉收缩压力值。
具体实施方式
现在，参考附图而描述本发明的特定的实施例。在本文中所使用的术语仅出于描述本发明的特定的实施例的目的，并不意在限制本发明的范围。另外，除非另有定义，在本文中使用的所有技术术语和科学术语具有与本发明所属的领域中的普通技术人员所普遍理解的意思相同的意思。
根据本发明的第一实施例，提供一种确定中心主动脉收缩压力值的方法10。在图1中，以流程图形式图解说明方法10。
如在图1中所显示地，方法10以获得代表动脉波形的血压(bp)测量(bp)值的集合开始(步骤12)。如在本领域中已知地，能够在动脉处测量动脉波形，并且，动脉波形可以直接使用或可以存储在数据库中。
将血压测量值的集合逆转，以便获得代表逆转动脉波形的逆转血压(rbp)测量值集的集合(步骤14)，该逆转动脉波形是参照有关血压轴的动脉波形的镜像。
血压测量值/逆转血压测量值的集合具有预定数量的测量值(n)。血压测量值和逆转血压测量值基于下面的公式而彼此相关：
rbp[1...n]＝bp[n...1]
预定数量(n)随后除以所需要的时间(以秒为单位)，以完成动脉波形，以便获得采样率(sr)。
在步骤16，采样率(sr)除以预定的除数(v)以确定整数区间值(i)。以该方式，整数区间值(i)可能服从于绝对或舍入数学函数。
然后，重复下面的步骤，从第一逆转血压测量值rbp[(q)]开始：
在步骤18，对逆转血压测量值{rbp[q]，rbp[q+1]，rbp[q+2]，rbp[q+3]，...，rbp[q+i-1]}进行求和，以形成求和值(s)。
在步骤20，求和值(s)除以整数区间值(i)，以形成平均逆转血压值(a)。
在步骤22，平均逆转血压值(a)存储在值的移动平均集中。
按顺序重复步骤18、20以及22，直到满足条件[q+i-1]等于预定数量(n)为止。
所获得的值的移动平均集代表移动平均波形。
步骤24识别逆转动脉波形上具有与逆转动脉波形和移动平均波形之间的交点处的血压值最接近的血压值的逆转血压值(CASPA)。
在步骤26中，通过在逆转动脉波形上定位与逆转血压值(CASPA)相邻的两个点(CASPA1和CASPA2)而确定中心主动脉收缩压力值。在步骤28，根据下面的公式而计算中心主动脉收缩压力值(CASP最终)：
(CASP最终)＝(CASPA1+CASPA2+CASPA)/3
根据本发明的第二实施例，提供一种分析动脉波形数据集以产生相应的中心主动脉收缩压力数据集的方法100。该实施例的一个常见的前身是医务人员根据如本领域技术人员所已知的任何技术而获得动脉的血压(bp)的动脉波形代表。然后，将动脉波形代表提供至中心处理站，以便确定相应的中心主动脉收缩压力值。在接收动脉波形代表之后，中心处理站如下地运算：
针对波形数据集中的每个动脉波形，获得逆转动脉波形，从而逆转动脉波形是参照有关血压轴的动脉波形的镜像。
在步骤110，将正在被处理的逆转动脉波形分成代表性的逆转血压(rbp)测量值的集合。血压测量值和逆转血压测量值的集合具有预定数量的测量值(n)。血压测量值和逆转血压测量值基于下面的公式而彼此相关：
rbp[1...n]＝bp[n...1]
每个预定数量(n)随后除以所需要的时间(以秒为单位)，以完成有关的动脉波形，以便获得采样率(sr)。
在步骤112，采样率(sr)除以预定的除数(v)以确定整数区间值(i)。
然后，重复下面的步骤，从第一逆转血压测量值rbp[(q)]开始：
在步骤118，对逆转血压测量值{rbp[q]，rbp[q+1]，rbp[q+2]，rbp[q+3]，...，rbp[q+i-1]}进行求和，以形成求和值(s)；
在步骤120，求和值(s)除以整数区间值(i)，以形成平均逆转血压值(a)。
在步骤122，平均逆转血压值(a)存储在移动平均集中。
按顺序重复步骤118、120以及122，直到满足条件[q+i-1]等于预定数量(n)为止。
所获得的值的移动平均集代表移动平均波形。
步骤124识别逆转动脉波形上血压值最接近逆转动脉波形和移动平均波形之间的交点的逆转血压值(CASPA)。
在步骤126中，通过在逆转动脉波形上定位与逆转血压值(CASPA)相邻的两个点(CASPA1和CASPA2)而获得中心主动脉收缩压力值。然后，在步骤128，根据下面的公式而计算中心主动脉收缩压力值(CASP最终)：
(CASP最终)＝(CASPA1+CASPA2+CASPA)/3
然后，按顺序重复步骤110、112、118、120、122、124、126、128，直到针对波形数据集中的每个动脉波形而计算相应的中心主动脉收缩压力值的时间为止。
根据本发明的第三实施例，其中，相似的数字引用第一实施例的相似的步骤，提供一种分析动脉波形数据集以产生相应的中心主动脉收缩压力值(未显示)的方法。该实施例基于申请人所进行的试验，其中，已发现，预定数量(n)能够是任何超过15的值。
在该实施例中，虽然优选预定的除数(v)为5以确定适当的区间，但已能够在整数区间值(i)处于一范围中的情况下获得基本上准确的主动脉收缩压力值，其中，该范围的极限如下地确定：
i范围＝n/(t×v)±(n/(t×30))
其中，在该实施例中，变量t代表单个波形的持续时间(以秒为单位)。
现在，将在以下示例的背景下更详细地描述该实施例：
根据步骤12，获得代表动脉波形的来自动脉的血压(bp)测量值的集合。在图3a中描绘该测量值。
基于如在步骤14中所描述地将血压测量值集的顺序逆转而获得逆转血压(rbp)测量值集。
完成该波形所需的持续时间是0.683秒。
在目前的实施例中，预定数量(n)设定在41。预定的除数(v)为5，整数区间值范围(i范围)如下地确定：
i范围＝n/(t×v)±(n/(t×30))，因此，
i范围＝41/(0.683×5)±(41/(0.683×30))，因此，
i范围＝12.01±(2.00)
由于区间不可能是除了整数值以外的值，因而i范围值被限定为10至14的范围。出于该示例的目的，应当使用整数区间值(i)12。
遵循步骤18至22的要求，获得值的移动平均集。
按照时间对值的移动平均集和逆转血压测量值进行绘图，如在图3b中所描绘地，获得平均波形和逆转动脉波形之间的交点。
CASPA是与逆转动脉波形和移动平均波形之间的交点最接近的血压值。如能够从绘图所看出地，CASPA值是116.8。CASPA1和CASPA2是逆转动脉波形上与CASPA相邻的两个血压值。如能够从绘图所看出地，CASPA1值是115.6，并且，CASPA2值是117.5。
所计算的CASP最终值是(115.6+117.5+116.8)/3＝116.6mmHg
在本发明的第四个且最优选的实施例中，提供一种确定如在本发明的第一实施例中所描述的中心主动脉收缩压力值的方法。然而，在该实施例中，根据下面的公式确定预定数量(n)：
n＝sr×t
其中，
sr＝用于记录该集合bp中的血压测量值的测量设备的采样率(以赫兹为单位)；并且
t＝完成一个动脉波形所需的时间(以秒为单位)。
在该实施例中，该集合bp中的血压测量值是由测量设备在采样率的每次重复所得到的测量值。此外，在本发明的该实施例中，区间值(i)处于极限如下地确定的范围中：
i范围＝sr/v±(sr/30)
在该情况下，sr再次代表用于记录该集合bp中的血压测量值的测量设备的采样率(以赫兹为单位)。另外，预定的除数(v)优选为5。该公式只有在值(sr)超过15时才可应用。最小极限应当为15，然而，出于图解说明的目的，在优选的实施例中使用值30。
使用如上面所阐明同一集合的血压测量值，但此时在该集合是具有60次每秒的采样率(sr)的采样设备所得到的血压值的顶点的情形下，如下地计算i范围值：
i范围＝sr/v±(sr/30)
i范围＝60/5±(60/30)
i范围＝12±2
这导致10和14之间的i范围值。然后，该方法的剩下的步骤能够如上面在第一实施例中所阐明地继续进行。
在所描述的所有实施例中，逆转动脉波形和移动平均波形之间的交点典型地出现在逆转动脉波形上血压测量密度值更高的区域处。这是图3a的区域A。如与区域B对比，在区域A中血压测量值密度更高。区域A和B描绘曲线的有关部分，但并非基于水平或垂直参考的测量值。
本领域技术人员应当意识到，上述的发明不限于所描述的实施例。特别地，在不背离本发明的范围的情况下，可以进行下面的修改和改进：
●可以从也适于执行本发明的方法的设备获得血压测量值的集合。可替代地，血压测量值的集合也可以从单独的设备获得并通信至适于执行本发明的方法的又一装置。
●可以修改本发明的第二实施例，从而波形数据集数据可以包含代表波形的数据集。
●可以将波形数据集与数据集中的每个波形的持续时间一起供给至执行该方法的实体。可替代地，执行该方法的实体可以通过可替代的手段而独立地确定每个波形的持续时间(例如通过接收具有x轴的预定长度的固定的时间值的图并基于该时间/距离关系而近似波形的持续时间，或者通过从诸如用于生成血压测量值的集合的采样率的其他组合值而推导持续时间)。
●为了成为最具代表性的动脉波形，形成该集合bp的血压值应当沿着动脉波形均匀分布。
●中心主动脉收缩压力值(CASP最终)可以识别为最接近逆转动脉波形和移动平均波形之间的交点的中心主动脉收缩压力集的平均值(CASPA)，即，CASP最终＝CASPA。
●可以使用获得交点的其他方式。这些方式包括使用导数方法、解决联立方程以及用于识别两个波形之间的交点的其他数学工具。
●可以通过例如获得中值的其他方法或具有使逆转动脉波形平滑化的效果的其他这样的方法而实现获得移动平均的手段。
本领域普通技术人员应当进一步意识到，上面所描述的并不互相排斥的特征能够结合以形成本发明的又一实施例。

Claims (13)

1.一种用于推导中心主动脉收缩压力的系统，包括：
动脉波形测量设备；以及
处理单元，
其特征在于，所述动脉波形测量设备以预定的区间得到血压测量值，直到所得到的血压测量值的集合代表至少一个动脉波形为止，代表一个动脉波形的血压测量值的集合然后被通信至所述处理单元，所述处理单元执行如下动作：
a)将所述血压测量值的集合的顺序逆转，以获得逆转血压集，所述逆转血压集代表根据下面的公式而确定的逆转动脉波形：
rbp[1...n]＝bp[n...1]，
其中：
bp代表所述血压测量值的集合，
rbp代表所述逆转血压集，
并且，n是所述预定的区间的数量；
b)确定整数区间值；
c)从逆转血压集中的第q个逆转血压测量值开始对所述逆转血压集中等于所述整数区间值的一系列连续的逆转血压测量读数求平均；
d)将平均值存储在平均集中；以及
e)识别所述逆转血压集和所述平均集上的交点；以及
f)将所述中心主动脉收缩压力CASP最终设定为所述逆转血压集中最接近所述交点的逆转血压值CASPA，
其中，以q值从1开始并且每次增加1来重复动作c至d，直到q值加上所述整数区间值等于所述集合中的血压测量值的数量为止。
2.如权利要求1所述的系统，其中，所述动作f还包括在所述逆转动脉波形上定位与所述逆转血压值CASPA所位于的点相邻的两个点的动作，其中，所述两个点分别具有逆转血压值CASPA1和CASPA2；根据以下公式而确定所述CASP最终：
CASP最终＝(CASPA1+CASPA2+CASPA)/3。
3.如权利要求1或2所述的系统，其中，血压测量值的所述集合基本上等于来自所述动脉波形的值的均匀分布。
4.如权利要求3所述的系统，其中，所述处理单元确定所述动脉波形的持续时间，并且其中，根据下面的公式而确定所述预定的区间的数量：
所述预定的区间的数量＝sr×t
其中，
sr＝用于记录所述集合中的所述血压测量值的测量设备的以Hz为单位的采样率；并且
t＝以秒为单位的所述动脉波形的持续时间。
5.如权利要求4所述的系统，其中，所述整数区间值是所述采样率的除法运算。
6.如权利要求5所述的系统，其中，所述整数区间值是所述采样率除以5。
7.如权利要求1或权利要求2所述的系统，其中，所述整数区间值是所述预定的区间的数量的除法运算。
8.如权利要求7所述的系统，其中，所述整数区间值在如下确定的边界的范围内：
i范围＝n/(t×v)±(n/(t×30))
其中，
n是所述预定的区间的数量；
t是以秒为单位的所述波形的持续时间；并且
v是预定的除数。
9.如权利要求8所述的系统，其中，所述预定的除数是5。
10.如权利要求1或权利要求2所述的系统，其中，所述整数区间值等于60除以预定的除数。
11.如权利要求10所述的系统，其中，所述整数区间值i是12。
12.如权利要求1或2所述的系统，其中，所述预定的区间的数量等于或大于12。
13.如权利要求12所述的系统，其中，所述预定的区间的数量为至少15。
CN200880132612.9A 2008-12-05 2008-12-05 推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法 Active CN102271574B (zh)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/SG2008/000468 WO2010064993A1 (en) 2008-12-05 2008-12-05 Method of deriving central aortic systolic pressure values and method for analysing an arterial dataset to derive the same

Publications (2)

Publication Number Publication Date
CN102271574A CN102271574A (zh) 2011-12-07
CN102271574B true CN102271574B (zh) 2015-09-09

Family

ID=42233476

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200880132612.9A Active CN102271574B (zh) 2008-12-05 2008-12-05 推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法

Country Status (15)

Country Link
US (1) US8657753B2 (zh)
EP (1) EP2373214B1 (zh)
JP (1) JP5448213B2 (zh)
KR (1) KR101523025B1 (zh)
CN (1) CN102271574B (zh)
AU (1) AU2008364802B2 (zh)
CA (1) CA2745026A1 (zh)
HK (1) HK1158480A1 (zh)
IL (1) IL213392A (zh)
MX (1) MX2011005854A (zh)
RU (1) RU2472428C1 (zh)
SG (1) SG171770A1 (zh)
TW (1) TWI484938B (zh)
WO (1) WO2010064993A1 (zh)
ZA (1) ZA201103953B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102014007519A1 (de) * 2014-05-22 2015-11-26 Fresenius Medical Care Deutschland Gmbh Verfahren und Vorrichtung zur Bestimmung des zentralen systolischen Blutdrucks
CN107072538B (zh) 2014-09-08 2021-07-13 苹果公司 将脉搏传导时间(ptt)测量系统电耦接到心脏以用于血压测量
WO2016040256A1 (en) 2014-09-08 2016-03-17 Braintree Analytics Llc Systems, devices, and methods for measuring blood pressure of a user
CN107106054B (zh) 2014-09-08 2021-11-02 苹果公司 使用多功能腕戴式设备进行血压监测
WO2016040263A1 (en) 2014-09-08 2016-03-17 Braintree Analytics Llc Wrist worn accelerometer for pulse transit time (ptt) measurements of blood pressure
US20220369936A1 (en) * 2021-05-23 2022-11-24 Biobeat Technologies Ltd. Blood pressure measuring device and method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5265011A (en) * 1989-04-03 1993-11-23 Eastern Medical Testing Services, Inc. Method for ascertaining the pressure pulse and related parameters in the ascending aorta from the contour of the pressure pulse in the peripheral arteries
US5746698A (en) * 1995-09-28 1998-05-05 Nederlandse Organisatie Voor Toegepast-Natuurwetenschappelijk Onderzoek Tno Method and device for determining brachial arterial pressure wave on the basis of nonivasively measured finger blood pressure wave

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU746552A1 (ru) * 1977-07-05 1980-07-07 Первый Московский Ордена Ленина И Ордена Трудового Красного Знамени Медицинский Институт Им. И.М. Сеченова Устройство дл обработки данных о давлении крови
US4307727A (en) * 1979-10-15 1981-12-29 Tech Engineering And Design, Inc. Wrist band transducer support and tensioning apparatus
IL86010A (en) * 1987-07-03 1993-02-21 Hi Bridge Kk Apparatus for measuring blood pressure
US5908027A (en) * 1994-08-22 1999-06-01 Alaris Medical Systems, Inc. Tonometry system for monitoring blood pressure
SG94349A1 (en) * 2000-10-09 2003-02-18 Healthstats Int Pte Ltd Method and device for monitoring blood pressure
WO2003090617A1 (fr) * 2002-04-24 2003-11-06 Colin Corporation Procede d'analyse d'une onde d'impulsion, logiciel d'analyse d'une onde d'impulsion, etc.
US7828739B2 (en) * 2002-08-27 2010-11-09 Precision Pulsus, Inc. Apnea detection system
US20070197924A1 (en) * 2004-03-05 2007-08-23 O'rourke Michael F Method and apparatus for determination of cardiac output from the arterial pressure pulse waveform
JP4668651B2 (ja) 2005-03-16 2011-04-13 フクダ電子株式会社 血圧計
RU2337613C1 (ru) * 2007-04-23 2008-11-10 Государственное образовательное учреждение высшего профессионального образования "Пензенский государственный университет"(ПГУ) Способ экспресс-измерения артериального давления
US8469895B2 (en) * 2007-06-07 2013-06-25 Healthstats International Pte Ltd Deriving central aortic systolic pressure and analyzing arterial waveform data to derive central aortic systolic pressure values

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5265011A (en) * 1989-04-03 1993-11-23 Eastern Medical Testing Services, Inc. Method for ascertaining the pressure pulse and related parameters in the ascending aorta from the contour of the pressure pulse in the peripheral arteries
US5746698A (en) * 1995-09-28 1998-05-05 Nederlandse Organisatie Voor Toegepast-Natuurwetenschappelijk Onderzoek Tno Method and device for determining brachial arterial pressure wave on the basis of nonivasively measured finger blood pressure wave

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Parametric Model Derivation of Transfer Function for Noninvasive Estimation of Aortic Pressure by Radial Tonometry;Barry Fetics等;《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》;19990630;第46卷(第6期);698-706 *

Also Published As

Publication number Publication date
IL213392A0 (en) 2011-07-31
US8657753B2 (en) 2014-02-25
AU2008364802B2 (en) 2014-12-11
CA2745026A1 (en) 2010-06-10
TW201026281A (en) 2010-07-16
AU2008364802A1 (en) 2010-06-10
IL213392A (en) 2017-10-31
TWI484938B (zh) 2015-05-21
JP5448213B2 (ja) 2014-03-19
EP2373214B1 (en) 2016-03-23
MX2011005854A (es) 2011-06-20
KR101523025B1 (ko) 2015-05-27
US20110230773A1 (en) 2011-09-22
KR20110094103A (ko) 2011-08-19
RU2472428C1 (ru) 2013-01-20
ZA201103953B (en) 2012-02-29
EP2373214A4 (en) 2014-03-05
CN102271574A (zh) 2011-12-07
EP2373214A1 (en) 2011-10-12
HK1158480A1 (zh) 2012-07-20
JP2012510850A (ja) 2012-05-17
WO2010064993A1 (en) 2010-06-10
SG171770A1 (en) 2011-07-28

Similar Documents

Publication Publication Date Title
CN102271574B (zh) 推导中心主动脉收缩压力值的方法和用于分析动脉数据集以推导中心主动脉收缩压力值的方法
US11172891B2 (en) Method and apparatus for deriving mean arterial pressure of a subject
US4676253A (en) Method and apparatus for noninvasive determination of cardiac output
EP2087836B1 (en) Apparatus and method for determining a physiological parameter
US11006885B2 (en) Apparatus for determining blood pressure
US8469895B2 (en) Deriving central aortic systolic pressure and analyzing arterial waveform data to derive central aortic systolic pressure values
US8257269B2 (en) Apparatus for analysing pulse using array of pressure sensors
Suzuki et al. Feature selection method for estimating systolic blood pressure using the Taguchi method
US7477934B2 (en) Method of monitoring human relaxation level, and user-operated heart rate monitor
CN109512412B (zh) 一种中心动脉血压测量装置
JP6043920B2 (ja) 胎児心電位信号抽出プログラム、胎児心電位信号判別装置及びこれを用いた妊婦見守りシステム
US20170055858A1 (en) Method and system for facilitating patient self-measuring and recording
KR20070056925A (ko) 혈압측정 방법 및 장치
CN109157204A (zh) 一种无袖带式腕部动脉血压测量方法及系统
WO2019153579A1 (zh) 一种基于心电信号的无气囊血压检测方法及系统
EP3435853B1 (en) A method of deriving systolic blood pressure and/or diastolic blood pressure of a subject
CN107920762A (zh) 利用脉搏传导时间来计算收缩压的装置及其方法
CN101420903A (zh) 获得中心动脉收缩压的方法和分析动脉波形数据以获得中心动脉收缩压值的方法
JP2798737B2 (ja) 心拍性同期信号採取装置
Poliński et al. A meta-analysis of pulse arrival time based blood pressure estimation
CN108309272A (zh) 一种基于服务器的云计算血压计
Setiawan et al. Development of non-occlusive blood pressure monitor with height correction for home used application
WO2012119665A1 (en) Method and device for long-term variability monitoring of cardiovascular parameters based on ambulatory registration of electrocardiogram and pulse wave signals
CN113598724A (zh) 用于评估心血管自主神经功能的装置
SU1262770A1 (ru) Устройство дл анализа тахоосциллограмм

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: JIANZI TECHNOLOGY (SUZHOU INDUSTRIAL PARK) CO., LT

Free format text: FORMER OWNER: HEALTHSTATS INTERNATIONAL PTE LTD.

Effective date: 20150813

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20150813

Address after: 3 floor, building 3, 35 Hongfeng Road, Suzhou Industrial Park, Jiangsu, Suzhou, 215021

Applicant after: Kin capital technology (Suzhou Industrial Park) Limited

Address before: Singapore Singapore

Applicant before: Healthstats International Pte Ltd.

C14 Grant of patent or utility model
GR01 Patent grant