CN113545753A - 波形检测方法、装置、可穿戴设备及计算机可读存储介质 - Google Patents
波形检测方法、装置、可穿戴设备及计算机可读存储介质 Download PDFInfo
- Publication number
- CN113545753A CN113545753A CN202110827269.7A CN202110827269A CN113545753A CN 113545753 A CN113545753 A CN 113545753A CN 202110827269 A CN202110827269 A CN 202110827269A CN 113545753 A CN113545753 A CN 113545753A
- Authority
- CN
- China
- Prior art keywords
- waveform
- time
- signal
- waveform signal
- real
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 141
- 238000000034 method Methods 0.000 claims abstract description 57
- 238000001914 filtration Methods 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 239000002131 composite material Substances 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 8
- 230000003139 buffering effect Effects 0.000 claims description 2
- 238000011897 real-time detection Methods 0.000 abstract description 33
- 238000012545 processing Methods 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 19
- 238000010586 diagram Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 238000007781 pre-processing Methods 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000009466 transformation Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003183 myoelectrical effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000006467 substitution reaction Methods 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- 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
-
- 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- Physiology (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biophysics (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Psychiatry (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Power Engineering (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种波形检测方法、装置、可穿戴设备及计算机可读存储介质。该方法先将采集到的波形信号缓存一定时间后,根据预设的基准参数对其进行检测,得到波形特征参数(盲估阶段),然后根据该波形特征参数对实时采集到的波形信号进行检测(近实时检测阶段)。检测得到的波形特征参数可避免个体差异对检测结果造成的误差影响,进而可提高检测结果的准确性。同时,在近实时检测阶段,每缓冲1.5个信号周期的波形信号即可进行一次检测,进而可以近实时的实现对持续性输入的波形信号的波峰波谷点定位。此外,本方案对波形信号进行滤波时采用的滤波算法的计算量较小,可降低对设备的运算能力的要求,适合在处理能力较弱的可穿戴式设备上运行。
Description
技术领域
本申请属于信号处理领域,具体涉及一种波形检测方法、装置、可穿戴设备及计算机可读存储介质。
背景技术
在监测领域中,可以通过传感器去采集一些数据,从而获取到由人体所产生的波信号。例如可以利用LED发出光信号,并使用对应的探测单元采集经过人体皮肤、组织和血管反射、吸收的衰减光,从而可以获得反应血管内血流变化的脉搏波信号;再例如可以通过探测单元去采集用于反应人体心率变化的心电图波信号。
在得到波形信号之后,需要对波形信号进行检测,以便根据检测结果分析出与人体有关的生理信息。
在现有技术中,一般是以健康人体的波形信号参数作为参考数据,并基于此来对待检测的波形信号进行检测。然而,由于年龄、种族、身体状况等因素所导致的个体差异,以健康人体的波形信号参数作为统一的参考数据将会导致待检测的波形信号的检测结果不准确。
发明内容
有鉴于此,本申请的目的在于提供一种波形检测方法、装置、可穿戴设备及计算机可读存储介质,可以提高波形检测结果准确度。
本申请的实施例是这样实现的:
第一方面,本申请实施例提供一种波形检测方法,所述方法包括:S1、对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;S2、根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
在本实施例中,先将采集到的波形信号缓存一定时间(第一时间段)后,根据预设的基准参数对波形信号进行检测,得到波形特征参数(盲估阶段),然后再根据波形特征参数来对实时采集到的波形信号进行近实时的检测(即近实时检测阶段),而不是如现有技术所示的根据大数据所得到的固定参数对实时采集到的波形信号进行检测,那么相较于现有技术,本申请实施例至少可以避免个体差异对检测结果造成的误差影响,进而可以提高检测结果的准确性。
结合第一方面实施例,在一种可能的实施方式中,在S2根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷之后,所述方法还包括:S3、根据与所述实时采集到的波形信号对应的波峰以及波谷,得到当前的检测参数;根据当前的所述检测参数,更新所述波形特征参数,并重复S2。
在本实施例中,在根据盲估阶段得到波形特征参数后,便进入近实时检测阶段。在近实时检测阶段时,每次检测出的波峰以及波谷均可以生成对应的检测参数,并以当前的检测参数实时更新当前的波形特征参数,以便下一次实时检测时,可以以更新后的波形特征参数为基准来进行检测。在本申请实施例中,考虑到同一个检测对象在不同时间所产生的波形信号有差异性,因此,可以通过实时更新波形特征参数的方式,来使得波形特征参数的准确性以及时效性不断得到优化,进而有利于提高后续检测结果的准确性。
结合第一方面实施例,在一种可能的实施方式中,所述S1对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数,包括:通过滤波器过滤所述缓存第一时间段后的波形信号的噪声,得到处理后的缓存波形信号;对所述处理后的缓存波形信号的信号值进行识别,得到与所述缓存第一时间段后的波形信号对应的信号极值序列、信号最大值、信号最小值、最大信号振幅、波峰以及波谷;根据确定出的所有所述波峰以及所有所述波谷,确定出所述波形特征参数所包括的信号周期;其中,与所述缓存第一时间段后的波形信号对应的信号极值序列中,第一个满足第一波峰条件的极大值为第一个波峰,其余满足第二波峰条件的极大值为其他波峰,第一个满足第一波谷条件的极小值为第一个波谷,其余满足第二波谷条件的极小值为其他波谷;所述第一波峰条件为:所述极大值与所述信号最大值的比值大于预设阈值;所述第二波峰条件为:所述极大值与前一个已确定出的相邻波谷的幅值之差与所述最大信号振幅的比值大于所述预设阈值;所述第一波谷条件为:所述极小值与所述信号最小值的比值大于预设阈值;所述第二波谷条件为:所述极小值与前一个已确定出的相邻波峰的幅值之差与所述最大信号振幅的比值大于所述预设阈值。
结合第一方面实施例,在一种可能的实施方式中,所述通过滤波器过滤所述缓存第一时间段后的波形信号的噪声,得到处理后的缓存波形信号,包括:通过高阶均值滤波器去除所述缓存第一时间段后的波形信号的基线漂移,得到校准缓存波形信号;通过低阶均值滤波器对所述校准缓存波形信号进行滤波,得到所述处理后的缓存波形信号;其中,所述高阶均值滤波器的阶数等于所述基准参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
在本实施例中,通过采用二次均值滤波的方式对波形信号进行预处理,可以减少计算量和对存储空间的需求,且不需要浮点支持,适合在计算能力和存储空间有限的设备上运行。此外,采用二次均值滤波的方式,能够较好地去除基线漂移。
结合第一方面实施例,在一种可能的实施方式中,在所述对实时采集到的波形信号进行检测之前,所述方法还包括:确定本次检测的波形信号对象,本次检测的所述波形信号对象的产生时间起点为前一次检测得到的波峰或波谷所在的时间点,其产生时间终点为距离所述当前时刻之后第二时间段时所对应的时间点;
相应的,所述对实时采集到的波形信号进行检测,包括:对实时采集到的波形信号缓存所述第二时间段后,得到本次检测的所述波形信号对象,且对其进行检测;其中,所述第一时间段对应的时长大于所述第二时间段对应的时长。
在本实施例中,在近实时检测阶段,通过对实时采集到的波形信号缓存第二时间段后再进行检测,可以保证每次近实时检测时的波形信号能够覆盖一个完整的真实信号周期。
结合第一方面实施例,在一种可能的实施方式中,所述对实时采集到的波形信号进行检测,得到与所述实时采集到的实时波形信号对应的波峰以及波谷,包括:通过滤波器过滤所述实时采集到的波形信号的噪声,得到处理后的实时波形信号;对所述处理后的实时波形信号的信号值进行识别,得到与所述实时采集到的波形信号对应的信号极值序列;根据所述波形特征参数所包括的信号周期以及前一次检测所得到的最后一组波峰波谷所出现的时间坐标,预测得到所述实时采集到的波形信号中的下一组波峰波谷所出现的时间坐标范围;获取与所述实时采集到的波形信号对应的信号极值序列中其时间坐标属于所述时间坐标范围内的极值子序列;将所述极值子序列所包括的极大值确定为所述波峰,所述极值子序列所包括的极小值确定为所述波谷。
结合第一方面实施例,在一种可能的实施方式中,所述将所述极值子序列所包括的极大值确定为所述波峰,包括:在确定所述极值子序列内包括多个极大值时,计算所述多个极大值中的每个极大值的波峰综合评分,并将所述波峰综合评分最小时所对应的极大值确定为所述波峰;
相应的,所述将极值子序列所包括的极小值确定为所述波谷,包括:在确定所述极值子序列内包括多个极小值时,计算所述多个极小值中的每个极小值的波谷综合评分,并将所述波谷综合评分最小时所对应的极小值确定为所述波谷;其中,所述波峰综合评分用于表征当前的计算对象作为波峰时所得到的检测参数与所述波形特征参数之间的差异性;所述波谷综合评分用于表征当前的计算对象作为波谷时所得到的检测参数与所述波形特征参数之间的差异性。
结合第一方面实施例,在一种可能的实施方式中,所述通过滤波器过滤所述实时采集到的波形信号的噪声,包括:通过高阶均值滤波器去除所述实时采集到的波形信号的基线漂移,得到校准实时波形信号;通过低阶均值滤波器对所述校准实时波形信号进行滤波,得到所述处理后的实时波形信号;其中,所述高阶均值滤波器的阶数等于所述波形特征参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
在本实施例中,通过采用二次均值滤波的方式对波形信号进行预处理,可以减少计算量和对存储空间的需求,且不需要浮点支持,适合在计算能力和存储空间有限的设备上运行。此外,采用二次均值滤波的方式,能够较好地去除基线漂移。
结合第一方面实施例,在一种可能的实施方式中,检测对象为人体,所述波形信号为光电容积脉搏波信号。
第二方面,本申请实施例提供一种波形检测装置,所述装置包括:第一检测模块以及第二检测模块。
第一检测模块,用于对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;
第二检测模块,用于根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述采集到的实时波形信号对应的波峰以及波谷。
结合第二方面实施例,在一种可能的实施方式中,所述装置还包括更新模块,用于根据与所述实时采集到的波形信号对应的波峰以及波谷,得到当前的检测参数;根据当前的所述检测参数,更新所述波形特征参数。
结合第二方面实施例,在一种可能的实施方式中,所述第一检测模块,用于通过滤波器过滤所述缓存第一时间段后的波形信号的噪声,得到处理后的缓存波形信号;对所述处理后的缓存波形信号的信号值进行识别,得到与所述缓存第一时间段后的波形信号对应的信号极值序列、信号最大值、信号最小值、最大信号振幅、波峰以及波谷;根据确定出的所有所述波峰以及所有所述波谷,确定出所述波形特征参数所包括的信号周期;其中,与所述缓存第一时间段后的波形信号对应的信号极值序列中,第一个满足第一波峰条件的极大值为第一个波峰,其余满足第二波峰条件的极大值为其他波峰,第一个满足第一波谷条件的极小值为第一个波谷,其余满足第二波谷条件的极小值为其他波谷;所述第一波峰条件为:所述极大值与所述信号最大值的比值大于预设阈值;所述第二波峰条件为:所述极大值与前一个已确定出的相邻波谷的幅值之差与所述最大信号振幅的比值大于所述预设阈值;所述第一波谷条件为:所述极小值与所述信号最小值的比值大于预设阈值;所述第二波谷条件为:所述极小值与前一个已确定出的相邻波峰的幅值之差与所述最大信号振幅的比值大于所述预设阈值。
结合第二方面实施例,在一种可能的实施方式中,所述第一检测模块,用于通过高阶均值滤波器去除所述缓存第一时间段后的波形信号的基线漂移,得到校准缓存波形信号;通过低阶均值滤波器对所述校准缓存波形信号进行滤波,得到所述处理后的缓存波形信号;其中,所述高阶均值滤波器的阶数等于所述基准参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
结合第二方面实施例,在一种可能的实施方式中,所述装置还包括确定模块,用于确定本次检测的波形信号对象,本次检测的所述波形信号对象的产生时间起点为前一次检测得到的波峰或波谷所在的时间点,其产生时间终点为距离所述当前时刻之后第二时间段时所对应的时间点;
相应的,所述第二检测模块,用于对实时采集到的波形信号缓存所述第二时间段后,得到本次检测的所述波形信号对象,且对其进行检测;其中,所述第一时间段对应的时长大于所述第二时间段对应的时长。
结合第二方面实施例,在一种可能的实施方式中,所述第二检测模块,用于通过滤波器过滤所述实时采集到的波形信号的噪声,得到处理后的实时波形信号;对所述处理后的实时波形信号的信号值进行识别,得到与所述实时采集到的波形信号对应的信号极值序列;根据所述波形特征参数所包括的信号周期以及前一次检测所得到的最后一组波峰波谷所出现的时间坐标,预测得到所述实时采集到的波形信号中的下一组波峰波谷所出现的时间坐标范围;获取与所述实时采集到的波形信号对应的信号极值序列中其时间坐标属于所述时间坐标范围内的极值子序列;将所述极值子序列所包括的极大值确定为所述波峰,所述极值子序列所包括的极小值确定为所述波谷。
结合第二方面实施例,在一种可能的实施方式中,所述第二检测模块,用于在确定所述极值子序列内包括多个极大值时,计算所述多个极大值中的每个极大值的波峰综合评分,并将所述波峰综合评分最小时所对应的极大值确定为所述波峰;
相应的,还用于在确定所述极值子序列内包括多个极小值时,计算所述多个极小值中的每个极小值的波谷综合评分,并将所述波谷综合评分最小时所对应的极小值确定为所述波谷;
其中,所述波峰综合评分用于表征当前的计算对象作为波峰时所得到的检测参数与所述波形特征参数之间的差异性;所述波谷综合评分用于表征当前的计算对象作为波谷时所得到的检测参数与所述波形特征参数之间的差异性。
结合第二方面实施例,在一种可能的实施方式中,所述第二检测模块,用于通过高阶均值滤波器去除所述实时波形信号的基线漂移,得到校准实时波形信号;通过低阶均值滤波器对所述校准实时波形信号进行滤波,得到所述处理后的实时波形信号;其中,所述高阶均值滤波器的阶数等于所述波形特征参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
结合第二方面实施例,在一种可能的实施方式中,检测对象为人体,所述波形信号为光电容积脉搏波信号。
第三方面,本申请实施例还提供一种可穿戴设备,包括:传感器、存储器和处理器,所述传感器、所述存储器分别和所述处理器连接;
所述传感器,用于采集波形信号;
所述处理器,用于根据预先保存的基准参数,对缓存第一时间段后的波形信号进行检测,得到波形特征参数;
所述处理器,还用于根据所述波形特征参数,对所述传感器所实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
结合第三方面实施例,在一种可能的实施方式中,所述处理器,还用于:根据所述波形特征参数,对所述传感器所实时采集到的且缓存第二时间段后的波形信号进行检测,得到与所述实时采集到的且缓存第二时间段后的波形信号对应的波峰以及波谷;所述第一时间段大于所述第二时间段。
结合第三方面实施例,在一种可能的实施方式中,所述可穿戴设备为手环。
第四方面,本申请实施例还提供一种非易失性计算机可读取存储介质(以下简称计算机可读存储介质),其上存储有计算机程序,所述计算机程序被计算机运行时执行上述第一方面实施例和/或结合第一方面实施例的任一种可能的实施方式提供的方法。
本申请的其他特征和优点将在随后的说明书阐述,并且,部分地从说明书中变得显而易见,或者通过实施本申请实施例而了解。本申请的目的和其他优点可通过在所写的说明书以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。通过附图所示,本申请的上述及其它目的、特征和优势将更加清晰。在全部附图中相同的附图标记指示相同的部分。并未刻意按实际尺寸等比例缩放绘制附图,重点在于示出本申请的主旨。
图1示出本申请实施例提供的一种波形检测方法的流程图之一。
图2示出本申请实施例提供的一种波形信号的示意图之一。
图3示出本申请实施例提供的一种波形信号的示意图之二。
图4示出本申请实施例提供的一种波形检测方法的流程图之二。
图5示出本申请实施例提供的一种波形检测装置的结构框图。
图6示出本申请实施例提供的一种可穿戴设备的结构示意图之一。
图7示出本申请实施例提供的一种可穿戴设备的结构示意图之二。
附图标记:100-可穿戴设备;110-处理器;120-存储器;130-传感器;140-通信模块;150-显示器;400-波形检测装置;410-第一检测模块;420-第二检测模块。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行描述。
应注意到:相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。同时,在本申请的描述中诸如“第一”、“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
再者,本申请中术语“和/或”,仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。
此外,针对现有技术中的波形检测方案所存在的缺陷是申请人在经过实践并仔细研究后得出的结果,因此,上述缺陷的发现过程以及在下文中本申请实施例针对上述缺陷所提出的解决方案,都应该被认定为是申请人对本申请做出的贡献。
为了解决上述问题,本申请实施例提供一种波形检测方法、装置、可穿戴设备及计算机可读存储介质,可以提高检测的准确性。
该技术可采用相应的软件、硬件以及软硬结合的方式实现。以下对本申请实施例进行详细介绍。
首先,本申请实施例提供一种波形检测方法。该方法可以应用于可穿戴设备,也可以应用于其他具备数据处理功能的电子设备。下面将结合图1对其所包含的步骤进行说明。
步骤S110:对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数。
步骤S120:根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
在本申请实施例中,当需要对检测对象所产生实时波形信号进行检测时,为了提高检测结果的准确性,可以先将采集到的波形信号缓存一定时间(第一时间段)后,根据预设的基准参数对其进行检测,得到检测对象的波形特征参数(即与步骤S110对应的盲估阶段),然后再根据该波形特征参数对实时采集到的波形信号进行近实时的检测(即与步骤S120对应的近实时检测阶段),而不是如现有技术所示的根据大数据所得到的固定参数对检测对象所产生的实时波形信号进行检测,那么相较于现有技术,本申请实施例至少可以避免个体差异对检测结果造成的误差影响,进而可以提高检测结果的准确性。
下面将针对各个步骤进行详细介绍。
步骤S110:对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数。
其中,产生波形信号的检测对象可以是人体,也可以是其他动物体。
在一些实施方式中,当每次需要执行波形检测任务时,均需要通过盲估阶段来计算本次检测对象的波形特征参数,以便后续基于该波形特征数据来对检测对象所产生的实时波形信号进行检测。
下面将针对盲估阶段计算本次检测对象的波形特征参数的过程进行介绍。
盲估阶段主要是将采集到的波形信号缓存一定时间(第一时间段)后,结合预先保存的基准参数来对缓存第一时间段后的波形信号(以下简称为缓存的波形信号)进行波峰点和/或波谷点检测,进而根据检测结果得到与检测对象相适配的,且可信度相对较高的波形特征参数。
其中,缓存的波形信号是属于一段时间内的持续性信号。
缓存的波形信号可以通过采集检测对象在当前所产生的波形信号得到。在一些实施方式中,可以实时地采集检测对象所产生的实时波形信号,并将采集第一时间段后所得到的属于这段时间内的持续性波形信号确定为缓存的波形信号。
至于波形信号需要缓存的最小长度,可以预先进行配置,例如在一些实施方式中,该缓存的波形信号的持续时间为8秒,以使得在缓存的波形信号中至少保存预设数值(例如5)个完整的信号周期。
至于预先保存的基准参数,可以包括信号周期。该信号周期一般是根据健康人体的波形信号大数据所得到的平均参数,对于检测对象而言,其可信度较低。
值得指出的是,该健康人体的波形信号与检测对象所产生的波形信号属于同一种类型的波形信号,例如都为光电容积脉搏波(PPG,Photo Plethysmo Graphy)信号。
在确定出基准参数以及缓存的波形信号后,即可根据基准参数,针对缓存的波形信号进行检测,进而得到可信度相对较高的波形特征参数。
其中,盲估阶段的检测过程主要包括三个阶段,分别为信号预处理、极值检测以及极值过滤。
由于在针对人体进行波形信号采集时,采集到的波形信号容易受到采集设备的位置、人体的运动状态、肌电干扰等众多噪声的影响,进而使得缓存波形信号中存在大量的噪声,导致缓存波形信号中存在基线漂移。
盲估阶段的信号预处理主要是通过滤波器去除缓存的波形信号中所存在的基线漂移。
在一些实施方式中,可以通过带通滤波、小波变换等现有方法去除基线漂移,进而得到处理后的缓存波形信号。但是此类方式难以去除频率在缓存波形信号的频率范围内的基线漂移,会严重影响后续的检测结果的准确性;此外,此类滤波算法的计算量较大,对设备的运算能力有一定的要求,不适合在处理能力较弱的可穿戴式设备上运行。
在另一些实施方式中,可以采用均值滤波的方式来去除基线漂移和高频噪声。均值滤波器的阶数由基准参数所确定。
具体的,可以先通过高阶均值滤波器(该滤波器为软件滤波器)去除缓存的波形信号的基线漂移,得到校准缓存波形信号。然后再通过低阶均值滤波器对校准缓存波形信号进行滤波,得到处理后的缓存波形信号。
其中,高阶均值滤波器的主要作用是对缓存的波形信号的基线进行预测,在预测出基线后,将缓存的波形信号的各个信号值与预测得到的基线信号相减,即可得到去除基线偏移的校准缓存波形信号。
其中,高阶均值滤波器的阶数等于基准参数所包括的信号周期,例如当针对PPG信号进行检测时,基准参数所包括的信号周期用于表征健康成人的平均脉搏周期(75次/分钟),那么可以将高阶均值滤波器的阶数确定为75(尽管这会导致波形失真,但是对波形信号的周期和振幅的影响在可接收范围内)。
针对得到的校准缓存波形信号,极有可能仍然存在大量的高频噪声,因此,再通过一个低阶均值滤波器对其进行二次滤波,以此来平滑校准缓存波形信号中的高频噪声,得到处理后的缓存波形信号。
其中,低阶均值滤波器的阶数等于预设值。
该预设值的确定主要与缓存的波形信号的采样频率以及其噪声频率等因素有关,需要根据实际情况进行调整,以便能够去除校准缓存波形信号中的高频噪声的同时,还能够尽可能避免波形失真。
在一实施例中,预设值的取值范围在3-5之间。
相较于传统的带通滤波、小波变换等方式,均值滤波的计算量更小,可以降低对设备的运算能力的要求,适合在处理能力较弱的可穿戴式设备上运行。此外,相较于传统的带通滤波、小波变换等方式,均值滤波具有更优的基线漂移去除效果,有利于提高后续检测结果的准确性。
其中,照图2以及图3分别示出原始的缓存的波形信号以及经过均值滤波后所得到的处理后的缓存波形信号。
在得到处理后的缓存波形信号之后,即可进入盲估阶段的极值检测阶段。
在本申请实施例中,可以针对处理后的缓存波形信号的信号值进行识别,得到与缓存的波形信号对应的信号极值序列、信号最大值、信号最小值以及最大信号振幅。
其中,如图3所示,经过处理后的缓存波形信号存在对应的信号值,其数字表现形式即为一段数字序列。
在本申请实施例中,可以针对数字序列中的每个数字进行求导,得到对应的导数。其中,导数为0的数值,即为处理后的缓存波形信号中的极值。
在得到处理后的缓存波形信号中的所有极值后,其所组成的序列即为信号极值序列。信号最大值为所有极值中的最大值,信号最小值为所有极值中的最小值。相邻的极大值与极小值之间的连线距离即为信号振幅,通过比较,即可得到最大信号振幅。
在得到与缓存的波形信号对应的信号极值序列、信号最大值、信号最小值以及最大信号振幅之后,即可进入盲估阶段的极值过滤阶段。
在盲估阶段中,由于缺少可信的波形特征参数,因此只能利用处理后的缓存波形样点的幅值、样点的极性(信号值的正负)以及振幅信息来检测与缓存的波形信号对应的信号极值序列中的极值点是否为波峰或波谷点来实现极值过滤。
具体的,可以依次检索与缓存的波形信号对应的信号极值序列中的每一个极值。
其中,针对极大值点,计算该极大值与信号最大值的比值;还将该极大值点与前一个已确定出的相邻波谷的信号值(幅值)之差确定为该极大值点的振幅,并计算其振幅与最大信号振幅的比值。若两个比值都大于预设阈值(例如为40%),那么该极大值点即为波峰。
当然,值得指出的是,在计算第一个波峰时,由于可能不存在前一个已确定出的相邻波谷,那么可以将极大值与信号最大值的比值大于预设阈值的极大值点确定为第一个波峰。
针对极小值点,计算该极小值与信号最小值的比值;还将该极小值点与前一个已确定出的相邻波峰的信号值(幅值)之差确定为该极小值点的振幅,并计算其振幅与最小信号幅值的比值。若两个比值都大于预设阈值(例如为40%),那么该极小值点即为波谷。
当然,值得指出的是,在计算第一个波峰时,由于可能不存在前一个已确定出的相邻波谷,那么可以将极大值与信号最大值的比值大于预设阈值的极大值点确定为第一个波峰。
在得到与缓存的波形信号对应的信号极值序列中的波峰波谷后,即可根据相邻波峰或相邻波谷之间的时间间距,确定出相对可信的波形特征参数,例如信号周期、信号振幅等参数。
步骤S120:根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
在通过盲估阶段获取到波形特征参数之后,即可基于波形特征参数对用户实时所产生的波形信号进行检测,即进入近实时检测阶段。
值得指出的是,在本申请实施例中,近实时检测阶段(即步骤S120)是一个随着时间的推移而持续性重复进行的过程。在每次近实时检测过程中,其检测的波形信号对象为用户在一段时间内所产生的波形信号。
可选的,在近实时检测阶段,为了保证每次近实时检测时的波形能够覆盖一个完整的真实信号周期,在一些实施方式中,在每次进行近实时检测之前,可以先确定本次检测的波形信号对象。
其中,本次检测的波形信号对象为检测对象所产生的一段波形信号,且该一段波形信号的产生时间起点为前一次检测(存在前一次实时检测过程,那么此处的前一次检测为前一次实时检测,否则,此处的前一次检测指的是盲估阶段对缓存波形信号的检测)得到的波峰或波谷所在的时间点,该一段波形信号的产生时间终点为距离当前时刻之后第二时间段时所对应的时间点。
其中,该第二时间段的时长尽量保证波形信号对象的波形可以覆盖1.5个信号周期(该信号周期为参考数据所包括的信号周期),以保证当人体所产生的实时波形信号的真实信号周期变长(人体所产生的波形信号存在波动,相应的,其信号周期也存在波动)时,所得到的波形信号对象依旧能够覆盖一个完整的真实信号周期。
相应的,在近实时检测阶段,当对接收到的由检测对象所产生的实时波形信号进行检测时,可以对实时采集到的波形信号缓存第二时间段后,得到本次检测的波形信号对象,然后再针对本次检测的波形信号对象进行检测。
当然,值得指出的是,盲估阶段时的第一时间段对应的时长大于近实时阶段时的第二时间段对应的时长。
其中,当波形信号为PPG信号,且该第二时间段的时长为1.5个信号周期时,相当于在近实时检测阶段,每缓冲1.5个信号周期的波形信号后即可进行一次近实时检测,对于用户而言,可以得到时延在2秒内的脉搏波数据,进而可以实现对持续性输入的波形信号的近实时的波峰波谷点定位。
当然,与盲估阶段类似的,每次近实时检测阶段主要包括三个阶段,分别为信号预处理、极值检测以及极值过滤。
其中,针对近实时检测阶段的信号预处理,其处理过程与盲估阶段的信号预处理相似,可以通过带通滤波、小波变换等现有方法去除基线漂移,进而得到处理后的实时波形信号;也可以通过均值滤波的方式来去除实时波形信号中的基线漂移和高频噪声,进而得到处理后的实时波形信号。均值滤波器的阶数由波形特征参数所确定。
值得指出的是,当采用均值滤波的方式来去除实时采集到的波形信号中的基线漂移和高频噪声时,可以先通过高阶均值滤波器去除采集到的实时波形信号的基线漂移,得到校准实时波形信号。然后再通过低阶均值滤波器对校准实时波形信号进行滤波,得到处理后的实时波形信号。
高阶均值滤波器的阶数等于相对可信的波形特征参数所包括的信号周期,低阶均值滤波器的阶数等于预设值。
该预设值的确定主要与实时采集到的波形信号的采样频率以及其噪声频率等因素有关,需要根据实际情况进行调整,以便能够去除校准实时波形信号中的高频噪声的同时,还能够尽可能避免波形失真。
一般而言,预设值的取值范围在3-5之间。
针对近实时检测阶段的极值检测阶段,与盲估阶段的极值检测阶段类似的,可以针对处理后的实时波形信号的信号值进行识别,得到与实时采集到的波形信号对应的信号极值序列。
针对近实时检测阶段的极值过滤阶段,可以先确定出前一次检测所得到的最后一组波峰波谷所出现的时间坐标。
其中,前文提及,近实时检测阶段是一个随着时间的推移而持续性重复进行的过程。
若本次近实时检测过程是该第一次近实时检测过程,那么上文中所涉及到的前一次检测所得到的最后一组波峰波谷是盲估阶段(若存在)所检测出的最后一组波峰波谷,或者默认为波形信号对象的起始位置为波峰波谷。
若本次近实时检测过程不是第一次近实时检测过程,那么上文中所涉及到的前一次检测所得到的最后一组波峰波谷,是相较于本次而言的前一次近实时检测阶段所检测出的最后一组波峰波谷。
在确定出前一次检测所得到的最后一组波峰波谷所出现的时间坐标之后,可以根据波形特征参数所包括的信号周期以及前一次检测所得到的最后一组波峰波谷所出现的时间坐标,预测得到实时波形信号中的下一组波峰波谷所出现的时间坐标范围。
假设波形特征参数所包括的信号周期为T,最后一组波峰所出现的时间坐标为S1,最后一组波谷所出现的时间坐标为S2,那么下一组波峰所出现的时间坐标范围为[(1-X)×(S1+T),(1+X)×(S1+T)],下一组波谷所出现的时间坐标范围为[(1-X)×(S2+T),(1+X)×(S2+T)]。
其中,X大于0小于1,例如可以为0.2。
在得到时间坐标范围后,可以获取与近实时波形信号对应的信号极值序列中其时间坐标属于时间坐标范围内的极值子序列,然后将极值子序列中所包括的极大值确定为下一组波峰,将极值子序列中所包括的极小值确定为下一组波谷。
当然,在一些实施方式中,在同一个极值子序列中可能存在多个极大值和/或多个极小值。此时,需要通过波峰综合评分来从多个极大值中筛选出一个极大值,并将其确定为下一组波峰;和/或,通过波谷综合评分来从多个极小值中筛选出一个极小值,并将其确定为下一组波谷。
可选的,若存在多个极大值,可以计算多个极大值中的每个极大值的波峰综合评分。其中,波峰综合评分用于表征当前的计算对象(极大值)作为波峰时所得到的检测参数(包括但不限于信号周期、信号振幅等)与波形特征参数(包括但不限于信号周期、信号振幅等)之间的差异性。
在得到多个波峰综合评分后,可以将波峰综合评分最小时所对应的极大值确定为下一组波峰。
在一种可选的实施方式中,波峰综合评分的计算规则为对A1、B1以及C1三项进行加权求和。
其中,A1用于表征:当前的计算对象(极大值)作为波峰时所得到的信号周期与波形特征参数所包括的信号周期之差;B1用于表征:当前的计算对象(极大值)作为波峰时所得到的信号振幅与波形特征参数所包括的信号振幅之差;C1用于表征当前的计算对象(极大值)的信号值与波形特征参数所包括的波峰的信号值之差。
可选的,若存在多个极小值,可以计算多个极小值中的每个极小值的波谷综合评分。其中,波谷综合评分用于表征当前的计算对象(极小值)作为波谷时所得到的检测参数(包括但不限于信号周期、信号振幅等)与波形特征参数(包括但不限于信号周期、信号振幅等)之间的差异性。
在得到多个波谷综合评分后,可以将波谷综合评分最小时所对应的极小值确定为下一组波谷。
在一种可选的实施方式中,波谷综合评分的计算规则为对A2、B2以及C2三项进行加权求和。
其中,A2用于表征:当前的计算对象(极小值)作为波谷时所得到的信号周期与波形特征参数所包括的信号周期之差;B2用于表征:当前的计算对象(极小值)作为波谷时所得到的信号振幅与波形特征参数所包括的信号振幅之差;C2用于表征当前的计算对象(极小值)的信号值与波形特征参数所包括的波谷的信号值之差。
在得到下一组波峰波谷之后,即可重复上述过程来执行下一次近实时检测过程。
此外,在一些实施方式中,在每次得到下一组波峰波谷之后,还可以根据当前的下一组波峰波谷,得到当前的检测参数(包括但不限于信号周期、信号振幅等),并根据当前的检测参数更新原有的波形特征参数,进而优化波形特征参数的准确性。
其中,更新的过程可以是计算两者之间的平均值。
具体的,请参照图4,可以在每执行一次步骤S120之后,再执行步骤S130:根据与当前实时采集到的波形信号对应的波峰以及波谷,得到当前的检测参数,然后根据当前的检测参数,更新当前的波形特征参数,并基于更新后的波形特征参数重复步骤S120。
如图5所示,本申请实施例还提供一种波形检测装置400,波形检测装置400可以包括:第一检测模块410以及第二检测模块420。
第一检测模块410,用于对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;
第二检测模块420,用于根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述采集到的实时波形信号对应的波峰以及波谷。
在一种可能的实施方式中,所述装置还包括更新模块,用于根据与所述实时采集到的波形信号对应的波峰以及波谷,得到当前的检测参数;根据当前的所述检测参数,更新所述波形特征参数。
在一种可能的实施方式中,所述第一检测模块410,用于通过滤波器过滤所述缓存第一时间段后的波形信号的噪声,得到处理后的缓存波形信号;对所述处理后的缓存波形信号的信号值进行识别,得到与所述缓存第一时间段后的波形信号对应的信号极值序列、信号最大值、信号最小值、最大信号振幅、波峰以及波谷;根据确定出的所有所述波峰以及所有所述波谷,确定出所述波形特征参数所包括的信号周期;其中,与所述缓存第一时间段后的波形信号对应的信号极值序列中,第一个满足第一波峰条件的极大值为第一个波峰,其余满足第二波峰条件的极大值为其他波峰,第一个满足第一波谷条件的极小值为第一个波谷,其余满足第二波谷条件的极小值为其他波谷;所述第一波峰条件为:所述极大值与所述信号最大值的比值大于预设阈值;所述第二波峰条件为:所述极大值与前一个已确定出的相邻波谷的幅值之差与所述最大信号振幅的比值大于所述预设阈值;所述第一波谷条件为:所述极小值与所述信号最小值的比值大于预设阈值;所述第二波谷条件为:所述极小值与前一个已确定出的相邻波峰的幅值之差与所述最大信号振幅的比值大于所述预设阈值。
在一种可能的实施方式中,所述第一检测模块410,用于通过高阶均值滤波器去除所述缓存第一时间段后的波形信号的基线漂移,得到校准缓存波形信号;通过低阶均值滤波器对所述校准缓存波形信号进行滤波,得到所述处理后的缓存波形信号;其中,所述高阶均值滤波器的阶数等于所述基准参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
在一种可能的实施方式中,所述装置还包括确定模块,用于确定本次检测的波形信号对象,本次检测的所述波形信号对象的产生时间起点为前一次检测得到的波峰或波谷所在的时间点,其产生时间终点为距离所述当前时刻之后第二时间段时所对应的时间点;
相应的,所述第二检测模块420,用于对实时采集到的波形信号缓存所述第二时间段后,得到本次检测的所述波形信号对象,且对其进行检测;其中,所述第一时间段对应的时长大于所述第二时间段对应的时长。
在一种可能的实施方式中,所述第二检测模块420,用于通过滤波器过滤所述实时采集到的波形信号的噪声,得到处理后的实时波形信号;对所述处理后的实时波形信号的信号值进行识别,得到与所述实时采集到的波形信号对应的信号极值序列;根据所述波形特征参数所包括的信号周期以及前一次检测所得到的最后一组波峰波谷所出现的时间坐标,预测得到所述实时采集到的波形信号中的下一组波峰波谷所出现的时间坐标范围;获取与所述实时采集到的波形信号对应的信号极值序列中其时间坐标属于所述时间坐标范围内的极值子序列;将所述极值子序列所包括的极大值确定为所述波峰,所述极值子序列所包括的极小值确定为所述波谷。
在一种可能的实施方式中,所述第二检测模块420,用于在确定所述极值子序列内包括多个极大值时,计算所述多个极大值中的每个极大值的波峰综合评分,并将所述波峰综合评分最小时所对应的极大值确定为所述波峰;
相应的,还用于在确定所述极值子序列内包括多个极小值时,计算所述多个极小值中的每个极小值的波谷综合评分,并将所述波谷综合评分最小时所对应的极小值确定为所述波谷;
其中,所述波峰综合评分用于表征当前的计算对象作为波峰时所得到的检测参数与所述波形特征参数之间的差异性;所述波谷综合评分用于表征当前的计算对象作为波谷时所得到的检测参数与所述波形特征参数之间的差异性。
在一种可能的实施方式中,所述第二检测模块420,用于通过高阶均值滤波器去除所述实时波形信号的基线漂移,得到校准实时波形信号;通过低阶均值滤波器对所述校准实时波形信号进行滤波,得到所述处理后的实时波形信号;其中,所述高阶均值滤波器的阶数等于所述波形特征参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
在一种可能的实施方式中,所述检测对象为人体,波形信号为光电容积脉搏波信号。
本申请实施例所提供的波形检测装置400,其实现原理及产生的技术效果和前述方法实施例相同,为简要描述,装置实施例部分未提及之处,可参考前述方法实施例中相应内容。
此外,本申请实施例还提供一种计算机可读存储介质,该计算机可读存储介质上存储有计算机程序,该计算机程序被计算机运行时,执行如上述的波形检测方法所包含的步骤。
此外,请参照图6,本申请实施例还提供一种用于实现本申请实施例的波形检测方法、装置的可穿戴设备100。
可选的,可穿戴设备100,可以是,但不限于手环、手表、眼镜、智能服饰等。
其中,可穿戴设备100可以包括:处理器110、存储器120、传感器130。
应当注意,图6所示的可穿戴设备100的组件和结构只是示例性的,而非限制性的,根据需要,可穿戴设备100也可以具有其他组件和结构。
例如,在一些情况下,可穿戴设备100还可以包括通信模块、显示模块等。
其中,通信模块用于使得可穿戴设备100可以与其他设备之间进行数据交互,例如可以是蓝牙模块、WIFI模块、Zigbee模块等。
显示模块可以用于显示检测结果。
处理器110、存储器120、传感器130以及其他可能出现于可穿戴设备100的组件相互之间直接或间接地电性连接,以实现数据的传输或交互。例如,处理器110、存储器120、传感器130以及其他可能出现的组件相互之间可通过一条或多条通讯总线或信号线实现电性连接。
传感器130用于采集波形信号。
存储器120用于存储程序,例如存储有前文出现的波形检测方法对应的程序或者前文出现的波形检测装置。可选的,当存储器120内存储有波形检测装置时,波形检测装置包括至少一个可以以软件或固件(firmware)的形式存储于存储器120中的软件功能模块。
可选的,波形检测装置所包括软件功能模块也可以固化在可穿戴设备100的操作系统(operating system,OS)中。
处理器110用于执行存储器120中存储的可执行模块,例如波形检测装置包括的软件功能模块或计算机程序。当处理器110在接收到执行指令后,可以执行计算机程序,例如执行:对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
当然,本申请任一实施例所揭示的方法都可以应用于处理器110中,或者由处理器110实现。
此外,请参照图7,在一些实施方式中,可穿戴设备可以包括:处理器110、存储器120、传感器130,通信模块140以及显示器150。
其中,传感器130,用于采集波形信号;
处理器110,用于根据预先保存的基准参数,对缓存第一时间段后的波形信号进行检测,得到波形特征参数;还用于根据波形特征参数,对所述传感器所实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
此外,在一些实施方式中,所述处理器110,还用于:根据所述波形特征参数,对所述传感器所实时采集到的且缓存第二时间段后的波形信号进行检测,得到与所述实时采集到的且缓存第二时间段后的波形信号对应的波峰以及波谷;所述第一时间段大于所述第二时间段。
具体的,处理器110在对缓存第一时间段后的波形信号以及实时采集到的波形信号进行检测时,可以借助于滤波器来进行滤波处理,而相应的滤波算法则存储于存储器120内。
至于上述过程的具体执行过程,请参照前述实施例部分,此处不再赘述。
综上所述,本发明实施例提出的波形检测方法、装置、可穿戴设备及计算机可读存储介质,该方法包括:S1、对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;S2、根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。在本实施例中,当需要对检测对象所产生实时波形信号进行检测时,为了提高检测结果的准确性,可以先将采集到的波形信号缓存一定时间(第一时间段)后,根据预设的基准参数对其进行检测,得到检测对象的波形特征参数(即与步骤S110对应的盲估阶段),然后再根据该波形特征参数对实时采集到的波形信号进行近实时的检测(即与步骤S120对应的近实时检测阶段),而不是如现有技术所示的根据大数据所得到的固定参数对检测对象所产生的实时波形信号进行检测,那么相较于现有技术,本申请实施例至少可以避免个体差异对检测结果造成的误差影响,进而可以提高检测结果的准确性。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和框图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或两个以上模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,笔记本电脑,服务器,或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。
Claims (14)
1.一种波形检测方法,其特征在于,所述方法包括:
S1、对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;
S2、根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
2.根据权利要求1所述的方法,其特征在于,在S2根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷之后,所述方法还包括:
S3、根据与所述实时采集到的波形信号对应的波峰以及波谷,得到当前的检测参数;根据当前的所述检测参数,更新所述波形特征参数,并重复S2。
3.根据权利要求1所述的方法,其特征在于,所述S1对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数,包括:
通过滤波器过滤所述缓存第一时间段后的波形信号的噪声,得到处理后的缓存波形信号;
对所述处理后的缓存波形信号的信号值进行识别,得到与所述缓存第一时间段后的波形信号对应的信号极值序列、信号最大值、信号最小值、最大信号振幅、波峰以及波谷;
根据确定出的所有所述波峰以及所有所述波谷,确定出所述波形特征参数所包括的信号周期;
其中,与所述缓存第一时间段后的波形信号对应的信号极值序列中,第一个满足第一波峰条件的极大值为第一个波峰,其余满足第二波峰条件的极大值为其他波峰,第一个满足第一波谷条件的极小值为第一个波谷,其余满足第二波谷条件的极小值为其他波谷;
所述第一波峰条件为:所述极大值与所述信号最大值的比值大于预设阈值;
所述第二波峰条件为:所述极大值与前一个已确定出的相邻波谷的幅值之差与所述最大信号振幅的比值大于所述预设阈值;
所述第一波谷条件为:所述极小值与所述信号最小值的比值大于预设阈值;
所述第二波谷条件为:所述极小值与前一个已确定出的相邻波峰的幅值之差与所述最大信号振幅的比值大于所述预设阈值。
4.根据权利要求3所述的方法,其特征在于,所述通过滤波器过滤所述缓存第一时间段后的波形信号的噪声,得到处理后的缓存波形信号,包括:
通过高阶均值滤波器去除所述缓存第一时间段后的波形信号的基线漂移,得到校准缓存波形信号;
通过低阶均值滤波器对所述校准缓存波形信号进行滤波,得到所述处理后的缓存波形信号;
其中,所述高阶均值滤波器的阶数等于所述基准参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
5.根据权利要求1-4中任一项所述的方法,其特征在于,在所述对实时采集到的波形信号进行检测之前,所述方法还包括:
确定本次检测的波形信号对象,本次检测的所述波形信号对象的产生时间起点为前一次检测得到的波峰或波谷所在的时间点,其产生时间终点为距离所述当前时刻之后第二时间段时所对应的时间点;
相应的,所述对实时采集到的波形信号进行检测,包括:
对实时采集到的波形信号缓存所述第二时间段后,得到本次检测的所述波形信号对象,且对其进行检测;
其中,所述第一时间段对应的时长大于所述第二时间段对应的时长。
6.根据权利要求1所述的方法,其特征在于,所述对实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷,包括:
通过滤波器过滤所述实时采集到的波形信号的噪声,得到处理后的实时波形信号;
对所述处理后的实时波形信号的信号值进行识别,得到与所述实时采集到的波形信号对应的信号极值序列;
根据所述波形特征参数所包括的信号周期以及前一次检测所得到的最后一组波峰波谷所出现的时间坐标,预测得到所述实时采集到的波形信号中的下一组波峰波谷所出现的时间坐标范围;
获取与所述实时采集到的波形信号对应的信号极值序列中其时间坐标属于所述时间坐标范围内的极值子序列;
将所述极值子序列所包括的极大值确定为所述波峰,所述极值子序列所包括的极小值确定为所述波谷。
7.根据权利要求6所述的方法,其特征在于,所述将所述极值子序列所包括的极大值确定为所述波峰,包括:
在确定所述极值子序列内包括多个极大值时,计算所述多个极大值中的每个极大值的波峰综合评分,并将所述波峰综合评分最小时所对应的极大值确定为所述波峰;
相应的,所述将极值子序列所包括的极小值确定为所述波谷,包括:
在确定所述极值子序列内包括多个极小值时,计算所述多个极小值中的每个极小值的波谷综合评分,并将所述波谷综合评分最小时所对应的极小值确定为所述波谷;
其中,所述波峰综合评分用于表征当前的计算对象作为波峰时所得到的检测参数与所述波形特征参数之间的差异性;所述波谷综合评分用于表征当前的计算对象作为波谷时所得到的检测参数与所述波形特征参数之间的差异性。
8.根据权利要求6所述的方法,其特征在于,所述通过滤波器过滤所述实时采集到的波形信号的噪声,包括:
通过高阶均值滤波器去除所述实时采集到的波形信号的基线漂移,得到校准实时波形信号;
通过低阶均值滤波器对所述校准实时波形信号进行滤波,得到所述处理后的实时波形信号;
其中,所述高阶均值滤波器的阶数等于所述波形特征参数所包括的周期;所述低阶均值滤波器的阶数等于预设值。
9.根据权利要求1所述的方法,其特征在于,检测对象为人体,所述波形信号为光电容积脉搏波信号。
10.一种波形检测装置,其特征在于,所述装置包括:
第一检测模块,用于对采集到的波形信号缓存第一时间段后,根据预设的基准参数对其进行检测,得到波形特征参数;
第二检测模块,用于根据所述波形特征参数,对实时采集到的波形信号进行检测,得到与所述采集到的实时波形信号对应的波峰以及波谷。
11.一种可穿戴设备,其特征在于,包括:传感器、存储器和处理器,所述传感器以及所述存储器和所述处理器连接;
所述传感器,用于采集波形信号;
所述处理器,用于根据预先保存的基准参数,对缓存第一时间段后的波形信号进行检测,得到波形特征参数;
所述处理器,还用于根据所述波形特征参数,对所述传感器所实时采集到的波形信号进行检测,得到与所述实时采集到的波形信号对应的波峰以及波谷。
12.根据权利要求11所述的可穿戴设备,其特征在于,所述处理器,还用于:根据所述波形特征参数,对所述传感器所实时采集到的且缓存第二时间段后的波形信号进行检测,得到与所述实时采集到的且缓存第二时间段后的波形信号对应的波峰以及波谷;
所述第一时间段大于所述第二时间段。
13.根据权利要求11所述的可穿戴设备,其特征在于,所述可穿戴设备为手环。
14.一种计算机可读存储介质,其特征在于,其上存储有计算机程序,所述计算机程序被计算机运行时执行如权利要求1-9中任一项所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110827269.7A CN113545753B (zh) | 2021-07-21 | 2021-07-21 | 波形检测方法、装置、可穿戴设备及计算机可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110827269.7A CN113545753B (zh) | 2021-07-21 | 2021-07-21 | 波形检测方法、装置、可穿戴设备及计算机可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113545753A true CN113545753A (zh) | 2021-10-26 |
CN113545753B CN113545753B (zh) | 2024-04-05 |
Family
ID=78132331
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110827269.7A Active CN113545753B (zh) | 2021-07-21 | 2021-07-21 | 波形检测方法、装置、可穿戴设备及计算机可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113545753B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114385012A (zh) * | 2022-01-17 | 2022-04-22 | 维沃移动通信有限公司 | 运动的识别方法、装置、电子设备和可读存储介质 |
CN114983372A (zh) * | 2022-05-30 | 2022-09-02 | 歌尔股份有限公司 | 一种可穿戴运动设备及其数据检测方法、装置及介质 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63212327A (ja) * | 1987-02-27 | 1988-09-05 | 松下電器産業株式会社 | 加速度脈波計 |
JP2004000316A (ja) * | 2002-04-26 | 2004-01-08 | U-Medica Inc | 脈波計測・解析装置及び脈波の平均波形算出・提示方法 |
JP2006055406A (ja) * | 2004-08-20 | 2006-03-02 | Denso Corp | 心拍検出装置 |
CN101991410A (zh) * | 2009-08-31 | 2011-03-30 | 深圳市理邦精密仪器股份有限公司 | 一种脉率搜索和计算方法 |
CN102247129A (zh) * | 2011-06-15 | 2011-11-23 | 西安电子科技大学 | 脉搏波非典型波峰波谷的识别方法 |
RU2500069C1 (ru) * | 2012-09-20 | 2013-11-27 | Общество с ограниченной ответственностью Научно-производственное предприятие "Новые Технологии Телекоммуникаций" | Способ генерирования кодов для формирования ансамблей сигналов в телекоммуникационных сетях |
CN104706337A (zh) * | 2015-02-11 | 2015-06-17 | 华东师范大学 | 脉搏波波峰波谷自动检测方法 |
CN105816163A (zh) * | 2016-05-09 | 2016-08-03 | 安徽华米信息科技有限公司 | 检测心率的方法、装置及可穿戴设备 |
CN109864705A (zh) * | 2019-01-07 | 2019-06-11 | 平安科技(深圳)有限公司 | 对脉搏波进行滤波的方法、装置和计算机设备 |
CN110825826A (zh) * | 2019-11-07 | 2020-02-21 | 深圳大学 | 聚类簇计算方法、装置、终端及存储介质 |
CN111665900A (zh) * | 2020-05-26 | 2020-09-15 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于填充函数算法的光伏阵列多峰最大功率跟踪方法 |
CN111657905A (zh) * | 2020-06-23 | 2020-09-15 | 中国医学科学院生物医学工程研究所 | 一种特征点检测方法、装置、设备及存储介质 |
CN112888062A (zh) * | 2021-03-16 | 2021-06-01 | 芯原微电子(成都)有限公司 | 数据同步方法、装置、电子设备及计算机可读存储介质 |
-
2021
- 2021-07-21 CN CN202110827269.7A patent/CN113545753B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63212327A (ja) * | 1987-02-27 | 1988-09-05 | 松下電器産業株式会社 | 加速度脈波計 |
JP2004000316A (ja) * | 2002-04-26 | 2004-01-08 | U-Medica Inc | 脈波計測・解析装置及び脈波の平均波形算出・提示方法 |
JP2006055406A (ja) * | 2004-08-20 | 2006-03-02 | Denso Corp | 心拍検出装置 |
CN101991410A (zh) * | 2009-08-31 | 2011-03-30 | 深圳市理邦精密仪器股份有限公司 | 一种脉率搜索和计算方法 |
CN102247129A (zh) * | 2011-06-15 | 2011-11-23 | 西安电子科技大学 | 脉搏波非典型波峰波谷的识别方法 |
RU2500069C1 (ru) * | 2012-09-20 | 2013-11-27 | Общество с ограниченной ответственностью Научно-производственное предприятие "Новые Технологии Телекоммуникаций" | Способ генерирования кодов для формирования ансамблей сигналов в телекоммуникационных сетях |
CN104706337A (zh) * | 2015-02-11 | 2015-06-17 | 华东师范大学 | 脉搏波波峰波谷自动检测方法 |
CN105816163A (zh) * | 2016-05-09 | 2016-08-03 | 安徽华米信息科技有限公司 | 检测心率的方法、装置及可穿戴设备 |
CN109864705A (zh) * | 2019-01-07 | 2019-06-11 | 平安科技(深圳)有限公司 | 对脉搏波进行滤波的方法、装置和计算机设备 |
CN110825826A (zh) * | 2019-11-07 | 2020-02-21 | 深圳大学 | 聚类簇计算方法、装置、终端及存储介质 |
CN111665900A (zh) * | 2020-05-26 | 2020-09-15 | 国网江苏省电力有限公司盐城供电分公司 | 一种基于填充函数算法的光伏阵列多峰最大功率跟踪方法 |
CN111657905A (zh) * | 2020-06-23 | 2020-09-15 | 中国医学科学院生物医学工程研究所 | 一种特征点检测方法、装置、设备及存储介质 |
CN112888062A (zh) * | 2021-03-16 | 2021-06-01 | 芯原微电子(成都)有限公司 | 数据同步方法、装置、电子设备及计算机可读存储介质 |
Non-Patent Citations (9)
Title |
---|
LOVRIC, M (LOVRIC, MILIVOJ) [1] ; KOMORSKY-LOVRIC, S (KOMORSKY-LOVRIC, SEBOJKA): "Theory of reverse scan square-wave voltammetry influenced by the kinetics of reactant adsorption", 《CENTRAL EUROPEAN JOURNAL OF CHEMISTRY 》, vol. 8, no. 3, 30 June 2010 (2010-06-30), pages 513 - 518 * |
WU, FX (WU, FEIXIANG) [1] ; XU, XW (XU, XUEWU) [2] ; MIAO, XR (MIAO, XUERONG) [1] ; CHEN, JY (CHEN, JIAYING) [1] ; SUN, YM (SUN, Y: "Effect of recombinant adenovirus coding for endomorphin-2 on neuropathic pain in rats", 《INTERNATIONAL JOURNAL OF CLINICAL AND EXPERIMENTAL PATHOLOGY 》, vol. 5, no. 9, 31 December 2012 (2012-12-31), pages 914 - 923 * |
刘圳鑫: "基于SIFT特征的适配区选取方法研究", 《CNKI》, 15 March 2020 (2020-03-15) * |
周莹, 曹美媛, 谷红霞, 等: "一种基于微分阈值法的脉搏波特征点提取优化算法", 《工业控制计算机》, vol. 32, no. 8, 31 December 2019 (2019-12-31), pages 59 - 61 * |
唐铭一;李凯;马小铁;: "脉搏波信号时域特征提取与算法的研究", 计算机与现代化, no. 04, 15 April 2010 (2010-04-15) * |
张翠娜: "心电信号去噪算法的研究", 万方, 3 September 2012 (2012-09-03) * |
胡中华, 虞卫峰, 余利锋: "基于精细化管理的县(区)级国土资源数据中心建设的实践与探索", 《江苏省土地学会2015年度学术年会 论文集》, 31 December 2015 (2015-12-31), pages 609 - 617 * |
虞卫峰, 唐慧明: "一种用CPLD实现视频信号运动检测的方法", 《电子技术应用》, no. 11, 31 December 2003 (2003-12-31), pages 64 - 67 * |
阿布都瓦斯提·吾拉木, 李召良, 秦其明等.: "全覆盖植被冠层水分遥感监测的一种方法:短波红外垂直失水指数", 《中国科学(D辑:地球科学)》, no. 07, 31 December 2007 (2007-12-31), pages 957 - 965 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114385012A (zh) * | 2022-01-17 | 2022-04-22 | 维沃移动通信有限公司 | 运动的识别方法、装置、电子设备和可读存储介质 |
CN114385012B (zh) * | 2022-01-17 | 2023-06-30 | 维沃移动通信有限公司 | 运动的识别方法、装置、电子设备和可读存储介质 |
CN114983372A (zh) * | 2022-05-30 | 2022-09-02 | 歌尔股份有限公司 | 一种可穿戴运动设备及其数据检测方法、装置及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113545753B (zh) | 2024-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Van Gent et al. | Analysing noisy driver physiology real-time using off-the-shelf sensors: Heart rate analysis software from the taking the fast lane project | |
JP3787336B2 (ja) | Ppgを用いた心拍検出装置及びその方法 | |
CN113545753B (zh) | 波形检测方法、装置、可穿戴设备及计算机可读存储介质 | |
US20220386929A1 (en) | Detection of atrial fibrillation | |
Sparacino et al. | Continuous glucose monitoring time series and hypo/hyperglycemia prevention: requirements, methods, open problems | |
US11129538B2 (en) | Heart rate estimation apparatus with state sequence optimization | |
WO2019036908A1 (zh) | 心率检测方法及装置、电子终端 | |
GB2488316A (en) | Method for determining respiration rate from uncorrupted signal segments | |
KR20220123376A (ko) | 심혈관 파라미터를 결정하기 위한 방법 및 시스템 | |
CN117530666B (zh) | 呼吸异常识别模型训练方法、呼吸异常识别方法及设备 | |
Shintomi et al. | Effectiveness of the heartbeat interval error and compensation method on heart rate variability analysis | |
Mesin | Heartbeat monitoring from adaptively down-sampled electrocardiogram | |
Lin et al. | A characteristic filtering method for pulse wave signal quality assessment | |
CN118201542A (zh) | 用于判断人是否睡着的装置、系统和方法 | |
Fioravanti et al. | Machine learning framework for Inter-Beat Interval estimation using wearable Photoplethysmography sensors | |
US20230057740A1 (en) | Sensing | |
CN113940637B (zh) | 脉搏波信号质量评估方法、装置及存储介质 | |
CN113693578B (zh) | 一种心率估计方法、装置、设备、系统及存储介质 | |
Shcherbakova et al. | Determination of characteristic points of electrocardiograms using multi-start optimization with a wavelet transform | |
JP2010213809A (ja) | 生体信号分析装置 | |
CN114098690A (zh) | 一种轻量级脉搏波信号噪声去除及心率检测方法 | |
CN114469040A (zh) | 心率检测方法及装置、存储介质及电子设备 | |
JP7252542B2 (ja) | 覚醒状態推定装置及び覚醒状態推定方法 | |
Madhan Mohan et al. | Spot measurement of heart rate based on morphology of PhotoPlethysmoGraphic (PPG) signals | |
Rodriguez et al. | Embedded system based on an ARM microcontroller to analyze heart rate variability in real time using wavelets |
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 |