CN110403582B - 一种用于分析脉波波形品质的方法 - Google Patents
一种用于分析脉波波形品质的方法 Download PDFInfo
- Publication number
- CN110403582B CN110403582B CN201910667483.3A CN201910667483A CN110403582B CN 110403582 B CN110403582 B CN 110403582B CN 201910667483 A CN201910667483 A CN 201910667483A CN 110403582 B CN110403582 B CN 110403582B
- Authority
- CN
- China
- Prior art keywords
- cluster
- pulse waveform
- samples
- clusters
- pulse
- 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
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
Abstract
本发明涉及医疗设备技术领域,尤其是指一种用于分析脉波波形品质的方法,包括依次进行的以下步骤:A:根据相似度对各个脉波波形进行群集划分;B:比较由步骤A获得的各个群集中的样本数量,选出样本数量最多的群集;C:若样本数量最多的群集数量为一个,则判定脉波波形品质良好;若样本数量最多的群集数量大于等于两个且小于四个,则判定脉波波形品质中等;若样本数量最多的群集数量大于等于四个,则判定脉波波形品质不良。本发明通过将获取的脉波波形进行群集划分,并根据样本数量最多的群集数量来判定脉波波形的品质,从而使得脉波波形的品质能够量化。
Description
技术领域
本发明涉及医疗设备技术领域,尤其是指一种用于分析脉波波形品质的方法。
背景技术
许多穿戴式装置都会通过动脉脉波的波形来侦测心跳、血压、心搏量、动脉弹性等心脏、血管相关参数。若脉波波形的品质不佳,即脉波波形的形状不正确,则计算出来的相关参数也会不正确,导致误导使用者。因此,有必要对脉波波形的品质进行分析,以便使用者知悉脉波波形品质,并以此参照各项参数的可靠性。
发明内容
本发明针对现有技术的问题提供一种用于分析脉波波形品质的方法。
本发明采用如下技术方案:一种用于分析脉波波形品质的方法,包括依次进行的以下步骤:A:根据相似度对各个脉波波形进行群集划分;B:比较由步骤A获得的各个群集中的样本数量,选出样本数量最多的群集;C:若样本数量最多的群集数量为一个,则判定脉波波形品质良好;若样本数量最多的群集数量大于等于两个且小于四个,则判定脉波波形品质中等;若样本数量最多的群集数量大于等于四个,则判定脉波波形品质不良。
作为优选,在步骤A中,进行群集划分的步骤包括:A1:获取脉波波形;A2:在A1步骤后,查看是否有已经建立的群集,若有,则将A1步骤中获取的脉波波形与最新建立的群集的群中心进行比较,并进入步骤A21;若无,则以A1步骤中获取的脉波波形建立新的群集并以该脉波波形作为所述群集的群中心,随后进入步骤A3; A21:若A1步骤中获取的脉波波形与最新建立的群集的群中心相似度大于等于75%,则将该脉波波形纳入所述群集并作为所述群集的群中心,随后进入步骤A3;若A1步骤中获取的脉波波形与最新建立的群集的群中心相似度小于75%,则以该脉波波形建立新的群集并以该脉波波形作为所述群集的群中心,随后进入步骤A3;A3:查看是否还有需要获取的脉波波形,若有,则进入A1;若无,则结束群集划分。
作为优选,在步骤A1后以及步骤A2前,比较步骤A1中获取的脉波波形与指数衰减函数的相关度,若相关度大于等于90%则进入步骤A2;若相关度小于90%则进入步骤A3。
作为优选,在步骤C中,若样本数量最多的群集数量为一个,则以样本数量最多的群集中所有脉波波形的平均脉波波形作为选定结果。
作为优选,在步骤C中,若样本数量最多的群集数量大于等于两个,则进行以下步骤:C1:若样本数量最多的群集中,每个群集的样本数量均为大于两个,将群集中所有脉波波形平均相似度最大的群集作为选定群集,并以选定群集的平均脉波波形作为选定结果;C2:若样本数量最多的群集中,每个群集的样本数量均小于等于两个,则判定群集划分失败并进行步骤A。
本发明的有益效果:通过将获取的脉波波形进行群集划分,并根据样本数量最多的群集数量来判定脉波波形的品质,从而使得脉波波形的品质能够量化。
具体实施方式
为了便于本领域技术人员的理解,下面结合实施例对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
一种用于分析脉波波形品质的方法,包括依次进行的以下步骤: A:根据相似度对各个脉波波形进行群集划分;B:比较由步骤A获得的各个群集中的样本数量,选出样本数量最多的群集;C:若样本数量最多的群集数量为一个,则判定脉波波形品质良好;若样本数量最多的群集数量大于等于两个且小于四个,则判定脉波波形品质中等;若样本数量最多的群集数量大于等于四个,则判定脉波波形品质不良。
脉波波形变形相对脉波波形正常来说是更小概率的事件,同时脉波波形变形后,与正常脉波波形的相似度会降低,因此在正常情况下,正常脉波波形数量较多自成一个群集,而且其他变形了脉波波形会形成其他的群集,此时穿戴式的医疗设备只要用正常脉波波形形成的群集来进行计算,就能获得较为可靠的参数,因此此时可以判定脉波波形品质良好。当各个变形的脉波波形之间的相似度较高时,变形的脉波波形也会形成样本数量较多的群集,此时穿戴式的医疗设备难以判定哪个群集是拥有正常脉波波形的群集,只能利用多个群集中的脉波波形来计算参数,当群集数量越多时,计算出来的参数的误差就越大。通过实验证明,当样本数量最多的群集数量大于等于两个且小于四个时,参数的误差仍在可接受的范围内,因此此时可认为脉波波形的品质为中等,但当样本数量最多的群集数量大于等于四个时,计算出来的参数误差较高,可靠性较弱,因此此时可认为脉波波形的品质为不良。
进一步地,在步骤A中,进行群集划分的步骤包括:A1:获取脉波波形;A2:在A1步骤后,查看是否有已经建立的群集,若有,则将A1步骤中获取的脉波波形与最新建立的群集的群中心进行比较,并进入步骤A21;若无,则以A1步骤中获取的脉波波形建立新的群集并以该脉波波形作为所述群集的群中心,随后进入步骤A3; A21:若A1步骤中获取的脉波波形与最新建立的群集的群中心相似度大于等于75%,则将该脉波波形纳入所述群集并作为所述群集的群中心,随后进入步骤A3;若A1步骤中获取的脉波波形与最新建立的群集的群中心相似度小于75%,则以该脉波波形建立新的群集并以该脉波波形作为所述群集的群中心,随后进入步骤A3;A3:查看是否还有需要获取的脉波波形,若有,则进入A1;若无,则结束群集划分。本发明中,将连续获得的且相似度较高的脉波波形划分入同一个群集。当一个群集中的样本数量较多时,可以表明穿戴式医疗设备在采集该群集中的脉波波形时,长期处于较为稳定的采集状态,亦即该群集中的脉波波形可靠性较强,利用该群集进行计算的参数更为可信。
进一步地,在步骤A1后以及步骤A2前,比较步骤A1中获取的脉波波形与指数衰减函数的相关度,若相关度大于等于90%则进入步骤A2;若相关度小于90%则进入步骤A3。在一个舒张期内,正常的脉波波形应与指数衰减函数有较强的相关度。相关度弱的脉波波形可直接定义为严重变形的脉波波形,因此本发明中,直接把相关度弱的脉波波形丢弃,即不把其纳入所有群集。
进一步地,在步骤C中,若样本数量最多的群集数量为一个,则以样本数量最多的群集中所有脉波波形的平均脉波波形作为选定结果。样本数量最多的群集只有一个时,可直接认为该群集为正常波形所在的群集,利用该群集的平均脉波波形作为选定结果以计算各项参数,即可获得可靠性强的参数。
进一步地,在步骤C中,若样本数量最多的群集数量大于等于两个,则进行以下步骤:C1:若样本数量最多的群集中,每个群集的样本数量均为大于两个,将群集中所有脉波波形平均相似度最大的群集作为选定群集,并以选定群集的平均脉波波形作为选定结果;C2:若样本数量最多的群集中,每个群集的样本数量均小于等于两个,则判定群集划分失败并进行步骤A。群集的样本数量均为大于两个时,可以表明穿戴式医疗设备在采集该群集中的脉波波形时,长期处于较为稳定的采集状态,亦即该群集中的脉波波形可靠性较强,利用该群集进行计算的参数更为可信。群集的样本数量均小于等于两个时,可以表明穿戴式医疗设备在采集该群集中的脉波波形时,没有处于较为稳定的采集状态,亦即该群集中的脉波波形可靠性较弱,利用该群集进行计算的参数不可信,因此此时应重新采集脉波波形并重新划分群集。
以上所述,仅是本发明较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明以较佳实施例公开如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当利用上述揭示的技术内容作出些许变更或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案内容,依据本发明技术是指对以上实施例所作的任何简单修改、等同变化与修饰,均属于本发明技术方案的范围内。
Claims (5)
1.一种用于分析脉波波形品质的方法,其特征在于:包括依次进行的以下步骤:
A:根据相似度对各个脉波波形进行群集划分;
B:比较由步骤A获得的各个群集中的样本数量,选出样本数量最多的群集;
C:若样本数量最多的群集数量为一个,则判定脉波波形品质良好;若样本数量最多的群集数量大于等于两个且小于四个,则判定脉波波形品质中等;若样本数量最多的群集数量大于等于四个,则判定脉波波形品质不良。
2.根据权利要求1所述的一种用于分析脉波波形品质的方法,其特征在于:在步骤A中,进行群集划分的步骤包括:
A1:获取脉波波形;
A2:在A1步骤后,查看是否有已经建立的群集,若有,则将A1步骤中获取的脉波波形与最新建立的群集的群中心进行比较,并进入步骤A21;若无,则以A1步骤中获取的脉波波形建立新的群集并以该脉波波形作为所述群集的群中心,随后进入步骤A3;
A21:若A1步骤中获取的脉波波形与最新建立的群集的群中心相似度大于等于75%,则将该脉波波形纳入所述群集并作为所述群集的群中心,随后进入步骤A3;若A1步骤中获取的脉波波形与最新建立的群集的群中心相似度小于75%,则以该脉波波形建立新的群集并以该脉波波形作为所述群集的群中心,随后进入步骤A3;
A3:查看是否还有需要获取的脉波波形,若有,则进入A1;若无,则结束群集划分。
3.根据权利要求2所述的一种用于分析脉波波形品质的方法,其特征在于:在步骤A1后以及步骤A2前,比较步骤A1中获取的脉波波形与指数衰减函数的相关度,若相关度大于等于90%则进入步骤A2;若相关度小于90%则进入步骤A3。
4.根据权利要求1所述的一种用于分析脉波波形品质的方法,其特征在于:在步骤C中,若样本数量最多的群集数量为一个,则以样本数量最多的群集中所有脉波波形的平均脉波波形作为选定结果。
5.根据权利要求1所述的一种用于分析脉波波形品质的方法,其特征在于:在步骤C中,若样本数量最多的群集数量大于等于两个,则进行以下步骤:
C1:若样本数量最多的群集中,每个群集的样本数量均为大于两个,将群集中所有脉波波形平均相似度最大的群集作为选定群集,并以选定群集的平均脉波波形作为选定结果;
C2:若样本数量最多的群集中,每个群集的样本数量均小于等于两个,则判定群集划分失败并进行步骤A。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910667483.3A CN110403582B (zh) | 2019-07-23 | 2019-07-23 | 一种用于分析脉波波形品质的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910667483.3A CN110403582B (zh) | 2019-07-23 | 2019-07-23 | 一种用于分析脉波波形品质的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110403582A CN110403582A (zh) | 2019-11-05 |
CN110403582B true CN110403582B (zh) | 2021-12-03 |
Family
ID=68362580
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910667483.3A Active CN110403582B (zh) | 2019-07-23 | 2019-07-23 | 一种用于分析脉波波形品质的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110403582B (zh) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1662921A (zh) * | 2002-06-24 | 2005-08-31 | 皇家飞利浦电子股份有限公司 | 比较各种初始群组集合以便确定用于群集一个电视节目集合的最佳初始集合的方法 |
JPWO2004097703A1 (ja) * | 2003-05-02 | 2006-07-13 | 富士通株式会社 | 製造工程解析支援方法、その方法をコンピュータに実行させるプログラム、プログラムプロダクトおよび記録媒体 |
CN101120615A (zh) * | 2005-02-22 | 2008-02-06 | 弗劳恩霍夫应用研究促进协会 | 近透明或透明的多声道编码器/解码器方案 |
CN102112984A (zh) * | 2008-07-29 | 2011-06-29 | 皇家飞利浦电子股份有限公司 | 用于生成图像选集的方法和设备 |
CN102460433A (zh) * | 2009-06-15 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | 选择代表图像的方法和设备 |
CN102547844A (zh) * | 2010-12-29 | 2012-07-04 | 通用电气公司 | 无线网络中的动态数据管理的系统和方法 |
CN104090861A (zh) * | 2013-02-20 | 2014-10-08 | 哈佛蒸汽锅炉检验和保险公司 | 动态离群值偏倚减少系统和方法 |
CN104586383A (zh) * | 2015-01-12 | 2015-05-06 | 北京世纪今科医疗器械有限公司 | 一种心电波形分类方法及装置 |
WO2015107268A1 (en) * | 2014-01-16 | 2015-07-23 | Aboa Legis Oy | Method and device for the detection of respiratory rate |
CN105411567A (zh) * | 2015-12-15 | 2016-03-23 | 北京世纪今科医疗器械有限公司 | 一种心电图数据分析方法及系统 |
CN105686818A (zh) * | 2016-03-30 | 2016-06-22 | 华为技术有限公司 | 测量装置及方法 |
CN106254175A (zh) * | 2016-07-26 | 2016-12-21 | 北京蓝海讯通科技股份有限公司 | 一种集群异常节点检测方法、应用和计算设备 |
CN106963353A (zh) * | 2015-12-30 | 2017-07-21 | 朝阳科技大学 | 评估动脉血流介导血管舒张功能的方法 |
CN107423757A (zh) * | 2017-07-14 | 2017-12-01 | 北京小米移动软件有限公司 | 聚类处理方法及装置 |
CN108897821A (zh) * | 2018-06-21 | 2018-11-27 | 江苏赛睿信息科技股份有限公司 | 自动生成数据结论的方法及装置 |
CN109344146A (zh) * | 2018-09-13 | 2019-02-15 | 海南电网有限责任公司电力科学研究院 | 一种大规模过电压数据智能自清洗方法 |
CN109490776A (zh) * | 2018-11-06 | 2019-03-19 | 杭州君谋科技有限公司 | 一种基于机器学习的手机振动马达良次品检测方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9020271B2 (en) * | 2012-07-31 | 2015-04-28 | Hewlett-Packard Development Company, L.P. | Adaptive hierarchical clustering algorithm |
EP3048972B1 (en) * | 2013-09-27 | 2022-11-09 | Koninklijke Philips N.V. | Processing apparatus, processing method and system for processing a physiological signal |
KR102257300B1 (ko) * | 2014-09-01 | 2021-05-27 | 삼성전자주식회사 | 심전도 신호를 이용한 사용자 인증 방법 및 장치 |
US11457872B2 (en) * | 2017-12-01 | 2022-10-04 | Samsung Electronics Co., Ltd. | Bio-signal quality assessment apparatus and bio-signal quality assessment method |
-
2019
- 2019-07-23 CN CN201910667483.3A patent/CN110403582B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1662921A (zh) * | 2002-06-24 | 2005-08-31 | 皇家飞利浦电子股份有限公司 | 比较各种初始群组集合以便确定用于群集一个电视节目集合的最佳初始集合的方法 |
JPWO2004097703A1 (ja) * | 2003-05-02 | 2006-07-13 | 富士通株式会社 | 製造工程解析支援方法、その方法をコンピュータに実行させるプログラム、プログラムプロダクトおよび記録媒体 |
CN101120615A (zh) * | 2005-02-22 | 2008-02-06 | 弗劳恩霍夫应用研究促进协会 | 近透明或透明的多声道编码器/解码器方案 |
CN102112984A (zh) * | 2008-07-29 | 2011-06-29 | 皇家飞利浦电子股份有限公司 | 用于生成图像选集的方法和设备 |
CN102460433A (zh) * | 2009-06-15 | 2012-05-16 | 皇家飞利浦电子股份有限公司 | 选择代表图像的方法和设备 |
CN102547844A (zh) * | 2010-12-29 | 2012-07-04 | 通用电气公司 | 无线网络中的动态数据管理的系统和方法 |
CN104090861A (zh) * | 2013-02-20 | 2014-10-08 | 哈佛蒸汽锅炉检验和保险公司 | 动态离群值偏倚减少系统和方法 |
WO2015107268A1 (en) * | 2014-01-16 | 2015-07-23 | Aboa Legis Oy | Method and device for the detection of respiratory rate |
CN104586383A (zh) * | 2015-01-12 | 2015-05-06 | 北京世纪今科医疗器械有限公司 | 一种心电波形分类方法及装置 |
CN105411567A (zh) * | 2015-12-15 | 2016-03-23 | 北京世纪今科医疗器械有限公司 | 一种心电图数据分析方法及系统 |
CN106963353A (zh) * | 2015-12-30 | 2017-07-21 | 朝阳科技大学 | 评估动脉血流介导血管舒张功能的方法 |
CN105686818A (zh) * | 2016-03-30 | 2016-06-22 | 华为技术有限公司 | 测量装置及方法 |
CN106254175A (zh) * | 2016-07-26 | 2016-12-21 | 北京蓝海讯通科技股份有限公司 | 一种集群异常节点检测方法、应用和计算设备 |
CN107423757A (zh) * | 2017-07-14 | 2017-12-01 | 北京小米移动软件有限公司 | 聚类处理方法及装置 |
CN108897821A (zh) * | 2018-06-21 | 2018-11-27 | 江苏赛睿信息科技股份有限公司 | 自动生成数据结论的方法及装置 |
CN109344146A (zh) * | 2018-09-13 | 2019-02-15 | 海南电网有限责任公司电力科学研究院 | 一种大规模过电压数据智能自清洗方法 |
CN109490776A (zh) * | 2018-11-06 | 2019-03-19 | 杭州君谋科技有限公司 | 一种基于机器学习的手机振动马达良次品检测方法 |
Non-Patent Citations (2)
Title |
---|
"A New Dissimilarity Measure for Clustering Seismic Signals";Benvegna;《IMAGE ANALYSIS AND PROCESSING》;20110916;第434-438页 * |
"穿戴式人体生理参数监测系统的研究与实现";郭维;《中国博士论文全文数据库医药卫生科技辑》;20130215;第E06-129页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110403582A (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105411567B (zh) | 一种心电图数据分析方法及系统 | |
KR101524226B1 (ko) | 신경망을 이용한 심장질환판별 방법 및 그 장치 | |
Meyer et al. | Combining algorithms in automatic detection of QRS complexes in ECG signals | |
CN101828918B (zh) | 基于波形特征匹配的心电信号r波峰检测方法 | |
CN104905785B (zh) | 心电图信号检测方法和系统 | |
CN108324271B (zh) | 心电信号识别方法、系统和心电监测设备 | |
WO2019234458A1 (en) | Detecting abnormalities in ecg signals | |
CA2861500A1 (en) | Methods and systems for atrial fibrillation detection | |
CN107622259B (zh) | 一种t波检测方法、心电数据分析方法及装置 | |
CN113057648A (zh) | 一种基于复合lstm结构的ecg信号分类方法 | |
CN103970975A (zh) | 心电数据处理方法和系统 | |
US11690526B2 (en) | Device for extraction of peaks from continuous intracranial pressure waveform and detection method for the same | |
CN110403582B (zh) | 一种用于分析脉波波形品质的方法 | |
CN104983415B (zh) | 心拍识别的方法及装置 | |
JP5688309B2 (ja) | 生体情報分析装置、生体情報分析方法及び生体情報分析プログラム | |
Khaustov et al. | An open-source standard T-wave alternans detector for benchmarking | |
CN116473569A (zh) | 一种心电图数据处理方法、系统及存储介质 | |
CN111345815B (zh) | 一种心电信号中qrs波的检测方法、装置、设备和存储介质 | |
EP2928363B1 (en) | Apparatus and method for determining the occurrence of a qrs complex in ecg data | |
CN114711780A (zh) | 多导联心电图信号处理方法、装置、设备及存储介质 | |
CN108836312B (zh) | 一种基于人工智能的进行杂波剔除的方法及系统 | |
RU2410023C1 (ru) | Способ выделения qrs-комплекса электрокардиосигнала | |
CN107170467B (zh) | 一种基于Lempel-Ziv复杂度的异常心音识别方法 | |
CN111345814A (zh) | 一种心电信号中心拍的分析方法、装置、设备和存储介质 | |
CN114652323B (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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220221 Address after: 523857 room 701, No. 284, Wusha Huannan Road, Chang'an Town, Dongguan City, Guangdong Province Patentee after: Dongguan Yingyu Technology Consulting Co.,Ltd. Address before: Room 404, 4 / F, building 10, Dongguan Taiwan Biotechnology cooperation center, No.1 Taoyuan Road, Taiwan hi tech park, Songshanhu hi tech Industrial Development Zone, Dongguan City, Guangdong Province, 523000 Patentee before: Hongrenyi medical equipment (Dongguan) Co.,Ltd. |
|
TR01 | Transfer of patent right |