CN104188637B - 一种基于波形匹配法的主动脉脉搏波传导时间获取方法 - Google Patents

一种基于波形匹配法的主动脉脉搏波传导时间获取方法 Download PDF

Info

Publication number
CN104188637B
CN104188637B CN201410493975.2A CN201410493975A CN104188637B CN 104188637 B CN104188637 B CN 104188637B CN 201410493975 A CN201410493975 A CN 201410493975A CN 104188637 B CN104188637 B CN 104188637B
Authority
CN
China
Prior art keywords
matched
pressure waveform
normalization
pulse wave
amplitude
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
CN201410493975.2A
Other languages
English (en)
Other versions
CN104188637A (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.)
Hefei Institutes of Physical Science of CAS
Anqing Normal University
Original Assignee
Hefei Institutes of Physical Science of CAS
Anqing Normal University
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 Hefei Institutes of Physical Science of CAS, Anqing Normal University filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201410493975.2A priority Critical patent/CN104188637B/zh
Publication of CN104188637A publication Critical patent/CN104188637A/zh
Application granted granted Critical
Publication of CN104188637B publication Critical patent/CN104188637B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明提供了一种基于波形匹配法的主动脉脉搏波传导时间获取方法,包括:采集颈动脉脉搏波和股动脉脉搏波;利用肱动脉血压值对颈动脉脉搏波和股动脉脉搏波进行标定,获得颈动脉压力波形和股动脉压力波形;选取颈、股动脉压力波形一段特定区域并进行归一化处理,获得两段待匹配波形;利用波形匹配法对两段待匹配波形进行分析,计算出主动脉脉搏波传导时间aoPTT。相对于传统的“foot-to-foot”法,本发明基于波形一段特定的区域而不是单个特征点计算aoPTT,提高了主动脉脉搏波速度aoPWV测量的准确性和抗干扰能力,具有更广的应用前景。

Description

一种基于波形匹配法的主动脉脉搏波传导时间获取方法
技术领域
本发明涉及生物医学工程技术领域,具体涉及一种基于波形匹配法的主动脉脉搏波传导时间获取方法。
背景技术
随着我国经济水平的发展,人民生活水平提高,不合理的饮食结构以及人口老年化导致心血管发病率和死亡率迅速上升。《中国心血管报告2012》显示,我国心血管病患者2.9亿,每5个成人中就有1人患有心血管疾病。我国每年约有350万人死于心血管疾病,占总死亡原因的41%,即每10秒钟心血管死亡一人。心血管病已经成为二十一世纪威胁人类健康的重大公共卫生问题。
血管动脉结构以及弹性功能的改变是造成心血管疾病的主要因素。动脉弹性的改变早于结构的改变,且结构的改变往往是不可逆的。因此,对心血管疾病人群进行早期检测,积极干预具有非常重要的意义。
依赖于“动脉硬化程度加剧导致由心脏输出的血液产生的波动(脉搏波)的传导速度会加快”这一原理,学者们提出使用主动脉脉搏波速度aoPWV作为无创动脉僵硬度检测方法。其中,主动脉脉搏波传导时间aoPTT是测量aoPWV的关键。目前,aoPTT一般是通过测量颈股动脉传导时间cfPTT获取,这种方法的关键在于识别同一心动周期内两路波形的起点,计算它们的时间差用于评估aoPTT,现有的aoPTT评估方法包括舒张期极小值法、一阶导数极大值法、二阶导数极大值法和切线相交法,这些方法都是基于波形单个特征点,称为“foot-to-foot”法。但是,由于动脉壁弹性的差异,脉搏波的形状在传播过程中会发生变化,用于识别波形起点位置的特征点会发生转移,导致“foot-to-foot”法评估的结果不够准确,影响aoPTT的测量。因此,准确测量aoPTT是目前急需解决的问题。
发明内容
本发明是为解决上述“foot-to-foot”法所存在的特征点识别不准导致aoPTT测量存在误差的问题,提供一种基于波形匹配法的主动脉脉搏波传导时间获取方法。
本发明为解决技术问题采用如下技术方案:
本发明基于波形匹配法的主动脉脉搏波传导时间获取方法,它是一种通过时移匹配技术分析实时采集的人体颈动脉和股动脉脉搏波序列获取主动脉脉搏波传导时间aoPTT的方法,其特点在于包括如下步骤:
(1)利用脉搏波传感器同步采样获取人体颈动脉脉搏波和股动脉脉搏波;
(2)利用人体肱动脉血压值对所述颈动脉脉搏波和股动脉脉搏波进行标定,分别获得具有血压特性的颈动脉脉搏波和具有血压特性的股动脉脉搏波;所述具有血压特性的颈动脉脉搏波即为颈动脉压力波形CP,所述具有血压特性的股动脉脉搏波即为股动脉压力波形FP;
(3)选择颈动脉压力波形CP上升阶段斜率最大值点,记为cMaxPoint;选择颈动脉压力波形CP舒张期极小值点,记为cMinPoint;选择股动脉压力波形FP上升阶段斜率最大值点,记为fMaxPoint;选择股动脉压力波形FP舒张期极小值点,记为fMinPoint;以位于cMaxPoint和cMinPoint之间的波形,作为颈动脉待匹配压力波形,以位于fMaxPoint和fMinPoint之间的波形,作为股动脉待匹配压力波形,对所述颈动脉待匹配压力波形和所述股动脉待匹配压力波形作如下归一化处理:
3a、比较cMinPoint处和fMinPoint处的幅值大小:
若cMinPoint处的幅值>fMinPoint处的幅值,则以cMinPoint处的幅值作为归一化最小幅值,以cMinPoint作为颈动脉待匹配压力波形的归一化起点,记为cStartPoint,以所述归一化最小幅值在股动脉待匹配压力波形上的位置,作为股动脉待匹配压力波形的归一化起点,记为fStartPoint;
若cMinPoint处的幅值≤fMinPoint处的幅值,则以fMinPoint处的幅值作为归一化最小幅值,以fMinPoint作为股动脉待匹配压力波形的归一化起点,记为fStartPoint,以所述归一化最小幅值在颈动脉待匹配压力波形上的位置,作为颈动脉待匹配压力波形的归一化起点,记为cStartPoint;
3b、比较cMaxPoint处和fMaxPoint处的幅值大小:
若cMaxPoint处的幅值≤fMaxPoint处的幅值,则以cMaxPoint处的幅值作为归一化最大幅值,以cMaxPoint作为颈动脉待匹配压力波形的归一化终点,记为cEndPoint,以所述归一化最大幅值在股动脉待匹配压力波形上的位置,作为股动脉待匹配压力波形的归一化终点,记为fEndPoint;
若cMaxPoint处的幅值>fMaxPoint处的幅值,则以fMaxPoint处的幅值作为归一化最大幅值,以fMaxPoint作为股动脉待匹配压力波形的归一化终点,记为fEndPoint,以所述归一化最大幅值在颈动脉待匹配压力波形上的位置,作为颈动脉待匹配压力波形的归一化终点,记为cEndPoint;
3c、以颈动脉待匹配压力波形中cStartPoint与cEndPoint之间的波形作为归一化后的颈动脉待匹配压力波形CPm;以股动脉待匹配压力波形中fStartPoint与fEndPoint之间的波形作为归一化后的股动脉待匹配压力波形FPm;
(4)对归一化后的颈动脉待匹配压力波形CPm和归一化后的股动脉待匹配压力波形FPm通过波形匹配法求取aoPTT:
将归一化后的颈动脉待匹配压力波形CPm沿着时间轴向归一化后的股动脉待匹配压力波形FPm移动,每次移动1个步骤(1)所用的采样周期,共移动的时间为T,所述T为cStartPoint距离fEndPoint之间的时间,共移动的次数m=T/采样周期;
在每次移动后计算CPm上各采样点与FPm中相应采样点的幅值差的绝对值之和SADi或幅值差的平方之和SSDi,其中i为移动次数,i=1,2,、、、,m;
令SADi或SSDi的最小值所对应的移动次数为p,则传导时间aoPTT=p×采样周期。
在每次移动后CPm上各采样点与FPm中相应采样点的幅值差的绝对值之和SADi是按如下方法进行计算:
(1)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于相同时间段上的A部分(即在时间轴上重叠的部分),计算获得位于A部分的CPm上各采样点与FPm上相同时间点处的采样点的幅值差的绝对值,然后计算由A部分CPm上所有采样点所获得的幅值差的绝对值之和,记为SADr;
(2)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于不同时间段上的B部分(即在时间轴上不重叠的部分),计算获得位于B部分的CPm上各采样点的幅值与cStartPoint处幅值差的绝对值之和SADpc,和位于B部分的的FPm上各采样点的幅值与fEndPoint处幅值差的绝对值之和SADpf,将SADpc和SADpf相加获得B部分上所有采样点的幅值差的绝对差值之和,记为SADp;
(3)将SADr和SADp相加,即得在每次移动后CPm与FPm中各采样点的幅值差的绝对值之和SADi
在每次移动后CPm上各采样点与FPm中相应采样点的幅值差的平方之和SSDi是按如下方法进行计算:
(1)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于相同时间段上的A部分(即在时间轴上重叠的部分),计算获得位于A部分的CPm上各采样点与FPm上相同时间点处的采样点的幅值差的平方,然后计算由A部分CPm上所有采样点所获得的幅值差的平方之和,记为SSDr;
(2)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于不同时间段上的B部分(即在时间轴上不重叠的部分),计算获得位于B部分的CPm上各采样点的幅值与cStartPoint处幅值差的平方之和SSDpc,和位于B部分的的FPm上各采样点的幅值与fEndPoint处幅值差的平方之和SSDpf,将SADpc和SADpf相加获得B部分上所有采样点的幅值差的平方之和,记为SSDp;
(3)将SSDr和SSDp相加,即得在每次移动后CPm与FPm中各采样点的幅值的差的绝对值之和SSDi
所述肱动脉血压值包括肱动脉收缩压SBP、肱动脉舒张压DBP和肱动脉平均压MAP;
所述肱动脉平均压的计算方法为MAP=DBP+0.4*(SBP-DBP);
步骤(2)所述利用人体肱动脉血压值对所述颈动脉脉搏波和股动脉脉搏波进行标定,是指利用肱动脉舒张压DBP和肱动脉平均压MAP对所述颈动脉脉搏波和股动脉脉搏波进行标定。
步骤(1)中并不局限于获取人体颈动脉脉搏波和股动脉脉搏波,可获取颈动脉脉搏波、股动脉脉搏波及踝动脉脉搏波中的任意两个或颈动脉脉搏波、桡动脉脉搏波及踝动脉脉搏波中的任意两个,所得相应归一化后的待匹配压力波形在进行步骤(4)时按如下规则进行(遵循人体由上向下的移动规律):
若获取的是颈动脉脉搏波和桡动脉脉搏波、或颈动脉脉搏波和踝动脉脉搏波,则在步骤(4)中由归一化后的颈动脉待匹配压力波形向归一化后的桡动脉待匹配压力波形或归一化后的踝动脉待匹配压力波形移动;
若获取的是桡动脉脉搏波和踝动脉脉搏波,则在步骤(4)中由归一化后的桡动脉待匹配压力波形向归一化后的踝动脉待匹配压力波形移动;
若获取的是股动脉脉搏波和踝动脉脉搏,则在步骤(4)中由归一化后的股动脉待匹配压力波形向归一化后的踝动脉待匹配压力波形移动。
与已有技术相比,本发明的有益效果体现在:
(1)本发明选取颈动脉和股动脉压力波形上一段区域而不是单个特征点计算主动脉脉搏波传导时间aoPTT,提高了aoPWV测量的准确性;
(2)本发明在波形噪声干扰极大的情况下也具有很好的鲁棒性,很好地解决了传统的“foot-to-foot”法存在的测量精度不足的问题,具有更广的应用前景。
附图说明
图1为本发明中获取aoPTT方法的总流程图;
图2为舒张期极小值点和上升段斜率最大值点获取示意图;
图3为待匹配压力波形归一化示意图;
图4为波形匹配法求取aoPTT示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明实施例作进一步详细的说明。
如图1所示,本实施例基于波形匹配法的主动脉脉搏波传导时间aoPTT获取方法是按如下步骤进行:
1、利用双路脉搏波传感器在人体颈动脉和股动脉处以1ms为采样周期同步采集颈动脉脉搏值和股动脉脉搏值,由颈动脉脉搏值和股动脉脉搏值的时间序列构成颈动脉脉搏波和股动脉脉搏波,利用示波法血压计测量人体肱动脉收缩压(SBP)和肱动脉舒张压(DBP);
2、用肱动脉舒张压(DBP)和肱动脉平均压(MAP,MAP=DBP+0.4*(SBP-DBP))对颈、股动脉脉搏波进行标定,获得具有血压特性的颈、股动脉脉搏波。具有血压特性的颈、股动脉脉搏波即为颈动脉压力波形CP和股动脉压力波形FP。标定即是以DBP和MAP将各颈、股动脉脉搏值转换为颈、股动脉压力值,由各压力值的时间序列分别构成了颈动脉压力波形CP和股动脉压力波形FP。
3、如图2所示,选择颈动脉压力波形CP上升阶段斜率最大值点,记为cMaxPoint;选择颈动脉压力波形CP舒张期极小值点,记为cMinPoint;选择股动脉压力波形FP上升阶段斜率最大值点,记为fMaxPoint;选择股动脉压力波形FP舒张期极小值点,记为fMinPoint;以位于cMaxPoint和cMinPoint之间的波形,作为颈动脉待匹配压力波形,以位于fMaxPoint和fMinPoint之间的波形,作为股动脉待匹配压力波形;所述上升段斜率最大值点通过对颈、股动脉压力波形求取一阶导数,寻找一阶导数极大值点,该极大值点对应的位置即为上升段斜率最大值点所在的位置。所述舒张期极小值点通过选取上升段斜率极大值点,与前一个脉搏波峰值点间的区域波形作为搜索区域,搜索区域最小值点,其最小值对应的位置即为舒张期极小值点所在的位置。
4、如图3所示,由于cMinPoint和fMinPoint幅值间以及cMaxPoint和fMaxPoint幅值间有一定差距,为了使两个匹配波形的高度一致,需要对颈动脉待匹配压力波形和股动脉待匹配压力波形四个端点进行归一化处理,使得两段波形的起点和终点分别在同一水平线上;具体实现为:
a、比较cMinPoint处和fMinPoint处的幅值大小:
若cMinPoint处的幅值>fMinPoint处的幅值,则以cMinPoint处的幅值作为归一化最小幅值,以cMinPoint作为颈动脉待匹配压力波形的归一化起点,记为cStartPoint,以所述归一化最小幅值在股动脉待匹配压力波形上的位置,作为股动脉待匹配压力波形的归一化起点,记为fStartPoint;
若cMinPoint处的幅值≤fMinPoint处的幅值,则以fMinPoint处的幅值作为归一化最小幅值,以fMinPoint作为股动脉待匹配压力波形的归一化起点,记为fStartPoint,以所述归一化最小幅值在颈动脉待匹配压力波形上的位置,作为颈动脉待匹配压力波形的归一化起点,记为cStartPoint;
b、比较cMaxPoint处和fMaxPoint处的幅值大小:
若cMaxPoint处的幅值≤fMaxPoint处的幅值,则以cMaxPoint处的幅值作为归一化最大幅值,以cMaxPoint作为颈动脉待匹配压力波形的归一化终点,记为cEndPoint,以所述归一化最大幅值在股动脉待匹配压力波形上的位置,作为股动脉待匹配压力波形的归一化终点,记为fEndPoint;
若cMaxPoint处的幅值>fMaxPoint处的幅值,则以fMaxPoint处的幅值作为归一化最大幅值,以fMaxPoint作为股动脉待匹配压力波形的归一化终点,记为fEndPoint,以所述归一化最大幅值在颈动脉待匹配压力波形上的位置,作为颈动脉待匹配压力波形的归一化终点,记为cEndPoint;
c、以颈动脉待匹配压力波形中cStartPoint与cEndPoint之间的波形作为归一化后的颈动脉待匹配压力波形CPm;以股动脉待匹配压力波形中fStartPoint与fEndPoint之间的波形作为归一化后的股动脉待匹配压力波形FPm。
5、如图4所示,对归一化后的颈动脉待匹配压力波形CPm和归一化后的股动脉待匹配压力波形FPm通过波形匹配法求取aoPTT:
将归一化后的颈动脉待匹配压力波形CPm沿着时间轴向归一化后的股动脉待匹配压力波形FPm移动,每次移动1个步骤(1)所用的采样周期,共移动的时间为T,所述T为cStartPoint距离fEndPoint之间的时间,共移动的次数m=T/采样周期;
在每次移动后计算CPm上各采样点与FPm中相应采样点的幅值差的绝对值之和SADi或幅值差的平方之和SSDi,其中i为移动次数,i=1,2,、、、,m;
令SADi或SSDi的最小值所对应的移动次数为p,则传导时间aoPTT=p×采样周期。
在每次移动后CPm上各采样点与FPm中相应采样点的幅值差的绝对值之和SADi是按如下方法进行计算:
(1)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于相同时间段上的A部分(即在时间轴上重叠的部分),计算获得位于A部分的CPm上各采样点与FPm上相同时间点处的采样点的幅值差的绝对值,然后计算由A部分CPm上所有采样点所获得的幅值差的绝对值之和,记为SADr;
(2)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于不同时间段上的B部分(即在时间轴上不重叠的部分),计算获得位于B部分的CPm上各采样点的幅值与cStartPoint处幅值差的绝对值之和SADpc,和位于B部分的的FPm上各采样点的幅值与fEndPoint处幅值差的绝对值之和SADpf,将SADpc和SADpf相加获得B部分上所有采样点的幅值差的绝对差值之和,记为SADp;
(3)将SADr和SADp相加,即得在每次移动后CPm与FPm中各采样点的幅值的差的绝对值之和SADi
在每次移动后CPm上各采样点与FPm中相应采样点的幅值差的平方之和SSDi是按如下方法进行计算:
(1)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于相同时间段上的A部分(即在时间轴上重叠的部分),计算获得位于A部分的CPm上各采样点与FPm上相同时间点处的采样点的幅值差的平方,然后计算由A部分CPm上所有采样点所获得的幅值差的平方之和,记为SSDr;
(2)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于不同时间段上的B部分(即在时间轴上不重叠的部分),计算获得位于B部分的CPm上各采样点的幅值与cStartPoint处幅值差的平方之和SSDpc,和位于B部分的的FPm上各采样点的幅值与fEndPoint处幅值差的平方之和SSDpf,将SADpc和SADpf相加获得B部分上所有采样点的幅值差的平方之和,记为SSDp;
(3)将所述的SSDr和SSDp相加,即得在每次移动后CPm与FPm中各采样点的幅值的差的绝对值之和SSDi
步骤1中并不局限于获取人体颈动脉脉搏波和股动脉脉搏波,可获取颈动脉脉搏波、股动脉脉搏波及踝动脉脉搏波中的任意两个或颈动脉脉搏波、桡动脉脉搏波及踝动脉脉搏波中的任意两个,所得相应归一化后的待匹配压力波形在进行步骤5时按如下规则进行(遵循人体由上向下的移动规律):
若获取的是颈动脉脉搏波和桡动脉脉搏波、或颈动脉脉搏波和踝动脉脉搏波,则在步骤(4)中由归一化后的颈动脉待匹配压力波形向归一化后的桡动脉待匹配压力波形或归一化后的踝动脉待匹配压力波形移动;
若获取的是桡动脉脉搏波和踝动脉脉搏波,则在步骤(4)中由归一化后的桡动脉待匹配压力波形向归一化后的踝动脉待匹配压力波形移动;
若获取的是股动脉脉搏波和踝动脉脉搏,则在步骤(4)中由归一化后的股动脉待匹配压力波形向归一化后的踝动脉待匹配压力波形移动。

Claims (4)

1.一种基于波形匹配法的主动脉脉搏波传导时间获取方法,其特征在于包括如下步骤:
(1)利用脉搏波传感器同步采样获取人体颈动脉脉搏波和股动脉脉搏波;
(2)利用人体肱动脉血压值对所述颈动脉脉搏波和股动脉脉搏波进行标定,分别获得具有血压特性的颈动脉脉搏波和具有血压特性的股动脉脉搏波;所述具有血压特性的颈动脉脉搏波即为颈动脉压力波形CP,所述具有血压特性的股动脉脉搏波即为股动脉压力波形FP;
(3)选择颈动脉压力波形CP上升阶段斜率最大值点,记为cMaxPoint;选择颈动脉压力波形CP舒张期极小值点,记为cMinPoint;选择股动脉压力波形FP上升阶段斜率最大值点,记为fMaxPoint;选择股动脉压力波形FP舒张期极小值点,记为fMinPoint;以位于cMaxPoint和cMinPoint之间的波形,作为颈动脉待匹配压力波形,以位于fMaxPoint和fMinPoint之间的波形,作为股动脉待匹配压力波形,对所述颈动脉待匹配压力波形和所述股动脉待匹配压力波形作如下归一化处理:
3a、比较cMinPoint处和fMinPoint处的幅值大小:
若cMinPoint处的幅值>fMinPoint处的幅值,则以cMinPoint处的幅值作为归一化最小幅值,以cMinPoint作为颈动脉待匹配压力波形的归一化起点,记为cStartPoint,以所述归一化最小幅值在股动脉待匹配压力波形上的位置,作为股动脉待匹配压力波形的归一化起点,记为fStartPoint;
若cMinPoint处的幅值≤fMinPoint处的幅值,则以fMinPoint处的幅值作为归一化最小幅值,以fMinPoint作为股动脉待匹配压力波形的归一化起点,记为fStartPoint,以所述归一化最小幅值在颈动脉待匹配压力波形上的位置,作为颈动脉待匹配压力波形的归一化起点,记为cStartPoint;
3b、比较cMaxPoint处和fMaxPoint处的幅值大小:
若cMaxPoint处的幅值≤fMaxPoint处的幅值,则以cMaxPoint处的幅值作为归一化最大幅值,以cMaxPoint作为颈动脉待匹配压力波形的归一化终点,记为cEndPoint,以所述归一化最大幅值在股动脉待匹配压力波形上的位置,作为股动脉待匹配压力波形的归一化终点,记为fEndPoint;
若cMaxPoint处的幅值>fMaxPoint处的幅值,则以fMaxPoint处的幅值作为归一化最大幅值,以fMaxPoint作为股动脉待匹配压力波形的归一化终点,记为fEndPoint,以所述归一化最大幅值在颈动脉待匹配压力波形上的位置,作为颈动脉待匹配压力波形的归一化终点,记为cEndPoint;
3c、以颈动脉待匹配压力波形中cStartPoint与cEndPoint之间的波形作为归一化后的颈动脉待匹配压力波形CPm;以股动脉待匹配压力波形中fStartPoint与fEndPoint之间的波形作为归一化后的股动脉待匹配压力波形FPm;
(4)对归一化后的颈动脉待匹配压力波形CPm和归一化后的股动脉待匹配压力波形FPm通过波形匹配法求取主动脉脉搏波传导时间aoPTT:
将归一化后的颈动脉待匹配压力波形CPm沿着时间轴向归一化后的股动脉待匹配压力波形FPm移动,每次移动1个步骤(1)所用的采样周期,共移动的时间为T,所述T为cStartPoint距离fEndPoint之间的时间,共移动的次数m=T/采样周期;
在每次移动后计算CPm上各采样点与FPm中相应采样点的幅值差的绝对值之和SADi或幅值差的平方之和SSDi,其中i为移动次数,i=1,2,…,m;
令SADi或SSDi的最小值所对应的移动次数为p,则主动脉脉搏波传导时间aoPTT=p×采样周期。
2.根据权利要求1所述的基于波形匹配法的主动脉脉搏波传导时间获取方法,其特征在于:在每次移动后CPm上各采样点与FPm中相应采样点的幅值差的绝对值之和SADi是按如下方法进行计算:
(1)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于相同时间段上的A部分,计算获得位于A部分的CPm上各采样点与FPm上相同时间点处的采样点的幅值差的绝对值,然后计算由A部分CPm上所有采样点所获得的幅值差的绝对值之和,记为SADr;
(2)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于不同时间段上的B部分,计算获得位于B部分的CPm上各采样点的幅值与cStartPoint处幅值差的绝对值之和SADpc,和位于B部分的的FPm上各采样点的幅值与fEndPoint处幅值差的绝对值之和SADpf,将SADpc和SADpf相加获得B部分上所有采样点的幅值差的绝对差值之和,记为SADp;
(3)将SADr和SADp相加,即得在每次移动后CPm与FPm中各采样点的幅值差的绝对值之和SADi
3.根据权利要求1所述的基于波形匹配法的主动脉脉搏波传导时间获取方法,其特征在于:在每次移动后CPm上各采样点与FPm中相应采样点的幅值差的平方之和SSDi是按如下方法进行计算:
(1)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于相同时间段上的A部分,计算获得位于A部分的CPm上各采样点与FPm上相同时间点处的采样点的幅值差的平方,然后计算由A部分CPm上所有采样点所获得的幅值差的平方之和,记为SSDr;
(2)对归一化后的颈动脉待匹配压力波形CPm与归一化后的股动脉待匹配压力波形FPm位于不同时间段上的B部分,计算获得位于B部分的CPm上各采样点的幅值与cStartPoint处幅值差的平方之和SSDpc,和位于B部分的的FPm上各采样点的幅值与fEndPoint处幅值差的平方之和SSDpf,将SADpc和SADpf相加获得B部分上所有采样点的幅值差的平方之和,记为SSDp;
(3)将SSDr和SSDp相加,即得在每次移动后CPm与FPm中各采样点的幅值的差的绝对值之和SSDi
4.根据权利要求1所述的基于波形匹配法的主动脉脉搏波传导时间获取方法,其特征在于:所述肱动脉血压值包括肱动脉收缩压SBP、肱动脉舒张压DBP和肱动脉平均压MAP;
所述肱动脉平均压的计算方法为MAP=DBP+0.4*(SBP-DBP);
步骤(2)所述利用人体肱动脉血压值对所述颈动脉脉搏波和股动脉脉搏波进行标定,是指利用肱动脉舒张压DBP和肱动脉平均压MAP对所述颈动脉脉搏波和股动脉脉搏波进行标定。
CN201410493975.2A 2014-09-24 2014-09-24 一种基于波形匹配法的主动脉脉搏波传导时间获取方法 Active CN104188637B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410493975.2A CN104188637B (zh) 2014-09-24 2014-09-24 一种基于波形匹配法的主动脉脉搏波传导时间获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410493975.2A CN104188637B (zh) 2014-09-24 2014-09-24 一种基于波形匹配法的主动脉脉搏波传导时间获取方法

Publications (2)

Publication Number Publication Date
CN104188637A CN104188637A (zh) 2014-12-10
CN104188637B true CN104188637B (zh) 2016-02-24

Family

ID=52074145

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410493975.2A Active CN104188637B (zh) 2014-09-24 2014-09-24 一种基于波形匹配法的主动脉脉搏波传导时间获取方法

Country Status (1)

Country Link
CN (1) CN104188637B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110279404A (zh) * 2019-06-17 2019-09-27 清华大学 一种脉搏波采集装置及方法

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104757959B (zh) * 2015-04-20 2017-05-03 上海工程技术大学 一种基于影像叠影的脉搏波传输速率检测方法和系统
JP7138363B2 (ja) * 2018-04-28 2022-09-16 深▲せん▼市大耳馬科技有限公司 脈波伝導パラメータの測定システム、測定方法及び測定装置
CN109480802B (zh) * 2018-12-29 2021-06-25 中国科学院合肥物质科学研究院 一种基于波形分析技术的血压参数估计系统及方法
CN111265196B (zh) * 2020-03-04 2021-03-09 云南大学 一种颈动脉脉搏波速估计方法及系统
CN115990002B (zh) * 2023-03-21 2024-03-19 首都医科大学宣武医院 一种生命体征监测系统及方法
CN115990001B (zh) * 2023-03-21 2024-04-05 首都医科大学宣武医院 一种可穿戴监护系统、可穿戴装置及存储介质

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101361652B (zh) * 2007-08-07 2012-06-27 深圳迈瑞生物医疗电子股份有限公司 自适应阈值的波形峰谷识别方法和装置
US8313439B2 (en) * 2009-03-20 2012-11-20 Massachusetts Institute Of Technology Calibration of pulse transit time measurements to arterial blood pressure using external arterial pressure applied along the pulse transit path
US9220437B2 (en) * 2010-05-14 2015-12-29 Centre For Development Of Advanced Computing Diagnostic classifications of pulse signal waveform data
CN102197996B (zh) * 2010-11-23 2013-04-17 中国科学院合肥物质科学研究院 利用桡动脉脉搏波获取主动脉脉搏波速度的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110279404A (zh) * 2019-06-17 2019-09-27 清华大学 一种脉搏波采集装置及方法
CN110279404B (zh) * 2019-06-17 2020-07-28 清华大学 一种脉搏波采集装置及方法

Also Published As

Publication number Publication date
CN104188637A (zh) 2014-12-10

Similar Documents

Publication Publication Date Title
CN104188637B (zh) 一种基于波形匹配法的主动脉脉搏波传导时间获取方法
CN107736880B (zh) 一种脉象分析方法及系统
CN107822608B (zh) 基于高斯混合模型的脉搏波特征提取方法
CN102270264B (zh) 生理信号质量评估系统及方法
WO2017024457A1 (zh) 连续血压测量装置、测量模型建立方法和系统
CN104116503B (zh) 一种无创连续血压的测量装置
CN107432736B (zh) 一种识别脉搏波形信号的方法
CN101732040B (zh) 一种无创多路脉搏波检测方法
CN109124610B (zh) 一种无创血压测量的抗干扰方法和装置
CN105286815A (zh) 一种基于波形时域特征的脉搏波信号特征点检测方法
CN103584854B (zh) 心电信号r波的提取方法
CN103284703B (zh) 基于上肢动脉信息的主动脉脉搏波传递时间测量方法
CN104757955A (zh) 一种基于脉搏波的人体血压预测方法
CN104068841B (zh) 一种测量心脏收缩时间参数的测量方法及装置
CN109288508A (zh) 一种基于crnn-bp的血压值智能测量方法
CN105942980B (zh) 一种听诊法血压计特征提取及分类的方法和系统
WO2023226293A1 (zh) 基于超声多普勒的连续血压测量装置和电子设备
CN105147269A (zh) 一种无创连续血压测量方法
CN104644151B (zh) 一种基于光电容积脉搏信号的压力脉搏波波形传播预测方法
CN107049289A (zh) 基于桡动脉波的无袖带血压测量计和测量方法
CN110840428B (zh) 基于一维U-Net网络的无创血压估计方法
CN111839488A (zh) 基于脉搏波的无创连续血压测量装置和方法
CN103393415A (zh) 连续血压变化的测量方法
CN208677373U (zh) 一种脉象分析系统
CN104983412A (zh) 一种中心动脉收缩期平均归一化血流波形模型及基于其获取主动脉脉搏波传递时间的方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant