CN108992054A - 一种脉搏信号峰值点检测方法及装置 - Google Patents
一种脉搏信号峰值点检测方法及装置 Download PDFInfo
- Publication number
- CN108992054A CN108992054A CN201810677236.7A CN201810677236A CN108992054A CN 108992054 A CN108992054 A CN 108992054A CN 201810677236 A CN201810677236 A CN 201810677236A CN 108992054 A CN108992054 A CN 108992054A
- Authority
- CN
- China
- Prior art keywords
- point
- curve
- peak
- amplitude versus
- frequency
- 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 30
- 238000001914 filtration Methods 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 31
- 230000009466 transformation Effects 0.000 claims abstract description 16
- 230000005622 photoelectricity Effects 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims description 18
- 238000009795 derivation Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims 1
- 239000000284 extract Substances 0.000 abstract description 2
- 210000000709 aorta Anatomy 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 4
- 238000010168 coupling process Methods 0.000 description 4
- 238000005859 coupling reaction Methods 0.000 description 4
- 230000010349 pulsation Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000002861 ventricular Effects 0.000 description 2
- 230000006399 behavior Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 210000003017 ductus arteriosus Anatomy 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000031700 light absorption Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000013186 photoplethysmography Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
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/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02416—Detecting, measuring or recording pulse rate or heart rate using photoplethysmograph signals, e.g. generated by infrared radiation
-
- 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/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- 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/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
Abstract
本申请公开了一种脉搏信号峰值点检测方法及装置。该方法包括：对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度‑频率曲线；从所述幅度‑频率曲线中选出幅度最大的点对应的频率f；对所述PPG原始信号进行低通滤波处理，得到峰值曲线；确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。还公开了相应的装置。通过傅里叶变换提取出PPG原始信号的频率值f，再以f为阈值对PPG原始信号进行低通滤波，并求出滤波后的信号的极大值点，最后对已得的频率f对极大值点中的错误峰值点进行剔除，即可精确的提取出PPG信号中的峰值点。
Description
技术领域
本申请涉及信息处理技术领域，尤其涉及一种脉搏信号峰值点检测方法及装置。
背景技术
光电容积脉搏波描记法(photoplethysmography，PPG)测量人体生理参数具有操作简单、性能稳定和适应性强等优点。随着光电技术在医疗应用中的扩大，人们对利用光电技术测量人体生理参数的研究也逐渐加强。在利用PPG信号进行心率变异性分析时，需要对采集的PPG信号进行峰值点提取。现有技术一般采用对信号进行滤波处理、差分处理及傅里叶变换来提取PPG信号中的峰值点，但这种方式提取的峰值点精确度较低，且容易提取出错误的峰值点。
发明内容
本申请提供了一种脉搏信号峰值点检测方法及装置，以实现对脉搏信号峰值点的精确检测。
第一方面，提供了一种脉搏信号峰值点检测方法，包括：对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线；从所述幅度-频率曲线中选出幅度最大的点对应的频率f；对所述PPG原始信号进行低通滤波处理，得到峰值曲线；确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。在一种可能的实现方式中，所述从所述幅度-频率曲线中选出幅度最大的点对应的频率f，包括：从所述幅度-频率曲线中选取预设频率之间的幅度-频率曲线，得到第一幅度-频率曲线；对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线；将所述第一幅度-频率曲线与所述第二幅度-频率曲线进行相减，并取绝对值，得到第三幅度-频率曲线；从所述第三幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到所述幅度值最大的点对应的频率值f。
在该实现方式中，所述预设频率为0.667HZ～2.667HZ。在另一种可能的实现方式中，所述对所述PPG原始信号进行低通滤波处理，得到峰值曲线，包括：将所述频率值f设为低通滤波的阈值；将所述PPG原始信号的频率值与f相比较，若所述PPG信号中的频率值小于或等于f，则保留所述频率值对应的信号，得到峰值曲线，若所述PPG信号中的频率值大于f，则去除所述频率值对应的信号。
在又一种可能的实现方式中，所述确定所述峰值曲线中的所有极大值点，得到峰值点，包括：根据所述频率值f，求得时间周期T，其中，T大于0；对所述第峰值曲线求导，找出所有极大值点；求出所述任意一个极大值点与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；若所述时间差t1和t2均大于T/1.5，则将所述任意一个极值点划分为峰值点，若所述时间差t1和t2均小于或等于T/1.5，则不将所述任意一个极值点划分为峰值点。
第二方面，提供了一种脉搏信号峰值点检测装置，包括：处理单元，用于对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线；选取单元，用于从所述幅度-频率曲线中选出幅度最大的点对应的频率f；滤波单元，用于对所述PPG原始信号进行低通滤波处理，得到峰值曲线；确定单元，用于确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
在一种可能的实现方式中，所述选取单元包括：第一选取子单元，用于从所述幅度-频率曲线中选取预设频率之间的幅度-频率曲线，得到第一幅度-频率曲线；第二选取子单元，用于对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线；第一计算子单元，用于将所述第一幅度-频率曲线与所述第二幅度-频率曲线进行相减，并取绝对值，得到第三幅度-频率曲线；第三选取子单元，用于从所述第三幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到所述幅度值最大的点对应的频率值f。
在该实现方式中，所述预设频率为0.667HZ～2.667HZ。
在另一种可能的实现方式中，所述滤波单元包括：设定子单元，用于将所述频率值f设为低通滤波的阈值；比较子单元，用于将所述PPG原始信号的频率值与f相比较，若所述PPG信号中的频率值小于或等于f，则保留所述频率值对应的信号，得到峰值曲线，若所述PPG信号中的频率值大于f，则去除所述频率值对应的信号。
在又一种可能的实现方式中，所述确定单元包括：第二计算子单元，用于根据所述频率值f，求得时间周期T，其中，T大于0；第三计算子单元，用于对所述第峰值曲线求导，找出所有极大值点；第四计算子单元，用于求出所述任意一个极大值点与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；划分子单元，用于若所述时间差t1和t2均大于T/1.5，则将所述任意一个极值点划分为峰值点，若所述时间差t1和t2均小于或等于T/1.5，则不将所述任意一个极值点划分为峰值点。
第三方面，提供了一种脉搏信号峰值点检测装置，包括：包括处理器、存储器；所述处理器被配置为支持所述装置执行上述第一方面及其任一种可能的实现方式的方法中相应的功能。存储器用于与处理器耦合，其保存所述装置必要的程序(指令)和数据。可选的，所述装置还可以包括输入/输出接口，用于支持所述装置与其他装置之间的通信。
第四方面，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有指令，当其在计算机上运行时，使得计算机执行上述第一方面及其任一种可能的实现方式的方法。
第五方面，提供了一种包含指令的计算机程序产品，当其在计算机上运行时，使得计算机执行上述第一方面及其任一种可能的实现方式的方法。
根据本申请提供的一种脉搏信号峰值点检测方法，通过傅里叶变换提取出PPG原始信号的频率值f，再以频率值f为阈值对PPG原始信号进行低通滤波，并求出滤波后的信号的极大值点，最后，利用已得的频率f对极大值点中的错误峰值点进行剔除，即可精确的提取出PPG信号中的峰值点。
附图说明
为了更清楚地说明本申请实施例或背景技术中的技术方案，下面将对本申请实施例或背景技术中所需要使用的附图进行说明。
图1为本申请实施例提供的一种脉搏信号峰值点检测方法的流程示意图；
图2为本申请实施例提供的另一种脉搏信号峰值点检测方法的流程示意图；
图3为本申请实施例提供的一种脉搏信号峰值点检测装置的结构示意图；
图4为本申请实施例提供的一种脉搏信号峰值点检测装置的硬件结构示意图。
具体实施方式
下面结合本申请实施例中的附图对本申请实施例进行描述。
请参阅图1，图1是本申请实施例提供的一种脉搏信号峰值点检测方法的流程示意图。
S101、对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线。
首先，通过心率传感器采集得到PPG原始信号，再对PPG原始信号进行傅里叶变换，将时间-频率曲线转换为幅度-频率曲线，得到PPG信号的幅度-频率曲线，以方便后续的处理。
S102、从所述幅度-频率曲线中选出幅度最大的点对应的频率f。
根据正常人的心率范围对幅度-频率曲线进行筛选，具体地，从PPG原始信号的幅度-频率曲线中选取心率范围内的信号，并令该信号的功率为P。再对选取出的信号进行滤波处理，得到PPG信号，令PPG信号的功率为Ps。将滤波处理前和滤波处理后的两条曲线的幅值相减，并取绝对值，得到新的幅度-频率曲线。从新的幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到该幅度值最大的点对应的频率值f。
S103、对所述PPG原始信号进行低通滤波处理，得到峰值曲线。
将频率值f设为低通滤波的阈值，并对PPG原始信号进行低通滤波，得到峰值曲线。具体地，PPG原始信号中频率大于或等于f的信号将被去除，频率小于f的信号将被保留，并得到峰值曲线。
S104、确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
通过对峰值曲线中的任意一点进行求导，得到峰值曲线中所有的极大值点。由于所求得的极大值点有部分极大值点并非峰值点，即极大值点中有误检测峰值点，因此，需要将这些误检测峰值点排除。
首先，根据上述所求得的频率f，求得与频率对应的周期T，即T＝1/f。然后，求出所述任意一个极大值点与相邻两个极大值点的时间差，并令两个时间差分别为t1和t2，其中，t1和t2均大于0。最后，结合t1、t2及T对所有极大值点进行判断，具体地，若某一极大值点的时间差t1和t2均大于T/1.5，则将该极大值点划分为峰值点，若某一极大值点的时间差t1和t2均小于或等于T/1.5，则不将该极值点划分为峰值点。
这样，就可将非峰值点的极大值点排除，得到PPG原始信号中的峰值点。
根据本申请提供的一种脉搏信号峰值点检测方法，通过傅里叶变换提取出PPG原始信号的频率值f，再以频率值f为阈值对PPG原始信号进行低通滤波，并求出滤波后的信号的极大值点，最后，利用已得的频率f对极大值点中的错误峰值点进行剔除，即可精确的提取出PPG信号中的峰值点。
请参阅图2，图2是本申请实施例提供的一种脉搏信号峰值点检测方法的流程示意图。
S201、通过心率传感器采集PPG原始信号。
利用心率传感器采集PPG原始信号。具体地，当心室收缩时，血液冲开主动脉瓣，并把血液射入主动脉中，主动脉内压突然增高，血管壁迅速膨大，当心室舒张时，主动脉压降低，主动脉壁因其具有弹性而回缩，这样，动脉管壁就随心室的收缩而出现周期性的起伏搏动，形成脉搏。由于皮肤、肌肉组织等对光的吸收在整个血液循环中是保持恒定不变的，而皮肤内的血液、容积在心脏作用下呈搏动性变化。一定波长的光束照射到指端皮肤表面时，光束将通过透射或反射方式传送到光电接收器，当心脏收缩时外周血容量最多光吸收量也最大，检测到的光强度最小，而在心脏舒张时，正好相反，检测到的光强度最大。将光源和光敏元件置于被测部位(如手指)的同一侧(或两侧)，光源发出的光照射在组织上，经反射(或透射)后被光敏元件接收，光敏元件再将脉动的光强度信号转变为脉动的电信号，从而可得到PPG原始信号。
S202、对所述PPG原始信号进行傅里叶变换，得到PPG信号的幅度-频率曲线。
通过对PPG原始信号进行傅里叶变换，将时间-频率曲线转换为幅度-频率曲线，得到PPG信号的幅度-频率曲线，以方便后续的处理。
S203、从所述幅度-频率曲线中选取频率为0.667HZ～2.667HZ之间的幅度-频率曲线，得到第一幅度-频率曲线。
通过上述脉搏波信号可得到相应的心率，心率是指正常人安静状态下每分钟心跳的次数，正常人在安静状态下的心率一般为60～100次/分，但由于存在年龄、性别或其他生理因素的影响，心率也会产生个体差异。本实施例是对正常人的脉搏信号中的峰值点进行检测，因此，综合各种因素考虑(包括正常人的运动状态)，取心率在40～160次/分钟之间的信号为正常信号，再将40～160次/分钟的心率转换为频率，即为0.667HZ～2.667HZ。
S204、对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线。
由于上述得到的第一幅度-频率曲线中有其他干扰信号，而这一部分干扰信号将对后续的处理产生很大的干扰，同时是曲线看起来很不平滑，且带来很大的误差，因此，在进行后续处理之前需要先对第一幅度-频率曲线进行滤波处理。
本实施例采用的是5点平滑滤波，本质就是去掉信号中不必要的快速起伏的成份，使信号看起来更平滑。采用简单平均法进行，求邻近5点信号的平均值，用以替代该点的信号值，以达到对信号平滑处理的效果。具体可由下式表达：
其中，x(·)为滤波处理前的信号值，y(·)为滤波处理后的信号值，n≥2，且n为整数。
通过对第一幅度-频率曲线进行5点平滑滤波处理，即可得第二幅度-频率曲线。
S205、将所述第一幅度-频率曲线与所述第二幅度-频率曲线做差，并取绝对值，得到第三幅度-频率曲线。
由上述可知，第二幅度-频率曲线是由第一幅度-频率曲线滤波处理后得到的，因此，两条曲线的频率范围一致，但相同频率对应的幅度值不同，通过将第一幅度-频率曲线与第二幅度-频率曲线同一频率下的幅度相减，得到幅度差值，再对差值取绝对值，并将该绝对值与原频率一一对应，即可得到第三幅度-频率曲线。
通过对第一幅度-频率曲线与第二幅度-频率曲线做差，并取绝对值，可减少心率传感器采集到的错误信号。
S206、从所述第三幅度-频率曲线中选取幅度最大的点所对应的频率f。
从第三幅度-频率曲线中找出幅度最大的点，并得到与该点相对应的频率f，用以后续的处理。
S207、以所述频率f为阈值，对所述PPG原始信号进行低通滤波处理，得到峰值曲线。
以上述得到的频率f为阈值，对PPG原始信号进行低通滤波处理，具体地，将PPG原始信号中频率大于f的信号全部去除，保留频率小于或等于f的信号，得到峰值曲线。
通过对PPG原始信号进行低通滤波处理，可将PPG原始信号中的高频噪声去除。
S208、确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
该步骤可参考前述实施例中的S104，在此不再赘述。
根据本申请提供的一种脉搏信号峰值点检测方法，通过傅里叶变换提取出PPG原始信号的频率值f，再以频率值f为阈值对PPG原始信号进行低通滤波，并求出滤波后的信号的极大值点，最后，利用已得的频率f对极大值点中的错误峰值点进行剔除，即可精确的提取出PPG信号中的峰值点。
请参阅图3，图3为本申请实施例提供的一种脉搏信号峰值点检测装置的结构示意图，该装置1000包括：处理单元11、选取单元12、滤波单元13、确定单元14。其中：
处理单元11，用于对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线；
选取单元12，用于从所述幅度-频率曲线中选出幅度最大的点对应的频率f；
滤波单元13，用于对所述PPG原始信号进行低通滤波处理，得到峰值曲线；
确定单元14，用于确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
进一步地，所述选取单元12包括：
第一选取子单元121，用于从所述幅度-频率曲线中选取预设频率之间的幅度-频率曲线，得到第一幅度-频率曲线；
第二选取子单元122，用于对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线；
第一计算子单元123，用于将所述第一幅度-频率曲线与所述第二幅度-频率曲线进行相减，并取绝对值，得到第三幅度-频率曲线；
第三选取子单元124，用于从所述第三幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到所述幅度值最大的点对应的频率值f。
其中，所述预设频率为0.667HZ～2.667HZ。
进一步地，所述滤波单元13包括：
设定子单元131，用于将所述频率值f设为低通滤波的阈值；
比较子单元132，用于将所述PPG原始信号的频率值与f相比较，若所述PPG信号中的频率值小于或等于f，则保留所述频率值对应的信号，得到峰值曲线，若所述PPG信号中的频率值大于f，则去除所述频率值对应的信号。
进一步地，所述确定单元包括：
第二计算子单元，用于根据所述频率值f，求得时间周期T，其中，T大于0；
第三计算子单元，用于对所述第峰值曲线求导，找出所有极大值点；
第四计算子单元，用于求出所述任意一个极大值点与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；
划分子单元，用于若所述时间差t1和t2均大于T/1.5，则将所述任意一个极值点划分为峰值点，若所述时间差t1和t2均小于或等于T/1.5，则不将所述任意一个极值点划分为峰值点。
根据本申请提供的一种脉搏信号峰值点检测装置，通过对PPG原始信号进行低通滤波，并求出滤波后的信号的极大值点，最后，从极大值点中筛选得到峰值点，可实现对PPG信号峰值点的精确提取。
图4为本申请实施例提供的一种脉搏信号峰值点检测装置的硬件结构示意图。该通信装置包括处理器21，还可以包括输入装置22、输出装置23和存储器34。该输入装置22、输出装置23、存储器24和处理器21之间通过总线相互连接。
存储器包括但不限于是随机存储记忆体(random access memory，RAM)、只读存储器(read-only memory，ROM)、可擦除可编程只读存储器(erasable programmable readonly memory，EPROM)、或便携式只读存储器(compact disc read-only memory，CD-ROM)，该存储器用于相关指令及数据。
输入装置用于输入数据和/或信号，以及输出装置用于输出数据和/或信号。输出装置和输入装置可以是独立的器件，也可以是一个整体的器件。
处理器可以包括是一个或多个处理器，例如包括一个或多个中央处理器(centralprocessing unit，CPU)，在处理器是一个CPU的情况下，该CPU可以是单核CPU，也可以是多核CPU。
存储器用于存储网络设备的程序代码和数据。
处理器用于调用该存储器中的程序代码和数据，执行如下步骤：
在一个实现方式中，所述处理器用于执行如下步骤：对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线；从所述幅度-频率曲线中选出幅度最大的点对应的频率f；对所述PPG原始信号进行低通滤波处理，得到峰值曲线；确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
在另一个实现方式中，所述处理器用于执行如下步骤：从所述幅度-频率曲线中选取预设频率之间的幅度-频率曲线，得到第一幅度-频率曲线；对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线；将所述第一幅度-频率曲线与所述第二幅度-频率曲线进行相减，并取绝对值，得到第三幅度-频率曲线；从所述第三幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到所述幅度值最大的点对应的频率值f。
在又一个实现方式中，所述处理器用于执行如下步骤：将所述频率值f设为低通滤波的阈值；将所述PPG原始信号的频率值与f相比较，若所述PPG信号中的频率值小于或等于f，则保留所述频率值对应的信号，得到峰值曲线，若所述PPG信号中的频率值大于f，则去除所述频率值对应的信号。
在又一个实现方式中，所述处理器用于执行如下步骤：根据所述频率值f，求得时间周期T，其中，T大于0；对所述第峰值曲线求导，找出所有极大值点；求出所述任意一个极大值点与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；
点与相邻两个极大值点的时间差和，其中，和均大于0；若所述时间差t1和t2均大于T/1.5，则将所述任意一个极值点划分为峰值点，若所述时间差t1和t2均小于或等于T/1.5，则不将所述任意一个极值点划分为峰值点。
可以理解的是，图4仅仅示出了一种脉搏信号峰值点检测装置的简化设计。在实际应用中，脉搏信号峰值点检测装置还可以分别包含必要的其他元件，包含但不限于任意数量的输入/输出装置、处理器、控制器、存储器等，而所有可以实现本申请实施例的脉搏信号峰值点检测装置都在本申请的保护范围之内。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统、装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。例如，该单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。所显示或讨论的相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
在上述实施例中，可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时，可以全部或部分地以计算机程序产品的形式实现。该计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行该计算机程序指令时，全部或部分地产生按照本申请实施例的流程或功能。该计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。该计算机指令可以存储在计算机可读存储介质中，或者通过该计算机可读存储介质进行传输。该计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(digital subscriber line，DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。该计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。该可用介质可以是只读存储器(read-onlymemory，ROM)，或随机存储存储器(random access memory，RAM)，或磁性介质，例如，软盘、硬盘、磁带、磁碟、或光介质，例如，数字通用光盘(digital versatile disc，DVD)、或者半导体介质，例如，固态硬盘(solid state disk，SSD)等。
Claims (10)
1.一种脉搏信号峰值点检测方法，其特征在于，包括：
对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线；
从所述幅度-频率曲线中选出幅度最大的点对应的频率f；
对所述PPG原始信号进行低通滤波处理，得到峰值曲线；
确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
2.如权利要求1所述方法，其特征在于，所述从所述幅度-频率曲线中选出幅度最大的点对应的频率f，包括：
从所述幅度-频率曲线中选取预设频率之间的幅度-频率曲线，得到第一幅度-频率曲线；
对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线；
将所述第一幅度-频率曲线与所述第二幅度-频率曲线进行相减，并取绝对值，得到第三幅度-频率曲线；
从所述第三幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到所述幅度值最大的点对应的频率值f。
3.如权利要求2所述方法，其特征在于，所述预设频率为0.667HZ～2.667HZ。
4.如权利要求1所述方法，其特征在于，所述对所述PPG原始信号进行低通滤波处理，得到峰值曲线，包括：
将所述频率值f设为低通滤波的阈值；
将所述PPG原始信号的频率值与f相比较，若所述PPG信号中的频率值小于或等于f，则保留所述频率值对应的信号，得到峰值曲线，若所述PPG信号中的频率值大于f，则去除所述频率值对应的信号。
5.如权利要求1所述方法，其特征在于，所述确定所述峰值曲线中的所有极大值点，得到峰值点，包括：
根据所述频率值f，求得时间周期T，其中，T大于0；
对所述第峰值曲线求导，找出所有极大值点；
求出所述任意一个极大值点与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；
若所述时间差t1和t2均大于T/1.5，则将所述任意一个极值点划分为峰值点，若所述时间差t1和t2均小于或等于T/1.5，则不将所述任意一个极值点划分为峰值点。
6.一种脉搏信号峰值点检测装置，其特征在于，包括：
处理单元，用于对采集的光电容积脉搏描记法PPG原始信号进行傅里叶变换，得到所述PPG原始信号的幅度-频率曲线；
选取单元，用于从所述幅度-频率曲线中选出幅度最大的点对应的频率f；
滤波单元，用于对所述PPG原始信号进行低通滤波处理，得到峰值曲线；
确定单元，用于确定所述峰值曲线中的所有极大值点，得到峰值点，其中，所述峰值点为所述极大值点去除误检测峰值点后的所有极大值点。
7.如权利要求6所述装置，其特征在于，所述选取单元包括：
第一选取子单元，用于从所述幅度-频率曲线中选取预设频率之间的幅度-频率曲线，得到第一幅度-频率曲线；
第二选取子单元，用于对所述第一幅度-频率曲线进行平滑滤波处理，得到第二幅度-频率曲线；
第一计算子单元，用于将所述第一幅度-频率曲线与所述第二幅度-频率曲线进行相减，并取绝对值，得到第三幅度-频率曲线；
第三选取子单元，用于从所述第三幅度-频率曲线中选出除左、右两个端点之外的幅度值最大的点，得到所述幅度值最大的点对应的频率值f。
8.如权利要求6所述装置，其特征在于，所述滤波单元包括：
设定子单元，用于将所述频率值f设为低通滤波的阈值；
比较子单元，用于将所述PPG原始信号的频率值与f相比较，若所述PPG信号中的频率值小于或等于f，则保留所述频率值对应的信号，得到峰值曲线，若所述PPG信号中的频率值大于f，则去除所述频率值对应的信号。
9.如权利要求6所述装置，其特征在于，所述确定单元包括：
第二计算子单元，用于根据所述频率值f，求得时间周期T，其中，T大于0；
第三计算子单元，用于对所述第峰值曲线求导，找出所有极大值点；
第四计算子单元，用于求出所述任意一个极大值点与相邻两个极大值点的时间差t1和t2，其中，t1和t2均大于0；
划分子单元，用于若所述时间差t1和t2均大于T/1.5，则将所述任意一个极值点划分为峰值点，若所述时间差t1和t2均小于或等于T/1.5，则不将所述任意一个极值点划分为峰值点。
10.一种可穿戴设备穿戴状态的检测装置，其特征在于，包括：处理器、存储器；所述计算机可读存储介质中存储有指令，当其在计算机上运行时，使得计算机执行如权利要求1～5任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810677236.7A CN108992054B (zh) | 2018-06-27 | 2018-06-27 | 一种脉搏信号峰值点检测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810677236.7A CN108992054B (zh) | 2018-06-27 | 2018-06-27 | 一种脉搏信号峰值点检测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108992054A true CN108992054A (zh) | 2018-12-14 |
| CN108992054B CN108992054B (zh) | 2021-11-12 |
Family
ID=64600326
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810677236.7A Active CN108992054B (zh) | 2018-06-27 | 2018-06-27 | 一种脉搏信号峰值点检测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108992054B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113032724A (zh) * | 2019-12-09 | 2021-06-25 | 中国石油化工股份有限公司 | 一种示踪剂产出曲线降噪处理方法、存储介质及电子设备 |
| CN113317757A (zh) * | 2021-04-30 | 2021-08-31 | 深圳麦格米特电气股份有限公司 | 一种获取生命体征数据的方法、光纤传感器及设备 |
| CN114041768A (zh) * | 2021-10-11 | 2022-02-15 | 广东宝莱特医用科技股份有限公司 | 脉率提取方法、设备及存储介质 |
| CN114521880A (zh) * | 2022-01-21 | 2022-05-24 | 中国人民解放军陆军军医大学 | 一种运动状态下心率计算方法、系统及计算机存储介质 |
| CN115349816A (zh) * | 2022-10-21 | 2022-11-18 | 万灵帮桥医疗器械(广州)有限责任公司 | 干涉信号峰值检测方法、装置、设备和存储介质 |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101732050A (zh) * | 2009-12-04 | 2010-06-16 | 西安交通大学 | 一种基于光电容积波的呼吸率监测方法 |
| CN103027667A (zh) * | 2011-09-30 | 2013-04-10 | Ge医疗系统环球技术有限公司 | 脉搏波的特征参数提取 |
| CN103153175A (zh) * | 2010-05-14 | 2013-06-12 | 先进计算发展中心 | 脉搏信号波形数据的诊断分类 |
| CN104349713A (zh) * | 2012-04-23 | 2015-02-11 | 图尔库大学 | 确定指示心脏功能障碍和异常的信息的方法和设备 |
| US20150105638A1 (en) * | 2009-07-30 | 2015-04-16 | Oxitone Medical Ltd. | Photoplethysmography Device and Method |
| CN104605838A (zh) * | 2015-01-19 | 2015-05-13 | 深圳市元征科技股份有限公司 | 一种心率检测的方法及装置 |
| US20160020798A1 (en) * | 2012-01-31 | 2016-01-21 | Broadcom Corporation | Systems and methods for enhancing audio quality of fm receivers |
| CN105644365A (zh) * | 2016-02-02 | 2016-06-08 | 深圳还是威健康科技有限公司 | 心率监测方法、装置及系统 |
| CN105832323A (zh) * | 2016-06-12 | 2016-08-10 | 电子科技大学 | 一种基于rc电路的光电血管容积图信号处理方法 |
| CN106994010A (zh) * | 2016-01-26 | 2017-08-01 | 深圳市新元素健康管理有限公司 | 一种基于ppg信号的心率检测方法及系统 |
| CN107137071A (zh) * | 2017-04-26 | 2017-09-08 | 可瑞尔科技(扬州)有限公司 | 一种分析心冲击信号用来计算短期心率值的方法 |
| CN107249449A (zh) * | 2015-03-16 | 2017-10-13 | 努沃集团有限公司 | 用于感测胎儿活动的系统、装置和方法 |
| US9814400B1 (en) * | 2015-05-26 | 2017-11-14 | Verily Life Sciences Llc | Method for improving accuracy of pulse rate estimation |
| CN107510462A (zh) * | 2017-09-06 | 2017-12-26 | 广州视源电子科技股份有限公司 | 血氧测量方法和装置 |
-
2018
- 2018-06-27 CN CN201810677236.7A patent/CN108992054B/zh active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150105638A1 (en) * | 2009-07-30 | 2015-04-16 | Oxitone Medical Ltd. | Photoplethysmography Device and Method |
| CN101732050A (zh) * | 2009-12-04 | 2010-06-16 | 西安交通大学 | 一种基于光电容积波的呼吸率监测方法 |
| CN103153175A (zh) * | 2010-05-14 | 2013-06-12 | 先进计算发展中心 | 脉搏信号波形数据的诊断分类 |
| CN103027667A (zh) * | 2011-09-30 | 2013-04-10 | Ge医疗系统环球技术有限公司 | 脉搏波的特征参数提取 |
| US20160020798A1 (en) * | 2012-01-31 | 2016-01-21 | Broadcom Corporation | Systems and methods for enhancing audio quality of fm receivers |
| CN104349713A (zh) * | 2012-04-23 | 2015-02-11 | 图尔库大学 | 确定指示心脏功能障碍和异常的信息的方法和设备 |
| CN104605838A (zh) * | 2015-01-19 | 2015-05-13 | 深圳市元征科技股份有限公司 | 一种心率检测的方法及装置 |
| CN107249449A (zh) * | 2015-03-16 | 2017-10-13 | 努沃集团有限公司 | 用于感测胎儿活动的系统、装置和方法 |
| US9814400B1 (en) * | 2015-05-26 | 2017-11-14 | Verily Life Sciences Llc | Method for improving accuracy of pulse rate estimation |
| CN106994010A (zh) * | 2016-01-26 | 2017-08-01 | 深圳市新元素健康管理有限公司 | 一种基于ppg信号的心率检测方法及系统 |
| CN105644365A (zh) * | 2016-02-02 | 2016-06-08 | 深圳还是威健康科技有限公司 | 心率监测方法、装置及系统 |
| CN105832323A (zh) * | 2016-06-12 | 2016-08-10 | 电子科技大学 | 一种基于rc电路的光电血管容积图信号处理方法 |
| CN107137071A (zh) * | 2017-04-26 | 2017-09-08 | 可瑞尔科技(扬州)有限公司 | 一种分析心冲击信号用来计算短期心率值的方法 |
| CN107510462A (zh) * | 2017-09-06 | 2017-12-26 | 广州视源电子科技股份有限公司 | 血氧测量方法和装置 |
Non-Patent Citations (5)
| Title |
|---|
| ELGENDI, MOHAMED; NORTON, IAN; BREARLEY, MATT; 等.: "Systolic Peak Detection in Acceleration", 《PLOS ONE》 * |
| SUN, BIAO; ZHANG, ZHILIN: "Photoplethysmography-Based Heart Rate Monitoring Using Asymmetric Least Squares Spectrum Subtraction and Bayesian Decision Theory", 《IEEE SENSORS JOURNAL》 * |
| 张爱华，王平，丑永新: "基于动态差分阈值的脉搏信号峰值检测算法", 《吉林大学学报（工学版）》 * |
| 赵雪峰: "三元复合驱采出液静电聚结特性研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 金晟来: "基于脉搏色素浓度法检测肝储备功能与心血流动力学参数的方法研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113032724A (zh) * | 2019-12-09 | 2021-06-25 | 中国石油化工股份有限公司 | 一种示踪剂产出曲线降噪处理方法、存储介质及电子设备 |
| CN113032724B (zh) * | 2019-12-09 | 2023-10-27 | 中国石油化工股份有限公司 | 一种示踪剂产出曲线降噪处理方法、存储介质及电子设备 |
| CN113317757A (zh) * | 2021-04-30 | 2021-08-31 | 深圳麦格米特电气股份有限公司 | 一种获取生命体征数据的方法、光纤传感器及设备 |
| CN114041768A (zh) * | 2021-10-11 | 2022-02-15 | 广东宝莱特医用科技股份有限公司 | 脉率提取方法、设备及存储介质 |
| CN114041768B (zh) * | 2021-10-11 | 2023-08-08 | 广东宝莱特医用科技股份有限公司 | 脉率提取方法、设备及存储介质 |
| CN114521880A (zh) * | 2022-01-21 | 2022-05-24 | 中国人民解放军陆军军医大学 | 一种运动状态下心率计算方法、系统及计算机存储介质 |
| CN114521880B (zh) * | 2022-01-21 | 2023-09-01 | 中国人民解放军陆军军医大学 | 一种运动状态下心率计算方法、系统及计算机存储介质 |
| CN115349816A (zh) * | 2022-10-21 | 2022-11-18 | 万灵帮桥医疗器械(广州)有限责任公司 | 干涉信号峰值检测方法、装置、设备和存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108992054B (zh) | 2021-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108992054A (zh) | 一种脉搏信号峰值点检测方法及装置 | |
| JP6761265B2 (ja) | 光電式容積脈波信号のノイズ除去方法とシステム | |
| JP6335199B2 (ja) | ハンドヘルド・デバイスを用いた非侵襲的な脈波伝播時間の判定 | |
| CN109189207B (zh) | 一种可穿戴设备穿戴状态的检测方法及装置 | |
| CN109222949A (zh) | 心率检测方法和心率检测装置 | |
| CN111091116A (zh) | 一种用于判断心律失常的信号处理方法及系统 | |
| CN108513665A (zh) | 用户身份识别的方法、装置和系统 | |
| JP6869691B2 (ja) | 信号特徴抽出方法及びその装置並びにコンピュータ読取可能記録媒体 | |
| KR20190011026A (ko) | 혈압 측정 장치 및 방법 | |
| Pangerc et al. | Robust detection of heart beats in multimodal records using slope-and peak-sensitive band-pass filters | |
| CN112568886A (zh) | 心律的检测方法、电子设备以及计算机可读存储介质 | |
| CN109077721B (zh) | 房颤检测装置及存储介质 | |
| CN111588367A (zh) | 心率检测方法、装置及计算机可读存储介质 | |
| Chin et al. | Bayesian real-time QRS complex detector for healthcare system | |
| CN110680307A (zh) | 一种运动环境下基于脉搏波传导时间的动态血压监测方法 | |
| Liu et al. | The accuracy on the common Pan-Tompkins based QRS detection methods through low-quality electrocardiogram database | |
| Sundar et al. | A comprehensive assessment of the performance of modern algorithms for enhancement of digital volume pulse signals | |
| CN114027804A (zh) | 一种脉象诊断的方法、装置及可读存储介质 | |
| Pal et al. | Improved heart rate detection using smart phone | |
| CN104739384B (zh) | 动脉波分析方法及其系统 | |
| WO2019146025A1 (ja) | 脈波算出装置、脈波算出方法及び脈波算出プログラム | |
| Lin et al. | A characteristic filtering method for pulse wave signal quality assessment | |
| Chou et al. | Modelling arterial blood pressure waveforms for extreme bradycardia and tachycardia by curve fitting with Gaussian functions | |
| CN106491118A (zh) | 基于安卓系统的心电图机的实时心率计算系统及其方法 | |
| Yan et al. | A resource-efficient, robust QRS detector using data compression and time-sharing architecture |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190514 Address after: 518000 Yuanzheng Industrial Park, North of Fifth Avenue, Bantian Street, Longgang District, Shenzhen City, Guangdong Province Applicant after: Yuanzheng Science and Technology Co., Ltd., Shenzhen City Address before: 518000 Room 201, building A, No. 1, Qian Wan Road, Qianhai Shenzhen Hong Kong cooperation zone, Shenzhen, Guangdong (Shenzhen Qianhai business secretary Co., Ltd.) Applicant before: HESVIT HEALTH TECH CO., LTD. |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |