CN110200612B - 电子血压计测量方法、系统和电子血压计 - Google Patents
电子血压计测量方法、系统和电子血压计 Download PDFInfo
- Publication number
- CN110200612B CN110200612B CN201810167529.0A CN201810167529A CN110200612B CN 110200612 B CN110200612 B CN 110200612B CN 201810167529 A CN201810167529 A CN 201810167529A CN 110200612 B CN110200612 B CN 110200612B
- Authority
- CN
- China
- Prior art keywords
- pressure
- air pressure
- value
- clinical data
- column vector
- 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 abstract description 28
- 230000035488 systolic blood pressure Effects 0.000 claims abstract description 78
- 230000035487 diastolic blood pressure Effects 0.000 claims abstract description 75
- 230000010349 pulsation Effects 0.000 claims abstract description 73
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 17
- 238000012937 correction Methods 0.000 claims description 33
- 230000036772 blood pressure Effects 0.000 claims description 20
- 238000000605 extraction Methods 0.000 claims description 9
- 238000003825 pressing Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 10
- 238000005070 sampling Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 238000001914 filtration Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 5
- 210000004204 blood vessel Anatomy 0.000 description 4
- 238000009530 blood pressure measurement Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 2
- 229910052753 mercury Inorganic materials 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 210000000707 wrist Anatomy 0.000 description 1
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/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
- A61B5/02225—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers using the oscillometric method
-
- 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/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
- A61B5/0225—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers the pressure being controlled by electric signals, e.g. derived from Korotkoff sounds
-
- 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/7278—Artificial waveform generation or derivation, e.g. synthesising signals from measured signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
-
- 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
- A61B5/02116—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave amplitude
-
- 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)
- Cardiology (AREA)
- Vascular Medicine (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Ophthalmology & Optometry (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提供了一种电子血压计测量方法、系统和电子血压计,涉及电子血压计测量的技术领域,包括:获取气压值序列;根据气压值序列,得到收缩压/舒张压;提取气压值序列的特征信息;根据特征信息对收缩压/舒张压进行修正。本发明可以通过气压值序列计算收缩压/舒张压,并通过提取气压脉动序列的特征信息,根据特征信息对收缩压/舒张压进行修正,提高了收缩压/舒张压的计算精度。
Description
技术领域
本发明涉及电子血压计测量技术领域,尤其是涉及一种电子血压计测量方法、系统和电子血压计。
背景技术
使用示波法原理的电子血压计,已在市场上出现多年。该类型的电子血压计通过将空气袋卷绕在人体的某一部分上,通过控制气泵与气阀等元件,可以对空气袋进行充/放气,使得空气袋内空气量发生变化,进而采集气压值。其中,对空气袋的充/放气,有两种方案,一种是快速充气后,平缓放气,并且气压脉动信息是在平缓放气的过程中采集的;另一种是平缓充气,快速放气,并且气压脉动信息是在平缓充气的过程中采集的。
无论是采用充气式还是放气式测量,其包络线的横坐标(基础气压)按从左至右为由大到小顺序,两种测量方式生成的包络线是一致的。都是通过气压值计算包络线,进而得到受测者的收缩压/舒张压。然而,这种方式计算出来的收缩压/舒张压仅仅使用包络线这一个参数进行计算,使用的参数过于单一,从而导致了血压计算结果精度不足。
发明内容
有鉴于此,本发明的目的在于提供电子血压计测量方法、系统和电子血压计,可以通过气压值序列计算收缩压/舒张压,并通过提取气压脉动序列的特征信息,根据特征信息对收缩压/舒张压进行修正,提高了收缩压/舒张压的计算精度。
第一方面,本发明实施例提供了一种电子血压计测量方法,包括:获取气压值序列;根据所述气压值序列,得到收缩压/舒张压;提取所述气压值序列的特征信息;根据所述特征信息对所述收缩压/舒张压进行修正。
结合第一方面,本发明实施例提供了第一方面的第一种可能的实施方式,其中,所述提取所述气压值序列的特征信息,包括:从所述气压值序列中,分别提取气压脉动序列,及基础气压值序列;计算气压脉动序列的一阶导数,选取所述气压脉动序列的一阶导数中的极小值;获取与所述极小值对应的气压脉动的基础气压值,以将所述极小值和极小值对应的基础气压值作为特征信息对所述收缩压/舒张压进行修正。
结合第一方面,本发明实施例提供了第一方面的第二种可能的实施方式,其中,所述根据所述特征信息对所述收缩压/舒张压进行修正,包括:根据算式(1)修正所述收缩压,
Sys=A+B*Sys0+C*Dpmin+D*Pdpmin(1)
其中,Sys为修正后的收缩压,Sys0为修正前的收缩压,Dpmin为气压脉动的一阶导数的极小值,Pdpmin为极小值对应的基础气压值,A、B、C、D均为系数,A的单位为mmHg,C的单位为s,B、D无量纲;
根据算式(2)修正所述舒张压,
Dia=E+F*Dia0+G*Dpmin+H*Pdpmin(2)
其中,Dia为修正后的舒张压,Dia0为修正前的舒张压,E、F、G、H均为系数,E的单位为mmHg,G的单位为s,F、H无量纲。
结合第一方面,本发明实施例提供了第一方面的第三种可能的实施方式,其中,还包括:根据多个临床数据的收缩压参比值、多个临床数据,计算算式(1)中的系数A、B、C、D的值;根据多个临床数据的舒张压参比值、多个临床数据,计算算式(2)中的系数E、F、G、H的值。
结合第一方面,本发明实施例提供了第一方面的第四种可能的实施方式,其中,所述根据多个临床数据的收缩压参比值、多个临床数据,计算算式(1)中的系数A、B、C、D的值,包括:获取多个临床数据的收缩压参比值组成收缩压参比值列向量Y1;根据临床数据的个数组成列向量x1;根据所述多个临床数据,得到每个临床数据对应的收缩压,组成列向量x2;根据所述多个临床数据,得到每个临床数据对应的气压脉动序列的一阶导数的极小值,组成列向量x3;根据所述每个临床数据对应的气压脉动序列的一阶导数的极小值,得到所述极小值对应的基础气压值,组成列向量x4;根据列向量x1、列向量x2、列向量x3、列向量x4、列向量Y1,得到算式(1)中的系数A、B、C、D。
结合第一方面,本发明实施例提供了第一方面的第五种可能的实施方式,其中,所述根据多个临床数据的舒张压参比值、多个临床数据,计算算式(2)中的系数E、F、G、H的值,包括:获取多个临床数据的舒张压参比值组成舒张压参比值列向量Y2;根据临床数据的个数组成列向量x5;根据所述多个临床数据,得到每个临床数据对应的舒张压,组成列向量x6;根据所述多个临床数据,得到每个临床数据对应的气压脉动序列的一阶导数的极小值,组成列向量x7;根据所述每个临床数据对应的气压脉动序列的一阶导数的极小值,得到所述极小值对应的基础气压值,组成列向量x8;根据列向量x5、列向量x6、列向量x7、列向量x8、列向量Y2,得到算式(2)中的系数E、F、G、H。
第二方面,本发明实施例还提供一种电子血压计测量系统,包括:获取模块,用于获取气压值序列;计算模块,所述计算模块与所述获取模块相连,用于根据所述气压值序列,得到收缩压/舒张压;提取模块,所述提取模块与所述获取模块相连,用于提取所述气压值序列的特征信息;修正模块,所述修正模块分别与所述计算模块与所述提取模块相连,根据所述特征信息对所述收缩压/舒张压进行修正。
结合第二方面,本发明实施例提供了第二方面的第一种可能的实施方式,其中,所述提取模块,具体用于:从所述气压值序列中,分别提取气压脉动序列,及基础气压值序列;计算气压脉动序列的一阶导数,选取所述气压脉动序列的一阶导数中的极小值;获取与所述极小值对应的气压脉动的基础气压值,以将所述极小值和极小值对应的基础气压值作为特征信息对所述收缩压/舒张压进行修正。
结合第二方面,本发明实施例提供了第二方面的第二种可能的实施方式,其中,所述修正模块,包括:收缩压修正子模块,用于根据算式(1)修正所述收缩压,
Sys=A+B*Sys0+C*Dpmin+D*Pdpmin(1)
其中,Sys为修正后的收缩压,Sys0为修正前的收缩压,Dpmin为气压脉动的一阶导数的极小值,Pdpmin为极小值对应的基础气压值,A、B、C、D为系数,A的单位为mmHg,C的单位为s,B、D无量纲;
舒张压修正子模块,用于根据算式(2)修正所述舒张压,
Dia=E+F*Dia0+G*Dpmin+H*Pdpmin(2)
其中,Dia为修正后的舒张压,Dia0为修正前的舒张压,E、F、G、H均为系数,E的单位为mmHg,G的单位为s,F、H无量纲。
第二方面,本发明实施例还提供一种电子血压计测量系统,包括:采集部件和根据上述实施例任一项所述的电子血压计测量系统;所述采集部件用于采集气压值。
本发明实施例带来了以下有益效果:可以通过气压值序列计算收缩压/舒张压,并通过提取气压脉动序列的特征信息,根据特征信息对收缩压/舒张压进行修正,提高了收缩压/舒张压的计算精度。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一个实施例提供的电子血压计测量方法的流程图;
图2为空气袋内的气压值的示意图;
图3为本发明另一个实施例提供的电子血压计测量方法的流程图;
图4为本发明实施例提供的气压脉动的示意图;
图5为本发明实施例提供的包络线的示意图;
图6为本发明实施例提供的收缩压/舒张压在包络线中的位置示意图;
图7为本发明实施例提供的气压脉动一阶导数的示意图;
图8为本发明实施例提供的电子血压计测量系统的结构图;
图9为本发明实施例提供的电子血压计的结构图;
图10为本发明实施例提供的提取气压值序列的特征信息的流程图。
图标:
200-电子血压计测量系统;210-获取模块;220-计算模块;230-提取模块;240-修正模块;300-电子血压计;310-测量部件。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,采用电子血压计测量的气压值,无论是充气式还是放气式测量,都是计算包络线,根据包络线的峰值计算受测者的收缩压/舒张压。然而,这种方式计算出来的收缩压/舒张压仅仅使用包络线这一个参数进行计算,使用的参数过于单一,从而导致了血压计算结果精度不足,基于此,本发明实施例提供的一种电子血压计测量方法、系统和电子血压计,可以通过气压值序列计算收缩压/舒张压,并通过提取气压脉动序列的特征信息,根据特征信息对收缩压/舒张压进行修正,提高了收缩压/舒张压的计算精度。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种电子血压计测量方法进行详细介绍,参见图1所述,包括:
S110:获取气压值序列。
具体来说,电子血压计包括采集部件和数据处理部件,采集部件以空气袋为例,将空气袋卷绕在人体的上臂或者手腕处,通过气泵和气阀等元件,对空气袋进行充/放气,采集空气袋内空气量的气压值变化,将采集到的气压值传输到数据处理部件,从而可以获取多个气压值,多个气压值组成气压值序列。
结合图2所示,不断采集充气过程中的空气袋的气压值,在图1中,横坐标为采样点,64点对应1秒,纵坐标为气压,单位为mmHg。
S120:根据气压值序列,得到收缩压/舒张压。
作为一个示例,结合图3所示,首先从采集到的气压值序列中分别提取气压脉动序列和基础气压值序列。因为,在测量气压值时,有两个因素的影响能够使得气压值形成如图2所示的曲线图,第一个因素是空气袋内气压变化,由充放气导致的气压变化称为基础气压变化;第二个是人体被卷绕部分的血管受到空气袋的压迫时,血管的搏动也会导致空气袋内气压变化,由血管搏动导致的气压变化称为气压脉动,当基础气压变化时,血管受空气袋压迫的程度也发生变化,气压脉动幅度也会发生变化,即,图2中的气压值是在不断平稳变化的基础气压之上,会有振荡幅度不断变化的气压脉动叠加,所以采集到气压值序列在后期处理时,可以提取气压脉动序列和基础气压值序列。
步骤S121:生成包络线。具体来说,先提取出所有的气压脉动值,得到如图4所示的气压脉动图,在图4中,横坐标为采样点,64点对应1秒,纵坐标为气压值,单位为mmHg。其中,从气压值中提取出气压脉动值的方法为通过高通滤波的方式进行提取,还可以通过带通滤波、滑动平均滤波等等方式进行提取。以上述所有的气压脉动值的振幅为纵坐标,以它们对应的基础气压为横坐标,横坐标和纵坐标的单位均mmHg,为生成包络线,如图5所示,其中,可以对包络线进行平滑处理,实线为平滑处理前的包络线,虚线为平滑处理后的包络线。
作为一个示例,对包络线进行平滑出的方法包括多点加权平均,傅里叶变换与反变换等。
步骤S122:识别出包络线峰值,即:识别出包络线的最大振幅。将包络线极大值乘以舒张压乘以舒张压系数后,得到舒张压对应的气压脉动振幅;在包络线左侧寻找等于其值的点,此点对应的基础气压即为舒张压;将前述包络线极大值乘以收缩压系数后,得到收缩压对应的气压脉动振幅;在包络线右侧寻找等于其值的点,此点对应的基础气压即为收缩压。结合图6所示,最大振幅、舒张压对应点和收缩压对应点。在图6中,横坐标为气压脉动对应的基础气压,单位为mmHg,纵坐标为气压脉动振幅,单位为mmHg,黑色曲线为平滑过的包络线。
另外,舒张压系数与收缩压系数均为根据临床数据分析得到,舒张压系数取值范围在[0.2,0.7]之间,收缩压系数取值范围在[0.4,0.9]之间。
值得注意的是,上述计算的收缩压/舒张压只是示例性的,对于其他计算方法本发明对此不受限制。
S130:提取气压值序列的特征信息。
步骤S130具体包括:
结合图10所示,S131:将气压值序列拆分为气压脉动序列和气压脉动序列对应的基础气压值序列。
根据步骤S120中已将说明从气压值序列中分别提取气压脉动序列。具体请参考步骤S120中的内容。对于从气压值序列中提取基础气压值序列的方法为通过低通滤波的方式进行提取,还可以通过带通滤波、滑动平均滤波等等方式进行提取。
S132:计算气压脉动序列的一阶导数,选取气压脉动序列的一阶导数中的极小值。
具体来说,包括:计算当前气压脉动值与前一次测量的气压脉动值之差,除以采样周期,得到当前气压脉动的一阶导数;从气压脉动序列的一阶导数中选取极小值作为所述气压脉动序列的一阶导数的极小值。
其中,气压脉动的一阶导数也可以折算为单位之间内气压脉动的变化的气压值。结合图7所示,为气压脉动的一阶导数示意图,横坐标为采样点,采样频率64Hz,即每1秒采样64次,纵坐标为气压脉动的一阶导数,单位为mmHg/s。
举例,如果有n个气压值组成一个气压序列,则可从中分别提取出n个气压脉动值和n个基础气压值,将气压脉动序列的相邻点两两相减并除以采样周期,可以计算出n-1个气压脉动的一阶导数,然后在其前面补充一个零值,可得到长度为n的一阶导数序列,在这个n个气压脉动的一阶导数中,通过比较选取一个数值最小的值为极小值。
S133:获取与极小值对应的气压脉动的基础气压值。
举例,长度为n的一阶导数值序列,与长度为n的基础气压值序列一一对应,如果第K个一阶导数值为极小值,则第K个一阶导数值对应的基础气压值就是想要得到的基础气压值。
S140:根据特征信息对收缩压/舒张压进行修正。
具体包括:根据算式(1)修正收缩压,
Sys=A+B*Sys0+C*Dpmin+D*Pdpmin(1)
其中,Sys为修正后的收缩压,Sys0为修正前的收缩压,Dpmin为气压脉动的一阶导数的极小值,Pdpmin为极小值对应的基础气压值,A、B、C、D为系数,A的单位为mmHg,C的单位为s,B、D无量纲;
根据算式(2)修正所述舒张压,
Dia=E+F*Dia0+G*Dpmin+H*Pdpmin(2)
其中,Dia为修正后的舒张压,Dia0为修正前的舒张压,E、F、G、H为系数,E的单位为mmHg,G的单位为s,F、G无量纲。
进一步的,还包括:根据多个临床数据的收缩压参比值、多个临床数据,计算算式(1)中的系数A、B、C、D的值;根据多个临床数据的舒张压参比值、多个临床数据,计算算式(2)中的系数E、F、G、H的值。
根据多个临床数据的收缩压参比值、多个临床数据,计算算式(1)中的系数A、B、C、D的值,具体包括:获取多个临床数据的收缩压参比值组成收缩压参比值列向量Y1;根据临床数据的个数组成列向量x1;根据多个临床数据,得到每个临床数据对应的收缩压,组成列向量x2;根据多个临床数据,得到每个临床数据对应的气压脉动序列的一阶导数的极小值,组成列向量x3;根据每个临床数据对应的气压脉动序列的一阶导数的极小值,得到极小值对应的基础气压值,组成列向量x4;根据列向量x1、列向量x2、列向量x3、列向量x4、列向量Y1,得到算式(1)中的系数A、B、C、D。
具体来说,算式(1)中的系数A、B、C、D的取值均是根据对有效临床测试数据的分析统计而得到的,有效临床数据要求临床数据量越大越好,且对受试人群的分布有要求,受试者类型/血压等级分布越广越好。如表1所示:
表1
表1为本发明所使用的临床数据数量与分布,临床数据的例数、分布范围、测试方法参考ISO81060:2013的要求。
其中,临床数据的收缩压参比值的意思是收缩压的参考比较数值,是已经通过临床数据得到的收缩压值。
在一些实施例中,根据列向量x1、列向量x2、列向量x3、列向量x4、列向量Y1,得到算式(1)中的系数A、B、C、D,包括:对Y=X*P等式中的P进行拟合,优选的,使用软件matlab中的regress函数;P是一个长度为4的列向量,其元素依次分别对应前述收缩压精确计算公式中的系数A、B、C、D。上述的拟合方法还可以使用其他方式,例如使用软件matlab中的ployfit函数、cftool工具箱等,以及其他软件如EXCEL的趋势预测/回归分析功能等。
根据多个临床数据的舒张压参比值、多个临床数据,计算算式(2)中的系数E、F、G、H的值,具体包括:获取多个临床数据的舒张压参比值组成舒张压参比值列向量Y2;根据临床数据的个数组成列向量x5;根据所述多个临床数据,得到每个临床数据对应的舒张压,组成列向量x6;根据所述多个临床数据,得到每个临床数据对应的气压脉动序列的一阶导数的极小值,组成列向量x7;根据每个临床数据对应的气压脉动序列的一阶导数的极小值,得到极小值对应的基础气压值,组成列向量x8;根据列向量x5、列向量x6、列向量x7、列向量x8、列向量Y2,得到算式(2)中的系数E、F、G、H。
本发明实施例所提供的根据多个临床数据的舒张压参比值、多个临床数据,计算算式(2)中的系数E、F、G、H的值,其和前述根据多个临床数据的收缩压参比值、多个临床数据,计算算式(1)中的系数A、B、C、D的值的实施例的实现原理相同,为简要描述,可参考前述部分中相应内容。
根据上述的方式得到的A、B、C、D、E、F、G、H的取值范围为:A的取值区间在[0,25],B的取值区间在[0,2],C的取值区间在[-0.1,0.1],D的取值区间在[-0.5,0.5],E的取值区间在[0,25],F的取值区间在[0,2],G的取值区间在[-0.1,0.1],H的取值区间在[-0.5,0.5]。
作为一个示例,随机选取10位成人受测者进行对比测量,每位受试者使用电子血压计(公知技术)、电子血压计(本发明方法)、水银血压计分别依次测试三次,并分别求均值并记录,如表2。测试精度的对比结果如表3。
表2本发明方法与公知技术的血压测试数据列表
如表2所示,受试者为10名,采用电子血压计为受试者在公知技术下计算出来收缩压/舒张压,使用本发明的方法计算出来的收缩压/舒张压,和采用水银血压计作为参考比较值。
表3本发明方法与公知技术的血压精度对比
表3中记载了采用公知技术计算的收缩压/舒张压的平均差和标准差,采用本发明的方式计算的收缩压/舒张压的平均差和标准差,可以看出采用本发明方法后,收缩压/舒张压的平均差与标准差均有所下降,精度提高。
参见图8所示,电子血压计测量系统200,包括:获取模块210、计算模块220、提取模块230、修正模块240。
其中,获取模块210,用于获取气压值序列。计算模块220与获取模块210相连,用于根据气压值序列,得到收缩压/舒张压。提取模块230与获取模块210相连,用于提取气压值序列的特征信息。修正模块240分别与计算模块220和提取模块230相连,根据特征信息对收缩压/舒张压进行修正。
在一些实施例中,提取模块230,具体用于:从气压值序列中,分别提取气压脉动序列,及基础气压值序列;计算气压脉动序列的一阶导数,选取气压脉动序列的一阶导数中的极小值;获取与极小值对应的气压脉动的基础气压值,以将极小值和极小值对应的基础气压值作为特征信息对收缩压/舒张压进行修正。
在一些实施例中,修正模块240,包括:收缩压修正子模块,用于根据算式(1)修正所述收缩压,
Sys=A+B*Sys0+C*Dpmin+D*Pdpmin(1)
其中,Sys为修正后的收缩压,Sys0为修正前的收缩压,Dpmin为气压脉动的一阶导数的极小值,Pdpmin为极小值对应的基础气压值,A、B、C、D为系数,A的单位为mmHg,C的单位为s,B、D无量纲;
舒张压修正子模块,用于根据算式(2)修正所述舒张压,
Dia=E+F*Dia0+G*Dpmin+H*Pdpmin(2)
其中,Dia为修正后的舒张压,Dia0为修正前的舒张压,E、F、G、H为系数,E的单位为mmHg,G的单位为s,F、H无量纲。
本发明实施例所提供的系统,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例中相应内容。
参见图9所示,电子血压计300,包括:测量部件310和根据上述实施例任一项所述的电子血压计测量系统200。其中,测量部件310用于向人体测量部位施加压力并采集气压值。
本发明实施例提供的电子血压计测量,与上述实施例提供的电子血压计测量系统具有相同的技术特征,所以也能解决相同的技术问题,达到相同的技术效果。
除非另外具体说明,否则在这些实施例中阐述的部件和步骤的相对步骤、数字表达式和数值并不限制本发明的范围。
在这里示出和描述的所有示例中,任何具体值应被解释为仅仅是示例性的,而不是作为限制,因此,示例性实施例的其他示例可以具有不同的值。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本发明实施例的描述中,除非另有明确的规定和限定,术语“相连”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中,需要说明的是,术语“上”、“下”、“左”、“右”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (6)
1.一种电子血压计测量系统,其特征在于,包括:
获取模块,用于获取气压值序列;
计算模块,所述计算模块与所述获取模块相连,用于根据所述气压值序列,得到收缩压/舒张压;
提取模块,所述提取模块与所述获取模块相连,用于提取所述气压值序列的特征信息;
修正模块,所述修正模块分别与所述计算模块与所述提取模块相连,根据所述特征信息对所述收缩压/舒张压进行修正;
所述提取模块,具体用于:从所述气压值序列中,分别提取气压脉动序列,及基础气压值序列;计算气压脉动序列的一阶导数,选取所述气压脉动序列的一阶导数中的极小值;获取与所述极小值对应的气压脉动的基础气压值,以将所述极小值和极小值对应的基础气压值作为特征信息对所述收缩压/舒张压进行修正。
2.根据权利要求1所述的电子血压计测量系统,其特征在于,所述修正模块,包括:
收缩压修正子模块,用于根据算式(1)修正所述收缩压,
Sys=A+B*Sys0+C*Dpmin+D*Pdpmin (1)
其中,Sys为修正后的收缩压,Sys0为修正前的收缩压,Dpmin为气压脉动序列的一阶导数的极小值,Pdpmin为极小值对应的基础气压值,A、B、C、D均为系数,A的单位为mmHg,C的单位为s;
舒张压修正子模块,用于根据算式(2)修正所述舒张压,
Dia=E+F*Dia0+G*Dpmin+H*Pdpmin (2)
其中,Dia为修正后的舒张压,Dia0为修正前的舒张压,E、F、G、H均为系数,E的单位为mmHg,G的单位为s。
3.根据权利要求2所述的电子血压计测量系统,其特征在于,所述修正模块还包括:
第一计算单元,用于根据多个临床数据的收缩压参比值、多个临床数据,计算算式(1)中的系数A、B、C、D的值;
第二计算单元,用于根据多个临床数据的舒张压参比值、多个临床数据,计算算式(2)中的系数E、F、G、H的值。
4.根据权利要求3所述的电子血压计测量系统,其特征在于,所述第一计算单元还用于:
获取多个临床数据的收缩压参比值组成收缩压参比值列向量Y1;
根据临床数据的个数组成列向量x1;
根据所述多个临床数据,得到每个临床数据对应的收缩压,组成列向量x2;
根据所述多个临床数据,得到每个临床数据对应的气压脉动序列的一阶导数的极小值,组成列向量x3;
根据所述每个临床数据对应的气压脉动序列的一阶导数的极小值,得到所述极小值对应的基础气压值,组成列向量x4;
根据列向量x1、列向量x2、列向量x3、列向量x4、列向量Y1,得到算式(1)中的系数A、B、C、D。
5.根据权利要求3所述的电子血压计测量系统,其特征在于,所述第二计算单元还用于:
获取多个临床数据的舒张压参比值组成舒张压参比值列向量Y2;
根据临床数据的个数组成列向量x5;
根据所述多个临床数据,得到每个临床数据对应的舒张压,组成列向量x6;
根据所述多个临床数据,得到每个临床数据对应的气压脉动序列的一阶导数的极小值,组成列向量x7;
根据所述每个临床数据对应的气压脉动序列的一阶导数的极小值,得到所述极小值对应的基础气压值,组成列向量x8;
根据列向量x5、列向量x6、列向量x7、列向量x8、列向量Y2,得到算式(2)中的系数E、F、G、H。
6.一种电子血压计,其特征在于,包括:测量部件和权利要求1-5任一项所述的电子血压计测量系统;
所述测量部件用于向人体受测部位施加压力并采集气压值。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810167529.0A CN110200612B (zh) | 2018-02-28 | 2018-02-28 | 电子血压计测量方法、系统和电子血压计 |
PCT/CN2018/088786 WO2019165717A1 (zh) | 2018-02-28 | 2018-05-29 | 电子血压计测量方法、系统、电子血压计及存储介质 |
EP18907855.3A EP3760110B1 (en) | 2018-02-28 | 2018-05-29 | Electronic sphygmomanometer measuring method, system, electronic sphygmomanometer, and storage medium |
US16/976,476 US11771331B2 (en) | 2018-02-28 | 2018-05-29 | Electronic sphygmomanometer measuring method, system, electronic sphygmomanometer, and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810167529.0A CN110200612B (zh) | 2018-02-28 | 2018-02-28 | 电子血压计测量方法、系统和电子血压计 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110200612A CN110200612A (zh) | 2019-09-06 |
CN110200612B true CN110200612B (zh) | 2021-02-26 |
Family
ID=67778696
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810167529.0A Active CN110200612B (zh) | 2018-02-28 | 2018-02-28 | 电子血压计测量方法、系统和电子血压计 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11771331B2 (zh) |
EP (1) | EP3760110B1 (zh) |
CN (1) | CN110200612B (zh) |
WO (1) | WO2019165717A1 (zh) |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4984577A (en) * | 1989-03-20 | 1991-01-15 | Hewlett-Packard Company | Oscillometric non-invasive method for measuring blood pressure and apparatus for automated oscillometric blood pressure measuring |
JPH05309073A (ja) * | 1992-05-13 | 1993-11-22 | Omron Corp | 電子血圧計 |
CN1394546A (zh) * | 2002-08-08 | 2003-02-05 | 天津市先石光学技术有限公司 | 一种血压测定装置及方法 |
CN1698535A (zh) * | 2004-05-18 | 2005-11-23 | 香港中文大学 | 测量血压变化率的方法 |
CN100577092C (zh) * | 2005-05-24 | 2010-01-06 | 北京新兴阳升科技有限公司 | 获取和处理与动脉血压相关信息柯氏音延迟时间的方法 |
JP4470876B2 (ja) * | 2005-12-20 | 2010-06-02 | オムロンヘルスケア株式会社 | 電子血圧計 |
JP5098721B2 (ja) * | 2008-03-14 | 2012-12-12 | オムロンヘルスケア株式会社 | 血圧測定装置、血圧導出プログラムおよび血圧導出方法 |
TWI383778B (zh) | 2009-09-03 | 2013-02-01 | Quanta Comp Inc | 具有消除動作雜訊功能之血壓量測裝置及方法 |
CN102202567A (zh) * | 2011-04-29 | 2011-09-28 | 华为技术有限公司 | 一种电子血压计的误差检测方法、装置及电子血压计 |
CN102283638A (zh) * | 2011-09-27 | 2011-12-21 | 江苏鱼跃信息系统有限公司 | 智能系数匹配的血压测量方法 |
JP6100705B2 (ja) * | 2014-01-22 | 2017-03-22 | 日本光電工業株式会社 | 血圧測定システム |
JP6615180B2 (ja) * | 2014-08-28 | 2019-12-04 | コーニンクレッカ フィリップス エヌ ヴェ | 振動非侵襲血圧(nibp)測定に対する方法及びnibp装置に対する制御ユニット |
CN104367308B (zh) * | 2014-11-10 | 2016-08-24 | 中国计量学院 | 一种可外部压力校准的电子血压计 |
KR102396248B1 (ko) * | 2015-06-12 | 2022-05-09 | 삼성전자주식회사 | 혈압 측정 장치 및 그의 동작 방법 |
CN105615845B (zh) | 2016-02-25 | 2020-05-19 | 广州视源电子科技股份有限公司 | 干扰脉搏信号检测方法及其系统 |
CN106510669B (zh) | 2016-11-14 | 2019-08-16 | 中国科学院电子学研究所 | 无袖带血压测量系统 |
CN107174230B (zh) * | 2017-06-22 | 2023-10-03 | 深圳市普瑞拓科技有限公司 | 智能血压计及其血压计量方法 |
CN107320088B (zh) * | 2017-06-23 | 2021-03-26 | 成都市欣康兴泰科技有限公司 | 无创连续脉血压测量方法 |
CN107320089A (zh) * | 2017-06-27 | 2017-11-07 | 西南大学 | 自校准的人体血压测量方法 |
-
2018
- 2018-02-28 CN CN201810167529.0A patent/CN110200612B/zh active Active
- 2018-05-29 EP EP18907855.3A patent/EP3760110B1/en active Active
- 2018-05-29 WO PCT/CN2018/088786 patent/WO2019165717A1/zh unknown
- 2018-05-29 US US16/976,476 patent/US11771331B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US11771331B2 (en) | 2023-10-03 |
WO2019165717A1 (zh) | 2019-09-06 |
EP3760110A4 (en) | 2021-05-12 |
CN110200612A (zh) | 2019-09-06 |
US20210000354A1 (en) | 2021-01-07 |
EP3760110A1 (en) | 2021-01-06 |
EP3760110B1 (en) | 2024-05-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108903929B (zh) | 心率检测修正的方法、装置、存储介质和系统 | |
KR101947676B1 (ko) | 생체 신호를 측정하는 장치 및 방법 | |
US20120157791A1 (en) | Adaptive time domain filtering for improved blood pressure estimation | |
EP3773157B1 (en) | Apparatus for use with a wearable cuff | |
US20080234589A1 (en) | Blood pressure algorithm | |
CN102579023B (zh) | 脉搏波信号处理方法和装置及电子血压测量装置 | |
CN110200612B (zh) | 电子血压计测量方法、系统和电子血压计 | |
CN112006669B (zh) | 一种基于血压测量模型的双通道计算方法的血压仪 | |
CN113413135A (zh) | 一种基于脉搏采集分析的方法、系统、装置及存储介质 | |
CN110801214A (zh) | 一种心率实时检测方法及系统 | |
KR101992763B1 (ko) | 비정상 심전도 신호 정보 출력 디바이스 및 방법 | |
CN114305359B (zh) | 血压数据采集设备、芯片 | |
JP3781956B2 (ja) | 脈波検出装置 | |
CN104902809B (zh) | 用于减少ecg信号中的运动伪影的方法和装置 | |
CN114748052A (zh) | 呼吸电阻抗成像激励参数选择方法、装置、设备与介质 | |
KR100650041B1 (ko) | 가속도 센서를 이용한 동잡음 제거수단을 구비하는오실로메트릭 방식 혈압계 | |
CN110876614B (zh) | 电子血压计 | |
CN109077711B (zh) | 动态心率数据获取方法、装置、穿戴设备及可读存储介质 | |
WO2020090763A1 (ja) | 処理装置、システム、処理方法、およびプログラム | |
KR20170037704A (ko) | 흉부임피던스를 측정을 위한 다중 전극 및 이를 이용한 흉부 임피던스 측정방법 | |
KR102041456B1 (ko) | 미네소타코드 출력 디바이스 및 방법 | |
EP3744238A1 (en) | Blood pressure measurement device and calculation method thereof | |
CN115721275B (zh) | 脉搏波数据采集和血压测量方法、装置、设备及介质 | |
US20210361951A1 (en) | Pulse discrimination device and electrocardiogram analyzer | |
CN106344022B (zh) | 一种呼吸率提取方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
GR01 | Patent grant | ||
GR01 | Patent grant |