CN114676738B - Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 - Google Patents
Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 Download PDFInfo
- Publication number
- CN114676738B CN114676738B CN202210585033.1A CN202210585033A CN114676738B CN 114676738 B CN114676738 B CN 114676738B CN 202210585033 A CN202210585033 A CN 202210585033A CN 114676738 B CN114676738 B CN 114676738B
- Authority
- CN
- China
- Prior art keywords
- signal
- ppg signal
- original
- frequency
- peak
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 75
- 238000013441 quality evaluation Methods 0.000 title claims abstract description 48
- 238000003672 processing method Methods 0.000 title claims abstract description 8
- 238000005311 autocorrelation function Methods 0.000 claims abstract description 82
- 238000001228 spectrum Methods 0.000 claims abstract description 45
- 238000001914 filtration Methods 0.000 claims description 59
- 238000012545 processing Methods 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 9
- 238000001303 quality assessment method Methods 0.000 claims description 8
- 238000013186 photoplethysmography Methods 0.000 claims description 7
- 238000012790 confirmation Methods 0.000 claims description 2
- 238000009532 heart rate measurement Methods 0.000 claims description 2
- 238000005314 correlation function Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 22
- 238000010586 diagram Methods 0.000 description 11
- 230000009471 action Effects 0.000 description 9
- 238000004364 calculation method Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 5
- 239000007789 gas Substances 0.000 description 5
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000012854 evaluation process Methods 0.000 description 4
- 230000036772 blood pressure Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000035772 mutation Effects 0.000 description 3
- 230000000737 periodic effect Effects 0.000 description 3
- 238000013139 quantization Methods 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 230000017531 blood circulation Effects 0.000 description 2
- 238000013135 deep learning Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000008667 sleep stage Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
-
- 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
-
- 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/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/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- 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
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/14—Classification; Matching by matching peak patterns
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Signal Processing (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Hematology (AREA)
- Vascular Medicine (AREA)
- Power Engineering (AREA)
- Mathematical Physics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明涉及PPG信号质量评估方法及装置以及PPG信号处理方法及系统，包括：S1、获取原始信号，S2、当其平均幅度差在第一预设范围时执行S3，否则执行S11；S3、获取有用信号；S4、当其高频和低频信噪比分别大于对应的第一门限值时执行S5，否则执行S11；S5、当其功率谱的额峰度和偏度分别大于第二门限值时执行S6，否则执行S11；S6、当相关函数的平均峰值周期在第二预设范围内时执行S7，否则执行S11；S7、当出现最大峰值大于第三门限值或峰值周期标准差小于第四门限值时执行S8，否则执行S11；S8、当自相关函数的峰峰差值标准差小于第五门限值时执行S9，否则执行S10；S9、判定为第一等级；S10、判定为第二等级；S11、判定为第三等级。
Description
技术领域
本发明涉及信号分析及处理技术领域，更具体地说，涉及PPG信号质量评估方法及装置以及PPG信号处理方法及系统。
背景技术
光电容积脉搏波PPG(Photoplethysmograph)技术是一种利用光电的方法测量一次心跳周期血流量变化的技术，其主要利用光电二极管射向皮肤，然后通过接收反射光的强度，获得能够反映血流量的波形。这些波形包含人体丰富的生理信息。由于PPG便携性、低成本以及非侵入性受到越来越多的科技公司、科研机构的欢迎，应用也非常广泛：心率检测、血氧计算、血压估计和睡眠分期等生理指标的计算。
然而PPG在实际采集的过程中难免受到工频噪声、运动伪影和环境光的干扰，当PPG波形严重失真或者存在严重噪声时，是无法通过滤波等预处理技术得到可用的PPG波形用于准确计算生理指标的。因此，在计算生理指标之前需要对PPG信号的质量进行评估，将严重失真的不可用PPG信号舍弃，将质量好的干净的PPG信号用于生理指标的计算，这样可以有效提高生理指标计算的准确性、降低可穿戴设备的功耗，并且在配合硬件电路选型以及嵌入式参数调节时，可以提供客观的PPG信号质量评估的量化指标。
目前已有的方法有模板匹配方法、特征提取加机器学习的方法和深度学习的方法。模板匹配的方法主要是将一段PPG信号分割成单个心跳周期，找到代表性的心跳周期模板然后和所有心跳周期进行匹配，将匹配比例高的评估为质量良好的PPG波形，然而这类方法太依赖周期分割的算法了，因为PPG波形的形态是十分复杂多样的，周期分割不准确会误将好的信号评价为差的信号，而且十分依赖模板的选取标准，模板选择不准确会严重影响算法的鲁棒性。基于机器学习和深度学习的方法也能取得很好的效果，但是需要依赖标签的准确性以及数据集的泛化性，而且给数据打标签的过程十分耗时。而且大部分方法的评估结果只是将信号分为可用或者不可用的二元分类，这样的分类结果会比较粗糙。
发明内容
本发明要解决的技术问题在于，提供PPG信号质量评估方法及装置以及PPG信号处理方法及系统。
本发明解决其技术问题所采用的技术方案是：构造一种PPG信号质量评估方法，包括步骤：
S1、获取预设时长的原始PPG信号，其中，所述PPG信号为光电容积脉搏波形信号；
S2、预处理所述原始PPG信号以得到若干预设长度的波形片段，获取若干所述波形片段的平均幅度差，判断所述平均幅度差是否在第一预设范围内，若是，则执行步骤S3，否则执行步骤S11；
S3、对所述原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取所述原始PPG信号的有用信号；
S4、基于所述有用信号和所述高通滤波结果获取所述原始PPG信号的高频信噪比，基于所述有用信号和所述低通滤波结果获取所述原始PPG信号的低频信噪比，判断所述高频信噪比和所述低频信噪比是否分别大于其对应的第一门限值，若是，执行步骤S5，否则执行步骤S11；
S5、获取所述有用信号的自相关函数，对所述自相关函数进行傅里叶变换以得到所述有用信号的功率谱，分别获取所述功率谱的峰度和偏度，判断所述功率谱的峰度和偏度是否分别大于其对应的第二门限值，若是，执行步骤S6，否则执行步骤S11；
S6、获取所述自相关函数的平均峰值周期，确认所述平均峰值周期是否在第二预设范围内，若是，则执行步骤S7，否则执行步骤S11；
S7、获取所述自相关函数的最大峰值和峰值周期标准差，确认是否出现所述最大峰值大于第三门限值或所述峰值周期标准差小于第四门限值，若是，则执行步骤S8，否则执行步骤S11；
S8、获取所述自相关函数的峰峰差值标准差，确认所述峰峰差值标准差是否小于第五门限值，若是，则执行步骤S9，否则，执行步骤S10；
S9、判定所述原始PPG信号为第一等级信号，并结束；
S10、判定所述原始PPG信号为第二等级信号，所述第二等级信号质量低于所述第一等级信号质量，并结束；
S11、判定所述原始PPG信号为第三等级信号，所述第三等级信号质量低于所述第二等级信号质量，并结束。
优选地，在本发明的PPG信号质量评估方法中，在所述步骤S1中，所述预设时长为大于或等于5s小于或等于15s。
优选地，在本发明的PPG信号质量评估方法中，在所述步骤S2中，
所述预处理所述原始PPG信号以得到若干预设长度的波形片段，包括：
对所述原始PPG信号进行中值滤波，对滤波后的信号以所述预设长度的滑动窗口依次进行无间隔截取以获取若干所述波形片段，其中所述波形片段之间不重叠；和/或，
所述获取所述若干波形片段的平均幅度差，包括：
获取每一波形片段的幅度差，并获取所有波形片段的幅度差的平均值为所述平均幅度差。
优选地，在本发明的PPG信号质量评估方法中，所述第一预设范围大于或等于150且小于或等于90000；和/或，
所述预设长度为1s。
优选地，在本发明的PPG信号质量评估方法中，在所述步骤S3中，所述对所述原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取所述原始PPG信号的有用信号；包括：
S31、以截止频率为第一频率的巴特沃斯高通滤波器对所述原始PPG信号进行高通滤波，以得到高频噪声信号；
S32、以截止频率为第二频率的巴特沃斯低通滤波器对所述原始PPG信号进行低通滤波，以得到低频噪声信号；
S33、对所述原始PPG信号分别剔除所述高频噪声信号和所述低频噪声信号，以剩余信号为所述有用信号。
优选地，在本发明的PPG信号质量评估方法中，在所述步骤S4中，
所述基于所述有用信号和所述高通滤波结果获取所述原始PPG信号的高频信噪比，包括：获取所述高频噪声信号与所述有用信号的比值为所述高频信噪比；
所述基于所述有用信号和所述低通滤波结果获取所述原始PPG信号的低频信噪比，包括：获取所述低频噪声信号与所述有用信号的比值为所述低频信噪比。
优选地，在本发明的PPG信号质量评估方法中，在所述步骤S4中，所述高频信噪比对应的第一门限值大于或等于6dB，所述低频信噪比对应的第一门限值大于或等于-15dB；和/或，
在所述步骤S5中，所述功率谱的偏度对应的第二门限值大于或等于12，所述功率谱的峰度对应的第二门限值大于或等于150。
优选地，在本发明的PPG信号质量评估方法中，在所述步骤S5中，所述获取所述有用信号的自相关函数，包括：对所述有用信号做单边归一化的自相关函数，其中所述自相关函数满足以下公式：
其中， 
为所述自相关函数，k为所述自相关函数的索引值，
为所述有用信
号，n为所述有用信号的索引值，N为所述有用信号的点数。
优选地，在本发明的PPG信号质量评估方法中，所述方法还包括：
基于二阶导数的波谷检测算法对自相关函数进行峰值检测以获取若干峰值坐标；
基于所述若干峰值坐标分别获取所述自相关函数的平均峰值周期、最大峰值、峰值周期标准差和峰峰差值标准差。
优选地，在本发明的PPG信号质量评估方法中，所述方法包括以下参数设置中的一个或多个；
所述第二预设范围为大于0.3s且小于2s；
所述第三门限值大于或等于0.6；
所述第四门限值小于或等于5；
所述第五门限值小于或等于0.07。
另，本发明还构造一种PPG信号处理方法，包括，通过上面任意一项所述的PPG信号质量评估方法获取原始PPG信号的质量评估结果；并在所述原始PPG信号为第三等级信号时，剔除所述原始PPG信号；
所述原始PPG信号为第二等级信号时，将所述原始PPG信号用于部分预设生理指标计算；
所述原始PPG信号为第一等级信号时，将原始PPG信号用于与所述PPG信号相关的所有生理指标计算。
另，本发明还构造一种PPG信号质量评估装置，包括：
原始信号获取单元，用于获取预设时长的原始PPG信号，其中，所述PPG信号为光电容积脉搏波形信号；
第一判断单元，用于预处理所述原始PPG信号以得到若干预设长度的波形片段，获取所述若干波形片段的平均幅度差，判断所述平均幅度差是否在第一预设范围内，若是，则输出肯定结果，否则输出否定结果；
有用信号获取单元，用于对所述原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取所述原始PPG信号的有用信号；
第二判断单元，用于基于所述有用信号和所述高通滤波结果获取所述原始PPG信号的高频信噪比，基于所述有用信号和所述低通滤波结果获取所述原始PPG信号的低频信噪比，判断所述高频信噪比和所述低频信噪比是否分别大于其对应的第一门限值，若是，则输出肯定结果，否则输出否定结果；
第三判断单元，用于获取所述有用信号的自相关函数，对所述自相关函数进行傅里叶变换以得到所述有用信号的功率谱，分别获取所述功率谱的峰度和偏度，判断所述功率谱的峰度和偏度是否分别大于其对应的第二门限值，若是，则输出肯定结果，否则输出否定结果；
第四判断单元，用于获取所述自相关函数的平均峰值周期，判断所述平均峰值周期是否在第二预设范围内，若是，则输出肯定结果，否则输出否定结果；
第五判断单元，用于获取所述自相关函数的最大峰值和峰值周期标准差，确认是否出现所述最大峰值大于第三门限值或所述峰值周期标准差小于第四门限值，若是，则输出肯定结果，否则输出否定结果；
第六判断单元，用于获取所述自相关函数的峰峰差值标准差，确认所述峰峰差值标准差是否小于第五门限值，若是，则输出肯定结果，否则输出否定结果；
结果确认单元，用于在所述第一判断单元、所述第二判断单元、所述第三判断单元、所述第四判断单元和所述第六判断单元中任意一个输出否定结果时，判定所述原始PPG信号为第三等级信号，在所述第六判断单元输出否定结果时，判定所述原始PPG信号为第二等级信号，在所述第六判断单元输出肯定结果时，判定所述原始PPG信号为第一等级信号，其中，所述第三等级信号质量低于所述第二等级信号质量，所述第二等级信号质量低于所述第一等级信号质量。
另，本发明还构造一种PPG信号处理系统，包括如上面所述的PPG信号质量评估装置，以及处理单元；
所述处理单元用于获取原始PPG信号的质量评估结果；并在所述原始PPG信号为第三等级信号时，剔除所述原始PPG信号；
所述原始PPG信号为第二等级信号时，将所述原始PPG信号用于部分预设生理指标计算；
所述原始PPG信号为第一等级信号时，将原始PPG信号用于与所述PPG信号相关的所有生理指标计算。
实施本发明的PPG信号质量评估方法及装置以及PPG信号处理方法及系统，具有以下有益效果：能够实现PPG信号的多层级量化指标的评价，以实现基于PPG信号的不同质量等级进行合理的利用。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明一实施例中的一种PPG信号质量评估方法的程序流程图；
图2是本发明一实施例中的原始PPG信号的示意图；
图3是本发明一实施例中的原始PPG信号的平均幅度差的示意图；
图4是本发明一实施例中的原始PPG信号的低频信噪比的示意图；
图5是本发明一实施例中的原始PPG信号的低频信噪比的示意图；
图6是本发明一实施例中的原始PPG信号的功率谱指标的示意图；
图7是本发明另一实施例中的原始PPG信号的功率谱指标的示意图；
图8是本发明一实施例中的原始PPG信号的自相关函数指标的示意图；
图9是本发明另一实施例中的原始PPG信号的自相关函数指标的示意图；
图10是本发明再一实施例中的原始PPG信号的自相关函数指标的示意图；
图11是本发明一实施例中的一种PPG信号质量评估装置的逻辑框图；
图12是本发明一实施例中的一种PPG信号质量处理系统的逻辑框图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
如图1所示，在本发明的一种PPG信号质量评估方法第一实施例中，包括：S1、获取
预设时长的原始PPG信号，其中，PPG信号为光电容积脉搏波形信号。具体的，可以通过采集
设备获取一段预设时长的PPG信息即对应原始PPG信号。其中，原始PPG信号的时间长度不能
设置太长，否则信号的波动十分剧烈，而且计算复杂度会随着时长的增加而升高。其原始
PPG信号的长度至少是5s，最多不超过15s。根据正常最低心率为30次每分钟，5s时间段的波
形至少包含两到三个脉冲周期。如图2所示，其可以选择原始PPG信号的长度为10s，其在对
应采样频率
时，信号的总点数为1000。
S2、预处理原始PPG信号以得到若干预设长度的波形片段，获取若干波形片段的平均幅度差，判断平均幅度差是否在第一预设范围内，若是，则执行步骤S3，否则执行步骤S11。具体的，当设备和皮肤接触不良时或者没有佩戴时，采集设备得到的原始PPG信号会很微弱甚至没有，此时原始PPG信号的平均幅度差很低。相反，当硬件电路中出现饱和噪声或者尖峰脉冲时，信号的平均幅度差会变得异常高。因此，可以对原始PPG信号进行分段，得到预设长度的若干波形片段，基于该若干波形片段得到其平均幅度差。对该平均幅度差进行判断。当平均幅度差满足第一预设范围的要求时，即此时可以判断该原始PPG信号的幅度正常，其可以继续执行步骤S3及其之后的动作。当平均幅度差不满足第一预设范围的要求时，则判断该原始PPG信号的幅度异常，此时可以直接执行后面的判定动作。即可以直接判定原始PPG信号为第三等级信号，并结束该原始PPG信号的质量评估过程。该过程可以定义为第一层级评估过程，正常PPG波形的平均幅度差的范围是几百到几万。但是当设备和皮肤接触不良或者没有佩戴时，PPG信号的平均幅度差低于150时，PPG信号的波形近似一条直线。如图3所示，当硬件电路中的电流调节得过大时，PPG信号的平均幅度差会大于90000，PPG信号的波形就会出现饱和失真。如果原始PPG信号的平均幅度差出现异常，将评估该原始PPG信号为差的信号，就没有必要进行下一层级的评估了。这里的150和9000对应的就是离散后信号的数字，其由设定的量化精度决定。
在一实施例中，在上述步骤S2中，预处理原始PPG信号以得到若干预设长度的波形片段，包括：对原始PPG信号进行中值滤波，对滤波后的信号以预设长度的滑动窗口依次进行无间隔截取以获取若干波形片段，其中波形片段之间不重叠。即，其可以首先对原始PPG信号进行窗口大小为5的中值滤波，去除原始PPG信号中离群点的影响，得到信号 x[n]。其中，在该中值滤波滤波过程中，五个数中取中间那个数，这样就可以去除异常离群点的影响。比如信号为【1，2，3，4，25】，最后输出的是3，避免了离群点25的影响了，以L点长（对应预设长度）的滑动窗口将信号x[n]分成M个不重叠的片段，每个片段记作：
其中，
为波形片段，下标
为该波形片段的索引值，即代表第几个小信号片
段，
为该波形片段的长度。如对信号x[n]，其中n取值是【0， N-1】，分别对应10s的信号。在M
为10时得到的波形片段
，其中
取值【0，M-1】，分别对应1s的信号片段。
在一实施例中，当信号的采样频率
信号x[n]长度为10s，L取100个点，
即对应的预设长度也可以理解为1s，所以M=10。
在一实施例中，在上述步骤S2中，获取若干波形片段的平均幅度差，包括：获取每
一波形片段的幅度差，并获取所有波形片段的幅度差的平均值为平均幅度差。即，对得到的
波形片段计算每个波形片段的幅度差
，
其中，
为
波形片段中的幅度最大值，
为
波形片段中
的幅度最小值，所以得到所有的波形片段的平均幅度差
为：
其中，对平均幅度差
的判断可以为，在其大于或等于150且小于或等于90000
判定其满足第一预设范围。
S3、对原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取原始PPG信号的有用信号。具体的，当在原始PPG信号的幅度指标满足要求时继续对原始PPG信号的噪声进行判断。虽然，PPG信号的主要频率是0.5Hz-5Hz，大于5Hz的频率是信号的高频毛刺噪声，能量占比少；小于0.5Hz的是基线漂移相关的低频噪声，占了绝大部分能量。但在信号处理过程中，需要对高频噪声和低频噪声的影响分别进行判断。假设一段信号没有基线漂移，即它的低频噪声很小，但是高频噪声非常大，那么这段信号就是高频噪声太大导致信号无法使用了。因此需要剔除其中的高频噪声和低频噪声，最终得到原始PPG信号的有用信号，并基于有用信号对原始PPG信号进行判定。
可选的，在步骤S3中，对原始PPG信号分别进行高通滤波和低通滤波，基于高通滤
波结果和低通滤波结果获取原始PPG信号的有用信号；包括：S31、以截止频率为第一频率的
巴特沃斯高通滤波器对原始PPG信号进行高通滤波，以得到高频噪声信号；S32、以截止频率
为第二频率的巴特沃斯低通滤波器对原始PPG信号进行低通滤波，以得到低频噪声信号；
S33、对原始PPG信号分别剔除高频噪声信号和低频噪声信号，以剩余信号为有用信号。具体
的，其可以对原始PPG信号分别通过高通滤波器进行滤波得到高频噪声信号和通过低通滤
波器进行滤波得到低频噪声信号。在这里高通滤波器和低通滤波器均可以采用巴特沃斯滤
波器。此外，可以设置高通滤波器的截止频率即第一频率
, 得到高频噪声信
号
，则高频噪声功率
如下：
可以设置低通滤波器的截止频率即第二频率
, 得到低频噪声信号
，则低频噪声功率
如下：
最后，将原始PPG信号分别减去高频和低频噪声信号，可以得到其对应的有用信号
，则有用信号功率
如下：
可以理解，对信号进行滤波后，得到高频噪声信号
和低频噪声信号
对应的长度与原始信号x[n]长度一样，例如，N=1000。
S4、基于有用信号和高通滤波结果获取原始PPG信号的高频信噪比，基于有用信号和低通滤波结果获取原始PPG信号的低频信噪比，判断高频信噪比和低频信噪比是否分别大于其对应的第一门限值，若是，执行步骤S5，否则执行步骤S11。具体的，由于高频噪声和低频噪声对原始PPG信号的影响结果不同，因此需要进行高频信噪比和低频信噪比的区分，其可以根据高通滤波结果和低通滤波结果分别获取原始PPG信号的高频信噪比和低频信噪比，并在其高频信噪比和低频信噪比均满足要求的情况下，判定原始PPG信号的噪声满足要求，可以继续执行后面步骤S5的动作。当出现高频噪声比和低频噪声比中任意一个不满足要求的情况，则可以判断该原始PPG信号的噪声异常，此时可以直接执行后面的判断动作。即可以直接判定原始PPG信号为第三等级信号，并结束该原始PPG信号的质量评估过程。该过程可以定义为第二层级评估过程。PPG信号的噪声主要来源于高频毛刺噪声和低频的漂移噪声。如图4所示，当低频信噪比小于-15dB时，PPG信号的整体漂移十分严重，波动十分剧烈，完全是不规则的状态。如图5所示，当高频信噪比小于6dB时，PPG波形形态完全被毛刺掩盖，即使去除掉高频噪声信号和低频噪声信号也无法观察到其形态和周期。总之，频低频信噪比太低时，信号十分微弱被噪声掩盖，甚至通带内噪声的能量也远比信号的能量高，或者根本没有检测到PPG信号。因此，在信噪比出现异常时可以将此原始PPG信号直接判定为差的信号，没必要进行第三层的评估。
可选的，在步骤S4中，基于有用信号和高通滤波结果获取原始PPG信号的高频信噪
比，包括：获取高频噪声信号与有用信号的比值为高频信噪比；基于有用信号和低通滤波结
果获取原始PPG信号的低频信噪比，包括：获取低频噪声信号与有用信号的比值为低频信噪
比。其具体的高频信噪比
获取过程如下：
其具体的低频信噪比
获取过程如下：
在一实施例中，高频信噪比对应的第一门限值大于或等于6dB，低频信噪比对应的
第一门限值大于或等于-15dB。即，在高频信噪比对应的第一门限值为6dB的时候，即在高频
信噪比
大于6dB时可以判定该高频信噪比满足要求。在低频信噪比对应的第一门
限值为-15dB时，即在低频信噪比
大于-15dB时可以判定低频信噪比满足要求。
S5、获取有用信号的自相关函数，对自相关函数进行傅里叶变换以得到有用信号的功率谱，分别获取功率谱的峰度和偏度，判断功率谱的峰度和偏度是否分别大于其对应的第二门限值，若是，执行步骤S6，否则执行步骤S11。具体的，在原始PPG信号的噪声指标满足要求时，对功率谱相关指标进行判断。其主要对原始PPG信号对应的频率普的偏度和峰度进行判断。当原始PPG信号中存在有间断点和尖脉冲时，这些突变会导致功率谱系数将分布整个频率轴上，无法得到好的稀疏表示。这类突变信号的高频信噪比和低频信噪比都可能很高，无法通过信噪比相关指标检测出来，但其可以通过计算有用信号的功率谱的偏度和峰度衡量其稀疏表示。其先获取有用信号的自相关函数，并基于该自相关函数进行傅里叶变化得到有用信号的功率谱。基于得到的功率谱获取其对应的偏度和峰度。该过程可以定义为对功率谱相关指标可进行评估的第三层级评估过程。可以理解，第一层级和第二层层级的评估过程为从整体指标的角度去评价信号，当信号中局部小区域出现间断点和尖脉冲时，这些突变会导致功率谱系数将分布整个频率轴上，无法得到很好的稀疏表示。如图6和图7所示，此时通过功率谱指标能很好区分信号中是否出现突变区域造成的异常。
在一具体实施例中，其可以对有用信号
做单边归一化的自相关函数，其具体
函数
表达式如下：
可以理解，由于N为有用信号的点数，对信号进行单边归一化的自相关函数的长度
也是1000点。为了区分原始信号x[n]，把自相关函数称为
，其中 k为所述自相关函数的
索引值，其取值范围是0到N-1，n代表着有用信号的索引值，其取值的范围也是0到N-1。
由维纳-辛钦定理可知，自相关函数的傅里叶变换对应为信号的功率谱。对自相关
函数
进行傅里叶变换得到有用信号的功率谱
满足如下表达式：
功率谱的峰度
和偏度
分别如下：
其中，μ为
的均值，σ为
方差。
对功率谱的偏度和峰度同时进行判断，在其均满足对应的门限要求时，即均大于其对应的第二门限值要求时，则判定功率谱的信号为正常，可以继续执行后面步骤S6的动作。当功率谱的偏度和峰度中任意一个不满足其对应的第二门限值，则可以判断该原始PPG信号的存在突变区域，此时直接执行后面的判断动作。即可以直接判定原始PPG信号为第三等级信号，并结束该原始PPG信号的质量评估过程。
在一实施例中，在步骤S5中，功率谱的偏度对应的第二门限值大于或等于12，功率谱的峰度对应的第二门限值大于或等于150。即，可以分别设置功率谱的偏度对应的第二门限值和功率谱的峰度对应的第二门限值。
S6、获取自相关函数的平均峰值周期，确认平均峰值周期是否在第二预设范围内，若是，则执行步骤S7，否则执行步骤S11；具体的，由于自相关函数计算的信号本身与其时间延迟之间的相关性，当信号存在周期性成分时，自相关函数会在周期的整数倍出现极大值。因此可以通过检测自相关函数的峰值检验其周期性。即对自相关函数的平均峰值周期进行判定。当其在第二预设范围内，则判定原始PPG信号在被采样信号例如心率信号的正常范围内，其可以继续步骤S7的动作。当自相关函数的平均峰值周期不在第二预设范围时，则判定原始PPG信号不在心率的正常范围内，此时可以直接执行后面的判断动作，即可以直接判定原始PPG信号为第三等级信号，并结束该原始PPG信号的质量评估过程。
在一实施例中，第二预设范围为大于0.3s且小于2s。即可以在一具体实施例中，设置当平均峰值周期大于0.3s小于2s时，判定其在第二预设范围内。
S7、获取自相关函数的最大峰值和峰值周期标准差，确认是否出现最大峰值大于第三门限值或峰值周期标准差小于第四门限值，若是，则执行步骤S8，否则执行步骤S11；具体的，在自相关函数的平均峰值周期满足要求即判定原始PPG信号在检测信号的正常范围内时，对自相关函数的最大峰值和峰值周期标准差进行判断，即对原始PPG信号周期性进行判断。当判定原始PPG信号存在周期性时，其可以继续步骤S8的动作。当判定原始PPG信号不存在周期性时，可以直接执行后面的判断动作，即可以直接判定原始PPG信号为第三等级信号，并结束该原始PPG信号的质量评估过程。其具体判定过程为在同时出现自相关函数的最大峰值大于第三门限值和峰值周期标准差小于第四门限值时，判定原始PPG信号不存在周期性，否则可以判定原始PPG信号存在周期性。
在一实施例中，第三门限值大于或等于0.6，第四门限值小于或等于5；具体的，其可以在相关函数的最大峰值大于0.6时判定原始PPG信号存在周期性，其也可以在自相关函数的峰值周期标准差小于5时判定原始PPG信号存在周期性。
S8、获取自相关函数的峰峰差值标准差，确认峰峰差值标准差是否小于第五门限值，若是，则执行步骤S9，否则，执行步骤S10；S9、判定原始PPG信号为第一等级信号，并结束；S10、判定原始PPG信号为第二等级信号，第二等级信号质量低于第一等级信号质量，并结束；S11、判定原始PPG信号为第三等级信号，第三等级信号质量低于第二等级信号质量，并结束。具体的，在根据自相关函数的最大峰值和峰值周期标准差即判定原始PPG信号在检测信号的存在周期性时，对自相关函数的峰峰差值标准差进行判断，根据峰峰差值标准差判定原始PPG信号的各个周期的信号形态是否相似，其在自相关函数的峰峰差值标准差小于第五门限值时，可以判定其信号形态相似，此时可以判定该判定原始PPG信号为第一等级信号，否则判定该原始PPG信号为第二等级信号。其中第一等级信号的质量优于第二等级信号，第二等级信号的质量优于第三等级信号。
可选的，本发明的PPG信号质量评估方法还包括：基于二阶导数的波谷检测算法对自相关函数进行峰值检测以获取若干峰值坐标；基于若干峰值坐标分别获取自相关函数的平均峰值周期、最大峰值、峰值周期标准差和峰峰差值标准差。具体的，其可以先使用基于二阶导数的波谷检测算法对自相关函数进行峰值检测，得到w个所有峰值坐标：
则自相关函数的最大峰值
为：
相邻的峰峰间隔
为：
则自相关函数的平均峰值周期
为：
则自相关函数的峰值周期标准差
为：
峰峰差值
为：
平均峰峰差值
为：
则自相关函数的峰峰差值的标准差
为：
对与自相关函数的相关指标判定的过程可以定义为第四层级评估。其通过相比直接对有用信号进行周期性检验，计算自相关函数，能够消除加性噪声对检验信号周期性的影响。自相关函数的平均峰值周期对应原始PPG信号的周期，根据正常心率30次~200次每分钟，平均峰值周期应该在0.3s和2s之间，超过正常心率范围的评估为差的信号（如图8所示）。当自相关函数的最大峰值大于0.6或者峰值周期标准差小于5时，可以判断出信号是存在周期性的，用具有周期性的信号计算心率等生理指标其准确性会更高，此时，当相邻周期PPG信号形态的一致性、相似性越高时，自相关函数应该是一个逐步衰减的正弦波信号，而个别存在周期性的信号其波形形态可能受到通带内的噪声影响，相邻周期PPG信号的波形形态各不相同，此时其自相关函数的峰值高低起伏不一，所以通过计算自相关函数的峰峰差值的标准差，能够判断自相关函数峰值的波动情况，当自相关函数的峰峰差值的标准差小于0.07时，评价为好的信号（如图9所示），否则评价为中等信号（如图10所示）。
根据以上过程中信号质量评估的层级关系，如果前一层的量化指标无法满足要求，就没必要进行下一层级的判断了，这样可以有效减少误检的概率，从而提升质量评估算法的准确性。
另，本发明的一种PPG信号处理方法，包括，通过上面任意一项的PPG信号质量评估方法获取原始PPG信号的质量评估结果；并在原始PPG信号为第三等级信号时，剔除原始PPG信号；原始PPG信号为第二等级信号时，将原始PPG信号用于部分预设生理指标计算；原始PPG信号为第一等级信号时，将原始PPG信号用于与所述PPG信号相关的所有生理指标计算。即，在PPG信号处理过程中，基于上面的PPG信号质量评估方法，对得到的不同等级的PPG信号进行不同的操作。对第三等级信号直接剔除，对第二等级信号用于心率等部分指标的计算过程，该部分指标为预设的部分指标。对第一等级信号则可以正常的使用，即可以用于心率和血压等所有的生理指标的计算，这里所有生理指标为与PPG相关的生理指标。
另，如图11所示，本发明的一种PPG信号质量评估装置100，包括：
原始信号获取单元110，用于获取预设时长的原始PPG信号，其中，PPG信号为光电容积脉搏波形信号；
第一判断单元121，用于预处理原始PPG信号以得到若干预设长度的波形片段，获取若干波形片段的平均幅度差，判断平均幅度差是否在第一预设范围内，若是，则输出肯定结果，否则输出否定结果；
有用信号获取单元130，用于对原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取原始PPG信号的有用信号；
第二判断单元122，用于基于有用信号和高通滤波结果获取原始PPG信号的高频信噪比，基于有用信号和低通滤波结果获取原始PPG信号的低频信噪比，判断高频信噪比和低频信噪比是否分别大于其对应的第一门限值，若是，则输出肯定结果，否则输出否定结果；
第三判断单元123，用于获取有用信号的自相关函数，对自相关函数进行傅里叶变换以得到有用信号的功率谱，分别获取功率谱的峰度和偏度，判断功率谱的峰度和偏度是否分别大于其对应的第二门限值，若是，则输出肯定结果，否则输出否定结果；
第四判断单元124，用于获取自相关函数的平均峰值周期，判断平均峰值周期是否在第二预设范围内，若是，则输出肯定结果，否则输出否定结果；
第五判断单元125，用于获取自相关函数的最大峰值和峰值周期标准差，确认是否出现最大峰值大于第三门限值或峰值周期标准差小于第四门限值，若是，则输出肯定结果，否则输出否定结果；
第六判断单元126，用于获取自相关函数的峰峰差值标准差，确认峰峰差值标准差是否小于第五门限值，若是，则输出肯定结果，否则输出否定结果；
结果确认单元140，用于在第一判断单元121、第二判断单元122、第三判断单元123、第四判断单元124和第五判断单元125中任意一个输出否定结果时，判定原始PPG信号为第三等级信号，在第六判断单元126输出否定结果时，判定原始PPG信号为第二等级信号，在第六判断单元127输出肯定结果时，判定原始PPG信号为第一等级信号，其中，第三等级信号质量低于第二等级信号质量，第二等级信号质量低于第一等级信号质量。具体的，这里的PPG信号质量评估装置各单元之间具体的配合操作过程具体可以参照上述PPG信号质量评估方法，这里不再赘述。
另，如图12所示，本发明的一种PPG信号处理系统，包括如上面的PPG信号质量评估装置100，以及处理单元200；处理单元用于获取原始PPG信号的质量评估结果；并在原始PPG信号为第三等级信号时，剔除原始PPG信号；原始PPG信号为第二等级信号时，将原始PPG信号用于部分预设生理指标计算；原始PPG信号为第一等级信号时，将原始PPG信号用于与PPG信号相关的所有生理指标计算。即， PPG信号处理系统处理过程中，通过上面的PPG信号质量评估装置得到不同的等级的PPG信号，并通过处理单元对得到的不同等级的PPG信号进行不同的操作。对第三等级信号直接剔除，对第二等级信号用于心率等部分指标的计算过程，该部分指标为预设的部分指标。对第一等级信号则可以正常的使用，即可以用于心率和血压等所有的生理指标的计算，这里所有生理指标为与PPG相关的生理指标。
可以理解的，以上实施例仅表达了本发明的优选实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制；应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，可以对上述技术特点进行自由组合，还可以做出若干变形和改进，这些都属于本发明的保护范围；因此，凡跟本发明权利要求范围所做的等同变换与修饰，均应属于本发明权利要求的涵盖范围。
Claims (13)
1.一种PPG信号质量评估方法，其特征在于，包括步骤：
S1、获取预设时长的原始PPG信号，其中，所述PPG信号为光电容积脉搏波形信号；
S2、预处理所述原始PPG信号以得到若干预设长度的波形片段，获取若干所述波形片段的平均幅度差，判断所述平均幅度差是否在第一预设范围内，若是，则执行步骤S3，否则执行步骤S11；
S3、对所述原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取所述原始PPG信号的有用信号；
S4、基于所述有用信号和所述高通滤波结果获取所述原始PPG信号的高频信噪比，基于所述有用信号和所述低通滤波结果获取所述原始PPG信号的低频信噪比，判断所述高频信噪比和所述低频信噪比是否分别大于其对应的第一门限值，若是，执行步骤S5，否则执行步骤S11；
S5、获取所述有用信号的自相关函数，对所述自相关函数进行傅里叶变换以得到所述有用信号的功率谱，分别获取所述功率谱的峰度和偏度，判断所述功率谱的峰度和偏度是否分别大于其对应的第二门限值，若是，执行步骤S6，否则执行步骤S11；
S6、获取所述自相关函数的平均峰值周期，确认所述平均峰值周期是否在第二预设范围内，若是，则执行步骤S7，否则执行步骤S11；
S7、获取所述自相关函数的最大峰值和峰值周期标准差，确认是否出现所述最大峰值大于第三门限值或所述峰值周期标准差小于第四门限值，若是，则执行步骤S8，否则执行步骤S11；
S8、获取所述自相关函数的峰峰差值标准差，确认所述峰峰差值标准差是否小于第五门限值，若是，则执行步骤S9，否则，执行步骤S10；
S9、判定所述原始PPG信号为第一等级信号，并结束；
S10、判定所述原始PPG信号为第二等级信号，所述第二等级信号质量低于所述第一等级信号质量，并结束；
S11、判定所述原始PPG信号为第三等级信号，所述第三等级信号质量低于所述第二等级信号质量，并结束。
2.根据权利要求1所述的PPG信号质量评估方法，其特征在于，在所述步骤S1中，所述预设时长为大于或等于5s小于或等于15s。
3.根据权利要求1所述的PPG信号质量评估方法，其特征在于，在所述步骤S2中，
所述预处理所述原始PPG信号以得到若干预设长度的波形片段，包括：
对所述原始PPG信号进行中值滤波，对滤波后的信号以所述预设长度的滑动窗口依次进行无间隔截取以获取若干所述波形片段，其中所述波形片段之间不重叠；和/或，
所述获取若干所述 波形片段的平均幅度差，包括：
获取每一波形片段的幅度差，并获取所有波形片段的幅度差的平均值为所述平均幅度差。
4.根据权利要求1所述的PPG信号质量评估方法，其特征在于，
所述第一预设范围大于或等于150且小于或等于90000；和/或，
所述预设长度为1s。
5.根据权利要求1所述的PPG信号质量评估方法，其特征在于，在所述步骤S3中，所述对所述原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取所述原始PPG信号的有用信号；包括：
S31、以截止频率为第一频率的巴特沃斯高通滤波器对所述原始PPG信号进行高通滤波，以得到高频噪声信号；
S32、以截止频率为第二频率的巴特沃斯低通滤波器对所述原始PPG信号进行低通滤波，以得到低频噪声信号；
S33、对所述原始PPG信号分别剔除所述高频噪声信号和所述低频噪声信号，以剩余信号为所述有用信号。
6.根据权利要求5所述的PPG信号质量评估方法，其特征在于，在所述步骤S4中，
所述基于所述有用信号和所述高通滤波结果获取所述原始PPG信号的高频信噪比，包括：获取所述高频噪声信号与所述有用信号的比值为所述高频信噪比；
所述基于所述有用信号和所述低通滤波结果获取所述原始PPG信号的低频信噪比，包括：获取所述低频噪声信号与所述有用信号的比值为所述低频信噪比。
7.根据权利要求1所述的PPG信号质量评估方法，其特征在于，在所述步骤S4中，所述高频信噪比对应的第一门限值大于或等于6dB，所述低频信噪比对应的第一门限值大于或等于-15dB；和/或，
在所述步骤S5中，所述功率谱的偏度对应的第二门限值大于或等于12，所述功率谱的峰度对应的第二门限值大于或等于150。
8.根据权利要求1所述的PPG信号质量评估方法，其特征在于，在所述步骤S5中，所述获取所述有用信号的自相关函数，包括：对所述有用信号做单边归一化的自相关函数，其中所述自相关函数满足以下公式：
其中，
为所述自相关函数，k为所述自相关函数的索引值，
为所述有用信号，n为所述有用信号的索引值，N为所述有用信号的点数。
9.根据权利要求1所述的PPG信号质量评估方法，其特征在于，所述方法还包括：
基于二阶导数的波谷检测算法对自相关函数进行峰值检测以获取若干峰值坐标；
基于所述若干峰值坐标分别获取所述自相关函数的平均峰值周期、最大峰值、峰值周期标准差和峰峰差值标准差。
10.根据权利要求1所述的PPG信号质量评估方法，其特征在于，所述方法包括以下参数设置中的一个或多个：
所述第二预设范围为大于0.3s且小于2s；
所述第三门限值大于或等于0.6；
所述第四门限值小于或等于5；
所述第五门限值小于或等于0.07。
11.一种PPG信号处理方法，其特征在于，包括，通过权利要求1至10任意一项所述的PPG信号质量评估方法获取原始PPG信号的质量评估结果；并在所述原始PPG信号为第三等级信号时，剔除所述原始PPG信号；
所述原始PPG信号为第二等级信号时，将所述原始PPG信号用于部分预设生理指标计算；
所述原始PPG信号为第一等级信号时，将原始PPG信号用于与所述PPG信号相关的所有生理指标计算。
12.一种PPG信号质量评估装置，其特征在于，包括：
原始信号获取单元，用于获取预设时长的原始PPG信号，其中，所述PPG信号为光电容积脉搏波形信号；
第一判断单元，用于预处理所述原始PPG信号以得到若干预设长度的波形片段，获取若干所述波形片段的平均幅度差，判断所述平均幅度差是否在第一预设范围内，若是，则输出肯定结果，否则输出否定结果；
有用信号获取单元，用于对所述原始PPG信号分别进行高通滤波和低通滤波，基于高通滤波结果和低通滤波结果获取所述原始PPG信号的有用信号；
第二判断单元，用于基于所述有用信号和所述高通滤波结果获取所述原始PPG信号的高频信噪比，基于所述有用信号和所述低通滤波结果获取所述原始PPG信号的低频信噪比，判断所述高频信噪比和所述低频信噪比是否分别大于其对应的第一门限值，若是，则输出肯定结果，否则输出否定结果；
第三判断单元，用于获取所述有用信号的自相关函数，对所述自相关函数进行傅里叶变换以得到所述有用信号的功率谱，分别获取所述功率谱的峰度和偏度，判断所述功率谱的峰度和偏度是否分别大于其对应的第二门限值，若是，则输出肯定结果，否则输出否定结果；
第四判断单元，用于获取所述自相关函数的平均峰值周期，判断所述平均峰值周期是否在第二预设范围内，若是，则输出肯定结果，否则输出否定结果；
第五判断单元，用于获取所述自相关函数的最大峰值和峰值周期标准差，确认是否出现所述最大峰值大于第三门限值或所述峰值周期标准差小于第四门限值，若是，则输出肯定结果，否则输出否定结果；
第六判断单元，用于获取所述自相关函数的峰峰差值标准差，确认所述峰峰差值标准差是否小于第五门限值，若是，则输出肯定结果，否则输出否定结果；
结果确认单元，用于在所述第一判断单元、所述第二判断单元、所述第三判断单元、所述第四判断单元和所述第六判断单元中任意一个输出否定结果时，判定所述原始PPG信号为第三等级信号，在所述第六判断单元输出否定结果时，判定所述原始PPG信号为第二等级信号，在所述第六判断单元输出肯定结果时，判定所述原始PPG信号为第一等级信号，其中，所述第三等级信号质量低于所述第二等级信号质量，所述第二等级信号质量低于所述第一等级信号质量。
13.一种PPG信号处理系统，其特征在于，包括如权利要求12所述的PPG信号质量评估装置，以及处理单元；
所述处理单元用于获取原始PPG信号的质量评估结果；并在
所述原始PPG信号为第三等级信号时，剔除所述原始PPG信号；
所述原始PPG信号为第二等级信号时，将所述原始PPG信号用于部分预设生理指标计算；
所述原始PPG信号为第一等级信号时，将原始PPG信号用于与所述PPG信号相关的所有生理指标计算。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210585033.1A CN114676738B (zh) | 2022-05-27 | 2022-05-27 | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 |
| PCT/CN2022/115490 WO2023226223A1 (zh) | 2022-05-27 | 2022-08-29 | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210585033.1A CN114676738B (zh) | 2022-05-27 | 2022-05-27 | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114676738A CN114676738A (zh) | 2022-06-28 |
| CN114676738B true CN114676738B (zh) | 2022-08-02 |
Family
ID=82079397
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210585033.1A Active CN114676738B (zh) | 2022-05-27 | 2022-05-27 | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN114676738B (zh) |
| WO (1) | WO2023226223A1 (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114676738B (zh) * | 2022-05-27 | 2022-08-02 | 广东玖智科技有限公司 | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 |
| CN115670397B (zh) * | 2022-11-17 | 2023-06-02 | 北京中科心研科技有限公司 | 一种ppg伪迹识别方法、装置、存储介质及电子设备 |
| CN118101088B (zh) * | 2024-04-23 | 2024-07-19 | 北京邮电大学 | 一种针对非完全重叠频谱的多信号检测方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103020472B (zh) * | 2012-12-27 | 2015-12-09 | 中国科学院深圳先进技术研究院 | 基于约束估计的生理信号质量评估方法和系统 |
| JP6811773B2 (ja) * | 2015-11-24 | 2021-01-13 | コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． | フォトプレチスモグラム（ｐｐｇ）信号品質を定量化する方法 |
| CN109598180A (zh) * | 2017-09-30 | 2019-04-09 | 深圳市岩尚科技有限公司 | 光电容积脉搏波的质量评估方法 |
| EP3479763B1 (en) * | 2017-11-06 | 2023-03-01 | Tata Consultancy Services Limited | System and method for photoplethysmogram (ppg) signal quality assessment |
| CN108523877B (zh) * | 2018-03-23 | 2021-02-09 | 南京中医药大学 | 一种心电信号质量辨识方法及其心电分析方法 |
| CN111493821B (zh) * | 2020-03-16 | 2023-03-14 | 浙江好络维医疗技术有限公司 | 一种基于modwt及中值滤波的ppg信号实时去噪方法 |
| CN112494001B (zh) * | 2020-11-25 | 2023-05-26 | 青岛歌尔智能传感器有限公司 | Ppg信号质量评估方法、装置、穿戴式设备及存储介质 |
| CN112545472B (zh) * | 2020-12-02 | 2022-11-01 | 成都心吉康科技有限公司 | 一种ppg信号质量评估的方法、装置、设备及储存介质 |
| CN114528888B (zh) * | 2022-04-25 | 2022-07-12 | 广东玖智科技有限公司 | 一种ppg信号聚类中心获取方法及装置和ppg信号处理方法及装置 |
| CN114676738B (zh) * | 2022-05-27 | 2022-08-02 | 广东玖智科技有限公司 | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 |
-
2022
- 2022-05-27 CN CN202210585033.1A patent/CN114676738B/zh active Active
- 2022-08-29 WO PCT/CN2022/115490 patent/WO2023226223A1/zh unknown
Also Published As
| Publication number | Publication date |
|---|---|
| CN114676738A (zh) | 2022-06-28 |
| WO2023226223A1 (zh) | 2023-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114676738B (zh) | Ppg信号质量评估方法及装置以及ppg信号处理方法及系统 | |
| Varghees et al. | A novel heart sound activity detection framework for automated heart sound analysis | |
| Vadrevu et al. | A robust pulse onset and peak detection method for automated PPG signal analysis system | |
| CN102319064B (zh) | 一种提高胎心率数据减速识别准确性的装置和方法 | |
| US20060094943A1 (en) | Use of time indexed plethysmographic spectral data in assessing saturation estimation validity | |
| CN111407315B (zh) | 一种针对超声多普勒胎心信号的质量评估方法 | |
| CN108888259A (zh) | 一种心电信号实时qrs波检测方法 | |
| CN104305992B (zh) | 一种交互式胎儿心电快速自动提取方法 | |
| CN112998690B (zh) | 一种基于脉搏波多特征融合的呼吸率提取方法 | |
| Shaik et al. | A method for QRS delineation based on STFT using adaptive threshold | |
| CN110974189A (zh) | 脉搏波的信号质量检测方法及装置、设备、系统 | |
| CN109350031A (zh) | 一种光电容积脉搏波信号特征点检测方法及装置 | |
| TW201438669A (zh) | 血氧飽和度檢測方法及裝置 | |
| Azami et al. | Automatic signal segmentation using the fractal dimension and weighted moving average filter | |
| CN112656392B (zh) | 一种低功耗ppg心率计算方法及电子设备 | |
| Lin et al. | A characteristic filtering method for pulse wave signal quality assessment | |
| CN115486849A (zh) | 心电信号质量评估方法、装置及设备 | |
| CN103211586B (zh) | 基于光学检测的无创压信号采集的方法和装置 | |
| CN111134634B (zh) | 一种基于聚类分析的光电容积脉搏波分析处理方法 | |
| CN106667474A (zh) | 一种适用于智能穿戴设备的心率测量方法及装置 | |
| Vadrevu et al. | Effective systolic peak detection algorithm using variational mode decomposition and center of gravity | |
| Haibing et al. | A method of continuous wavelet transform for qrs wave detection in ecg signal | |
| Sapoznikov et al. | Detection of regularities in heart rate variations by linear and non-linear analysis: power spectrum versus approximate entropy | |
| CN110432885A (zh) | 一种光电容积脉搏波噪声去除方法 | |
| Fedotov et al. | A QRS-complex detector of the electrocardiogram signal for the long-term monitoring of the patient’s condition |
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 | ||
| CP03 | Change of name, title or address | ||
| CP03 | Change of name, title or address |
Address after: 518000, Building 3, 1301, Sunshine Yuehai Garden Phase II, No. 3818 Baishi Road, Yuegui Community, Yuehai Street, Nanshan District, Shenzhen, Guangdong Province Patentee after: Shenzhen Jiuzhi Technology Co.,Ltd. Country or region after: China Address before: 528400 room 411-18, floor 4, building C, innovation center plant, No. 34, Xiangshan Avenue, Cuiheng new area, Zhongshan City, Guangdong Province Patentee before: Guangdong Jiuzhi Technology Co.,Ltd. Country or region before: China |