一种非接触式血压测量方法及系统与流程

本发明涉及血压测量领域，特别是涉及一种非接触式血压测量方法及系统。
背景技术：
：血压(bloodpressure，bp)是指血液在血管内流动时作用于单位面积血管壁的侧压力，它是推动血液在血管内流动的动力。血压是一项很重要的生理指标，会随着人的生理活动变化而变化，例如当人情绪激动时会引起血压升高。目前普遍使用的血压测量设备主要基于示波法和柯式音法，基于这两种原理方法的血压测量设备，需要使用袖带对人体血管进行压迫然后才能测量到人体血压值，依赖于设备对血管的接触和压迫。技术实现要素：本发明的目的是提供一种非接触式血压测量方法及系统，避免对人体的直接接触和压迫。为实现上述目的，本发明提供了如下方案：一种非接触式血压测量方法，包括：获取用户的面部视频；从所述面部视频中筛选第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧；所述第一感兴趣区域与所述第二感兴趣区域为每个用户的面部中血管分布量超过预设值且互不重叠的两个区域；分别对所述第一感兴趣区域的多个图像帧和所述第二感兴趣区域的多个图像帧提取三基色的视频迹曲线；所述视频迹曲线为按图像帧在视频中的时间顺序经过空间滤波形成的曲线；分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行盲源分离，得到心率的估计值；以所述心率的估计值为先验知识，通过盲源信号提取确定第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线；根据所述第一感兴趣区域的光电容积脉搏波曲线和所述第二感兴趣区域的光电容积脉搏波曲线求取ppg信号时滞；将所述ppg信号时滞输入时滞与血压的关系方程得到血压值。可选的，所述时滞与血压的关系方程的构建方法包括：获取单个用户的多个面部视频以及每个面部视频对应的实际血压值，得到多个样本视频和多个实际血压值；从每个所述样本视频中筛选第一感兴趣样本区域的多个图像帧和第二感兴趣样本区域的多个图像帧；所述第一感兴趣样本区域和所述第二感兴趣样本区域为提供样本视频的用户面部中血管分布量超过预设值且互不重叠的两个区域；分别对所述第一感兴趣样本区域的多个图像帧和所述第二感兴趣样本区域的多个图像帧提取三基色的样本视频迹曲线；所述样本视频迹曲线为按图像帧在样本视频中的时间顺序经过空间滤波形成的曲线；分别对所述第一感兴趣样本区域的三基色的样本视频迹曲线和所述第二感兴趣样本区域的三基色的样本视频迹曲线进行盲源分离，得到样本的心率估计值；以所述样本的心率估计值为先验知识，通过盲源信号提取确定第一感兴趣样本区域的光电容积脉搏波曲线和第二感兴趣样本区域的光电容积脉搏波曲线；根据所述第一感兴趣样本区域的光电容积脉搏波曲线和所述第二感兴趣样本区域的光电容积脉搏波曲线求取样本ppg信号时滞；建立所述样本ppg信号时滞与所述ppg信号时滞对应的实际血压值之间的线性回归方程，得到时滞与血压的关系方程。可选的，所述以所述心率的估计值为先验知识，通过盲源信号提取确定第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线，具体包括：分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行带通滤波和三次样条插值，得到第一感兴趣区域的三基色插值后曲线和第二感兴趣区域的三基色插值后曲线；以所述心率的估计值为先验知识，分别对所述第一感兴趣区域的三基色插值后曲线和所述第二感兴趣区域的三基色插值后曲线进行盲源信号提取，得到第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线。可选的，所述从所述面部视频中筛选第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧，具体包括：根据所述面试视频的视频图像分析面部的血管分布，确定血管分布量超过预设值且中心距离超过预设距离的两个互不重叠区域，得到第一感兴趣区域和第二感兴趣区域；选取包含所述第一感兴趣区域和所述第二感兴趣区域的一帧图像作为基准帧；以所述基准帧为起始帧，选取所述基准帧之后的连续多帧图像；以所述基准帧为模板，利用模板匹配算法筛选所述连续多帧图像中的第一感兴趣区域和第二感兴趣区域，得到第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧。本发明还公开一种非接触式血压测量系统，包括：面部视频获取模块，用于获取用户的面部视频；图像帧筛选模块，用于从所述面部视频中筛选第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧；所述第一感兴趣区域与所述第二感兴趣区域为每个用户的面部中血管分布量超过预设值且互不重叠的两个区域；三基色视频迹曲线提取模块，用于分别对所述第一感兴趣区域的多个图像帧和所述第二感兴趣区域的多个图像帧提取三基色的视频迹曲线；所述视频迹曲线为按图像帧在视频中的时间顺序经过空间滤波形成的曲线；盲源分离模块，用于分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行盲源分离，得到心率的估计值；脉搏波曲线确定模块，用于以所述心率的估计值为先验知识，通过盲源信号提取确定第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线；时滞计算模块，用于根据所述第一感兴趣区域的光电容积脉搏波曲线和所述第二感兴趣区域的光电容积脉搏波曲线求取ppg信号时滞；血压确定模块，用于将所述ppg信号时滞输入时滞与血压的关系方程得到血压值。可选的，该系统还包括方程构建模块，用于构建所述时滞与血压的关系方程；所述方程构建模块包括：视频与血压获取单元，用于获取单个用户的多个面部视频以及每个面部视频对应的实际血压值，得到多个样本视频和多个实际血压值；样本图像帧筛选单元，用于从每个所述样本视频中筛选第一感兴趣样本区域的多个图像帧和第二感兴趣样本区域的多个图像帧；所述第一感兴趣样本区域和所述第二感兴趣样本区域为提供样本视频的用户面部中血管分布量超过预设值且互不重叠的两个区域；样本三基色视频迹曲线提取单元，用于分别对所述第一感兴趣样本区域的多个图像帧和所述第二感兴趣样本区域的多个图像帧提取三基色的样本视频迹曲线；所述样本视频迹曲线为按图像帧在样本视频中的时间顺序经过空间滤波形成的曲线；样本盲源分离单元，用于分别对所述第一感兴趣样本区域的三基色的样本视频迹曲线和所述第二感兴趣样本区域的三基色的样本视频迹曲线进行盲源分离，得到样本的心率估计值；样本脉搏波曲线确定单元，用于以所述心率样本的估计值为先验知识，通过盲源信号提取确定第一感兴趣样本区域的光电容积脉搏波曲线和第二感兴趣样本区域的光电容积脉搏波曲线；样本ppg信号时滞计算单元，用于根据所述第一感兴趣样本区域的光电容积脉搏波曲线和所述第二感兴趣样本区域的光电容积脉搏波曲线求取样本ppg信号时滞；线性回归方程建立单元，用于建立所述样本ppg信号时滞与所述ppg信号时滞对应的实际血压值之间的线性回归方程，得到时滞与血压的关系方程。可选的，所述脉搏波曲线确定模块包括：滤波与插值单元，用于分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行带通滤波和三次样条插值，得到第一感兴趣区域的三基色插值后曲线和第二感兴趣区域的三基色插值后曲线；盲源信号提取单元，用于以所述心率的估计值为先验知识，分别对所述第一感兴趣区域的三基色插值后曲线和所述第二感兴趣区域的三基色插值后曲线进行盲源信号提取，得到第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线。可选的，所述图像帧筛选模块包括：感兴趣区域确定单元，用于根据所述面试视频的视频图像分析面部的血管分布，确定血管分布量超过预设值且中心距离超过预设距离的两个互不重叠区域，得到第一感兴趣区域和第二感兴趣区域；基准帧确定单元，用于选取包含所述第一感兴趣区域和所述第二感兴趣区域的一帧图像作为基准帧；多帧图像提取单元，用于以所述基准帧为起始帧，选取所述基准帧之后的连续多帧图像；模板匹配单元，用于以所述基准帧为模板，利用模板匹配算法筛选所述连续多帧图像中的第一感兴趣区域和第二感兴趣区域，得到第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧。根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明所公开的非接触式血压测量方法及系统，通过拍摄用户的面部视频来实现血压的测量，从而避免设备对人体的直接接触和压迫，提高血压测量的便捷性。同时，避免设备在人体的安装和拆卸过程，用户只需正对镜头即可，提高了测量效率。附图说明为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本发明实施例1非接触式血压测量方法的方法流程图；图2为本发明实施例2非接触式血压测量方法中感兴趣区域示意图；图3为本发明实施例3非接触式血压测量系统的系统结构图。具体实施方式下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。本发明的原理是基于脉搏波传播速度与血压收缩压之间的相关关系，由于人脸皮下组织的毛细血管含有丰富的光容积变化信号，血管光容积的变化通过环境光反射被摄像机所记录，通过算法计算可以提取出血管的随时间变化的光容积变化信号，即光电容积脉搏波ppg信号，然后利用脉搏波传播速度与血压收缩压之间的相关关系进行血压收缩压的估计。实施例1：图1为本发明实施例1非接触式血压测量方法的方法流程图。参见图1，该非接触式血压测量方法，包括：步骤101：获取用户的面部视频；步骤102：从所述面部视频中筛选第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧；所述第一感兴趣区域与所述第二感兴趣区域为每个用户的面部中血管分布量超过预设值且互不重叠的两个区域；步骤103：分别对所述第一感兴趣区域的多个图像帧和所述第二感兴趣区域的多个图像帧提取三基色的视频迹曲线；所述视频迹曲线为按图像帧在视频中的时间顺序经过空间滤波形成的曲线；步骤104：分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行盲源分离，得到心率的估计值；步骤105：以所述心率的估计值为先验知识，通过盲源信号提取确定第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线；步骤106：根据所述第一感兴趣区域的光电容积脉搏波曲线和所述第二感兴趣区域的光电容积脉搏波曲线求取ppg信号时滞；步骤107：将所述ppg信号时滞输入时滞与血压的关系方程得到血压值。作为一种可选的实施方式，步骤107中时滞与血压的关系方程的构建方法包括：获取单用户的多个面部视频以及每个面部视频对应的实际血压值，得到多个样本视频和多个实际血压值；从每个所述样本视频中筛选第一感兴趣样本区域的多个图像帧和第二感兴趣样本区域的多个图像帧；所述第一感兴趣样本区域和所述第二感兴趣样本区域为提供样本视频的用户面部中血管分布量超过预设值且互不重叠的两个区域；分别对所述第一感兴趣样本区域的多个图像帧和所述第二感兴趣样本区域的多个图像帧提取三基色的样本视频迹曲线；所述样本视频迹曲线为按图像帧在样本视频中的时间顺序经过空间滤波形成的曲线；分别对所述第一感兴趣样本区域的三基色的样本视频迹曲线和所述第二感兴趣样本区域的三基色的样本视频迹曲线进行盲源分离，得到样本的心率估计值；以所述样本的心率估计值为先验知识，通过盲源信号提取确定第一感兴趣样本区域的光电容积脉搏波曲线和第二感兴趣样本区域的光电容积脉搏波曲线；根据所述第一感兴趣样本区域的光电容积脉搏波曲线和所述第二感兴趣样本区域的光电容积脉搏波曲线求取样本ppg信号时滞；建立所述样本ppg信号时滞与所述ppg信号时滞对应的实际血压值之间的线性回归方程，得到时滞与血压的关系方程。作为一种可选的实施方式，步骤105具体包括：分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行带通滤波和三次样条插值，得到第一感兴趣区域的三基色插值后曲线和第二感兴趣区域的三基色插值后曲线；以所述心率的估计值为先验知识，分别对所述第一感兴趣区域的三基色插值后曲线和所述第二感兴趣区域的三基色插值后曲线进行盲源信号提取，得到第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线。作为一种可选的实施方式，步骤102具体包括：根据所述面试视频的视频图像分析面部的血管分布，确定血管分布量超过预设值且中心距离超过预设距离的两个互不重叠区域，得到第一感兴趣区域和第二感兴趣区域；选取包含所述第一感兴趣区域和所述第二感兴趣区域的一帧图像作为基准帧；以所述基准帧为起始帧，选取所述基准帧之后的连续多帧图像；以所述基准帧为模板，利用模板匹配算法筛选所述连续多帧图像中的第一感兴趣区域和第二感兴趣区域，得到第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧。实施例2：该实施例2中，由于拍摄的面部视频含有人脸抖动、环境光变化以及成像传感器的噪声，盲源分离算法(blindsourceseparation，bss)或者经验模态分解(empiricalmodedecomposition，emd)不足以分离出精度足够的光电容积脉搏波ppg信号用于血压的预估，因而引入了盲源提取算法(blindsourceextraction，bse)，用于提取高精度的光电容积脉搏波ppg信号，用于血压值的预估。步骤一、首先需要构建脉搏波时滞与血压的关系方程。拍摄人脸面部视频，得到多个样本视频；测量每个面部视频所属用户的实际血压值，得到多个实际血压值。拍摄过程中保持环境光没有明显变化、人脸保持不动以及摄像机不要发生移动，从而在数据采集阶段减少部分误差。步骤二、以矩形框形状，使矩形框中包含有对比度较高的部分，以提高后续模板匹配算法的准确性。图2为本发明实施例2非接触式血压测量方法中感兴趣区域示意图。参见图2，以虚线矩形框选中的区域作为感兴趣区域roi的第一帧，两个感兴趣区域分别记为第一感兴趣样本区域roi1和第二感兴趣样本区域roi2；两个roi的选择依据的是血管在面部的分布情况，脉搏流动在这两个roi引发的光学信号更易于测量。步骤三、以roi1和roi2为基准帧，使用模板匹配算法选择后续的1s-10s(优选为3s)时间长度的视频帧，以矩形框的形状，参考图2中的实线矩形框，选取提取ppg信号的区域的视频数据，记为第一感兴趣样本区域的多个图像帧video1和第二感兴趣样本区域的多个图像帧video2；第一感兴趣样本区域的多个图像帧video1为第一感兴趣样本区域内的某个目标区域的连续多帧图像。第二感兴趣样本区域的多个图像帧video2为第二感兴趣样本区域内的某个目标区域的连续多帧图像。模板匹配算法使用归一化相关系数匹配法，相似程度计算公式如下：上述公式中，t表示模板图像，i表示下一帧要匹配的图像，w是模板图像的宽，h是模板图像的高，r表示匹配的结果分数，x和y是下一帧图要匹配的图像位置的坐标，x'和y'为当前要匹配的模版和待匹配的模板要遍历的图像像素坐标。步骤四、以video1为原始数据，求取样本视频迹曲线，得到第一感兴趣样本区域的三基色的样本视频迹曲线r1、g1、b1，以video2为原始数据，求取视频迹曲线，得到第二感兴趣样本区域的三基色的样本视频迹曲线r2、g2、b2；其中三基色视频迹样本曲线的计算公式如下：上述公式中，分子代表t时刻时感兴趣区域的像素值总和，分母代表t时刻时感兴趣区域的像素个数。roin为感兴趣区域的编号，和表示感兴趣区域中要遍历的坐标。步骤五、对步骤四中得到的三路样本视频迹曲线rroin、groin、broin进行盲源分离(bss)，进而计算获得样本的心率估计值：s501：首先对三路样本视频迹曲线做白化处理；其中v是白化矩阵，wi1、wi2、wi3是三路样本视频迹曲线白化后的表示。其中d是cx的特征组成的对角矩阵，p是cx的单位特征向量，e表示求期望。cx表示信号的协方差矩阵。s502：对s501中计算得到的wi1、wi2、wi3进行盲源分离，使用基于互信息极小化准则，其目标函数为：其中线性变换公式为seperate＝dm·sourceh(en)＝-e[logp(en)]上述公式中，info(seperate)衡量分离信号的互信息与原始信号之间互信息值的差异，dm表示解混矩阵，seperate表示解混后的三路信号，seperatej表示解混后的三路信号其中一路信号，source表示解混前的信号即wi1、wi2、wi3，en是待求信息熵的输入信号，即seperate和source。h(·)为信息熵，p(·)为待求信息熵的输入信号的概率分布。s503：经过s502中盲源分离后得到的三路信号为独立成分，从r1、g1、b1分别得到独立成分c11、c12、c13，从r2、g2、b2分别得到独立成分c21、c22、c23。在c11、c12、c13里筛选出与g1的皮尔逊相关系数最高的信号，记为c1；在c21、c22、c23里筛选出与g2的皮尔逊相关系数最高的信号，记为c2；其中相关系数计算公式为：上述公式中，ρc,g表示相关系数；c为所计算的独立成分，g为对应的绿色通道数据，e代表求期望，μx表示c的均值，μy表示g的均值。s504：使用巴特沃斯带通滤波器(butterworth带通滤波器)，对c1和c2进行滤波，分别求取c1和c2的频谱峰值，记为m1和m2，计算作为心率样本的估计值。步骤六、使用butterworth带通滤波器对步骤四中得到的r1、g1、b1、r2、g2、b2进行滤波，对滤波结果使用三次样条算法进行插值，得到第一感兴趣样本区域的三基色插值后样本曲线r1f、g1f、b1f和第二感兴趣样本区域的三基色插值后样本曲线r2f、g2f、b2f；其中三次样条插值求解公式为：在每个子区间zk≤z≤zk+1中，插值方程为gk(z)＝ak+bk(z-zk)+ck(z-zk)2+dk(z-zk)3ak＝yk上述公式中，ak、bk、ck、dk为插值方程的四个系数，hk表示步长，mk为二次微分值，通过求解矩阵方程获得。z表示要进行插值的时间坐标，zk表示要进行插值的坐标所在子区间的左端点，zk+1表示要进行插值的坐标所在子区间的右端点，yk表示要进行插值的坐标所在子区间的左端点的数据值，字母k表示插值数据所分成子区间的区间编号。步骤七、以步骤s503求得的样本心率的估计值为先验知识，分别对r1f、g1f、b1f和r2f、g2f、b2f利用盲源提取算法进行盲源信号提取，得到的信号提取结果记为第一感兴趣区域的样本脉搏波曲线res1和第二感兴趣区域的样本脉搏波曲线res2；其中该步骤中采用的盲源提取算法为改进后的盲源提取算法，引入自动更新权重算子，有机结合了基于信号间隔的盲源提取算法和基于负熵的盲源提取算法。改进后的盲源提取算法的目标函数和迭代过程如下：j(w)＝(1-ρ)e{q(wt(x(t)-bx(t-τ)))}-ρ[e{q(wtx(t))}-e{g(v)}]2w←w-μw((1-ρ)e{q′(wt(x(t)-bx(t-τ)))}(x(t)-bx(t-τ))-2ρ(e{q(wtx(t))}-e{g(v)})e{q′(wtx(t))}x(t))b←b+μb((1-ρ)e{q′(wt(x(t)-bx(t-τ)))}x(t-τ))ρ←ρ-μρ(-e{q(wt(x(t)-bx(t-τ)))}-[e{q(wtx(t))}-e{g(v)}]2)上述公式中，j(w)是目标函数值，用于度量提取信号的好坏程度；q(·)和g(·)均为可微函数，g(v)是标准正态分布的数据的非二次函数值，v是零均值和单位方差的高斯变量，q′(·)是函数q(·)的导数，w是解混向量，e{·}是求期望，ρ是自动更新权重算子，μw、μb、μρ为各个参数的学习率。b是信号间隔衰减系数，x(t)是原始信号，τ是信号时延(间隔)，τ可由先验知识心率估计值得到，转换公式为：τ＝60/hr上述公式hr代表心率估计值。步骤八、计算步骤七中的res1和res2的相关图，求取一个时延周期之内的峰值，作为ppg信号之间的时延。步骤九、由于选取的roi区域不变，所以两个roi之间的距离不会发生改变，信号之间的时滞可以等同反映信号的传播速度。因为脉搏波波速与人体收缩压有线性相关的关系，而时滞可以反映ppg信号的传播速度，所以可以通过信号之间的时滞来对血压收缩压进行回归。收集数据样本，通过线性回归构建时滞与血压关系方程，如表1为样本数据中的几个样本ppg信号时滞与实际血压值的对应关系表。表1时滞血压τ1y1τ2y2τ3y3τ4y4回归方程为：hθ(τ)＝θ0+θ1τ求解的损失函数为：上述公式中，θ0和θ1是回归方程的两个系数，τ表示输入的时滞，nn为样本集里样本个数，τc是第c个样本的时滞，hθ(τc)表示第c个样本预测的血压值，c为遍历样本集是当前计算的样本标号，利用梯度下降或者极值法求取最优参数θ0和θ1。步骤一至步骤九为用于利用样本数据构建时滞与血压关系方程的过程。从步骤十开始为实际血压的测量过程。步骤十、在要测取血压的视频中，以矩形框形状，参见图2，使矩形框中包含有对比度较高的部分，分别截取两个感兴趣区域的第一帧，分别记为第一感兴趣区域roi3和第二感兴趣区域roi4；步骤十一、以roi3和roi4为基准帧，使用与步骤三相同的模板匹配算法选择后续的的1s-10s(优选为3s)时间长度的视频帧，以矩形框的形状，参考图2中的实线矩形框，选取提取ppg信号的区域的视频数据，记为第一感兴趣区域的多个图像帧video3和第二感兴趣区域的多个图像帧video4；步骤十二、以video3为原始数据，求取视频迹曲线，得到第一感兴趣区域的三基色视频迹曲线r3、g3、b3，以video4为原始数据，求取视频迹曲线，得到第二感兴趣区域的三基色视频迹曲线r4、g4、b4。其中r、g、b的计算公式与步骤四相同；步骤十三、对步骤十二中得到的三路视频迹曲线ri、gi、bi进行盲源分离(bss)，进而计算获得心率的估计值，计算过程与步骤五相同；步骤十四、使用butterworth带通滤波器对步骤十二中得到的r3、g3、b3、r4、g4、b4进行滤波，对滤波结果使用三次样条算法进行插值，得到得到第一感兴趣区域的三基色插值后曲线r3f、g3f、b3f和第二感兴趣区域的三基色插值后曲线r4f、g4f、b4f；步骤十五、以步骤十三求得的心率的估计值为先验知识，分别对r3f、g3f、b3f和r4f、g4f、b4f利用盲源提取算法进行盲源信号提取，得到的信号提取结果记为第一感兴趣区域的光电容积脉搏波曲线res3和第二感兴趣区域的光电容积脉搏波曲线res4；步骤十六、计算步骤十五中的res3和res4相关图，求取一个时延周期之内的峰值，作为两路ppg信号之间的时滞τ。步骤十七、将步骤十六得到的信号时滞τ输入步骤九获得的时滞与血压关系方程，求得血压的收缩压预估值。血压回归方程：hθ(τ)＝θ0+θ1τ其中θ0、θ1为步骤九中所求解结果。本发明的非接触式血压测量方法采用非接触式的测量方式，能够使血压的测量过程免去压迫血管的不适，避免由于操作经验不足而导致的测量误差，以及可以实现血压的长时间连续监控。并且，引入盲源提取算法实现非接触式的光电容积脉搏波ppg信号的提取，计算量少，提取结果精度高，算法鲁棒性高；本申请的盲源提取算法有机结合了基于信号时延的盲源提取算法和基于负熵的盲源提取算法，并且引入一个自动权衡算子，自动权衡两个目标函数的权重，进一步提高分离结果的精度，提高算法的鲁棒性；即使拍摄的面部视频有环境光变化的噪声和人脸运动噪声，使用盲源提取算法也可以分离出精度较高的光电容积脉搏波ppg信号。实施例3：图3为本发明实施例3非接触式血压测量系统的系统结构图。参见图3，该非接触式血压测量系统，包括：面部视频获取模块301，用于获取用户的面部视频；图像帧筛选模块302，用于从所述面部视频中筛选第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧；所述第一感兴趣区域与所述第二感兴趣区域为每个用户的面部中血管分布量超过预设值且互不重叠的两个区域；三基色视频迹曲线提取模块303，用于分别对所述第一感兴趣区域的多个图像帧和所述第二感兴趣区域的多个图像帧提取三基色的视频迹曲线；所述视频迹曲线为按图像帧在视频中的时间顺序经过空间滤波形成的曲线；盲源分离模块304，用于分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行盲源分离，得到心率的估计值；脉搏波曲线确定模块305，用于以所述心率的估计值为先验知识，通过盲源信号提取确定第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线；时滞计算模块306，用于根据所述第一感兴趣区域的光电容积脉搏波曲线和所述第二感兴趣区域的光电容积脉搏波曲线求取ppg信号时滞；血压确定模块307，用于将所述ppg信号时滞输入时滞与血压的关系方程得到血压值。可选的，该血压测量系统还包括方程构建模块308，用于构建所述时滞与血压的关系方程；所述方程构建模块308包括：视频与血压获取单元，用于获取单用户的多个面部视频以及每个面部视频对应的实际血压值，得到多个样本视频和多个实际血压值；样本图像帧筛选单元，用于从每个所述样本视频中筛选第一感兴趣样本区域的多个图像帧和第二感兴趣样本区域的多个图像帧；所述第一感兴趣样本区域和所述第二感兴趣样本区域为提供样本视频的用户面部中血管分布量超过预设值且互不重叠的两个区域；样本三基色曲线提取单元，用于分别对所述第一感兴趣样本区域的多个图像帧和所述第二感兴趣样本区域的多个图像帧提取三基色的样本视频迹曲线；所述样本视频迹曲线为按图像帧在样本视频中的时间顺序经过空间滤波形成的曲线；样本盲源分离单元，用于分别对所述第一感兴趣样本区域的三基色的视频迹样本曲线和所述第二感兴趣样本区域的三基色的视频迹样本曲线进行盲源分离，得到样本的心率估计值；样本脉搏波曲线确定单元，用于以所述样本的心率估计值为先验知识，通过盲源信号提取确定第一感兴趣样本区域的光电容积脉搏波曲线和第二感兴趣样本区域的光电容积脉搏波曲线；样本ppg信号时滞计算单元，用于根据所述第一感兴趣样本区域的光电容积脉搏波曲线和所述第二感兴趣样本区域的光电容积脉搏波曲线求取样本ppg信号时滞；线性回归方程建立单元，用于建立所述样本ppg信号时滞与所述ppg信号时滞对应的实际血压值之间的线性回归方程，得到时滞与血压的关系方程。可选的，所述脉搏波曲线确定模块305包括：滤波与插值单元，用于分别对所述第一感兴趣区域的三基色的视频迹曲线和所述第二感兴趣区域的三基色的视频迹曲线进行带通滤波和三次样条插值，得到第一感兴趣区域的三基色插值后曲线和第二感兴趣区域的三基色插值后曲线；盲源信号提取单元，用于以所述心率的估计值为先验知识，分别对所述第一感兴趣区域的三基色插值后曲线和所述第二感兴趣区域的三基色插值后曲线进行盲源信号提取，得到第一感兴趣区域的光电容积脉搏波曲线和第二感兴趣区域的光电容积脉搏波曲线。可选的，所述图像帧筛选模块302包括：感兴趣区域确定单元，用于根据所述面试视频的视频图像分析面部的血管分布，确定血管分布量超过预设值且中心距离超过预设距离的两个互不重叠区域，得到第一感兴趣区域和第二感兴趣区域；基准帧确定单元，用于选取包含所述第一感兴趣区域和所述第二感兴趣区域的一帧图像作为基准帧；多帧图像提取单元，用于以所述基准帧为起始帧，选取所述基准帧之后的连续多帧图像；模板匹配单元，用于以所述基准帧为模板，利用模板匹配算法筛选所述连续多帧图像中的第一感兴趣区域和第二感兴趣区域，得到第一感兴趣区域的多个图像帧和第二感兴趣区域的多个图像帧。根据本发明提供的具体实施例，本发明公开了以下技术效果：本发明所公开的非接触式血压测量方法及系统，通过拍摄用户的面部视频来实现血压的测量，从而避免设备对人体的直接接触和压迫，提高血压测量的便捷性。同时，避免设备在人体的安装和拆卸过程，用户只需正对镜头即可，提高了测量效率。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。当前第1页12