CN108652610A - 一种多人心跳的非接触式检测方法 - Google Patents
一种多人心跳的非接触式检测方法 Download PDFInfo
- Publication number
- CN108652610A CN108652610A CN201810561513.8A CN201810561513A CN108652610A CN 108652610 A CN108652610 A CN 108652610A CN 201810561513 A CN201810561513 A CN 201810561513A CN 108652610 A CN108652610 A CN 108652610A
- Authority
- CN
- China
- Prior art keywords
- value
- jumped
- queue
- detection method
- concurrent queue
- 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
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- Psychiatry (AREA)
- Cardiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明公开一种多人心跳的非接触式检测方法,应用于机器智能化技术领域,为解决现有技术中心跳检测方法的预处理过程较慢问题,本发明通过在GPU缓冲里构造并发队列,基于非接触式采集的皮肤反射信号,高效、实时地计算心率值。
Description
技术领域
本发明属于机器智能化技术领域,特别涉及一种基于GPU的多人心跳检测技术。
背景技术
心率,指人心脏跳动的频率,通常采用每分钟心跳的次数作为单位。心跳与人的健康状况和情绪状态密切相关。据报道,心跳状况与多种疾病,如心血管疾病及心脏疾病等密切相关,心率异常是多种心血管病的高危因素,通过测量它,可预知不少心血管病的发生,这样就能尽早开展诊治,以获得更好的疗效。人的兴奋、激动、吸烟、饮酒、喝浓茶或咖啡后,或见于感染、发热、休克、贫血、缺氧、甲亢、心力衰竭等病理状态下,心率都会发生相应变化。因而及时检测并或者心跳的速率及其变化趋势,在医学及对人的情绪行为判断方面具有重要意义。
检测心跳主要有两类方法,接触式和非接触式的。采用听诊器等方法测心跳,属于接触式的方法,特点是测试比较准确,但要求检测心跳时必须接触被测人员的身体,因而限制比较多。非接触式的检测方法,则要求在检测心跳时必须与被测人员的身体保持一定距离,通常基于人体皮肤表面对微波、红外线、可见光的反射频谱进行分析,采用信号分析的方法获取心率估计值。当心跳检测涉及多人时,需要对多个人并行地收集相应的心率样本数据,进行预处理,以便进一步进行计算,但现有预处理方法基于CPU并发性能较差。
发明内容
为解决现有技术中心跳检测方法的预处理过程较慢问题,本发明提出一种多人心跳的非接触式检测方法,通过在GPU缓存里构造并发队列,能支持上百人的心跳检测过程的并发预处理。
本发明采用的技术方案:一种多人心跳的非接触式检测方法,包括:
S1、构建一个支持M个并发队列的GPU缓存;
S2、在GPU里为每一个并发队列分配一个线程;
S3、每一个GPU线程内持续采样同一检测对象的同一部位裸露皮肤的色度变化值,并将采样值加入并发队列;
S4、在GPU线程内,计算对应并发队列的中采样值的最大值与最小值,若根据最大值与最小值得到的值范围超出设定值范围,则执行步骤S5,否则执行步骤S6;
S5、从并发队列的起始位置逐一删除采样值,重新计算值范围,直到值范围小于或等于设定门限范围;
S6、按照时间等间隔对并发队列中的采样值数据进行插值处理,并进行汉明处理,完成与处理;
S7、根据步骤S6的预处理结果计算每个检测对象的心率。
进一步地,步骤S1所述每个并发队列为一个循环数组,包含三个域:ID、Start、N;ID为人的标识,由CPU传入;N是并发队列内采样数据的个数;Start为并发队列中第一个采样数据的位置。
进一步地,当步骤S3所述采样个数少于N时,则继续缓存,直至采样个数超过N。
更进一步地,当采样个数超过N时,逐一删除并发队列Start位置处的采样值,然后将新的采样数据加入并发队列。
进一步地,步骤S6所述插值处理后,该并发队列中数据的总数为N。
进一步地,步骤S7具体包括以下分步骤:
S71、根据与处理结果,按并发队列进行一维DFT变换,从而计算每个心率对应的频谱及其增益;
S72、将步骤S71的计算结果中的频谱进行过滤,剩下满足条件的频谱;通过对比满足条件频谱所对应的增益,找出最大增益值对应的频谱;
S73、根据步骤S72找出的频谱得到对应检测对象的估计心率值。
进一步地,步骤S72所述满足条件的频谱具体为在设定频谱范围内的频谱。
本发明的有益效果:一种多人心跳的非接触式检测方法,通过在GPU缓冲里构造并发队列,基于非接触式采集的皮肤反射信号,高效、实时地计算心率值,本发明的方法能支持上百人的心跳检测过程的并发预处理,由于减少了CPU的负担,能够高效、实时地计算多人心率值,可工作于集成了GPU的计算机或手机等设备上,为高效、实时地计算人群心率值创造条件。
附图说明
图1为本发明实施例提供的方案流程图。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
为了克服现有技术中心跳检测方法的预处理过程较慢缺点,本法提出了基于GPU的一种多人心跳的非接触式检测方法,主要内容如下:
首先,基于GPU的裸露皮肤的色度变化值预处理。持续采样同一人的同一部位裸露皮肤的色度变化值,并按每N个采样值为一组进行分组,N为整数。若采样个数少于N,则持续缓存,直到采样个数超过N。每人对应一个线程,在每个线程里对采样后的数据按照时间等间隔对组内的数据进行插值处理。插值处理后,采样个数不变,然后进行汉明处理,完成预处理过程。
其次,在GPU中对每个并发队列里的采样数据进行一维DFT变换。基于上一步的预处理结果,分别计算每个频率对应的频谱f(k),及其增益g(k),如下:
g(k)=Fx×Fx+Fy×Fy
再者,基于GPU频谱增益分析,估算心率值并做过滤处理。基于上一步的DFT计算结果,将心率合理频谱范围外的频率值过滤掉,然后比较合理频谱范围内的各个频率分量的增益。设心率范围为[cl=40次/分钟,ch=160次/分钟],则心率范围对应的频谱下标范围为[il,ih]。
然后,对频谱下标范围内的增益值进行对比,找到增益值最大的频率下标iMax,则估计的心率值为B。
最后,当只有一个队列时,直接输出B,作为心跳检测结果;否则当多于一个并发队列时,则在CPU中将各组的心率计算结果,再进一步采用取平均或kalman过滤,输出最终计算值,作为心跳检测结果。
如图1所示,本发明的一种多人心跳检测方法,包括:
S1、构建一个支持M个并发队列的GPU缓存;每个并发队列是一个循环数组,包含三个域:ID、Start、Size。ID对应人的标识,由CPU传入。Size是队列内采样数据的个数。队列的最大长度设为QueueSize,可配置。Start对应队列中第一个采样数据的位置,(Start+Size-1)%QueueSize对应队列中最后一个采样数据的位置。
S2、在GPU里为每一个并发队列分配一个线程,共M个线程;
S3、每一个GPU线程内持续采样同一检测对象的同一部位裸露皮肤的色度变化值,并将采样值加入并发队列;具体过程为:
基于GPU的裸露皮肤的色度变化值预处理。持续采样同一人的同一部位裸露皮肤的色度变化值,并按每N个采样值为一组进行分组(N为整数,本实施例中取值为64,即一个并发队列中可缓存64个采样数据)。若采样个数少于64,则持续缓存,直到采样个数超过64,则删除Start位置处的采样数据,然后将新的采样数据加入并发队列;
裸露皮肤不限于面部,只要是同一块皮肤即可。将分组处理后的数据,提交给GPU进行计算。
S4、在GPU线程内,计算对应并发队列的中采样值的最大值与最小值,若根据最大值与最小值得到的值范围超出设定值范围SpikeLimit,则执行步骤S5,否则执行步骤S6;这一步骤主要用于过滤掉跳动比较大的采样值,SpikeLimit是一个可配置参数,默认设置为30,在实际应用中根据需要进行取值。
S5、从并发队列的起始位置逐一删除采样值,重新计算值范围,直到值范围小于或等于设定门限SpikeLimit;若是值范围超出SpikeLimit,则从并发队列的起始位置开始每次删除一个采样数据,重新计算值范围,直到值范围小于或等于设定门限SpikeLimit。
S6、按照时间等间隔对并发队列中的采样值数据进行插值处理,插值处理后采样个数不变,还是64个,并进行汉明处理,完成与处理;
S7、根据步骤S6的预处理结果计算每个检测对象的心率;包括以下分步骤:
S71、GPU中按队列进行一维DFT变换,方法是基于上一步的预处理结果,分别计算每个频率对应的频谱f(k),及其增益g(k),如下:
g(k)=Fx×Fx+Fy×Fy,N为并发队列中的采样数据的个数,Fx、Fy分别为f(k)的实部和虚部,i为复数单位。
S72、GPU中,基于上一步的DFT计算结果,将心率合理频谱范围外的频率值过滤掉,然后比较合理频谱范围内的各个频率分量的增益。本实施例中合理频谱范围为[cl=40次/分钟,ch=160次/分钟],则心率范围对应的频谱下标范围为[il,ih]。
设N=64,fps=30,则频谱下标范围为[il=2,ih=6]。
然后对频谱下标范围内的增益值进行对比,找到增益值最大的频率下标iMax,则估计的心率值为B。
比如,假设iMax=3,则B≈85次/分钟。
S73、将不做S72估计的心率值作为每个人的心率检测结果。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (7)
1.一种多人心跳的非接触式检测方法,其特征在于,包括:
S1、构建一个支持M个并发队列的GPU缓存;
S2、在GPU里为每一个并发队列分配一个线程;
S3、每一个GPU线程内持续采样同一检测对象的同一部位裸露皮肤的色度变化值,并将采样值加入并发队列;
S4、在GPU线程内,计算对应并发队列的中采样值的最大值与最小值,若根据最大值与最小值得到的值范围超出设定值范围,则执行步骤S5,否则执行步骤S6;
S5、从并发队列的起始位置逐一删除采样值,重新计算值范围,直到值范围小于或等于设定门限范围;
S6、按照时间等间隔对并发队列中的采样值数据进行插值处理,并进行汉明处理,完成与处理;
S7、根据步骤S6的预处理结果计算每个检测对象的心率。
2.根据权利要求1所述的一种多人心跳的非接触式检测方法,其特征在于,步骤S1所述每个并发队列为一个循环数组,包含三个域:ID、Start、N;ID为人的标识,由CPU传入;N是并发队列内采样数据的个数;Start为并发队列中第一个采样数据的位置。
3.根据权利要求2所述的一种多人心跳的非接触式检测方法,其特征在于,当步骤S3所述采样个数少于N时,则继续缓存,直至采样个数超过N。
4.根据权利要求3所述的一种多人心跳的非接触式检测方法,其特征在于,当采样个数超过N时,逐一删除并发队列Start位置处的采样值,然后将新的采样数据加入并发队列。
5.根据权利要求2所述的一种多人心跳的非接触式检测方法,其特征在于,步骤S6所述插值处理后,该并发队列中数据的总数为N。
6.根据权利要求1所述的一种多人心跳的非接触式检测方法,其特征在于,步骤S7具体包括以下分步骤:
S71、根据与处理结果,按并发队列进行一维DFT变换,从而计算每个心率对应的频谱及其增益;
S72、将步骤S71的计算结果中的频谱进行过滤,剩下满足条件的频谱;通过对比满足条件频谱所对应的增益,找出最大增益值对应的频谱;
S73、根据步骤S72找出的频谱得到对应检测对象的估计心率值。
7.根据权利要求6所述的一种多人心跳的非接触式检测方法,其特征在于,步骤S72所述满足条件的频谱具体为在设定频谱范围内的频谱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810561513.8A CN108652610B (zh) | 2018-06-04 | 2018-06-04 | 一种多人心跳的非接触式检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810561513.8A CN108652610B (zh) | 2018-06-04 | 2018-06-04 | 一种多人心跳的非接触式检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108652610A true CN108652610A (zh) | 2018-10-16 |
CN108652610B CN108652610B (zh) | 2020-12-22 |
Family
ID=63774928
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810561513.8A Active CN108652610B (zh) | 2018-06-04 | 2018-06-04 | 一种多人心跳的非接触式检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108652610B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080079739A1 (en) * | 2006-09-29 | 2008-04-03 | Abhay Gupta | Graphics processor and method for controlling a display panel in self-refresh and low-response-time modes |
US20100184518A1 (en) * | 2006-08-14 | 2010-07-22 | Wms Gaming Inc. | Applying graphical characteristics to graphical objects in a wagering game machine |
CN106361316A (zh) * | 2016-08-30 | 2017-02-01 | 苏州品诺维新医疗科技有限公司 | 一种多人心跳检测系统及获取多人心跳变化曲线的方法 |
CN106725411A (zh) * | 2017-02-15 | 2017-05-31 | 成都皓图智能科技有限责任公司 | 一种基于gpu加速的非接触式心跳检测方法 |
CN107066378A (zh) * | 2017-01-24 | 2017-08-18 | 成都皓图智能科技有限责任公司 | 一种实时在线gpu通用计算模拟及调试方法 |
CN107168782A (zh) * | 2017-04-24 | 2017-09-15 | 复旦大学 | 一种基于Spark与GPU的并行计算系统 |
-
2018
- 2018-06-04 CN CN201810561513.8A patent/CN108652610B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100184518A1 (en) * | 2006-08-14 | 2010-07-22 | Wms Gaming Inc. | Applying graphical characteristics to graphical objects in a wagering game machine |
US20080079739A1 (en) * | 2006-09-29 | 2008-04-03 | Abhay Gupta | Graphics processor and method for controlling a display panel in self-refresh and low-response-time modes |
CN106361316A (zh) * | 2016-08-30 | 2017-02-01 | 苏州品诺维新医疗科技有限公司 | 一种多人心跳检测系统及获取多人心跳变化曲线的方法 |
CN107066378A (zh) * | 2017-01-24 | 2017-08-18 | 成都皓图智能科技有限责任公司 | 一种实时在线gpu通用计算模拟及调试方法 |
CN106725411A (zh) * | 2017-02-15 | 2017-05-31 | 成都皓图智能科技有限责任公司 | 一种基于gpu加速的非接触式心跳检测方法 |
CN107168782A (zh) * | 2017-04-24 | 2017-09-15 | 复旦大学 | 一种基于Spark与GPU的并行计算系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108652610B (zh) | 2020-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105266772B (zh) | 一种生理参数的测量方法 | |
CN102499694B (zh) | 用于消除血氧饱和度监测干扰的方法 | |
Couceiro et al. | Detection of motion artifact patterns in photoplethysmographic signals based on time and period domain analysis | |
CN106725411B (zh) | 一种基于gpu加速的非接触式心跳检测方法 | |
Li et al. | A 410-nW efficient QRS processor for mobile ECG monitoring in 0.18-μm CMOS | |
Lee et al. | Effective PPG sensor placement for reflected red and green light, and infrared wristband-type photoplethysmography | |
CN108113649A (zh) | 生物信号处理设备和生物测量学信息检测设备及方法 | |
CN107908928A (zh) | 一种基于深度学习技术的血红蛋白动态光谱分析预测方法 | |
CN110141249B (zh) | 基于ppg信号的无创血糖监测方法、系统、设备及介质 | |
Tekeste et al. | A nano-watt ECG feature extraction engine in 65-nm technology | |
CN107106080A (zh) | 呼吸状态估计装置、便携式设备、佩戴型仪器、程序、介质、呼吸状态估计方法以及呼吸状态估计器 | |
CN106798547A (zh) | 一种基于统计的臂踝脉搏波延时的计算方法 | |
CN109189207A (zh) | 一种可穿戴设备穿戴状态的检测方法及装置 | |
CN108992054A (zh) | 一种脉搏信号峰值点检测方法及装置 | |
Rinkevičius et al. | Photoplethysmogram signal morphology-based stress assessment | |
CN106073800B (zh) | 基于绝对差值和提取的动态光谱数据处理方法及其装置 | |
CN108652610A (zh) | 一种多人心跳的非接触式检测方法 | |
Sabbadini et al. | Atrial fibrillation detection by means of edge computing on wearable device: a feasibility assessment | |
Zheng et al. | Adaptive notch-filtration to effectively recover photoplethysmographic signals during physical activity | |
CN104182956B (zh) | 一种血管提取方法 | |
Ferdinando et al. | Photoplethysmography signal analysis to assess obesity, age group and hypertension | |
Wang et al. | Analysis of pulse waveforms preprocessing | |
CN107837083A (zh) | 基于最小二乘支持向量机的j波自动检测方法 | |
CN108664956A (zh) | 人体静脉生物特征信息采集过程中的活体识别方法及装置 | |
CN107981862A (zh) | 一种核磁共振环境下的脑信号在线去噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |