一种非接触式心率变异性测量方法技术

本发明专利技术公开了一种非接触式心率变异性测量方法，属于人体心理压力的评价技术领域。首先使用高速彩色摄像机非接触式的采集检测者面部的区域视频图像，根据每帧视频图像，提取人体面部特征点，并根据特征点选取各图像的目标区域；然后对每帧视频图像的目标区域内的所有像素点的红色、绿色和蓝色三个通道的信号进行独立成分分析，并选取心率相关分量k(t)，对其进行频谱分析和滤波处理，提取脉搏波特征点，得到逐拍心动周期，计算心率变异性的时域和频域参数。最后根据心率变异性的时域和频域参数对心理压力进行评价和分析。本发明专利技术可以无接触式的测量心率变异性，其测量精度高，减少了对检测对象的影响，获得心理压力的评价较为客观。

【技术实现步骤摘要】
一种非接触式心率变异性测量方法
本专利技术涉及一种非接触式心率变异性测量方法，应用在人体心理压力的评价

技术介绍
目前心率变异性的测量主要基于心电信号，即通过心电计算心动周期，并根据逐拍心动周期计算心率变异性。由于心电测量需要与人体接触，其应用受到很多局限，例如被试者可能会抗拒心电测量。目前非接触式的心率测量，通过摄像机拍摄人面部视频获得；如申请号：201310172275.9公开的非接触式心率测量系统及测量方法，提出当拍摄视频后，通过对面部图像的多个颜色分量进行盲源分离，筛选出与心率相关的分量，然后通过频谱分析得到心率。但是在对人体的疲劳程度、心理状态、压力大小等进行评价和分析时，仅有心率是不够的，还需要进行心率变异性的测量。上述方法在测量心率变异性时存在缺陷，主要体现在：1、心率分析需要依赖频谱分析，但是频谱分析只能反映较长时间内的平均心率，而不能反映这一段时间内心率的变化，即心率变异性。2、使用盲源分离后筛选出的与心率相关的信号分量来计算心率时，存在较大的噪声和干扰；由于计算心率变异性比心率更易受噪声的干扰，盲源分离后分量中的噪声虽然对计算心率带来的误差很小，但是会对计算心率变异性引入巨大的误差，甚至错误。3、普通的摄像机采集图像的速率为每秒钟25-30帧，即采样率为25-30Hz；而心率变异性测量要求的时间精度为20ms，根据奈奎斯特采样定律，采样率不能低于100Hz；因此，普通摄像机难以满足心率变异性测量的精度要求。
技术实现思路
本专利技术针对现有的非接触式测量方法难以满足心率变异性测量要求的缺陷，提出一种非接触式心率变异性测量方法。克服了上述三方面的缺陷：在心率分析时不完全依赖于频谱分析，实现了心率变异性的计算；对盲源分离后的分量进行了进一步的滤波，提高了心率变异性计算的准确性；使用高速摄像机，采样频率可以满足心率变异性测量精度要求。具体步骤如下：步骤一、针对某检测者，使用高速彩色摄像机非接触式的采集该检测者面部的区域视频图像；采集过程在明亮的环境中进行，摄像机采集速率不低于每秒钟100帧，采集过程不短于3分钟。步骤二、根据每帧视频图像，分别提取人体面部特征点，并根据特征点选取各图像的目标区域；对每一帧图像进行人脸识别，并提取特征点；提取过程为：用视频图像中的每个点分别与人脸特征Dlib跨平台计算机视觉库中的特征点进行对比，衡量这两个点周围的领域图像的相似度，最后，提取视频图像中相似度最高的点，即为人脸图像中的特征点。提取的特征点包括每侧脸颊边缘各2个、鼻翼每侧各1个以及鼻梁顶端1个；将每侧面部鼻梁顶端、鼻翼一侧及脸颊边缘的特征点依次相连，得到一个四边形；人脸两侧的四边形区域定义为目标区域。步骤三、对每帧视频图像的目标区域内的所有像素点的红色、绿色和蓝色三个通道的信号进行独立成分分析，并选取心率相关分量k(t)；具体步骤为：步骤301，针对每帧视频图像，对目标区域内所有像素点的红色、绿色和蓝色的通道信号值分别求取平均值。即：r(i)为目标区域内的第i个像素点的红色通道信号值；g(i)为目标区域内的第i个像素点的绿色通道信号值；b(i)为目标区域内的第i个像素点的蓝色通道信号值。为目标区域内的所有像素点的红色通道信号值的平均值；为目标区域内的所有像素点的绿色通道信号值的平均值；为目标区域内的所有像素点的蓝色通道信号值的平均值；步骤302，将所有视频图像计算所得的三个通道信号强度平均值分别按视频采集的时间排列，得到三个时间序列r(t),g(t),b(t)。r(t)是所有视频图像的红色通道信号的平均值按视频采集的时间排列得到的时间序列；g(t)是所有视频图像的绿色通道信号的平均值按视频采集的时间排列得到的时间序列；b(t)是所有视频图像的蓝色通道信号的平均值按视频采集的时间排列得到的时间序列；步骤303，对三个时间序列进行独立成分分析ICA，分解得到三个独立成分分量，记为i(t),j(t),k(t)；计算公式如下：其中A是一个3*3的矩阵，通过使i(t),j(t),k(t)相关性最小的最优化算法获得。步骤304，将三个独立成分分量分别与g(t)分量进行相关性分析，把相关性最大的独立成分分量选出，并命名为心率相关分量k(t)。相关性最大的分量包含脉搏波信号。步骤四、对心率相关分量k(t)进行频谱分析和滤波处理，提取脉搏波特征点，得到逐拍心动周期；首先，对心率相关分量k(t)的频谱分析，提取0.5-3Hz的频段中能量最大点所对应的频率，得到心率频率fc；然后，根据心率频率fc，选取截止频率为(fc-0.3)～(fc+0.3)Hz的带通滤波器，对心率相关分量k(t)进行频域滤波，结果记为kc(t)。进而，对滤波结果kc(t)进行脉搏波特征点提取，并依次记录脉搏波峰值点出现的时刻p1,p2,……，pn。提取出现频率约为fc的局部极大值点作为脉搏波峰值点。最后，根据相邻极大值点出现的时刻计算逐拍心动周期y＝{y1,y2,……,yn-1}。计算公式为：y1＝p2-p1,y2＝p3-p2,……,yn-1＝pn-pn-1；步骤五、根据逐拍心动周期计算心率变异性的时域和频域参数。首先，对逐拍心动周期包含的随机噪声进行平滑滤波；针对当前时刻的逐拍心动周期yn，平滑滤波公式如下：j表示采样间隔的数量；然后，将极大值出现的时刻p1,p2,……，pn作为横坐标，将对应的平滑滤波后的逐拍心动周期y1,y2,……,yn-1作为纵坐标，生成新的时间序列y(p)。进而，对时间序列y(p)进行三次样条曲线插值，使其采样频率等于视频采集的频率fs，插值后的时间序列记为y(t)，称为心率变异性信号。最后，对心率变异性信号y(t)进行时域和频域分析，得到心率变异性参数。对心率变异性信号进行时域分析，得到参数：MeanNN、SDNN、RMSSD和SDSD；其中，N是参与心率变异性计算的心动周期的个数；μ为相邻y(n)差值的均值。y(n)是心率变异性信号y(t)中n时刻的心动周期值，y(n+1)表示心率变异性信号y(t)中n+1时刻的心动周期值。对心率变异性信号进行频域分析，得到pLF、pHF、LFHFRatio、VLF、LF和HF。其中，对心率变异性信号y(t)进行傅里叶变换，频谱图当中的0.003-0.04Hz的能量综合为VLF；0.04-0.15Hz的能量总和为LF；0.15-0.4Hz的能量总和为HF；其它参数计算方法如下：pLF＝LF/(LF+HF)pHF＝HF/(LF+HF)LFHFRatio＝LF/HF步骤六、根据心率变异性的时域和频域参数对心理压力进行评价和分析。本专利技术的优点在于，可以无接触式的测量心率变异性，其测量精度可与心电测量所得的心率变异性相当；由于非接触式的特点，可以减少对检测对象的影响，获得较为客观的心理压力的评价。附图说明图1为本专利技术一种非接触式心率变异性测量方法的整体流程示意图；图2为本专利技术提取的人脸特征点和目标区域示意图；图3为本专利技术的心率相关分量及其频谱示意图；图4为本专利技术的心率相关分量滤波后脉搏波特征点提取结果。具体实施方式下面将结合附图和实施例对本专利技术作进一步的详细说明。本专利技术提供了一种非接触式心率变异性测量方法，可在受试者无察觉的情况下，对其心率变异性进行测量，从而较为真实的反映其心理本文档来自技高网...

【技术保护点】
1.一种非接触式心率变异性测量方法，其特征在于，具体步骤如下：步骤一、针对某检测者，使用高速彩色摄像机非接触式的采集该检测者面部的区域视频图像；步骤二、根据每帧视频图像，分别提取人体面部特征点，并根据特征点选取各图像的目标区域；步骤三、对每帧视频图像的目标区域内的所有像素点的红色、绿色和蓝色三个通道的信号进行独立成分分析，并选取心率相关分量k(t)；具体步骤为：步骤301，针对每帧视频图像，对目标区域内所有像素点的红色、绿色和蓝色的通道信号值分别求取平均值；即：

【技术特征摘要】
1.一种非接触式心率变异性测量方法，其特征在于，具体步骤如下：步骤一、针对某检测者，使用高速彩色摄像机非接触式的采集该检测者面部的区域视频图像；步骤二、根据每帧视频图像，分别提取人体面部特征点，并根据特征点选取各图像的目标区域；步骤三、对每帧视频图像的目标区域内的所有像素点的红色、绿色和蓝色三个通道的信号进行独立成分分析，并选取心率相关分量k(t)；具体步骤为：步骤301，针对每帧视频图像，对目标区域内所有像素点的红色、绿色和蓝色的通道信号值分别求取平均值；即：r(i)为目标区域内的第i个像素点的红色通道信号值；g(i)为目标区域内的第i个像素点的绿色通道信号值；b(i)为目标区域内的第i个像素点的蓝色通道信号值；为目标区域内的所有像素点的红色通道信号值的平均值；为目标区域内的所有像素点的绿色通道信号值的平均值；为目标区域内的所有像素点的蓝色通道信号值的平均值；步骤302，将所有视频图像计算所得的三个通道信号强度平均值分别按视频采集的时间排列，得到三个时间序列r(t),g(t),b(t)；r(t)是所有视频图像的红色通道信号的平均值按视频采集的时间排列得到的时间序列；g(t)是所有视频图像的绿色通道信号的平均值按视频采集的时间排列得到的时间序列；b(t)是所有视频图像的蓝色通道信号的平均值按视频采集的时间排列得到的时间序列；步骤303，对三个时间序列进行独立成分分析ICA，分解得到三个独立成分分量，记为i(t),j(t),k(t)；计算公式如下：其中A是一个3*3的矩阵，通过使i(t),j(t),k(t)相关性最小的最优化算法获得；步骤304，将三个独立成分分量分别与g(t)分量进行相关性分析，把相关性最大的独立成分分量选出，并命名为心率相关分量k(t)；相关性最大的分量包含脉搏波信号；步骤四、对心率相关分量k(t)进行频谱分析和滤波处理，提取脉搏波特征点，得到逐拍心动周期；步骤五、根据逐拍心动周期计算心率变异性的时域和频域参数；步骤六、根据心率变异性的时域和频域参数对心理压力进行评价和分析。2.如权利要求1所述的一种非接触式心率变异性测量方法，其特征在于，所述的步骤二中对每一帧图像进行人脸识别，并提取特征点；提取过程为：用视频图像中的每个点分别与人脸特征Dlib跨平台计算机视觉库中的特征点进行对比，衡量这两个点周围的领域图像的相似度，最后，提取视频图像中相似度最高的点，即为人脸图像中的特征点；提取的特征点...

【专利技术属性】
技术研发人员：张弛，齐子杰，王维，罗凯，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11