CN111297339B - 一种光体积变化描记图法信号的生成方法和装置 - Google Patents
一种光体积变化描记图法信号的生成方法和装置 Download PDFInfo
- Publication number
- CN111297339B CN111297339B CN202010110280.7A CN202010110280A CN111297339B CN 111297339 B CN111297339 B CN 111297339B CN 202010110280 A CN202010110280 A CN 202010110280A CN 111297339 B CN111297339 B CN 111297339B
- Authority
- CN
- China
- Prior art keywords
- signal
- red
- digital signal
- identifier
- red light
- 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
- 238000013186 photoplethysmography Methods 0.000 title claims abstract description 114
- 238000000034 method Methods 0.000 title claims abstract description 89
- 238000001914 filtration Methods 0.000 claims abstract description 132
- 238000012545 processing Methods 0.000 claims abstract description 92
- 238000004364 calculation method Methods 0.000 claims abstract description 23
- 238000001514 detection method Methods 0.000 claims description 56
- 230000008569 process Effects 0.000 claims description 28
- 230000001360 synchronised effect Effects 0.000 claims description 22
- 238000000605 extraction Methods 0.000 claims description 15
- 238000005286 illumination Methods 0.000 claims description 12
- 238000007689 inspection Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 8
- 238000000926 separation method Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000012163 sequencing technique Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 description 13
- 230000003287 optical effect Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 230000009467 reduction Effects 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 4
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000010223 real-time analysis Methods 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000032683 aging Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000013500 data storage Methods 0.000 description 2
- 210000000624 ear auricle Anatomy 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 210000001061 forehead Anatomy 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002996 emotional effect Effects 0.000 description 1
- 230000003511 endothelial effect Effects 0.000 description 1
- 210000003811 finger Anatomy 0.000 description 1
- 230000033083 heart process Effects 0.000 description 1
- 229940088597 hormone Drugs 0.000 description 1
- 239000005556 hormone Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000007619 statistical method 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/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
-
- 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/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0077—Devices for viewing the surface of the body, e.g. camera, magnifying lens
-
- 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
Abstract
本发明实施例涉及一种光体积变化描记图法信号的生成方法和装置,所述方法包括:根据预置的红光像素阈值范围对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧红、绿色通道数据计算生成第一红、绿光数字信号;对第一红、绿光数字信号进行信号带通滤波生成第二红、绿光数字信号;获取信号判断标识符;根据信号判断标识符判断是否对第二红光数字信号和第二绿光数字信号进行信号最大频差判断处理并生成第一判断结果;当第一判断结果为达标信号标识时对第二红和或绿光数字信号进行信噪比判断生成第二判断结果;当第二判断结果为达标信号标识时生成PPG信号。
Description
技术领域
本发明涉及电生理信号处理技术领域,特别涉及一种光体积变化描记图法信号的生成方法和装置。
背景技术
人体当前的健康状态、激素水平、情绪状态乃至生活方式等因素都会影响心脏的工作状态,反之,通过对心脏运动的轨迹进行分析也可以反向获悉人体的健康与疾病状态信息。要获得心脏运动的轨迹数据,就需要对心脏运动进行监测和数据采集,传统的方法是通过使用心电图(electrocardiogram,ECG)监测技术来完成的。该方式需要将导联电极连接到身体多个部位来进行心电信号的采集,具体操作时,被监测对象需要静卧或者再日常生活中时刻佩戴多种接触式设备(心贴、心带等等),对监测对象的行动与生活造成诸多不便。
光体积变化描记图法(Photoplethysmography,PPG),是借助光电手段在活体组织中检测血液容积变化的一种无创检测方法。在心脏搏动时,对血管内单位面积的血流量形成周期性变化,与之对应的血液体积也相应发生变化,从而导致反映血液吸收光量的PPG信号也呈现周期性变化趋势。PPG信号可以通过手指、耳朵等区域测量得到。通过PPG信号同样可以对心脏活动做一系列定向分析,且PPG方式相对ECG方式可以提高监测对象的体验舒适度。但在实际应用中,我们发现,当前的PPG信号采集装置和方式比较刻板:照射光源只采用单一的红光或者红外光居多,采集前端常见最多的是指夹、耳夹之类的定制机械设备,采集的数据需要经过特定传感器进行分析等等,诸如此类的问题都加大了上层医疗应用对PPG信号数据的获取难度。
发明内容
本发明的目的,就是针对现有技术的缺陷,提供一种光体积变化描记图法信号的生成方法和装置,不规定特定照射光源类型采用普通灯光即可,简化了PPG信号前端采集装置使用一般的移动终端也可实现;通过对录取的视频数据进行分析并生成PPG信号,即能满足本地设备实时分析与数据备份的要求,又能满足远程传输与分析的需求;对有效信号的识别过程不依赖特定传感器,降低了应用层面对PPG数据的获取难度;通过设置信号判断标识符,提高对多种拍摄设备的兼容性。
为实现上述目的,本发明实施例第一方面提供了一种光体积变化描记图法信号的生成方法,其特征在于,所述方法包括:
根据预置的红光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧红色通道数据计算,生成第一红光数字信号;并根据预置的绿光像素阈值范围,对所述皮表视频数据帧图像序列的所有所述皮表视频数据帧图像依次进行单帧绿色通道数据计算,生成第一绿光数字信号;所述皮表视频数据帧图像序列包括多个所述皮表视频数据帧图像;
根据预置的带通滤波频率阈值范围对所述第一红光数字信号进行信号带通滤波预处理生成第二红光数字信号,并根据所述带通滤波频率阈值范围对所述第一绿光数字信号进行信号带通滤波预处理生成第二绿光数字信号;
获取预置的信号判断标识符;所述信号判断标识符包括红绿光同检标识和红光单检标识两种标识符;
根据所述信号判断标识符判断是否对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理,当所述信号判断标识符为所述红绿光同检标识时对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理生成第一判断结果,当所述信号判断标识符为所述红光单检标识时直接设置所述第一判断结果为达标信号标识;
当所述第一判断结果为所述达标信号标识时,根据所述信号判断标识符对所述第二红光数字信号和或所述第二绿光数字信号进行信号信噪比判断处理,生成第二判断结果;
当所述第二判断结果为所述达标信号标识时,根据所述信号判断标识符、所述第二红光数字信号和或所述第二绿光数字信号进行光体积变化描记图法PPG信号生成处理,生成PPG信号。
优选的,所述皮表视频数据帧图像序列是,调用视频拍摄设备以预置的采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据;对所述皮表视频数据进行帧图像提取处理生成的。
进一步的,所述调用视频拍摄设备以预置的采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据,具体包括:
调用所述视频拍摄设备的灯光照明模块对所述生物体局部皮肤表面提供照明;所述视频拍摄设备包括灯光照明模块和视频拍摄模块;所述视频拍摄模块至少包括镜头模块;
当所述灯光照明模块与所述镜头模块的物理间距小于或等于预置的设备间距阈值时,将所述生物体局部皮肤表面同时覆盖所述灯光照明模块与所述镜头模块;当所述灯光照明模块与所述镜头模块的物理间距大于所述设备间距阈值时,将所述生物体局部皮肤表面覆盖所述镜头模块;
调用所述视频拍摄模块以所述采集时间阈值为拍摄时间长度对所述生物体局部皮肤表面进行所述连续采集拍摄操作生成所述皮表视频数据。
进一步的,所述方法还包括:
在调用所述视频拍摄设备以所述采集时间阈值为拍摄时间长度对所述生物体局部皮肤表面进行所述连续采集拍摄操作的过程中,以预置的缓存片段时间阈值为片段提取长度对所述视频拍摄设备采集拍摄的视频数据进行连续视频片段提取处理生成多个缓存片段;
提取所述缓存片段的帧图像并将提取出的所有帧图像按时间顺序进行排序,生成缓存片段帧图像序列;所述缓存片段包括多个所述帧图像;
对所述缓存片段帧图像序列的所有缓存片段帧图像,根据所述红光像素阈值范围进行红色像素点占比统计生成红色占比参数,并根据预置的杂色像素阈值范围进行杂色像素点占比统计生成杂色占比参数;所述缓存片段帧图像序列包括多个所述缓存片段帧图像;
根据所述红色占比参数对所述缓存片段进行皮肤表面脱离拍摄错误识别判断,当所述红色占比参数小于预置的红色占比下限阈值时停止所述连续采集拍摄操作并生成皮肤表面脱离拍摄错误信息;
根据所述杂色占比参数对所述缓存片段进行皮肤表面脱离拍摄错误识别判断,当所述杂色占比参数大于预置的杂色占比上限阈值时停止所述连续采集拍摄操作并生成所述皮肤表面脱离拍摄错误信息。
优选的,
所述第一判断结果包括所述达标信号标识和不达标信号标识两种判断结果;
所述第二判断结果包括所述达标信号标识和所述不达标信号标识两种判断结果。
优选的,所述根据预置的红光像素阈值范围,对所述皮表视频数据帧图像序列的所有所述皮表视频数据帧图像依次进行单帧红色通道数据计算,生成第一红光数字信号;并根据预置的绿光像素阈值范围,对所述皮表视频数据帧图像序列的所有所述皮表视频数据帧图像依次进行单帧绿色通道数据计算,生成第一绿光数字信号,具体包括:
步骤61,初始化临时第一红光数字信号为空,初始化临时第一绿光数字信号为空,初始化第一索引的值为1,初始化第一总数为所述皮表视频数据帧图像序列的皮表视频数据帧图像总数;
步骤62,从所述皮表视频数据帧图像序列中提取与所述第一索引对应的所述皮表视频数据帧图像生成第一索引皮表视频数据帧图像;
步骤63,根据所述红光像素阈值范围,在所述第一索引皮表视频数据帧图像中统计所有满足所述红光像素阈值范围的像素点生成红色像素点集合,并对所述红色像素点集合的所有像素点的像素值按指定加权平均方式进行加权平均计算,生成第一索引帧红色通道数据;将所述第一索引帧红色通道数据作为红光信号点数据向所述临时第一红光数字信号进行信号点添加操作;
步骤64,根据所述绿光像素阈值范围,在所述第一索引皮表视频数据帧图像中统计所有满足所述绿光像素阈值范围的像素点生成绿色像素点集合,并对所述绿色像素点集合的所有像素点的像素值按指定加权平均方式进行加权平均计算,生成第一索引帧绿色通道数据;将所述第一索引帧绿色通道数据作为绿光信号点数据向所述临时第一绿光数字信号进行信号点添加操作;
步骤65,将所述第一索引的值加1;
步骤66,判断所述第一索引是否大于所述第一总数,如果所述第一索引小于或等于所述第一总数则转至步骤62,如果所述第一索引大于所述第一总数则转至步骤67;
步骤67,设置所述第一红光数字信号为所述临时第一红光数字信号;设置所述第一绿光数字信号为所述临时第一绿光数字信号。
优选的,所述根据预置的带通滤波频率阈值范围对所述第一红光数字信号进行信号带通滤波预处理生成第二红光数字信号,并根据所述带通滤波频率阈值范围对所述第一绿光数字信号进行信号带通滤波预处理生成第二绿光数字信号,具体包括:
根据所述带通滤波频率阈值范围,在所述第一红光数字信号中,将信号频率低于所述带通滤波频率阈值范围的低频噪声信号点和信号频率高于所述带通滤波频率阈值范围的高频噪声信号点从所述第一红光数字信号中进行数字信号滤波处理,生成所述第二红光数字信号;
根据所述带通滤波频率阈值范围,在所述第一绿光数字信号中,将信号频率低于所述带通滤波频率阈值范围的低频噪声信号点和信号频率高于所述带通滤波频率阈值范围的高频噪声信号点从所述第一绿光数字信号中进行数字信号滤波处理,生成所述第二绿光数字信号。
优选的,所述根据所述信号判断标识符判断是否对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理,当所述信号判断标识符为所述红绿光同检标识时对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理生成第一判断结果,当所述信号判断标识符为所述红光单检标识时直接设置所述第一判断结果为达标信号标识,具体包括:
步骤81,初始化临时第一判断结果为空;
步骤82,判断所述信号判断标识符是否为所述红绿光同检标识,当所述信号判断标识符为所述红绿光同检标识时转至步骤83,当所述信号判断标识符为所述红光单检标识时转至步骤85;
步骤83,对所述第二红光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成红光频域信号,对所述第二绿光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成绿光频域信号;从所述红光频域信号中提取能量最高频率生成红光最大频率,从所述绿光频域信号中提取能量最高频率生成绿光最大频率;根据所述红光最大频率与所述绿光最大频率的频率差生成红绿最大频差;
步骤84,根据预置的最大频差阈值范围对所述红绿最大频差进行信号最大频差判断处理,当所述红绿最大频差未超过所述最大频差阈值范围时设置所述临时第一判断结果为所述达标信号标识,当所述红绿最大频差超过所述最大频差阈值范围时设置所述临时第一判断结果为所述不达标信号标识;转至步骤86;
步骤85,设置所述临时第一判断结果为所述达标信号标识;
步骤86,设置所述第一判断结果为所述临时第一判断结果。
优选的,所述当所述第一判断结果为所述达标信号标识时,根据所述信号判断标识符对所述第二红光数字信号和或所述第二绿光数字信号进行信号信噪比判断处理,生成第二判断结果,具体包括:
步骤91,初始化临时第二判断结果为空;
步骤91,当所述第一判断结果为所述达标信号标识时,根据所述信号判断标识符进行判断,当所述信号判断标识符为所述红绿光同检标识时转至步骤92,当所述信号判断标识符为所述红光单检标识时转至步骤95;
步骤92,根据预置的带阻滤波频率阈值范围,对所述第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足所述带阻滤波频率阈值范围的有效信号从所述第二红光数字信号中滤除生成红光噪声信号;并根据所述带阻滤波频率阈值范围,对所述第二绿光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足所述带阻滤波频率阈值范围的有效信号从所述第二绿光数字信号中滤除生成绿光噪声信号;
步骤93,计算所述第二红光数字信号的信号能量生成红光信号能量,计算所述红光噪声信号的信号能量生成红光噪声能量,根据所述红光信号能量减去所述红光噪声能量的差生成有效红光信号能量,根据所述有效红光信号能量与所述红光噪声能量的比值生成红光信噪比;计算所述第二绿光数字信号的信号能量生成绿光信号能量,计算所述绿光噪声信号的信号能量生成绿光噪声能量,根据所述绿光信号能量减去所述绿光噪声能量的差生成有效绿光信号能量,根据所述有效绿光信号能量与所述绿光噪声能量的比值生成绿光信噪比;
步骤94,如果所述红光信噪比与所述绿光信噪比均小于预置的信噪比阈值则设置所述临时第二判断结果为所述不达标信号标识;如果所述红光信噪比与所述绿光信噪比中任一个大于或等于所述信噪比阈值则设置所述临时第二判断结果为所述达标信号标识;转至步骤97;
步骤95,根据所述带阻滤波频率阈值范围,对所述第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足所述带阻滤波频率阈值范围的有效信号从所述第二红光数字信号中滤除生成所述红光噪声信号;计算所述第二红光数字信号的信号能量生成所述红光信号能量,计算所述红光噪声信号的信号能量生成所述红光噪声能量,根据所述红光信号能量减去所述红光噪声能量的差生成所述有效红光信号能量,根据所述有效红光信号能量与所述红光噪声能量的比值生成所述红光信噪比;
步骤96,如果所述红光信噪比小于所述信噪比阈值则设置所述临时第二判断结果为所述不达标信号标识;如果所述红光信噪比大于或等于所述信噪比阈值则设置所述临时第二判断结果为所述达标信号标识;
步骤97,设置所述第二判断结果为所述临时第二判断结果。
优选的,所述当所述第二判断结果为所述达标信号标识时,根据所述信号判断标识符、所述第二红光数字信号和或所述第二绿光数字信号进行光体积变化描记图法PPG信号生成处理,生成PPG信号,具体包括:
当所述第二判断结果为所述达标信号标识时,设置所述PPG信号为空;初始化所述PPG信号的通道标识符为空,初始化所述PPG信号的红光数字信号为空,初始化所述PPG信号的绿光数字信号为空;
当所述信号判断标识符为所述红绿光同检标识时,设置所述PPG信号的所述通道标识符为红绿通道标识,设置所述PPG信号的所述红光数字信号为所述第二红光数字信号,设置所述PPG信号的所述绿光数字信号为所述第二绿光数字信号;
当所述信号判断标识符为所述红光单检标识时,设置所述通道标识符为红通道标识,设置所述PPG信号的所述红光数字信号为所述第二红光数字信号。
优选的,所述方法还包括:
当所述第一判断结果为所述不达标信号标识时,停止PPG信号生成处理流程,并生成信号质量错误信息;
当所述第二判断结果为所述不达标信号标识时,停止PPG信号生成处理流程,并生成信号质量错误信息。
本发明实施例第一方面提供的一种光体积变化描记图法信号的生成方法,不规定特定照射光源类型采用普通灯光即可,简化了PPG信号前端采集装置使用一般的移动终端也可实现;通过对录取的视频数据进行分析并生成PPG信号,即能满足本地设备实时分析与数据备份的要求,又能满足远程传输与分析的需求;对有效信号的识别过程不依赖特定传感器,降低了应用层面对PPG数据的获取难度;通过设置信号判断标识符,提高对多种拍摄设备的兼容性。
本发明实施例第二方面提供了一种设备,该设备包括存储器和处理器,存储器用于存储程序,处理器用于执行第一方面及第一方面的各实现方式中的方法。
本发明实施例第三方面提供了一种包含指令的计算机程序产品,当计算机程序产品在计算机上运行时,使得计算机执行第一方面及第一方面的各实现方式中的方法。
本发明实施例第四方面提供了一种计算机可读存储介质,计算机可读存储介质上存储有计算机程序,计算机程序被处理器执行时实现第一方面及第一方面的各实现方式中的方法。
附图说明
图1为本发明实施例一提供的一种光体积变化描记图法信号的生成方法示意图;
图2为本发明实施例二提供的一种光体积变化描记图法信号的生成方法示意图;
图3为本发明实施例三提供的一种光体积变化描记图法信号的生成装置的设备结构示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部份实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
在进行实例描述之前,首先对本发明实施例的实现技术进行简要描述如下。
常规PPG信号采集,是采用特定光源对皮肤表面进行照射,然后使用特定光源接收传感器对由该皮肤表面反射或者透射的光信号进行接收;该传感器的分析模块会按时间顺序对光强变化进行统计分析最终以归一化的PPG信号方式对动态脉搏波进行呈现。
参考PPG信号分析原理,在连续的心动过程中,我们又发现测试者皮肤表面对自然光中的红光与绿光的光能吸收与常规PPG采集设备的特定光源的光能吸收趋势是一样的。那么借助另外的介质,只要能够体现时间周期内皮肤表面红、绿光变化趋势,也就同样可以获得归一化的PPG信号。
基于上述发现,本发明实施例相对常规PPG信号采集设备,不设置特定光源要求,不使用特定光源对皮肤表明进行照射并使用特定传感器对光强进行采集。本发明实施例采用自然光(即使使用闪光灯之类的强光源只是为了提高光照强度而非设定光源性质),使用常规的拍摄装置(比如移动终端的摄像头之类常见的摄像装置),对测试者的指定皮肤表面(手指头、耳垂、脖子、额头等等)进行拍摄,并将获得的视频数据进行帧图像提取就生成一个帧图像的序列。本发明实施例就是以这些连续的帧图像作为另外的介质对红、绿光的变化趋势进行捕捉。具体实现方法是,将帧图像中的红、绿色像素值进行加权统计生成该帧图像对应时间点的原始红、光信号点,将原始红、绿光信号进行相应的滤波降噪处理生成PPG信号的红、绿数据通道值。本发明实施例提供一种对摄像设备的兼容方式:信号判断标识符(包括红绿光同检标识和红光单检标识两个取值),当摄像设备的镜头对绿光有明显屈光缺陷时,就选取红光单检标识只使用红光形成PPG信号;当摄像设备的镜头没有明显屈光缺陷时,选取红绿光同检标识同时参考红、绿光信号形成PPG信号。
如图1为本发明实施例一提供的一种光体积变化描记图法信号的生成方法示意图所示,本方法主要包括如下步骤:
步骤1,调用视频拍摄设备以预置的采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据;对皮表视频数据进行帧图像提取处理生成皮表视频数据帧图像序列;
其中,皮表视频数据帧图像序列包括多个皮表视频数据帧图像;
具体包括:步骤11,调用视频拍摄设备的灯光照明模块对生物体局部皮肤表面提供照明;
其中,视频拍摄设备包括灯光照明模块和视频拍摄模块;视频拍摄模块至少包括镜头模块;
此处,在日常所见的,只要带有摄像头、具备视频拍摄功能的移动终端均可充当本发明实施例的视频拍摄设备;在拍摄之前,需要打开移动终端的闪光灯并保持常亮,以此满足稳定的照射光源需求;
步骤12,当灯光照明模块与镜头模块的物理间距小于或等于预置的设备间距阈值时,将生物体局部皮肤表面同时覆盖灯光照明模块与镜头模块;当灯光照明模块与镜头模块的物理间距大于设备间距阈值时,将生物体局部皮肤表面覆盖镜头模块;
此处,在拍摄时,需要保持被拍摄生物体的局部皮肤表面轻微按压在移动终端的摄像头上;对于闪光灯和摄像头比较接近的移动终端,需要将皮肤表面同时覆盖住闪光灯和摄像头;对于闪光灯和摄像头比较远的移动终端,只需要将皮肤表面覆盖住摄像头即可;拍摄过程中皮肤表面需保持静止、用力均匀;
步骤13,调用视频拍摄模块以采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据。
此处,采集时间阈值根据实际应用的数据存储能力和分析能力可多次进行预定义。
此处,生物体局部皮肤表面可以是指尖,可以是耳朵的耳垂部分,还可以是额头部位等等,原理上只要具备较好透射或者反射效果的皮肤表面都可适用。
此处,在调用视频拍摄设备以采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作的过程中,实施例一以预置的缓存片段时间阈值为片段提取长度对视频拍摄设备采集拍摄的视频数据进行连续视频片段提取处理生成多个缓存片段;例如,假设采集时间阈值为5分钟、缓存片段时间阈值为1秒,则意味着在这5分钟的拍摄过程中,每完成1秒的拍摄就会将这一秒的视频作为缓存片段进行存储,并同时启动对缓存数据的进一步分析;
进一步的,实施例一提取缓存片段的帧图像并将提取出的所有帧图像按时间顺序进行排序,生成缓存片段帧图像序列;对缓存片段帧图像序列的所有缓存片段帧图像,根据红光像素阈值范围进行红色像素点占比统计生成红色占比参数,并根据预置的杂色像素阈值范围进行杂色像素点占比统计生成杂色占比参数;例如,假设1秒的视频包括24帧,则意味着缓每个个缓存片段帧图像序列包括24帧缓存片段帧图像;此处,假设拍摄过程中皮肤表面是完全覆盖在摄像镜头上的话,则视频中帧图像里大多数都是红色像素点,反之如果皮肤表面未压紧或者已经移开摄像头位置则都可能造成红色像素点减少、杂色像素点增多,因此,需要对每一个缓存片段帧图像进行红色占比参数与杂色占比参数的计算和比对;
最后,对每个缓存片段进行皮肤表面脱离拍摄错误识别判断,当红色占比参数小于预置的红色占比下限阈值时停止连续采集拍摄操作并生成皮肤表面脱离拍摄错误信息,当杂色占比参数大于预置的杂色占比上限阈值时停止连续采集拍摄操作并生成皮肤表面脱离拍摄错误信息;发生该类错误时,基于上位应用的实际应用需求,可根据这个错误信息要求移动终端对用户进行警示提醒、手动重拍或者自动重拍等等操作。
步骤2,根据预置的红光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧红色通道数据计算,生成第一红光数字信号;并根据预置的绿光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧绿色通道数据计算,生成第一绿光数字信号;
具体包括:步骤21,初始化临时第一红光数字信号为空,初始化临时第一绿光数字信号为空,初始化第一索引的值为1,初始化第一总数为皮表视频数据帧图像序列的皮表视频数据帧图像总数;
步骤22,从皮表视频数据帧图像序列中提取与第一索引对应的皮表视频数据帧图像生成第一索引皮表视频数据帧图像;
步骤23,根据红光像素阈值范围,在第一索引皮表视频数据帧图像中统计所有满足红光像素阈值范围的像素点生成红色像素点集合,并对红色像素点集合的所有像素点的像素值按指定加权平均方式进行加权平均计算,生成第一索引帧红色通道数据;将第一索引帧红色通道数据作为红光信号点数据向临时第一红光数字信号进行信号点添加操作;
此处的指定加权平均方式的常规采用总和平均方式,例如,提取第1皮表视频数据帧图像中所有满足红光像素阈值范围(因为不同位置受到光源的照射因内部结构或者血管的反射和透射的程度不一样,会导致光通过率由差异,进一步会导致拍摄下来的视频中红色也会产生颜色深浅的差别,因此采用像素阈值范围)的像素点总数生成红色像素点总数,提取第1皮表视频数据帧图像中所有满足红光像素阈值范围的像素点的像素值进行总和计算生成红色像素值总和,那么第一索引帧红色通道数据就等于红色像素值总和/红色像素点总数;
步骤24,根据绿光像素阈值范围,在第一索引皮表视频数据帧图像中统计所有满足绿光像素阈值范围的像素点生成绿色像素点集合,并对绿色像素点集合的所有像素点的像素值按指定加权平均方式进行加权平均计算,生成第一索引帧绿色通道数据;将第一索引帧绿色通道数据作为绿光信号点数据向临时第一绿光数字信号进行信号点添加操作;
此处的指定加权平均方式的常规采用总和平均方式,那么第一索引帧绿色通道数据就等于绿色像素值总和/绿色像素点总数;
步骤25,将第一索引的值加1;
步骤26,判断第一索引是否大于第一总数,如果第一索引小于或等于第一总数则转至步骤22,如果第一索引大于第一总数则转至步骤27;
步骤27,设置第一红光数字信号为临时第一红光数字信号;设置第一绿光数字信号为临时第一绿光数字信号。
步骤21-27是对步骤2的进一步详细阐述,实施例一通过步骤2将由皮表视频数据转化而来的帧图像序列中的所有帧图像进行多种光源信息的提取操作,本实施例中暂时针对红光和绿光两种光源;对光信号的提取方式,就是通过对帧图像中特定像素的加权平均计算得到一个像素均值,并以此代表该光源在所在帧图像中的颜色通道数据;按时间先后顺序,对视频中的每一帧都做同样的处理,可以得到两组一维数字信号:第一红光数字信号和第一绿光数字信号。
步骤3,根据预置的带通滤波频率阈值范围对第一红光数字信号进行信号带通滤波预处理生成第二红光数字信号,并根据带通滤波频率阈值范围对第一绿光数字信号进行信号带通滤波预处理生成第二绿光数字信号;
具体包括:步骤31,根据带通滤波频率阈值范围,在第一红光数字信号中,将信号频率低于带通滤波频率阈值范围的低频噪声信号点和信号频率高于带通滤波频率阈值范围的高频噪声信号点从第一红光数字信号中进行数字信号滤波处理,生成第二红光数字信号;
步骤32,根据带通滤波频率阈值范围,在第一绿光数字信号中,将信号频率低于带通滤波频率阈值范围的低频噪声信号点和信号频率高于带通滤波频率阈值范围的高频噪声信号点从第一绿光数字信号中进行数字信号滤波处理,生成第二绿光数字信号。
此处,步骤31-32是对步骤3的进一步阐述,步骤3是对通过视频数据提取的两种光源的数字信号进行信号信号预处理,即降噪处理;此处,实施例一使用的降噪手段是带通滤波方式,即预置一个带通滤波频率阈值范围,基于带通滤波原理对低于或高于该频段的信号、干扰和噪声进行信号抑制处理;一般此处的带通滤波频率阈值范围常见的0.5赫兹到10赫兹;在某些移动终端上进行带通滤波处理时,考虑到移动终端的处理能力,使用的是有限长单位冲激响应((Finite lmpulse Response,FIR)滤波模块。
步骤4,获取预置的信号判断标识符;
其中,信号判断标识符包括红绿光同检标识和红光单检标识两种标识符。
此处,之所以要设置信号判断标识符进行选择判断是为了提高对拍摄设备的兼容度,因为在具体应用场景中,有些移动终端因为镜头质量、设备老化、偏光算法落后等多种原因,会导致其对绿光的处理能力较差,通过这些移动终端所采集的视频数据乃至视频的帧图像数据中,基本很难获得足够数量及像素值达标的像素信息,也就无法产生绿光数字信号;在这种情况下,就需要针对不同情况通过设置信号判断标识符进行灵活处理:红绿光同检标识意味着红光、绿光信号都需要做识别操作,红光单检标识则标识只需对红光信号做识别操作。
步骤5,根据信号判断标识符判断是否对第二红光数字信号和第二绿光数字信号进行信号最大频差判断处理,当信号判断标识符为红绿光同检标识时对第二红光数字信号和第二绿光数字信号进行信号最大频差判断处理生成第一判断结果,当信号判断标识符为红光单检标识时直接设置第一判断结果为达标信号标识;
具体包括:步骤51,初始化临时第一判断结果为空;
步骤52,判断信号判断标识符是否为红绿光同检标识,当信号判断标识符为红绿光同检标识时转至步骤53,当信号判断标识符为红光单检标识时转至步骤55;
此处,如果信号判断标识符具体被设置为红绿光同检标识时,说明当前用于拍摄的拍摄设备拍摄对红绿光的处理能力正常,可以转至步骤53通过对第二红光数字信号和第二绿光数字信号进行频率差的关联检测来判断二者的信号质量;如果信号判断标识符具体被设置为红光单检标识时,说明当前用于拍摄的拍摄设备拍摄对绿光的处理能力较差,上述的第二绿光数字信号将会在后续的光信号处理过程中被抛弃,也即无法被引用到最大频差判断处理中,为保证流程能够继续执行,此处会转至步骤55,不对红绿光进行频差比较而直接将第一判断结果设置为达标信号标识;
步骤53,对第二红光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成红光频域信号,对第二绿光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成绿光频域信号;从红光频域信号中提取能量最高频率生成红光最大频率,从绿光频域信号中提取能量最高频率生成绿光最大频率;根据红光最大频率与绿光最大频率的频率差生成红绿最大频差;
步骤54,根据预置的最大频差阈值范围对红绿最大频差进行信号最大频差判断处理,当红绿最大频差未超过最大频差阈值范围时设置临时第一判断结果为达标信号标识,当红绿最大频差超过最大频差阈值范围时设置临时第一判断结果为不达标信号标识;转至步骤56;
此处,首先通过离散傅里叶变换得到第二红色数字信号和第二绿色数字信号的频域信号,通过频域信号得到能量最高的频率(一般这个频率通常对应着心率),此处基本原理是检查这两个数字信号的能量最高的频率是否一致,如果误差在允许误差范围之内,则继续分析,如果误差较大,则认为信号质量不好,需要重新拍摄;
步骤55,设置临时第一判断结果为达标信号标识;
此处,如前文所述,部分移动终端型号的皮表视频数据中绿色数字信号质量较差,为兼容多种移动终端,设置信号判断标识符为红光单检标识,将频差检测信号质量的处理流程直接跳过;
步骤56,设置第一判断结果为临时第一判断结果。
步骤6,当第一判断结果为达标信号标识时,根据信号判断标识符对第二红光数字信号和或第二绿光数字信号进行信号信噪比判断处理,生成第二判断结果;
具体包括:步骤61,初始化临时第二判断结果为空;
步骤61,当第一判断结果为达标信号标识时,根据信号判断标识符进行判断,当信号判断标识符为红绿光同检标识时转至步骤62,当信号判断标识符为红光单检标识时转至步骤65;
此处,如果信号判断标识符具体被设置为红绿光同检标识时,说明当前用于拍摄的拍摄设备拍摄对红绿光的处理能力正常,可以转至步骤62使用第二红光数字信号和第二绿光数字信号的信噪比与预定阈值进行比对进一步判断二者的信号质量;如果信号判断标识符具体被设置为红光单检标识时,说明当前用于拍摄的拍摄设备拍摄对绿光的处理能力较差,上述的第二绿光数字信号将会在后续的光信号处理过程中被抛弃,此处会转至步骤65只使用第二红光数字信号的信噪比与预定阈值进行比对进一步判断第二红光数字信号的信号质量;
步骤62,根据预置的带阻滤波频率阈值范围,对第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足带阻滤波频率阈值范围的有效信号从第二红光数字信号中滤除生成红光噪声信号;并根据带阻滤波频率阈值范围,对第二绿光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足带阻滤波频率阈值范围的有效信号从第二绿光数字信号中滤除生成绿光噪声信号;
步骤63,计算第二红光数字信号的信号能量生成红光信号能量,计算红光噪声信号的信号能量生成红光噪声能量,根据红光信号能量减去红光噪声能量的差生成有效红光信号能量,根据有效红光信号能量与红光噪声能量的比值生成红光信噪比;计算第二绿光数字信号的信号能量生成绿光信号能量,计算绿光噪声信号的信号能量生成绿光噪声能量,根据绿光信号能量减去绿光噪声能量的差生成有效绿光信号能量,根据有效绿光信号能量与绿光噪声能量的比值生成绿光信噪比;
步骤64,如果红光信噪比与绿光信噪比均小于预置的信噪比阈值则设置临时第二判断结果为不达标信号标识;如果红光信噪比与绿光信噪比中任一个大于或等于信噪比阈值则设置临时第二判断结果为达标信号标识;转至步骤67;
此处,步骤62-64,首先对红绿光做二次滤波处理:本次滤波是一种带阻滤波方式,即将属于带阻滤波频率阈值范围内的信号进行抑制,具体的采用的是多阶巴特沃斯带阻滤波方式(例如,4阶巴特沃斯带阻滤波、3阶巴特沃斯带阻滤波、2阶巴特沃斯带阻滤波、1阶巴特沃斯带阻滤波等),通过带阻滤波将有效信号滤掉、对噪声及干扰信号进行保留生成噪声信号;接着对有效信号和噪声信号的能量进行计算生成信噪比;最后根据预置的信噪比阈值对红、绿光数字信号进行信号质量达标与否的识别操作;
步骤65,根据带阻滤波频率阈值范围,对第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足带阻滤波频率阈值范围的有效信号从第二红光数字信号中滤除生成红光噪声信号;计算第二红光数字信号的信号能量生成红光信号能量,计算红光噪声信号的信号能量生成红光噪声能量,根据红光信号能量减去红光噪声能量的差生成有效红光信号能量,根据有效红光信号能量与红光噪声能量的比值生成红光信噪比;
步骤66,如果红光信噪比小于信噪比阈值则设置临时第二判断结果为不达标信号标识;如果红光信噪比大于或等于信噪比阈值则设置临时第二判断结果为达标信号标识;
此处,步骤65-66,是因部分移动终端型号的皮表视频数据中绿色数字信号质量较差,信号判断标识符具体为红光单检标识,所以只对红光做二次滤波处理,通过带阻滤波将有效信号滤掉、对噪声及干扰信号进行保留生成红光噪声信号;接着对有效信号和噪声信号的能量进行计算生成红光信噪比;最后根据预置的信噪比阈值对红光数字信号进行信号质量达标与否的识别操作;
步骤67,设置第二判断结果为临时第二判断结果。
步骤7,当第二判断结果为达标信号标识时,根据信号判断标识符、第二红光数字信号和或第二绿光数字信号进行光体积变化描记图法PPG信号生成处理,生成PPG信号;
具体包括:步骤71,当第二判断结果为达标信号标识时,设置PPG信号为空;初始化PPG信号的通道标识符为空,初始化PPG信号的红光数字信号为空,初始化PPG信号的绿光数字信号为空;
步骤72,当信号判断标识符为红绿光同检标识时,设置PPG信号的通道标识符为红绿通道标识,设置PPG信号的红光数字信号为第二红光数字信号,设置PPG信号的绿光数字信号为第二绿光数字信号;
步骤73,当信号判断标识符为红光单检标识时,设置通道标识符为红通道标识,设置PPG信号的红光数字信号为第二红光数字信号。
此处,PPG信号包括:通道标识符、红光数字信号和绿光数字信号。
此处,如果信号判断标识符为红绿光同检标识,说明用于采集拍摄的移动终端的红、绿光处理能力是合格的,那么通道标识符会被设置为红绿通道标识表示PPG信号包括两种光数字信号(或者说两种通道数据);如果信号判断标识符为红光单检标识,说明用于采集拍摄的移动终端的绿光处理能力是不合格的,那么通道标识符会被设置为红通道标识表示PPG信号只包括一种光数字信号(或者说一种通道数据)。
如图2为本发明实施例二提供的一种光体积变化描记图法信号的生成方法示意图所示,本方法主要包括如下步骤:
步骤101,调用视频拍摄设备以预置的采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据;对皮表视频数据进行帧图像提取处理生成皮表视频数据帧图像序列;
其中,皮表视频数据帧图像序列包括多个皮表视频数据帧图像。
此处,在日常所见的,只要带有摄像头、具备视频拍摄功能的移动终端均可充当本发明实施例的视频拍摄设备;在拍摄之前,需要打开移动终端的闪光灯并保持常亮,以此满足稳定的照射光源需求;在拍摄时,需要保持被拍摄生物体的局部皮肤表面轻微按压在移动终端的摄像头上;对于闪光灯和摄像头比较接近的移动终端,需要将皮肤表面同时覆盖住闪光灯和摄像头;对于闪光灯和摄像头比较远的移动终端,只需要将皮肤表面覆盖住摄像头即可;拍摄过程中皮肤表面需保持静止、用力均匀;采集时间阈值根据实际应用的数据存储能力和分析能力可多次进行预定义。
步骤102,根据预置的红光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧红色通道数据计算,生成第一红光数字信号;并根据预置的绿光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧绿色通道数据计算,生成第一绿光数字信号。
此处,实施例二通过步骤102将由皮表视频数据转化而来的帧图像序列中的所有帧图像进行多种光源信息的提取操作,本实施例中暂时针对红光和绿光两种光源;对光信号的提取方式,就是通过对帧图像中特定像素的加权平均计算得到一个像素均值,并以此代表该光源在所在帧图像中的颜色通道数据;按时间先后顺序,对视频中的每一帧都做同样的处理,可以得到两段一维数字信号:第一红光数字信号和第一绿光数字信号。
步骤103,根据预置的带通滤波频率阈值范围对第一红光数字信号进行信号带通滤波预处理生成第二红光数字信号,并根据带通滤波频率阈值范围对第一绿光数字信号进行信号带通滤波预处理生成第二绿光数字信号。
此处,步骤103,是对通过视频数据提取的两种光源的数字信号进行信号信号预处理,即降噪处理;此处,实施例一使用的降噪手段是带通滤波方式,即预置一个带通滤波频率阈值范围,基于带通滤波原理对低于或高于该频段的信号、干扰和噪声进行信号抑制处理;一般此处的带通滤波频率阈值范围常见的0.5赫兹到10赫兹;在某些移动终端上进行带通滤波处理时,考虑到移动终端的处理能力,使用的是有限长单位冲激响应((Finitelmpulse Response,FIR)滤波模块。
步骤104,获取预置的信号判断标识符;
其中,信号判断标识符包括红绿光同检标识和红光单检标识两种标识符。
此处,之所以要设置信号判断标识符进行选择判断是为了提高对拍摄设备的兼容度,因为在具体应用场景中,有些移动终端因为镜头质量、设备老化、偏光算法落后等多种原因,会导致其对绿光的处理能力较差,通过这些移动终端所采集的视频数据乃至视频的帧图像数据中,基本很难获得足够数量及像素值达标的像素信息,也就无法产生绿光数字信号;在这种情况下,就需要针对不同情况通过设置信号判断标识符进行灵活处理:红绿光同检标识意味着红光、绿光信号都需要做识别操作,红光单检标识则标识只需对红光信号做识别操作。
步骤105,根据信号判断标识符判断是否对第二红光数字信号和第二绿光数字信号进行信号最大频差判断处理,当信号判断标识符为红绿光同检标识时对第二红光数字信号和第二绿光数字信号进行信号最大频差判断处理生成第一判断结果,当信号判断标识符为红光单检标识时直接设置第一判断结果为达标信号标识;
具体包括:步骤1051,初始化临时第一判断结果为空;
步骤1052,判断信号判断标识符是否为红绿光同检标识,当信号判断标识符为红绿光同检标识时转至步骤1053,当信号判断标识符为红光单检标识时转至步骤1055;
此处,如果信号判断标识符具体被设置为红绿光同检标识时,说明当前用于拍摄的拍摄设备拍摄对红绿光的处理能力正常,可以转至步骤1053通过对第二红光数字信号和第二绿光数字信号进行频率差的关联检测来判断二者的信号质量;如果信号判断标识符具体被设置为红光单检标识时,说明当前用于拍摄的拍摄设备拍摄对绿光的处理能力较差,上述的第二绿光数字信号将会在后续的光信号处理过程中被抛弃,也即无法被引用到最大频差判断处理中,为保证流程能够继续执行,此处会转至步骤1055,不对红绿光进行频差比较而直接将第一判断结果设置为达标信号标识;
步骤1053,对第二红光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成红光频域信号,对第二绿光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成绿光频域信号;从红光频域信号中提取能量最高频率生成红光最大频率,从绿光频域信号中提取能量最高频率生成绿光最大频率;根据红光最大频率与绿光最大频率的频率差生成红绿最大频差;
步骤1054,根据预置的最大频差阈值范围对红绿最大频差进行信号最大频差判断处理,当红绿最大频差未超过最大频差阈值范围时设置临时第一判断结果为达标信号标识,当红绿最大频差超过最大频差阈值范围时设置临时第一判断结果为不达标信号标识;转至步骤1056;
此处,首先通过离散傅里叶变换得到第二红色数字信号和第二绿色数字信号的频域信号,通过频域信号得到能量最高的频率(一般这个频率通常对应着心率),此处基本原理是检查这两个数字信号的能量最高的频率是否一致,如果误差在允许误差范围之内,则继续分析,如果误差较大,则认为信号质量不好,需要重新拍摄;
步骤1055,设置临时第一判断结果为达标信号标识;
此处,如前文所述,部分移动终端型号的皮表视频数据中绿色数字信号质量较差,为兼容多种移动终端,设置信号判断标识符为红光单检标识,将频差检测信号质量的处理流程直接跳过;
步骤1056,设置第一判断结果为临时第一判断结果。
步骤106,判断第一判断结果是否为达标信号标识,如果第一判断结果为不达标信号标识则转至步骤107,如果第一判断结果为达标信号标识则转至步骤108。
此处,当第一判断结果为不达标信号标识时,应转至步骤107立即停止PPG信号分析,并生成信号质量错误信息,后续处理流程会将信号质量错误信息向上位应用发送,由上位应用决定是否向被拍摄用户示警、提醒主动重新拍摄还是自动重新拍摄等;当第一判断结果为达标信号标识时,转至步骤108对光信号进行进一步的质量判定操作。
步骤107,停止PPG信号生成处理流程,并生成信号质量错误信息;转至步骤111。
步骤108,当第一判断结果为达标信号标识时,根据信号判断标识符对第二红光数字信号和或第二绿光数字信号进行信号信噪比判断处理,生成第二判断结果;
具体包括:步骤1081,初始化临时第二判断结果为空;
步骤1081,当第一判断结果为达标信号标识时,根据信号判断标识符进行判断,当信号判断标识符为红绿光同检标识时转至步骤1082,当信号判断标识符为红光单检标识时转至步骤1085;
此处,如果信号判断标识符具体被设置为红绿光同检标识时,说明当前用于拍摄的拍摄设备拍摄对红绿光的处理能力正常,可以转至步骤1082使用第二红光数字信号和第二绿光数字信号的信噪比与预定阈值进行比对进一步判断二者的信号质量;如果信号判断标识符具体被设置为红光单检标识时,说明当前用于拍摄的拍摄设备拍摄对绿光的处理能力较差,上述的第二绿光数字信号将会在后续的光信号处理过程中被抛弃,此处会转至步骤1085只使用第二红光数字信号的信噪比与预定阈值进行比对进一步判断第二红光数字信号的信号质量;
步骤1082,根据预置的带阻滤波频率阈值范围,对第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足带阻滤波频率阈值范围的有效信号从第二红光数字信号中滤除生成红光噪声信号;并根据带阻滤波频率阈值范围,对第二绿光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足带阻滤波频率阈值范围的有效信号从第二绿光数字信号中滤除生成绿光噪声信号;
步骤1083,计算第二红光数字信号的信号能量生成红光信号能量,计算红光噪声信号的信号能量生成红光噪声能量,根据红光信号能量减去红光噪声能量的差生成有效红光信号能量,根据有效红光信号能量与红光噪声能量的比值生成红光信噪比;计算第二绿光数字信号的信号能量生成绿光信号能量,计算绿光噪声信号的信号能量生成绿光噪声能量,根据绿光信号能量减去绿光噪声能量的差生成有效绿光信号能量,根据有效绿光信号能量与绿光噪声能量的比值生成绿光信噪比;
步骤1084,如果红光信噪比与绿光信噪比均小于预置的信噪比阈值则设置临时第二判断结果为不达标信号标识;如果红光信噪比与绿光信噪比中任一个大于或等于信噪比阈值则设置临时第二判断结果为达标信号标识;转至步骤1087;
此处,步骤1082-1084,首先对红绿光做二次滤波处理:本次滤波是一种带阻滤波方式,即将属于带阻滤波频率阈值范围内的信号进行抑制,具体的采用的是多阶巴特沃斯带阻滤波方式(例如,4阶巴特沃斯带阻滤波、3阶巴特沃斯带阻滤波、2阶巴特沃斯带阻滤波、1阶巴特沃斯带阻滤波等),通过带阻滤波将有效信号滤掉、对噪声及干扰信号进行保留生成噪声信号;接着对有效信号和噪声信号的能量进行计算生成信噪比;最后根据预置的信噪比阈值对红、绿光数字信号进行信号质量达标与否的识别操作;
步骤1085,根据带阻滤波频率阈值范围,对第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足带阻滤波频率阈值范围的有效信号从第二红光数字信号中滤除生成红光噪声信号;计算第二红光数字信号的信号能量生成红光信号能量,计算红光噪声信号的信号能量生成红光噪声能量,根据红光信号能量减去红光噪声能量的差生成有效红光信号能量,根据有效红光信号能量与红光噪声能量的比值生成红光信噪比;
步骤1086,如果红光信噪比小于信噪比阈值则设置临时第二判断结果为不达标信号标识;如果红光信噪比大于或等于信噪比阈值则设置临时第二判断结果为达标信号标识;
此处,步骤1085-1086,是因部分移动终端型号的皮表视频数据中绿色数字信号质量较差,信号判断标识符具体为红光单检标识,所以只对红光做二次滤波处理,通过带阻滤波将有效信号滤掉、对噪声及干扰信号进行保留生成红光噪声信号;接着对有效信号和噪声信号的能量进行计算生成红光信噪比;最后根据预置的信噪比阈值对红光数字信号进行信号质量达标与否的识别操作;
步骤1087,设置第二判断结果为临时第二判断结果。
步骤109,判断第二判断结果是否为达标信号标识,如果第二判断结果为不达标信号标识则转至步骤107,如果第二判断结果为达标信号标识则转至步骤110。
此处,当第二判断结果为不达标信号标识时,应转至步骤107立即停止PPG信号分析,并生成信号质量错误信息,后续处理流程会将信号质量错误信息向上位应用发送,由上位应用决定是否向被拍摄用户示警、提醒主动重新拍摄还是自动重新拍摄等;;当第二判断结果为达标信号标识时,转至步骤110对最终生成的PPG信号进行设置。
步骤110,当第二判断结果为达标信号标识时,根据信号判断标识符、第二红光数字信号和或第二绿光数字信号进行光体积变化描记图法PPG信号生成处理,生成PPG信号;
具体包括:步骤1101,当第二判断结果为达标信号标识时,设置PPG信号为空;初始化PPG信号的通道标识符为空,初始化PPG信号的红光数字信号为空,初始化PPG信号的绿光数字信号为空;
步骤1102,当信号判断标识符为红绿光同检标识时,设置PPG信号的通道标识符为红绿通道标识,设置PPG信号的红光数字信号为第二红光数字信号,设置PPG信号的绿光数字信号为第二绿光数字信号;
步骤1103,当信号判断标识符为红光单检标识时,设置通道标识符为红通道标识,设置PPG信号的红光数字信号为第二红光数字信号。
此处,PPG信号包括:通道标识符、红光数字信号和绿光数字信号。
此处,如果信号判断标识符为红绿光同检标识,说明用于采集拍摄的移动终端的红、绿光处理能力是合格的,那么通道标识符会被设置为红绿通道标识表示PPG信号包括两种光数字信号(或者说两种通道数据);如果信号判断标识符为红光单检标识,说明用于采集拍摄的移动终端的绿光处理能力是不合格的,那么通道标识符会被设置为红通道标识表示PPG信号只包括一种光数字信号(或者说一种通道数据)。
步骤111,将生成的数据,向上位应用发送。
此处,如果是从步骤107跳转过来的处理流程,那么生成的数据具体为信号质量错误信息;如果是从步骤110执行下来的处理流程,那么生成的数据具体为PPG信号。
如图3为本发明实施例三提供的一种光体积变化描记图法信号的生成装置的设备结构示意图所示,该设备包括:处理器和存储器。存储器可通过总线与处理器连接。存储器可以是非易失存储器,例如硬盘驱动器和闪存,存储器中存储有软件程序和设备驱动程序。软件程序能够执行本发明实施例提供的上述方法的各种功能;设备驱动程序可以是网络和接口驱动程序。处理器用于执行软件程序,该软件程序被执行时,能够实现本发明实施例提供的方法。
需要说明的是,本发明实施例还提供了一种计算机可读存储介质。该计算机可读存储介质上存储有计算机程序,该计算机程序被处理器执行时,能够实现本发明实施例提供的方法。
本发明实施例还提供了一种包含指令的计算机程序产品。当该计算机程序产品在计算机上运行时,使得处理器执行上述方法。
本发明实施例提供的一种光体积变化描记图法信号的生成方法和装置,通过使用本发明实施例的方法,不规定特定照射光源类型采用普通灯光即可,简化了PPG信号前端采集装置使用一般的移动终端也可实现;通过对录取的视频数据进行分析并生成PPG信号,即能满足本地设备实时分析与数据备份的要求,又能满足远程传输与分析的需求;对有效信号的识别过程不依赖特定传感器,降低了应用层面对PPG数据的获取难度;通过设置信号判断标识符,提高对多种拍摄设备的兼容性。
专业人员应该还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以用硬件、处理器执行的软件模块,或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (13)
1.一种光体积变化描记图法信号的生成方法,其特征在于,所述方法包括:
根据预置的红光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧红色通道数据计算,生成第一红光数字信号;并根据预置的绿光像素阈值范围,对所述皮表视频数据帧图像序列的所有所述皮表视频数据帧图像依次进行单帧绿色通道数据计算,生成第一绿光数字信号;所述皮表视频数据帧图像序列包括多个所述皮表视频数据帧图像;
根据预置的带通滤波频率阈值范围对所述第一红光数字信号进行信号带通滤波预处理生成第二红光数字信号,并根据所述带通滤波频率阈值范围对所述第一绿光数字信号进行信号带通滤波预处理生成第二绿光数字信号;
获取预置的信号判断标识符;所述信号判断标识符包括红绿光同检标识和红光单检标识两种标识符;
根据所述信号判断标识符判断是否对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理,当所述信号判断标识符为所述红绿光同检标识时对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理生成第一判断结果,当所述信号判断标识符为所述红光单检标识时直接设置所述第一判断结果为达标信号标识;
当所述第一判断结果为所述达标信号标识时,根据所述信号判断标识符对所述第二红光数字信号和或所述第二绿光数字信号进行信号信噪比判断处理,生成第二判断结果;
当所述第二判断结果为所述达标信号标识时,根据所述信号判断标识符、所述第二红光数字信号和或所述第二绿光数字信号进行光体积变化描记图法PPG信号生成处理,生成PPG信号;
其中,所述根据预置的红光像素阈值范围,对皮表视频数据帧图像序列的所有皮表视频数据帧图像依次进行单帧红色通道数据计算,生成第一红光数字信号;并根据预置的绿光像素阈值范围,对所述皮表视频数据帧图像序列的所有所述皮表视频数据帧图像依次进行单帧绿色通道数据计算,生成第一绿光数字信号,具体包括:
步骤61,初始化临时第一红光数字信号为空,初始化临时第一绿光数字信号为空,初始化第一索引的值为1,初始化第一总数为所述皮表视频数据帧图像序列的皮表视频数据帧图像总数;
步骤62,从所述皮表视频数据帧图像序列中提取与所述第一索引对应的所述皮表视频数据帧图像生成第一索引皮表视频数据帧图像;
步骤63,根据所述红光像素阈值范围,在所述第一索引皮表视频数据帧图像中统计所有满足所述红光像素阈值范围的像素点生成红色像素点集合,并对所述红色像素点集合的所有像素点的像素值按指定加权平均方式进行加权平均计算,生成第一索引帧红色通道数据;将所述第一索引帧红色通道数据作为红光信号点数据向所述临时第一红光数字信号进行信号点添加操作;
步骤64,根据所述绿光像素阈值范围,在所述第一索引皮表视频数据帧图像中统计所有满足所述绿光像素阈值范围的像素点生成绿色像素点集合,并对所述绿色像素点集合的所有像素点的像素值按指定加权平均方式进行加权平均计算,生成第一索引帧绿色通道数据;将所述第一索引帧绿色通道数据作为绿光信号点数据向所述临时第一绿光数字信号进行信号点添加操作;
步骤65,将所述第一索引的值加1;
步骤66,判断所述第一索引是否大于所述第一总数,如果所述第一索引小于或等于所述第一总数则转至步骤62,如果所述第一索引大于所述第一总数则转至步骤67;
步骤67,设置所述第一红光数字信号为所述临时第一红光数字信号;设置所述第一绿光数字信号为所述临时第一绿光数字信号。
2.根据权利要求1所述的光体积变化描记图法信号的生成方法,其特征在于,
所述皮表视频数据帧图像序列是,调用视频拍摄设备以预置的采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据;对所述皮表视频数据进行帧图像提取处理生成的。
3.根据权利要求2所述的光体积变化描记图法信号的生成方法,其特征在于,所述调用视频拍摄设备以预置的采集时间阈值为拍摄时间长度对生物体局部皮肤表面进行连续采集拍摄操作生成皮表视频数据,具体包括:
调用所述视频拍摄设备的灯光照明模块对所述生物体局部皮肤表面提供照明;所述视频拍摄设备包括灯光照明模块和视频拍摄模块;所述视频拍摄模块至少包括镜头模块;
当所述灯光照明模块与所述镜头模块的物理间距小于或等于预置的设备间距阈值时,将所述生物体局部皮肤表面同时覆盖所述灯光照明模块与所述镜头模块;当所述灯光照明模块与所述镜头模块的物理间距大于所述设备间距阈值时,将所述生物体局部皮肤表面覆盖所述镜头模块;
调用所述视频拍摄模块以所述采集时间阈值为拍摄时间长度对所述生物体局部皮肤表面进行所述连续采集拍摄操作生成所述皮表视频数据。
4.根据权利要求2所述的光体积变化描记图法信号的生成方法,其特征在于,所述方法还包括:
在调用所述视频拍摄设备以所述采集时间阈值为拍摄时间长度对所述生物体局部皮肤表面进行所述连续采集拍摄操作的过程中,以预置的缓存片段时间阈值为片段提取长度对所述视频拍摄设备采集拍摄的视频数据进行连续视频片段提取处理生成多个缓存片段;
提取所述缓存片段的帧图像并将提取出的所有帧图像按时间顺序进行排序,生成缓存片段帧图像序列;所述缓存片段包括多个所述帧图像;
对所述缓存片段帧图像序列的所有缓存片段帧图像,根据所述红光像素阈值范围进行红色像素点占比统计生成红色占比参数,并根据预置的杂色像素阈值范围进行杂色像素点占比统计生成杂色占比参数;所述缓存片段帧图像序列包括多个所述缓存片段帧图像;
根据所述红色占比参数对所述缓存片段进行皮肤表面脱离拍摄错误识别判断,当所述红色占比参数小于预置的红色占比下限阈值时停止所述连续采集拍摄操作并生成皮肤表面脱离拍摄错误信息;
根据所述杂色占比参数对所述缓存片段进行皮肤表面脱离拍摄错误识别判断,当所述杂色占比参数大于预置的杂色占比上限阈值时停止所述连续采集拍摄操作并生成所述皮肤表面脱离拍摄错误信息。
5.根据权利要求1所述的光体积变化描记图法信号的生成方法,其特征在于,
所述第一判断结果包括所述达标信号标识和不达标信号标识两种判断结果;
所述第二判断结果包括所述达标信号标识和所述不达标信号标识两种判断结果。
6.根据权利要求1所述的光体积变化描记图法信号的生成方法,其特征在于,所述根据预置的带通滤波频率阈值范围对所述第一红光数字信号进行信号带通滤波预处理生成第二红光数字信号,并根据所述带通滤波频率阈值范围对所述第一绿光数字信号进行信号带通滤波预处理生成第二绿光数字信号,具体包括:
根据所述带通滤波频率阈值范围,在所述第一红光数字信号中,将信号频率低于所述带通滤波频率阈值范围的低频噪声信号点和信号频率高于所述带通滤波频率阈值范围的高频噪声信号点从所述第一红光数字信号中进行数字信号滤波处理,生成所述第二红光数字信号;
根据所述带通滤波频率阈值范围,在所述第一绿光数字信号中,将信号频率低于所述带通滤波频率阈值范围的低频噪声信号点和信号频率高于所述带通滤波频率阈值范围的高频噪声信号点从所述第一绿光数字信号中进行数字信号滤波处理,生成所述第二绿光数字信号。
7.根据权利要求5所述的光体积变化描记图法信号的生成方法,其特征在于,所述根据所述信号判断标识符判断是否对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理,当所述信号判断标识符为所述红绿光同检标识时对所述第二红光数字信号和所述第二绿光数字信号进行信号最大频差判断处理生成第一判断结果,当所述信号判断标识符为所述红光单检标识时直接设置所述第一判断结果为达标信号标识,具体包括:
步骤81,初始化临时第一判断结果为空;
步骤82,判断所述信号判断标识符是否为所述红绿光同检标识,当所述信号判断标识符为所述红绿光同检标识时转至步骤83,当所述信号判断标识符为所述红光单检标识时转至步骤85;
步骤83,对所述第二红光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成红光频域信号,对所述第二绿光数字信号使用离散傅里叶变换进行数字信号时域频域转换生成绿光频域信号;从所述红光频域信号中提取能量最高频率生成红光最大频率,从所述绿光频域信号中提取能量最高频率生成绿光最大频率;根据所述红光最大频率与所述绿光最大频率的频率差生成红绿最大频差;
步骤84,根据预置的最大频差阈值范围对所述红绿最大频差进行信号最大频差判断处理,当所述红绿最大频差未超过所述最大频差阈值范围时设置所述临时第一判断结果为所述达标信号标识,当所述红绿最大频差超过所述最大频差阈值范围时设置所述临时第一判断结果为所述不达标信号标识;转至步骤86;
步骤85,设置所述临时第一判断结果为所述达标信号标识;
步骤86,设置所述第一判断结果为所述临时第一判断结果。
8.根据权利要求5所述的光体积变化描记图法信号的生成方法,其特征在于,所述当所述第一判断结果为所述达标信号标识时,根据所述信号判断标识符对所述第二红光数字信号和或所述第二绿光数字信号进行信号信噪比判断处理,生成第二判断结果,具体包括:
步骤91,初始化临时第二判断结果为空;
步骤91,当所述第一判断结果为所述达标信号标识时,根据所述信号判断标识符进行判断,当所述信号判断标识符为所述红绿光同检标识时转至步骤92,当所述信号判断标识符为所述红光单检标识时转至步骤95;
步骤92,根据预置的带阻滤波频率阈值范围,对所述第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足所述带阻滤波频率阈值范围的有效信号从所述第二红光数字信号中滤除生成红光噪声信号;并根据所述带阻滤波频率阈值范围,对所述第二绿光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足所述带阻滤波频率阈值范围的有效信号从所述第二绿光数字信号中滤除生成绿光噪声信号;
步骤93,计算所述第二红光数字信号的信号能量生成红光信号能量,计算所述红光噪声信号的信号能量生成红光噪声能量,根据所述红光信号能量减去所述红光噪声能量的差生成有效红光信号能量,根据所述有效红光信号能量与所述红光噪声能量的比值生成红光信噪比;计算所述第二绿光数字信号的信号能量生成绿光信号能量,计算所述绿光噪声信号的信号能量生成绿光噪声能量,根据所述绿光信号能量减去所述绿光噪声能量的差生成有效绿光信号能量,根据所述有效绿光信号能量与所述绿光噪声能量的比值生成绿光信噪比;
步骤94,如果所述红光信噪比与所述绿光信噪比均小于预置的信噪比阈值则设置所述临时第二判断结果为所述不达标信号标识;如果所述红光信噪比与所述绿光信噪比中任一个大于或等于所述信噪比阈值则设置所述临时第二判断结果为所述达标信号标识;转至步骤97;
步骤95,根据所述带阻滤波频率阈值范围,对所述第二红光数字信号进行多阶巴特沃斯带阻滤波处理,将信号频率满足所述带阻滤波频率阈值范围的有效信号从所述第二红光数字信号中滤除生成所述红光噪声信号;计算所述第二红光数字信号的信号能量生成所述红光信号能量,计算所述红光噪声信号的信号能量生成所述红光噪声能量,根据所述红光信号能量减去所述红光噪声能量的差生成所述有效红光信号能量,根据所述有效红光信号能量与所述红光噪声能量的比值生成所述红光信噪比;
步骤96,如果所述红光信噪比小于所述信噪比阈值则设置所述临时第二判断结果为所述不达标信号标识;如果所述红光信噪比大于或等于所述信噪比阈值则设置所述临时第二判断结果为所述达标信号标识;
步骤97,设置所述第二判断结果为所述临时第二判断结果。
9.根据权利要求1所述的光体积变化描记图法信号的生成方法,其特征在于,所述当所述第二判断结果为所述达标信号标识时,根据所述信号判断标识符、所述第二红光数字信号和或所述第二绿光数字信号进行光体积变化描记图法PPG信号生成处理,生成PPG信号,具体包括:
当所述第二判断结果为所述达标信号标识时,设置所述PPG信号为空;初始化所述PPG信号的通道标识符为空,初始化所述PPG信号的红光数字信号为空,初始化所述PPG信号的绿光数字信号为空;
当所述信号判断标识符为所述红绿光同检标识时,设置所述PPG信号的所述通道标识符为红绿通道标识,设置所述PPG信号的所述红光数字信号为所述第二红光数字信号,设置所述PPG信号的所述绿光数字信号为所述第二绿光数字信号;
当所述信号判断标识符为所述红光单检标识时,设置所述通道标识符为红通道标识,设置所述PPG信号的所述红光数字信号为所述第二红光数字信号。
10.根据权利要求1所述的光体积变化描记图法信号的生成方法,其特征在于,所述方法还包括:
当所述第一判断结果为不达标信号标识时,停止PPG信号生成处理流程,并生成信号质量错误信息;
当所述第二判断结果为不达标信号标识时,停止PPG信号生成处理流程,并生成信号质量错误信息。
11.一种设备,包括存储器和处理器,其特征在于,所述存储器用于存储程序,所述处理器用于执行如权利要求1至10任一项所述的方法。
12.一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行如权利要求1至10任一项所述的方法。
13.一种计算机可读存储介质,包括指令,当所述指令在计算机上运行时,使所述计算机执行根据权利要求1至10任一项所述的方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010110280.7A CN111297339B (zh) | 2020-02-21 | 2020-02-21 | 一种光体积变化描记图法信号的生成方法和装置 |
EP20920454.4A EP4108168A4 (en) | 2020-02-21 | 2020-11-18 | METHOD AND DEVICE FOR GENERATION OF PHOTOPLETHYSMOGRAPHY SIGNALS |
PCT/CN2020/129634 WO2021164348A1 (zh) | 2020-02-21 | 2020-11-18 | 一种光体积变化描记图法信号的生成方法和装置 |
US17/904,688 US20230071469A1 (en) | 2020-02-21 | 2020-11-18 | Method and device for generating photoplethysmography signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010110280.7A CN111297339B (zh) | 2020-02-21 | 2020-02-21 | 一种光体积变化描记图法信号的生成方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111297339A CN111297339A (zh) | 2020-06-19 |
CN111297339B true CN111297339B (zh) | 2022-07-29 |
Family
ID=71156736
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010110280.7A Active CN111297339B (zh) | 2020-02-21 | 2020-02-21 | 一种光体积变化描记图法信号的生成方法和装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20230071469A1 (zh) |
EP (1) | EP4108168A4 (zh) |
CN (1) | CN111297339B (zh) |
WO (1) | WO2021164348A1 (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111297339B (zh) * | 2020-02-21 | 2022-07-29 | 乐普(北京)医疗器械股份有限公司 | 一种光体积变化描记图法信号的生成方法和装置 |
CN114073520A (zh) * | 2020-08-12 | 2022-02-22 | 华为技术有限公司 | 基于绿光的血氧检测设备及其血氧检测方法和介质 |
CN112315437A (zh) * | 2020-10-12 | 2021-02-05 | 乐普(北京)医疗器械股份有限公司 | 一种基于视频数据进行血压预测的系统 |
CN114176583A (zh) * | 2021-12-03 | 2022-03-15 | Oppo广东移动通信有限公司 | 血氧测量方法及相关装置 |
CN117084644B (zh) * | 2023-10-10 | 2024-03-29 | 荣耀终端有限公司 | 一种检测用户生理特征指标的方法、装置及可穿戴设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108471989A (zh) * | 2016-01-15 | 2018-08-31 | 皇家飞利浦有限公司 | 用于生成承载对象的生命体征信息的光体积描记图像的设备、系统和方法 |
CN109820499A (zh) * | 2018-12-24 | 2019-05-31 | 杨爽 | 基于视频的高抗干扰心率检测方法、电子设备及存储介质 |
CN110251115A (zh) * | 2019-06-20 | 2019-09-20 | 张远 | 基于体表视频的ppg信号提取方法、系统、设备及介质 |
CN110558966A (zh) * | 2019-08-12 | 2019-12-13 | 平安科技(深圳)有限公司 | 一种脉搏检测装置 |
Family Cites Families (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10244949B2 (en) * | 2012-10-07 | 2019-04-02 | Rhythm Diagnostic Systems, Inc. | Health monitoring systems and methods |
TWI546052B (zh) * | 2013-11-14 | 2016-08-21 | 財團法人工業技術研究院 | 影像式心率活動偵測裝置及其方法 |
JP6349075B2 (ja) * | 2013-11-22 | 2018-06-27 | 三星電子株式会社Samsung Electronics Co.,Ltd. | 心拍数測定装置及び心拍数測定方法 |
TWI556793B (zh) * | 2015-07-09 | 2016-11-11 | 國立臺灣科技大學 | 非接觸式即時生理訊號及動作偵測方法 |
US10888280B2 (en) * | 2016-09-24 | 2021-01-12 | Sanmina Corporation | System and method for obtaining health data using a neural network |
CN105395184B (zh) * | 2015-12-04 | 2018-05-08 | 华中科技大学 | 生物组织血流、血氧和血容量的多参数检测方法及装置 |
GB201522406D0 (en) * | 2015-12-18 | 2016-02-03 | Xim Ltd | A method, information processing apparatus and server for determining a physiological parameter of an individual |
CA2958010C (en) * | 2016-02-19 | 2021-09-07 | Covidien Lp | System and methods for video-based monitoring of vital signs |
CN111297339B (zh) * | 2020-02-21 | 2022-07-29 | 乐普(北京)医疗器械股份有限公司 | 一种光体积变化描记图法信号的生成方法和装置 |
-
2020
- 2020-02-21 CN CN202010110280.7A patent/CN111297339B/zh active Active
- 2020-11-18 US US17/904,688 patent/US20230071469A1/en active Pending
- 2020-11-18 WO PCT/CN2020/129634 patent/WO2021164348A1/zh unknown
- 2020-11-18 EP EP20920454.4A patent/EP4108168A4/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108471989A (zh) * | 2016-01-15 | 2018-08-31 | 皇家飞利浦有限公司 | 用于生成承载对象的生命体征信息的光体积描记图像的设备、系统和方法 |
CN109820499A (zh) * | 2018-12-24 | 2019-05-31 | 杨爽 | 基于视频的高抗干扰心率检测方法、电子设备及存储介质 |
CN110251115A (zh) * | 2019-06-20 | 2019-09-20 | 张远 | 基于体表视频的ppg信号提取方法、系统、设备及介质 |
CN110558966A (zh) * | 2019-08-12 | 2019-12-13 | 平安科技(深圳)有限公司 | 一种脉搏检测装置 |
Also Published As
Publication number | Publication date |
---|---|
WO2021164348A1 (zh) | 2021-08-26 |
CN111297339A (zh) | 2020-06-19 |
EP4108168A4 (en) | 2023-08-02 |
US20230071469A1 (en) | 2023-03-09 |
EP4108168A1 (en) | 2022-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111297339B (zh) | 一种光体积变化描记图法信号的生成方法和装置 | |
US10004410B2 (en) | System and methods for measuring physiological parameters | |
CN110269600B (zh) | 基于多元经验模态分解与联合盲源分离的非接触式视频心率检测方法 | |
CN111297347B (zh) | 一种生成光体积变化描记图法信号的方法和装置 | |
RU2669616C2 (ru) | Устройство и способ определения показателей жизнедеятельности субъекта | |
US20120190947A1 (en) | Physiological parameter monitoring with a mobile communication device | |
CN106793962A (zh) | 用于使用视频图像来连续估计人体血压的方法和装置 | |
CN111243739A (zh) | 抗干扰的生理参数遥测方法及系统 | |
CN112233813A (zh) | 一种基于ppg的非接触式无创心率呼吸测量方法及系统 | |
CN112890792A (zh) | 一种基于网络摄像头的云计算心血管健康监护系统及方法 | |
CN110251115B (zh) | 基于体表视频的ppg信号提取方法、系统、设备及介质 | |
CN111387959A (zh) | 一种基于ippg的非接触式生理参数检测方法 | |
US20200015688A1 (en) | Blood pressure measurement method, device and storage medium | |
CN112315437A (zh) | 一种基于视频数据进行血压预测的系统 | |
CN109620195B (zh) | 基于可穿戴设备的多波长信号融合心率检测方法及系统 | |
CN107863156A (zh) | 一种高血压筛查与管理方法与装置 | |
Kyal et al. | A method to detect cardiac arrhythmias with a webcam | |
CN114387479A (zh) | 一种基于人脸视频的非接触式心率测量方法及系统 | |
CN109770885A (zh) | 一种基于预览帧的快速心率检测方法 | |
Corino et al. | Identification of atrial fibrillation episodes using a camera as contactless sensor | |
Djeldjli et al. | Robust heart activity measurement using webcam | |
CN115399745A (zh) | 基于红外成像的心率监测方法、系统、装置和存储介质 | |
CN116889389A (zh) | 基于视频的远程心率检测方法和系统 | |
Sabatini et al. | Heart rate analysis through smartphone camera | |
Djeldjli et al. | A robust photoplethysmographic imaging for contactless heart and respiratory rates measurement using a simple webcam |
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 |