CN115990007B - 中心动脉压波形拟合方法、监控装置和手表设备 - Google Patents
中心动脉压波形拟合方法、监控装置和手表设备 Download PDFInfo
- Publication number
- CN115990007B CN115990007B CN202211529875.1A CN202211529875A CN115990007B CN 115990007 B CN115990007 B CN 115990007B CN 202211529875 A CN202211529875 A CN 202211529875A CN 115990007 B CN115990007 B CN 115990007B
- Authority
- CN
- China
- Prior art keywords
- pulse
- waveform
- pulse pressure
- mapping relation
- pressure
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000004872 arterial blood pressure Effects 0.000 title claims abstract description 48
- 238000012806 monitoring device Methods 0.000 title claims abstract description 14
- 230000035485 pulse pressure Effects 0.000 claims abstract description 168
- 238000013507 mapping Methods 0.000 claims abstract description 132
- 238000001228 spectrum Methods 0.000 claims abstract description 59
- 230000003595 spectral effect Effects 0.000 claims abstract description 30
- 230000000747 cardiac effect Effects 0.000 claims abstract description 24
- 210000002321 radial artery Anatomy 0.000 claims description 43
- 238000001514 detection method Methods 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 11
- 238000012216 screening Methods 0.000 claims description 9
- 238000003860 storage Methods 0.000 claims description 8
- 210000000707 wrist Anatomy 0.000 claims description 7
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 18
- 230000036772 blood pressure Effects 0.000 abstract description 11
- 238000005284 basis set Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 9
- 230000035487 diastolic blood pressure Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 8
- 230000035488 systolic blood pressure Effects 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 239000000758 substrate Substances 0.000 description 5
- 206010021143 Hypoxia Diseases 0.000 description 4
- 230000007954 hypoxia Effects 0.000 description 4
- 230000017074 necrotic cell death Effects 0.000 description 4
- 210000001367 artery Anatomy 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 210000002302 brachial artery Anatomy 0.000 description 2
- 210000001715 carotid artery Anatomy 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 210000000709 aorta Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000036581 peripheral resistance Effects 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000010992 reflux Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
- 201000002282 venous insufficiency Diseases 0.000 description 1
Classifications
-
- 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
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提供一种中心动脉压波形拟合方法、监控装置和手表设备,属于血压监测技术领域,方法包括:获取待测用户的第一脉搏波形;基于第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集,第一脉压映射关系根据获取的第一脉搏样本波形与获取的第二脉搏样本波形的脉压映射关系确定;基于第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;第二脉压映射关系根据中心动脉压波形与第二脉搏波形的脉压映射关系确定。本发明根据第一脉压映射关系和第二脉压映射关系,只需采集待测用户的第一脉搏波形,即可拟合得到中心动脉压波形,无需持续施加压力,不会阻断静脉回流,适合持续监控,方便制成可穿戴设备。
Description
技术领域
本发明涉及血压监测技术领域,尤其涉及一种中心动脉压波形拟合方法、监控装置和手表设备。
背景技术
中心动脉压指的是与心脏靠近的动脉血压,是心血管风险的一个重要预测指标。评估中心动脉压的方法包括侵入式方法和非侵入式方法,侵入式方法包括心导管插入手术和使用压力传感导管记录升主动脉的血压。然而侵入式方法的技术要求高,不适用于普通人群的常规筛查。非侵入式方法指的是从主动脉远端部位推导中心动脉压,比如颈动脉、桡动脉或肱动脉,其中,桡动脉血压采集的是单动脉压力波,是非侵入式量测方法中比较常用的方法。
目前中心动脉压的检测装置大多设置成手持式压平式压力传感设备,主要以一次性量测为主,长期使用会阻断静脉回流,导致组织缺氧、坏死,无法实现持续监控。
发明内容
本发明提供一种中心动脉压波形拟合方法、监控装置和手表设备,用以解决现有中心动脉压的检测装置大多设置成手持式压平式压力传感设备,主要以一次性量测为主,长期使用会阻断静脉回流,导致组织缺氧、坏死,无法实现持续监控的缺陷。
本发明提供一种中心动脉压波形拟合方法,包括:
通过第一检测法获取待测用户的第一脉搏波形;
基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集;其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
根据本发明提供的一种中心动脉压波形拟合方法,所述第一脉压映射关系是采用以下方法确定的:
对所述第一脉搏样本波形进行频域变换,确定第三频谱特征基底集;
对所述第二脉搏样本波形进行频域变换,确定第四频谱特征基底集;
根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系。
根据本发明提供的一种中心动脉压波形拟合方法,所述根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系,包括:
获取所述第三频谱特征基底集中在基频和/或多个倍频上的频谱特征基底的第一线性组合值,以及所述第四频谱特征基底集中在所述基频和/或所述多个倍频上的频谱特征基底的第二线性组合值;
基于所述第一线性组合值与所述第二线性组合值,确定所述第一脉压映射关系。
根据本发明提供的一种中心动脉压波形拟合方法,在基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形之前,还包括:
采集所述第二脉搏样本波形;
根据所述第二脉搏样本波形,获取动脉特征参数集;
利用所述动脉特征参数集对所述第二脉压映射关系进行校准。
根据本发明提供的一种中心动脉压波形拟合方法,所述利用所述动脉特征参数集对所述第二脉压映射关系进行校准,包括:
基于所述待测用户的与历史中心脉压波形,从所述动脉特征参数集中筛选出有效动脉特征参数;
利用所述有效动脉特征参数对所述第二脉压映射关系进行校准。
根据本发明提供的一种中心动脉压波形拟合方法,在所述第二脉搏波形为桡动脉波形的情况下,在从所述动脉特征参数集中筛选出有效动脉特征参数之前,还包括:
基于所述待测用户的历史桡动脉参数特征数据库,和/或所述待测用户的参考桡动脉参数特征,对所述动脉特征参数集进行优化,以剔除异常动脉特征参数;
所述历史桡动脉参数特征数据库是基于所述待测用户的多个历史桡动脉波形构建的,所述参考桡动脉参数特征是基于血压计所测量到的所述待测用户的任一历史桡动脉波形构建的。
根据本发明提供的一种中心动脉压波形拟合方法,还包括:
根据所述第二频谱特征基底集,得到第二脉搏波形。
根据本发明提供的一种中心动脉压波形拟合方法,所述第一检测法为光电容脉搏波描记法,所述第二脉搏波形用于表征采用示波法测得的脉搏波形。
本发明还提供一种中心动脉压的监控装置,包括:
脉搏波形采集模块,用于运行第一检测法获取待测用户的第一脉搏波形;
特征映射模块,用于基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集;其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
波形拟合模块,用于基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
本发明还提供一种手表设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述程序时实现所述中心动脉压波形拟合方法。
本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述中心动脉压波形拟合方法。
本发明提供的中心动脉压波形拟合方法、监控装置和手表设备,根据第一脉压映射关系和第二脉压映射关系,只需采集待测用户的第一脉搏波形,即可拟合得到中心动脉压波形,不会对待测用户持续施加压力,不会阻断静脉回流,适合持续监控,方便制成可穿戴设备。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的中心动脉压波形拟合方法的流程示意图;
图2是本发明实施例提供的第一脉压映射关系的确定流程示意图;
图3是本发明实施例提供的第一脉搏样本波形的波形示意图;
图4是本发明实施例提供的第二脉搏样本波形的波形示意图;
图5是本发明实施例提供的对第二脉压映射关系进行校准的流程示意图之一;
图6是本发明实施例提供的对第二脉压映射关系进行校准的流程示意图之二;
图7是本发明实施例提供的中心动脉压的监控装置的原理框图;
图8是本发明实施例提供的电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本申请实施例中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”可是第一特征在第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”可以是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度小于第二特征。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本申请实施例的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
图1是本发明实施例提供的中心动脉压波形拟合方法的流程示意图,参照图1,本发明提供一种中心动脉压波形拟合方法,包括:
S110,通过第一检测法获取待测用户的第一脉搏波形;
S120,基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集;其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
S130,基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
步骤S110~步骤S130用于实现中心动脉压的连续监控,通过获取待测用户的第一脉搏波形,拟合得到中心动脉压波形,以实现对中心动脉压的监控。
可选的,所述第一脉压映射关系和所述第二脉压映射关系是预先确定的。所述第一脉压映射关系是第一脉搏样本波形与第二脉搏样本波形在频域上的脉压映射关系,所述第二脉压映射关系是所述中心动脉压波形与第二脉搏波形在频域上的脉压映射关系。
可选的,第一脉搏波形和第一脉搏样本波形为待测用户手腕处的脉搏波信号;第二脉搏波形和第二脉搏样本波形可以为颈动脉、桡动脉或肱动脉处的脉搏波信号。
可选的,第一检测法可以为光电检测法或者磁感应检测法。
可以理解的是,本发明实施例提供的中心动脉压波形拟合方法,根据第一脉压映射关系和第二脉压映射关系,只需采集待测用户的第一脉搏波形,即可拟合得到中心动脉压波形,不会对待测用户持续施加压力,不会阻断静脉回流,适合持续监控,方便制成可穿戴设备。
在上述实施例的基础上,作为一个可选的实施例,所述第一脉压映射关系是采用以下方法确定的:
对所述第一脉搏样本波形进行频域变换,确定第三频谱特征基底集;
对所述第二脉搏样本波形进行频域变换,确定第四频谱特征基底集;
根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系。
可选的,可通过傅里叶变换或者小叶变换,对第一脉搏样本波形、第二脉搏样本波形进行频域变换,得到表征频谱特征的垂直正交基底,根据垂直正交基底得到对应的第三频谱特征基底集和第四频谱特征基底集。
可以理解的是,本发明实施例通过第三频谱特征基底集和第四频谱特征基底集确定第一脉压映射关系,因第三频谱特征基底集和第四频谱特征基底集为频域特征参数,能够涵盖第一脉搏样本波形与第二脉搏样本波形90%以上的波形能量,提高了第一脉压映射关系的准确性。
图2是本发明实施例提供的第一脉压映射关系的确定流程示意图;参照图2,在上述实施例的基础上,作为一个可选的实施例,所述根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系,包括:
S210,获取所述第三频谱特征基底集中在基频和/或多个倍频上的频谱特征基底的第一线性组合值,以及所述第四频谱特征基底集中在基频和/或多个倍频上的频谱特征基底的第二线性组合值;
S220,基于所述第一线性组合值与所述第二线性组合值,确定所述第一脉压映射关系。
图3是本发明实施例提供的第一脉搏样本波形的波形示意图;参照图3,可选的,预先通过第一检测法获取待测用户手腕处10~60秒的脉搏波信号,从脉搏波信号中选取具有代表性的第一脉搏样本波形,记为X[k],k值为人体脉搏波形参数,是诊断心血管外周阻力和血管壁硬化程度的一个重要指标。设傅里叶变换或小叶变换的基底函数为ψωn(t),对第一脉搏样本波形进行傅里叶变换或小叶变换,得到一系列垂直正交基底与系数An,第三频谱特征基底集为/>其中n=1,2,...,α,α为正整数,需满足涵盖第一脉搏样本波形能量的90%及以上。
图4是本发明实施例提供的第二脉搏样本波形的波形示意图;参照图4,可选的,预先通过第二检测法获取待测用户的桡动脉波形,从中选取具有代表性的第二脉搏样本波形,记为Y[s],s为桡动脉的参数指标。设傅里叶变换或小叶变换的基底函数为ξωn(t),对第二脉搏样本波形进行傅里叶变换或小叶变换,得到一系列χn垂直正交基底与系数Bn,第四频谱特征基底集合为{χn,Bn},其中n=1,2,...,α,α为正整数,需满足涵盖桡动脉脉搏波形能量的90%及以上。
可选的,可通过气囊在待测用户手腕处从0mmHg开始加压,示波法获取待测用户的桡动脉波形,加压到桡动脉波形峰-峰值为最大脉搏波的1/3以下为止。现有的连续监控设备的施加压力介于舒张压和收缩压之间,约100mmHg,会阻断静脉回流,长时间会造成组织缺氧、坏死。而本发明实施例施加的最大压力为最大脉搏波值的1/3,且仅施加一次压力,不会造成组织缺氧、坏死,更适合做穿戴设备,长时间穿戴,持续监控使用者的中心动脉压。
可选的,通过压力传感器检测对待监控人员目标位置施加的压力数值,压力传感器的压力响应范围为0~300mmHg,非线性度<0.7%。
可选的,代表性的第二脉搏样本波形与平均血压的差值为-10~12mmHg的脉搏波形,其中,平均血压基于示波法确定。
可选的,下面对所述第一脉压映射关系的计算公式进行举例说明,但并不局限于以下实施例。
其中,U、U1、U2、W、W1、W2为组合系数,上述组合系数可以为复变的参数。n为母函数的n倍倍频或谐频。
可选的,所述第一脉压映射关系可以简化为如下公式:
可选的,根据待测用户的个性化特征(年龄、性格、有无病史等),可从公式1~3/3.1中选择最适配的第一脉压映射关系,用以提高中心动脉压波形拟合的准确度。
可以理解的是,本发明实施例通过对所述第一线性组合值与所述第二线性组合值,确定所述第一脉压映射关系在所述基频和/或预设倍频上的映射关系,提供了一种第一脉压映射关系的确定方案,能够提高第一脉压映射关系的准确性,为后续中心动脉压波形的拟合提供便利。
图5是本发明实施例提供的对第二脉压映射关系进行校准的流程示意图之一;参照图5,在上述实施例的基础上,作为一个可选的实施例,在基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形之前,还包括:
S510,采集所述第二脉搏样本波形;
S520,根据所述第二脉搏样本波形,获取动脉特征参数集;
S530,利用所述动脉特征参数集对所述第二脉压映射关系进行校准。
可选的,第二脉压映射关系可采用一般化转换函数(Generalized Transferfunction,GTF(n))表示,为了提高第二脉压映射关系的准确性,本发明可根据所述动脉特征参数集对所述第二脉压映射关系进行校准。
可选的,基于示波法获取待测用户桡动脉的舒张压、收缩压、平均血压等动脉特征参数集。其中,设SBPi,DBPi,PPi分别表示基于示波法第i次计算出的桡动脉收缩压、桡动脉舒张压、桡动脉脉压的参数值。
结合待测用户输入的商用血压计取得的桡动脉收缩压SBPnew、桡动脉舒张压DBPnew、桡动脉脉压PPnew,得到滚动的代表性的桡动脉舒张压SBP’、收缩压DBP’、平均血压MAP’、动脉脉压PP’,公式如下:
SBP′=fr*f(SBPi,DBPi,PPi)+(1-fr)SBPnew (4)
DBP′=fr*f(SBPi,DBPi,PPi)+(1-fr)DBPnew (5)
PP′=fr*f(SBPi,DBPi,PPi)+(1-fr)PPnew (6)
其中,fr为设定的计算系数。
将SBP’、DBP’、PP’进行加减、比例或者线性组合得到第二脉压映射关系的校正系数ModifedFactor。
可以理解的是,本发明实施例通过对所述第二脉压映射关系进行校准,能够提高中心动脉压波形拟合的准确度。
图6是本发明实施例提供的对第二脉压映射关系进行校准的流程示意图之二;参照图6,在上述实施例的基础上,作为一个可选的实施例,所述利用所述动脉特征参数集对所述第二脉压映射关系进行校准,还包括:
S610,基于所述待测用户的历史第二脉搏波形与历史中心脉压波形,从所述动脉特征参数集中筛选出有效动脉特征参数;
S620,利用所述有效动脉特征参数对所述第二脉压映射关系进行校准。
可选的,根据待测用户的历史第二脉搏波形与历史中心脉压波形以及波形特征参数的取值,从所述动脉特征参数集中筛选出符合取值的有效动脉特征参数;比如,代表性的第二脉搏样本波形与平均血压的差值为-10~12mmHg的脉搏波形。
可以理解的是,本发明实施例通过设置筛选的步骤,能够有效提高动脉特征参数的准确性,进而提高了校准参数的准确性。
在上述实施例的基础上,作为一个可选的实施例,在所述第二脉搏波形为桡动脉波形的情况下,在从所述动脉特征参数集中筛选出有效动脉特征参数之前,还包括:
基于所述待测用户的历史桡动脉参数特征数据库,和/或所述待测用户的参考桡动脉参数特征,对所述动脉特征参数集进行优化,以剔除异常动脉特征参数;
所述历史桡动脉参数特征数据库是基于所述待测用户的多个历史桡动脉波形构建的,所述参考桡动脉参数特征是基于血压计所测量到的所述待测用户的任一历史桡动脉波形构建的。
可选的,将示波法采集到的桡动脉的舒张压、收缩压、平均血压等多参数特征,基于窗口的移动平均法或权重平均法与所述待测用户的历史桡动脉参数特征数据库进行比对,计算公式如下:
其中,SBPi、DBPi、PPi分别表示第i次计算出的桡动脉收缩压、舒张压、脉压。
可选的,校准后的第二脉压映射关系可用个人化中心动脉压转换函数(Personalized Transfer function,简称PTF(n)表示,计算公式如下:
PTF(n)=ff(Trans(n),GTF(n),ModifedFactor) (8)
其中,ff为设定的映射关系。
可选的,校准后的第二脉压映射关系的计算公式还可简化为如下公式:
PTF(n)=Trans(n)*GTF(n)*ModifedFactor (9)
可以理解的是,本发明实施例通过对所述动脉特征参数集进行优化,有利于提高第二映射关系的准确度。
在上述实施例的基础上,作为一个可选的实施例,在获取第二频谱特征基底集之后,还包括:
根据所述第二频谱特征基底集,得到第二脉搏波形。
可选的,可对所述第二频谱特征基底集进行傅里叶反变换或者小波反变换,得到第二脉搏波形。
可以理解的是,本发明实施例通过得到第二脉搏波形,以便确定计算的准确性,若是第二脉搏波形有误,则可省略后续的拟合步骤,减少计算成本。
在上述实施例的基础上,作为一个可选的实施例,所述第一检测法为光电容脉搏波描记法,所述第二脉搏波形用于表征采用示波法测得的脉搏波形。
可选的,光电容积法使用的LED光波长范围在660~1500nm范围内,LED光功耗电流范围在10~520mA。
可选的,第一脉搏波形的采样频率大于200Hz,采样时长大于10秒。
设采集到的第一脉搏波形为Xnew[k],对第一脉搏波形进行傅里叶变换或者小波变换,生成所述第一脉搏波形对应的第一频谱特征基底集
可以只从第一频谱特征基底集中选择几个重要特征,降低计算量,提高计算速度。比如n=1,2,4,7,透过公式3.1,公式4~7,计算出当下的TransS(n),ModifedFactor。
根据TransS(n),ModifedFactor,计算得到校准后的第二脉压映射关系,拟合得到中心动脉压波形CentralBP[k],计算公式如下:
再从其中的最高最低点得到中心动脉收缩压CentralSBP、中心动脉舒张压CentralDBP,中心动脉平均血压CentralMBP。
与现有的手持式压平式压力传感设备需要专业操作,只能一次性诊断和测量不同,基于本发明实施例提供的方法制成的监控设备,能够实现长时间持续监控,只需一次校准,在持续监控步骤中无需再对人体施压,仅以光电法获取脉搏波波形即可获取准确的,具有个人化特征的中心动脉压波形和中心动脉压特征参数。
采用上述任一实施例制成的穿戴设备,均可以实现一次校准,长时间持续监控的作用,在持续监控过程中,仅需光电容积法获取待监控人员的脉搏波波形,功耗比较小,可以达到10~100mW。另外,监控过程中无需再对人体施加压力,不会阻断静脉回流,对人体健康更加友好,长时间佩戴无压迫负担。
下面对本发明提供的中心动脉压的监控装置进行描述,下文描述的中心动脉压的监控装置与上文描述的中心动脉压波形拟合方法可相互对应参照。
图7是本发明提供的中心动脉压的监控装置的原理框图,参照图7,本发明还提供一种中心动脉压的监控装置,包括:
脉搏波形采集模块710,用于运行第一检测法获取待测用户的第一脉搏波形;
特征映射模块720,用于基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集;其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
波形拟合模块730,用于基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
在一个实施例中,还包括:
脉压确定模块,用于对所述第一脉搏样本波形进行频域变换,确定第三频谱特征基底集;对所述第二脉搏样本波形进行频域变换,确定第四频谱特征基底集;根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系。
在一个实施例中,所述脉压确定模块还用于:
获取所述第三频谱特征基底集中在基频和/或多个倍频上的频谱特征基底的第一线性组合值,以及所述第四频谱特征基底集中在基频和/或多个倍频上的频谱特征基底的第二线性组合值;
基于所述第一线性组合值与所述第二线性组合值,确定所述第一脉压映射关系。
在一个实施例中,还包括:
校准模块,用于采集待测用户的第二脉搏样本波形;根据所述第二脉搏样本波形,获取动脉特征参数集;利用所述动脉特征参数集对所述第二脉压映射关系进行校准。
在一个实施例中,所述校准模块还用于:
基于所述待测用户的历史第二脉搏波形与历史中心脉压波形,从所述动脉特征参数集中筛选出有效动脉特征参数;
利用所述有效动脉特征参数对所述第二脉压映射关系进行校准。
在一个实施例中,所述校准模块还用于:
基于所述待测用户的历史桡动脉参数特征数据库,和/或所述待测用户的参考桡动脉参数特征,对所述动脉特征参数集进行优化,以剔除异常动脉特征参数;
所述历史桡动脉参数特征数据库是基于所述待测用户的多个历史桡动脉波形构建的,所述参考桡动脉参数特征是基于血压计所测量到的所述待测用户的任一历史桡动脉波形构建的。
在一个实施例中,所述波形拟合模块730还用于:
根据所述第二频谱特征基底集,得到第二脉搏波形。
在一个实施例中,所述第一检测法为光电容脉搏波描记法,所述第二脉搏波形用于表征采用示波法测得的脉搏波形。
本发明还提供一种手表设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上述任一种所述中心动脉压的监控方法。
图8示例了一种手表设备的实体结构示意图,如图8所示,该手表设备可以包括:处理器(processor)810、通信接口(Communications Interface)820、存储器(memory)830和通信总线840,其中,处理器88,通信接口820,存储器830通过通信总线840完成相互间的通信。处理器88可以调用存储器830中的逻辑指令,以执行中心动脉压波形拟合方法,包括:
通过第一检测法获取待测用户的第一脉搏波形;
基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集,其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
此外,上述的存储器830中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面,本发明还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,计算机程序可存储在非暂态计算机可读存储介质上,所述计算机程序被处理器执行时,计算机能够执行上述各方法所提供的中心动脉压波形拟合方法,包括:
通过第一检测法获取待测用户的第一脉搏波形;
基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集;其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
又一方面,本发明还提供一种非暂态计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现以执行上述各方法提供的中心动脉压波形拟合方法,包括:
通过第一检测法获取待测用户的第一脉搏波形;
基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集,其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形,其中,所述第二脉压映射关系是根据中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种中心动脉压波形拟合方法,其特征在于,包括:
通过第一检测法获取待测用户的第一脉搏波形;
基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集,其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形,其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的;
在基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形之前,还包括:
采集所述第二脉搏样本波形;
根据所述第二脉搏样本波形,获取动脉特征参数集;
利用所述动脉特征参数集对所述第二脉压映射关系进行校准,对动脉特征参数进行加减、比例或者线性组合得到第二脉压映射关系的校正系数;
预先通过第一检测法获取待测用户手腕处的脉搏波信号,从脉搏波信号中选取具有代表性的第一脉搏样本波形,其中,通过气囊在待测用户手腕处从零开始加压,获取待测用户的桡动脉波形,加压到桡动脉波形峰-峰值为最大脉搏波的1/3以下为止。
2.根据权利要求1所述的中心动脉压波形拟合方法,其特征在于,所述第一脉压映射关系是采用以下方法确定的:
对所述第一脉搏样本波形进行频域变换,确定第三频谱特征基底集;
对所述第二脉搏样本波形进行频域变换,确定第四频谱特征基底集;
根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系。
3.根据权利要求2所述的中心动脉压波形拟合方法,其特征在于,所述根据所述第三频谱特征基底集和所述第四频谱特征基底集,确定所述第一脉压映射关系,包括:
获取所述第三频谱特征基底集中在基频和/或多个倍频上的频谱特征基底的第一线性组合值,以及所述第四频谱特征基底集中在所述基频和/或所述多个倍频上的频谱特征基底的第二线性组合值;
基于所述第一线性组合值与所述第二线性组合值,确定所述第一脉压映射关系。
4.根据权利要求1所述的中心动脉压波形拟合方法,其特征在于,所述利用所述动脉特征参数集对所述第二脉压映射关系进行校准,包括:
基于所述待测用户的历史第二脉搏波形与历史中心脉压波形,从所述动脉特征参数集中筛选出有效动脉特征参数;
利用所述有效动脉特征参数对所述第二脉压映射关系进行校准。
5.根据权利要求4所述的中心动脉压波形拟合方法,其特征在于,在所述第二脉搏波形为桡动脉波形的情况下,在从所述动脉特征参数集中筛选出有效动脉特征参数之前,还包括:
基于所述待测用户的历史桡动脉参数特征数据库,和/或所述待测用户的参考桡动脉参数特征,对所述动脉特征参数集进行优化,以剔除异常动脉特征参数;
所述历史桡动脉参数特征数据库是基于所述待测用户的多个历史桡动脉波形构建的,所述参考桡动脉参数特征是基于血压计所测量到的所述待测用户的任一历史桡动脉波形构建的。
6.根据权利要求1所述的中心动脉压波形拟合方法,其特征在于,还包括:
根据所述第二频谱特征基底集,得到第二脉搏波形。
7.根据权利要求6所述的中心动脉压波形拟合方法,其特征在于,所述第一检测法为光电容脉搏波描记法,所述第二脉搏波形用于表征采用示波法测得的脉搏波形。
8.一种中心动脉压的监控装置,其特征在于,包括:
脉搏波形采集模块,用于运行第一检测法获取待测用户的第一脉搏波形;
特征映射模块,用于基于所述第一脉搏波形对应的第一频谱特征基底集与第一脉压映射关系,获取第二频谱特征基底集;其中,所述第一脉压映射关系是根据获取的第一脉搏样本波形与获取的第二脉搏样本波形之间的脉压映射关系确定的;
波形拟合模块,用于基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形;其中,所述第二脉压映射关系是根据所述中心动脉压波形与第二脉搏波形之间的脉压映射关系确定的;
在基于所述第二频谱特征基底集与第二脉压映射关系,得到中心动脉压波形之前,还包括:
采集所述第二脉搏样本波形;
根据所述第二脉搏样本波形,获取动脉特征参数集;
利用所述动脉特征参数集对所述第二脉压映射关系进行校准,对动脉特征参数进行加减、比例或者线性组合得到第二脉压映射关系的校正系数;
预先通过第一检测法获取待测用户手腕处的脉搏波信号,从脉搏波信号中选取具有代表性的第一脉搏样本波形,其中,通过气囊在待测用户手腕处从零开始加压,获取待测用户的桡动脉波形,加压到桡动脉波形峰-峰值为最大脉搏波的1/3以下为止。
9.一种手表设备,包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1至7任一项所述中心动脉压波形拟合方法。
10.一种非暂态计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1至7任一项所述中心动脉压波形拟合方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211529875.1A CN115990007B (zh) | 2022-11-30 | 2022-11-30 | 中心动脉压波形拟合方法、监控装置和手表设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211529875.1A CN115990007B (zh) | 2022-11-30 | 2022-11-30 | 中心动脉压波形拟合方法、监控装置和手表设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115990007A CN115990007A (zh) | 2023-04-21 |
CN115990007B true CN115990007B (zh) | 2024-01-23 |
Family
ID=85991423
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211529875.1A Active CN115990007B (zh) | 2022-11-30 | 2022-11-30 | 中心动脉压波形拟合方法、监控装置和手表设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115990007B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1990011043A1 (en) * | 1989-03-24 | 1990-10-04 | Eastern Medical Testing Services, Inc. | A method and apparatus for ascertaining the contour of the pressure pulse in the central arteries from the contour of the pressure pulse in the peripheral arteries |
KR20080044224A (ko) * | 2008-04-30 | 2008-05-20 | (주)한별메디텍 | 비침습적 연속 혈압, 동맥탄성도 측정방법 |
JP2013230342A (ja) * | 2012-04-27 | 2013-11-14 | Microlife Corp | 中心動脈血圧の推定方法及びその装置 |
CN103908236A (zh) * | 2013-05-13 | 2014-07-09 | 天津点康科技有限公司 | 一种自动血压测量系统 |
CN104188642A (zh) * | 2014-09-24 | 2014-12-10 | 中国科学院合肥物质科学研究院 | 一种无创获取中心动脉压力波形及其相关参数的装置与获取方法 |
CN107785081A (zh) * | 2017-12-12 | 2018-03-09 | 北京动亮健康科技有限公司 | 计算中心血流动力学指标的方法、装置、存储介质及设备 |
CN109069031A (zh) * | 2016-10-20 | 2018-12-21 | 京东方科技集团股份有限公司 | 用于确定对象血压的设备和方法 |
CN109512412A (zh) * | 2018-11-21 | 2019-03-26 | 重庆理工大学 | 一种中心动脉血压测量装置 |
CN113499048A (zh) * | 2021-07-22 | 2021-10-15 | 重庆理工大学 | 一种基于CNN-BiLSTM的中心动脉压波形重构的系统及方法 |
CN113925474A (zh) * | 2020-07-10 | 2022-01-14 | 樱简(苏州)科技有限责任公司 | 一种基于pwtt原理的无袖带式血压监测系统 |
CN114947782A (zh) * | 2022-06-09 | 2022-08-30 | 重庆理工大学 | 一种基于ppg信号的中心动脉压波形重构的系统及方法 |
CN115281636A (zh) * | 2022-09-01 | 2022-11-04 | 广东乐心医疗电子股份有限公司 | 血压数据的处理方法、装置以及电子设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10849508B2 (en) * | 2016-06-03 | 2020-12-01 | Fourth Frontier Technologies Pvt. Ltd. | System and method for continuous monitoring of blood pressure |
WO2018167728A1 (en) * | 2017-03-17 | 2018-09-20 | Atcor Medical Pty Ltd | Central aortic blood pressure and waveform calibration method |
-
2022
- 2022-11-30 CN CN202211529875.1A patent/CN115990007B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1990011043A1 (en) * | 1989-03-24 | 1990-10-04 | Eastern Medical Testing Services, Inc. | A method and apparatus for ascertaining the contour of the pressure pulse in the central arteries from the contour of the pressure pulse in the peripheral arteries |
KR20080044224A (ko) * | 2008-04-30 | 2008-05-20 | (주)한별메디텍 | 비침습적 연속 혈압, 동맥탄성도 측정방법 |
JP2013230342A (ja) * | 2012-04-27 | 2013-11-14 | Microlife Corp | 中心動脈血圧の推定方法及びその装置 |
CN103908236A (zh) * | 2013-05-13 | 2014-07-09 | 天津点康科技有限公司 | 一种自动血压测量系统 |
CN104188642A (zh) * | 2014-09-24 | 2014-12-10 | 中国科学院合肥物质科学研究院 | 一种无创获取中心动脉压力波形及其相关参数的装置与获取方法 |
CN109069031A (zh) * | 2016-10-20 | 2018-12-21 | 京东方科技集团股份有限公司 | 用于确定对象血压的设备和方法 |
CN107785081A (zh) * | 2017-12-12 | 2018-03-09 | 北京动亮健康科技有限公司 | 计算中心血流动力学指标的方法、装置、存储介质及设备 |
CN109512412A (zh) * | 2018-11-21 | 2019-03-26 | 重庆理工大学 | 一种中心动脉血压测量装置 |
CN113925474A (zh) * | 2020-07-10 | 2022-01-14 | 樱简(苏州)科技有限责任公司 | 一种基于pwtt原理的无袖带式血压监测系统 |
CN113499048A (zh) * | 2021-07-22 | 2021-10-15 | 重庆理工大学 | 一种基于CNN-BiLSTM的中心动脉压波形重构的系统及方法 |
CN114947782A (zh) * | 2022-06-09 | 2022-08-30 | 重庆理工大学 | 一种基于ppg信号的中心动脉压波形重构的系统及方法 |
CN115281636A (zh) * | 2022-09-01 | 2022-11-04 | 广东乐心医疗电子股份有限公司 | 血压数据的处理方法、装置以及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN115990007A (zh) | 2023-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7479111B2 (en) | Methods for measuring blood pressure with automatic compensations | |
US20200268254A1 (en) | Continuous Blood Pressure Measurement | |
US20070149884A1 (en) | Pressure-based system and method for determining cardiac stroke volume | |
KR20170115550A (ko) | 대상자의 평균 동맥압을 도출하는 방법 및 장치 | |
Natarajan et al. | Photoplethysmography fast upstroke time intervals can be useful features for cuff-less measurement of blood pressure changes in humans | |
CN107785081B (zh) | 计算中心血流动力学指标的方法、装置、存储介质及设备 | |
US20200253564A1 (en) | Continuous Blood Pressure Measurement | |
CN109872820B (zh) | 一种无袖带血压测量方法、装置、设备及存储介质 | |
TW201400086A (zh) | 一種藉壓脈帶壓力脈波震盪訊號以預估中央動脈脈搏壓的方法及其裝置 | |
US8864678B2 (en) | Blood pressure measuring method and blood pressure manometer | |
Barbé et al. | Logistic ordinal regression for the calibration of oscillometric blood pressure monitors | |
Tanaka et al. | Accuracy assessment of a noninvasive device for monitoring beat-by-beat blood pressure in the radial artery using the volume-compensation method | |
CN115990007B (zh) | 中心动脉压波形拟合方法、监控装置和手表设备 | |
CN113470805A (zh) | 建立血压模型的方法 | |
Podaru et al. | Blood pressure estimation based on synchronous ECG and PPG recording | |
Sidhu et al. | Performance analysis of oscillometric blood pressure estimation techniques in cardiac patients | |
CN115734745A (zh) | 具有体积描记传感器的可穿戴设备 | |
US11154206B2 (en) | Viscoelasticity characteristics acquisition device, viscoelasticity characteristics acquisition method, viscoelasticity characteristics acquisition program, and recording medium recording said program | |
TWI228981B (en) | System and method for measurement and analysis of pulses from finger | |
Tang et al. | Continuous cuffless blood pressure estimation based on pulse transit time: a new evaluation with invasive BP reference | |
Hassan et al. | A portable continuous blood pressure monitoring kit | |
López-Lozada et al. | Blood Pressure Estimation Algorithm by a Cuff-Based Monitoring Unit | |
Wang et al. | An arterial compliance classification of the upper arm based on the cuff pressure signal | |
KEERTHANA et al. | NON-INVASIVE CUFFLESS BLOOD PRESSURE MEASUREMENT | |
CN113598724A (zh) | 用于评估心血管自主神经功能的装置 |
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 |