CN113297978B - 活体检测方法、装置及电子设备 - Google Patents

Abstract

本申请适用于多光谱检测技术领域，特别涉及一种活体检测方法、装置及电子设备，该活体检测方法包括：获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素；根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一反射率值，在第二特征波段的第二反射率值；将所述第一反射率值，所述第二反射率值，以及所述第一反射率值与第二反射率值的比值输入活体检测模型获得活体检测结果。本申请实施例可以提高活体检测的精度。

Description

活体检测方法、装置及电子设备

技术领域

本申请涉及多光谱检测技术领域，特别涉及一种活体检测方法、装置及电子设备。

背景技术

活体检测是在一些身份验证场景确定对象真实生理特征的方法。

以人脸识别应用场景为例进行说明。例如，活体检测能通过眨眼、张嘴、摇头、点头等组合动作，使用人脸关键点定位和人脸追踪等技术，验证用户是否为真实活体本人操作。可有效抵御照片、换脸、面具、遮挡以及屏幕翻拍等常见的攻击手段，从而帮助用户甄别欺诈行为，保障用户的利益。

又如，由于多光谱图像中含有更丰富的场景信息，可以利用物体表面不同的反射率进行活体检测，降低系统误检率，防御非人脸材质的假体。

但是随着身份验证场景贯穿人们生活的诸多方面，又与人们的利益息息相关，因此，如何提高活体检测精度是亟需解决的技术问题。

发明内容

有鉴于此，本申请实施例提供了一种活体检测方法、装置及电子设备，能够获得精度较高的活体检测结果。

第一方面，本申请一实施例提供了一种活体检测方法，包括：

获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素；

根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一反射率值，在第二特征波段的第二反射率值；

将所述第一反射率值，所述第二反射率值，以及所述第一反射率值与所述第二反射率值的比值输入活体检测模型获得活体检测结果。

本实施例将两个波段各自的反射率值以及两个波段之间反射率比值组合成三维特征，输入到活体检测模型，通过三维特征获得更高精度的检测结果，能够满足产品的高精度需求。

作为第一方面的一种实现方式，所述根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一反射率值，在第二特征波段的第二反射率值，包括：

根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一多光谱响应值Dw1以及在第二特征波段的第二多光谱响应值Dw2；

根据所述多光谱图像获取第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2；

计算所述至少一个像素在第一特征波段的第一反射率值Dw1/Sw1，计算所述至少一个像素在第二特征波段的第二反射率值Dw2/Sw2。

作为第一方面的一种实现方式，所述根据所述多光谱图像获取第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2，包括：

确定所述多光谱图像中每个像素的多光谱响应值；

根据所述多光谱图像重构RGB图像；

将所述RGB图像转换成灰度图像；

确定所述灰度图像中灰度值小于阈值或灰度值小于或等于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2。

本实现方式将灰度值小于阈值，或灰度值小于或等于阈值的目标区域找出来，基于目标区域计算第一特征波段和第二特征波段的光源光谱响应值，可以提高获取到的光源光谱的准确度。

作为第一方面的一种实现方式，所述灰度图像中各像素对应的灰度值根据该像素在所述RGB图像中的三通道数值计算得到。

作为第一方面的一种实现方式，根据公式deta＝abs(1-G/B)+abs(1-R/B)计算各像素对应的灰度值，其中，R、G和B表示RGB图像中各像素的三通道数值，即R值、G值和B值，abs表示绝对值函数。

作为第一方面的一种实现方式，所述将所述RGB图像转换成灰度图像之后，还包括：

根据所述灰度图像确定阈值。

作为第一方面的一种实现方式，根据所述灰度图像确定阈值包括：对所述灰度图像进行直方图统计，根据直方图统计结果中最小数值区间的区间参数确定阈值。

作为第一方面的一种实现方式，所述根据直方图统计结果中最小数值区间的区间参数确定阈值，包括：

根据直方图统计结果中最小数值区间的区间边界数值和像素占比确定阈值。

作为第一方面的一种实现方式，所述根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2，包括：

计算所述目标区域在所述多光谱图像中对应的各像素的第一特征波段的多光谱响应值的平均值，得到第一光源光谱响应值Sw1；计算所述目标区域在所述多光谱图像中对应的各像素的第二特征波段的多光谱响应值的平均值，得到第二光源光谱响应值Sw2。

作为第一方面的另一种实现方式，所述根据所述多光谱图像获取第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2，包括：

确定所述多光谱图像中每个像素的多光谱响应值；

获取RGB图像，将RGB图像与多光谱图像进行匹配，获得匹配后的RGB图像；

将匹配后的RGB图像转换成灰度图像；

确定所述灰度图像中灰度值小于阈值或灰度值小于或等于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2。

本实现方式将灰度值小于阈值，或灰度值小于或等于阈值的目标区域找出来，基于目标区域计算第一特征波段和第二特征波段的光源光谱响应值，可以提高获取到的光源光谱的准确度。

作为第一方面的一种实现方式，所述灰度图像中各像素对应的灰度值根据该像素在匹配后的RGB图像中的三通道数值计算得到。

作为第一方面的一种实现方式，根据公式deta＝abs(1-G/B)+abs(1-R/B)计算各像素对应的灰度值，其中，R、G和B表示匹配后的RGB图像中各像素的三通道数值，即R值、G值和B值，abs表示绝对值函数。

作为第一方面的一种实现方式，所述将匹配后的RGB图像转换成灰度图像之后，还包括：

根据所述灰度图像确定阈值。

作为第一方面的一种实现方式，根据所述灰度图像确定阈值包括：对所述灰度图像进行直方图统计，根据直方图统计结果中最小数值区间的区间参数确定阈值。

作为第一方面的一种实现方式，所述根据直方图统计结果中最小数值区间的区间参数确定阈值，包括：

根据直方图统计结果中最小数值区间的区间边界数值和像素占比确定阈值。

作为第一方面的一种实现方式，所述根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2，包括：

计算所述目标区域在所述多光谱图像中对应的各像素的第一特征波段的多光谱响应值的平均值，得到第一光源光谱响应值Sw1；计算所述目标区域在所述多光谱图像中对应的各像素的第二特征波段的多光谱响应值的平均值，得到第二光源光谱响应值Sw2。

作为第一方面的一种实现方式，所述第一特征波段包括真实人体皮肤的吸收峰波段；所述第二特征波段包括真实人体皮肤的非吸收峰波段。

第二方面，本申请实施例提供一种活体检测装置，包括：

获取模块，用于获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素；

确定模块，用于根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一反射率值，在第二特征波段的第二反射率值；

检测模块，用于将所述第一反射率值，所述第二反射率值，以及所述第一反射率值与所述第二反射率值的比值输入活体检测模型获得活体检测结果。

第三方面，本申请实施例提供了一种电子设备，包括：存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如第一方面或第一方面任一实现方式所述的活体检测方法。

第四方面，本申请实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如第一方面或第一方面任一实现方式所述的活体检测方法。

第五方面，本申请实施例提供了一种计算机程序产品，当计算机程序产品在电子设备上运行时，使得电子设备执行如第一方面或第一方面任一实现方式所述的活体检测方法。

可以理解的是，上述第二方面至第五方面的有益效果可以参见上述第一方面中的相关描述，在此不再赘述。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1是本申请一实施例提供的一种多光谱反射率图像获取方法的实现流程示意图；

图2是本申请一实施例提供的一种对灰度图像进行直方图统计的统计结果示意图；

图3是本申请一实施例提供的另一种多光谱反射率图像获取方法的实现流程示意图；

图4是本申请一实施例提供的一种活体检测方法的实现流程示意图；

图5是本申请一实施例提供的另一种活体检测方法的实现流程示意图；

图6是本申请一实施例提供的一种活体检测装置的示意图；

图7是本申请一实施例提供的另一种活体检测装置的示意图；

图8是本申请一实施例提供的一种活体检测装置中确定模块的示意图；

图9是本申请一实施例提供的一种活体检测装置中光源光谱响应值确定子模块的示意图；

图10是本申请另一实施例提供的一种活体检测装置中光源光谱响应值确定子模块的示意图；

图11是本申请另一实施例提供的一种活体检测装置中光源光谱响应值确定子模块的示意图；

图12是本申请一实施例提供的一种电子设备的示意图。

具体实施方式

以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本发明实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本发明的描述。

在本申请说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合，并且包括这些组合。

在本申请说明书中描述的“一个实施例”或“一些实施例”等意味着在本申请的一个或多个实施例中包括结合该实施例描述的特定特征、结构或特点。由此，在本说明书中的不同之处出现的语句“在一个实施例中”、“在一些实施例中”、“在其他一些实施例中”、“在另外一些实施例中”等不是必然都参考相同的实施例，而是意味着“一个或多个但不是所有的实施例”，除非是以其他方式另外特别强调。术语“包括”、“包含”、“具有”及它们的变形都意味着“包括但不限于”，除非是以其他方式另外特别强调。

此外，在本申请的描述中，“多个”的含义是两个或两个以上。术语“第一”、“第二”、“第三”和“第四”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

为了说明本发明所述的技术方案，下面通过具体实施例来进行说明。

光源估计方法一般有如下两种：第一种，基于白世界的光源光谱估计方法。该方法找到多光谱图像最亮的区域，求取其平均光谱作为光源光谱。该方法对最亮区域是白色区域时，还原效果较好。第二种，基于灰世界的光源光谱估计方法。该方法求取整个多光谱图像的平均光谱作为光源光谱。该方法对颜色丰富的场景，还原效果较好。

这两种方法都是基于对整个多光谱图像的模糊预估，来估计光源光谱。例如，基于白世界的光源光谱预估，取多光谱图像中最亮的区域作为光源光谱，如果最亮的区域不是白色，则预估的误差较大。又如，基于灰世界的光源光谱预估，取多光谱图像中所有像素的平均值作为光源光谱，如果图像中白色区域很少且有大面积单一色时，则预估的误差较大。

以上两种方法在使用不同光源的应用场景下，适应性不强，误差较大。为了解决如何更准确预估环境光或光源光谱(或称为环境光或光源近似光谱)的技术问题，本申请实施例提供了一种获取多光谱反射率图像的方法，可以通过获取多光谱图像，在该多光谱图像中定位光源区域，从而根据该光源区域的多光谱信息确定光源光谱。

实施例一

图1是本申请一实施例提供的一种多光谱反射率图像获取方法的实现流程示意图，本实施例中的多光谱反射率图像获取方法可由电子设备执行。电子设备包括但不限于计算机、平板电脑、服务器、手机或多光谱相机等。服务器包括但不限于独立服务器或云服务器等。本实施例中的多光谱反射率图像获取方法适用于需要预估当前环境中光源光谱(或光源近似光谱)的情形。如图1所示，多光谱反射率图像获取方法可以包括步骤S110至步骤S150。

S110，获取多光谱图像，确定所述多光谱图像中每个像素的多光谱响应值。

其中，多光谱图像为单张多光谱图像。通过多光谱相机采集任一场景(场景中存在环境光或光源)的多光谱图像。单张多光谱图像包含的信息包括各个像素的响应值信息，响应值信息代表反射到多光谱相机的光在多光谱相机上的响应。该响应值信息是随着光源的强度、光源光谱的形状、光源的打光方向而变化的。

多光谱相机的通道数量可以为几个至十几个，例如八通道、九通道或十六通道等。本实施例对多光谱相机的通道数量以及各通道的波段不作具体限定。为了更好的理解本实施例，后续以九通道的多光谱相机作为多光谱相机的示例，应理解，示例性描述不能解释为对本实施例的具体限制。

作为一非限制性示例，多光谱相机为九通道的多光谱相机，九通道的多光谱相机的每个像素能获得x1，x2，x3，x4，x5，x6，x7，x8，x9这九个响应值。也就是说，每个像素的多光谱响应值为对应九个通道的九个响应值。其中，x1代表具有q1响应曲线特性的第一通道的响应值；x2代表具有q2响应曲线特性的第二通道的响应值；x3代表具有q3响应曲线特性的第三通道的响应值；......；x9代表具有q9响应曲线特性的第三通道的响应值。也就是说，xi代表具有qi响应曲线特性的第i通道的响应值，i取值为1至9的整数。

S120，根据所述多光谱图像重构红绿蓝(red green blue，RGB)图像。

RGB图像中每个像素有三个通道的响应值，即R通道的R值、G通道的G值和B通道的B值。根据多光谱图像重构RGB图像，就是根据多光谱图像中每个像素的多光谱响应值计算该像素的R值、G值和B值。

作为一实现方式，步骤S120，根据所述多光谱图像重构RGB图像，包括如下步骤S121至S124。

S121，获取多光谱相机的九个通道的量子效率(quantum efficiency，QE)响应曲线。

具体地，获取多光谱相机的九个通道的QE响应曲线矩阵，QE响应曲线矩阵可以记为q1，q2，q3，q4，q5，q6，q7，q8，q9。其中，矩阵q1为第一通道的响应曲线，矩阵q2为第二通道的响应曲线，......，矩阵q9为第九通道的响应曲线。也就是说，矩阵qj为第j通道的响应曲线，j取值为1至9的整数。需要说明的是，针对一款固定的多光谱相机(或多光谱硬件)，这些响应曲线可以通过测试得到。测试得到这些曲线后，可以预先存储在电子设备的存储器中，在需要时调用即可。

S122，获取三刺激值曲线，即r曲线，g曲线和b曲线。

获取真实三原色表色系统(CIE 1931RGB系统)的光谱三刺激值曲线，包括r曲线，g曲线和b曲线。需要说明的是，这些曲线为已知，可以从CIE标准查到。三个曲线预先存储在电子设备的存储器中，在需要时调用即可。

S123，利用九通道的QE响应曲线对三刺激值曲线进行线性拟合，得到拟合参数。

具体地，用线性拟合的方法，分别将r曲线，g曲线和b曲线用九通道的响应曲线，即q1，q2，q3，q4，q5，q6，q7，q8，q9曲线进行线性拟合。线性拟合的公式如下：

r＝a1*q1+a2*q2+a3*q3+a4*q4+a5*q5+a6*q6+a7*q7+a8*q8+a9*q9；

g＝b1*q1+b2*q2+b3*q3+b4*q4+b5*q5+b6*q6+b7*q7+b8*q8+b9*q9；

b＝c1*q1+c2*q2+c3*q3+c4*q4+c5*q5+c6*q6+c7*q7+c8*q8+c9*q9。

用偏最小二乘法解以上方程的解，求得拟合参数的值，即以下各参数的值：

a1，a2，a3，a4，a5，a6，a7，a8，a9；

b1，b2，b3，b4，b5，b6，b7，b8，b9；

c1，c2，c3，c4，c5，c6，c7，c8，c9。

S124，根据所述拟合参数和每个像素的所述多光谱响应值进行拟合计算，得到每个像素的R值、G值和B值。

具体地，根据步骤S110确定多光谱图像中某一像素的九通道响应值为：x1，x2，x3，x4，x5，x6，x7，x8，x9，根据步骤S123计算得到拟合参数，在步骤S124中根据拟合参数和该像素的九通道响应值进行拟合计算，得到该像素的R值、G值和B值。公式如下：

R＝a1*x1+a2*x2+a3*x3+a4*x4+a5*x5+a6*x6+a7*x7+a8*x8+a9*x9；

G＝b1*x1+b2*x2+b3*x3+b4*x4+b5*x5+b6*x6+b7*x7+b8*x8+b9*x9；

B＝c1*x1+c2*x2+c3*x3+c4*x4+c5*x5+c6*x6+c7*x7+c8*x8+c9*x9。

经过拟合计算得到多光谱图像中每个像素的R值、G值和B值，就得出了整幅多光谱图像对应的RGB图像，即根据多光谱图像重构了RGB图像。

在其他一些实施例中，在重构RGB图像之后，还可以对RGB图像进行白平衡，得到经白平衡后的RGB图像，可以记为RGB_wb图像。在这些实施例中，在后续的步骤S130中，是将RGB_wb图像转换成灰度图像。

在一些实现方式中，可以直接借用现有的白平衡方法，例如灰世界、白世界或自动阈值等方法，对RGB图像做白平衡处理，得到经白平衡的RGB图像RGB_wb。借用白平衡这个步骤，可以使得在后续步骤S140中得到的deta值接近0的区域，与灰色或白色区域更好地对应上，可以更准确得到区域选择结果，从而得到更准确的光源光谱。

S130，将所述RGB图像转换成灰度图像。

其中，灰度图像可称为deta图像。灰度图像中各像素对应的灰度值根据该像素在所述RGB图像中的多通道数值计算得到。

根据RGB图像中每个像素的R、G和B三个通道的R值、G值和B值，计算该像素的灰度值(或deta值)，根据各个像素的灰度值(或deta值)就得到与RGB图像对应的灰度图像(或deta图像)。也就是说，灰度图像(或deta图像)各像素对应的灰度值(或deta值)根据该像素在所述RGB图像中的多通道数值，即R值、G值和B值计算得到。

作为一非限制性示例，提取RGB图像的R、G和B三个通道，针对每一个像素，根据公式deta＝abs(1-G/B)+abs(1-R/B)求取该像素对应的deta值，并将该deta值作为灰度值赋值给灰度图像的该像素，根据各个像素的灰度值得到deta图像，其中，公式中的abs代表绝对值函数。

S140，确定所述灰度图像中灰度值小于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值计算光源光谱响应值。

具体地，确定所述灰度图像(或deta图像)中灰度值(或deta值)小于阈值的目标区域。其中，阈值t可以取接近0的某个数值。根据多光谱图像中对应该目标区域的各像素的多光谱响应值计算光源的光谱响应值。

在本实施例中，找到deta图像中deta值接近0的区域，是为了找到R值、G值、B值这三个值均接近的区域。当R值、G值、B值这三个值接近时，可能为白色区域，也可能为不同灰度的灰色区域。由于白色区域和/或灰色区域的反射率是一条直线，而入射的光源光谱曲线和反射的光源光谱曲线一致，仅仅是亮度上的差异。因而，白色区域和/或灰色区域的光谱能更加准确的反映光源光谱。

作为本实施例一实现方式，对deta图像进行直方图统计，即用直方图统计deta图像中deta数据的分布，根据deta图像的直方图统计结果确定阈值t。具体地，对deta图像进行直方图统计后，根据直方图统计结果中最小数值区间的区间参数确定阈值t。区间参数包括但不限于像素个数，像素占比，区间边界数值等中的一个或多个。

作为一非限制性示例，对灰度图像(或deta图像)进行直方图统计的统计过程如下：首先找到灰度值(或deta值)的最小值M0和最大值M10；然后将最小值M0和最大值M10之间划为10个范围(或称为数值区间)，10个数值区间从小到大依次为：[M0，M1)，[M1，M2)，[M2，M3)，[M3，M4)，[M4，M5)，[M5，M6)，[M6，M7)，[M7，M8)，[M8，M9)，[M9，M10]，其中，M0，M1，M2，M3，M4，M5，M6，M7，M8，M9，M10可称为区间数值M。统计灰度值大于等于M0且小于M1的像素个数，即第一数值区间或最小数值区间的像素个数，该像素个数占总像素个数的比例为h1，即第一数值区间的像素占比h为h1。以同样方法得到第二数值区间至第十数值区间的像素占比h依次为：h2，h3，h4，h5，h6，h7，h8，h9，h10。对deta图像进行直方图统计的统计结果的示意图如图2所示。针对第一数值区间或最小数值区间，t＝M0+(M1-M0)*h1。每个数值区间对应的t值不一样，与该数值区间的区间数值M以及h值相关。在本实施例中，只需要找到第一数值区间的t值，即确定使deta接近0的t值。

作为另一非限制性示例，首先找到灰度值(或deta值)的最小值M0和最大值M10；然后将最小值M0和最大值M10之间划为10个范围(或称为数值区间)。确定第一个数值区间，即最小数值区间的区间参数，具体地，统计deta值大于等于M0且小于M1的像素个数，即最小数值区间的像素个数，该像素个数占总像素个数的比例为h1，即第一数值区间的像素占比为h1。最后，根据M0，M1和h1确定预设数值t。例如，t＝M0+(M1-M0)*h1。这样就确定了使deta接近0的t值。

在阈值t确定以后，统计deta<t的目标区域，即找到灰度图像中deta值接近0的目标区域，计算该目标区域在多光谱图像中对应的各像素的九通道中每个通道的平均值。该目标区域的平均多光谱数据，即为近似的光源光谱。例如，灰度图像中deta<t的目标区域包括N个像素，N为正整数，获取目标区域在多光谱图像中对应的N个像素各自的九通道的多光谱响应值，针对九通道中每个通道，计算N个像素的多光谱响应值的平均值，平均值作为光源光谱响应值。N个像素中每个像素对应九通道的多光谱响应值，因而平均值是对应九个通道的九个数值。

在本实施例其他实现方式中，阈值t确定以后，统计deta<＝t的目标区域。

需要说明的是，本实施例直方图统计时数值区间的划分数量可以取经验值，例如可以根据现有拍摄数据的经验获得。区间划分得越多，得到的目标区域的deta值更接近0，得到的光源光谱在理论上会越准确，但是当区间划分得足够多导致deta值接近0的目标区域只包括几个像素时，可能会由于噪声过大，反而使得到的光源光谱噪声过大，因此区间的划分数量需要折中考虑，既不能太大，也不能太小。本申请对此不予具体限制。

还需要说明的是，直方图统计时划分好的多个数值区间，可以包括左开右闭区间，左闭右开区间，左开右开区间，左闭右闭区间等中的一种或多种的组合。本申请对此不予具体限制。

S150，根据所述多光谱图像中每个像素的多光谱响应值和所述光源光谱响应值，获取多光谱反射率图像。

作为本实施例一实现方式，根据步骤S110确定了多光谱图像中每个像素的多光谱响应值，根据步骤S140确定了光源光谱响应值，因此，在步骤S150中，将多光谱图像中每个像素的多光谱响应值除以光源光谱响应值获得多光谱反射率图像。

作为一非限制性示例，多光谱图像中某一像素的九通道多光谱响应值为x1，x2，x3，x4，x5，x6，x7，x8，x9。光源光谱响应值，即九个通道的多光谱响应值的平均值为y1，y2，y3，y4，y5，y6，y7，y8，y9。计算x1/y1，x2/y2，x3/y3，x4/y4，x5/y5，x6/y6，x7/y7，x8/y8，x9/y9得到该像素的反射率，求出每个像素的反射率后，这样就得到了多光谱图像对应的多光谱反射率图。

本实施例借用多光谱图像可以还原RGB图像的优势，从还原的RGB图像中寻找到白色或灰色区域，由于多光谱图像中白色或灰色区域的光谱为最接近光源的光谱，本方案中增加了区域选择这一步骤，将该区域的平均光谱作为光源的近似光谱，预估的光源光谱精度较高，可适用于采用不同光源的场景，基于该光源光谱计算的多光谱反射率图像也更精准。

实施例二

图3是本申请另一实施例提供的一种多光谱反射率图像获取方法的实现流程示意图，本实施例中的多光谱反射率图像获取方法可由电子设备执行。如图3所示，多光谱反射率图像获取方法可以包括步骤S210至步骤S250。应理解，实施例二与实施例一的相同之处，请详见前述方法的描述，此处不再赘述。

S210，获取多光谱图像，确定所述多光谱图像中每个像素的多光谱响应值。

S220，获取RGB图像，将RGB图像与多光谱图像进行匹配，获得匹配后的RGB图像。

在实施例一中RGB图像是根据多光谱图像重构获得的，因此RGB图像和多光谱图像具有相同视角。而在实施例二中，通过另一个相机，即彩色相机获取同一场景的RGB图像，因此RGB图像与通过多光谱相机获取的多光谱图像视角不同，需要进行匹配操作。

作为本实施例一实现方式，将RGB图像中的像素点和多光谱图像的像素点一一对应，例如将RGB图像中某物体和多光谱图像中该物体的像素点对应。当通过RGB图像找到灰色白色区域时，再通过对应关系找到多光谱图像中的灰色白色区域，计算该区域的多通道响应的平均值，作为近似的光源光谱响应值。

在本实施例中，彩色相机和多光谱相机相邻设置，两者位置越近，两者的接收端或成像端拍摄的视场越接近，这样在匹配过程中RGB图像和多光谱图像具备更多的对应像素点，从而可以增加光源光谱预估结果的精度。

S230，将匹配后的RGB图像转换成灰度图像。

其中，灰度图像中各像素对应的灰度值根据该像素在匹配后的RGB图像中的多通道数值计算得到。

S240，确定所述灰度图像中灰度值小于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值计算光源光谱响应值。

S250，根据所述多光谱图像中每个像素的多光谱响应值和所述光源光谱响应值，获取多光谱反射率图像。

实施例二与实施例一的不同之处在于步骤S120与S220不同，其他步骤相同或相似。在实施例一中RGB图像是根据多光谱图像重构获得的，因而RGB图像和多光谱图像是由同一相机获得的，两者具有相同视角，因而实施例一预估的光源光谱精度较实施例二更高。

基于实施例一和实施例二的方法获得的多光谱反射率图像，进行活体检测或应用在其他的模型，在不同光源下应用时均可以有更好的鲁棒性。例如，基于多光谱反射率图像做活体检测，其分析结果不随光源的变化而变化，具有较好的鲁棒性。接下来介绍一种活体检测方法。

由于真实人体皮肤和假体(例如假手指或假面具等)在几个特征波段的光谱特性有较大的区别，用本申请提供的波段比方法能应用皮肤本身的光谱特性，来排除掉大部分的假体，足够满足普通产品的精度需要。例如真实人体皮肤的特性包括：在420至440nm(单位：纳米)波段，皮肤特有黑色素吸收；在550至590nm波段，皮肤特有血红蛋白吸收；在960至980nm波段，皮肤特有水分吸收；在800至850nm波段，皮肤吸收较弱(即反射较高)等。对于精度要求较高的支付类消费场景，该波段比方法可以作为多光谱活检的第一步判断，将大部分假体排除掉，当遇到精度特别高的假体，再用机器学习或深度学习等精度更高的模型来判断。该波段比方法计算过程简单，且受环境光和暗噪声等因素的影响较小。

实施例三

图4是本申请另一实施例提供的一种活体检测方法的实现流程示意图，本实施例中的活体检测方法可由电子设备执行。如图4所示，活体检测方法可以包括步骤S310至步骤S340。

S310，获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素。

其中，人体皮肤包括但不限于人体中未被覆盖的某个部位或某个区域的皮肤，例如人脸皮肤、或人脸的某个区域的皮肤、或手指的皮肤等。

通过多光谱相机获取包含人体皮肤的多光谱图像。所述多光谱图像至少包含一个像素。需要说明的是，该至少一个像素是针对人体皮肤成像的像素。

S320，确定所述至少一个像素分别在第一特征波段和第二特征波段的第一多光谱响应值Dw1及第二多光谱响应值Dw2。

其中，根据多光谱图像，确定至少一个像素在第一特征波段的第一多光谱响应值Dw1，以及在第二特征波段的第二多光谱响应值Dw2。

根据实施一的描述可知，多光谱图像包括每个像素的多通道的多光谱响应值。在实施例一中，并不限定通道的数量和波段，而在实施例三中，多通道至少包括第一特征波段和第二特征波段这两个通道，并不限定其他通道的数量和波段。也就是说，在实施例三中，多光谱相机的通道数量至少为2个，至少包括第一特征波段和第二特征波段这两个通道。多光谱图像包括每个像素的至少两个通道的多光谱响应值，即包括第一特征波段的第一多光谱响应值Dw1，以及第二特征波段的第二多光谱响应值Dw2。因而，可以根据多光谱图像确定人体皮肤对应的至少一个像素的第一特征波段的第一多光谱响应值Dw1及第二特征波段的第二多光谱响应值Dw2。

在实施例三中，可以根据真实人体皮肤的反射光谱特性，选择两个有代表性的波段，即第一特征波段w1和第二特征波段w2。

在一些实现方式中，第一特征波段w1选择为真实人体皮肤特有的吸收峰波段，在该波段假体和真实人体皮肤的反射率有较大差异。例如420至440nm波段或在该波段内的某个波段，该波段为真实人体皮肤特有黑色素吸收波段；又如550至590nm波段或在该波段内的某个波段，该波段为真实人体皮肤特有血红蛋白吸收波段；再如960至980nm波段或在该波段内的某个波段，该波段为真实人体皮肤特有水分吸收波段。

在一些实现方式中，第二特征波段w2选择为真实人体皮肤的非吸收峰波段，即真实人体皮肤吸收较弱(或反射较高)的波段，例如800至850nm波段或在该波段内的某个波段。

S330，根据所述多光谱图像分别获取第一特征波段和第二特征波段的第一光源光谱响应值Sw1以及第二光源光谱响应值Sw2。

其中，根据多光谱图像获取第一特征波段的第一光源光谱响应值Sw1，以及第二特征波段的第二光源光谱响应值Sw2。

在实施例三的一些实现方式中，可以用现有技术获取多光谱图像在第一特征波段的第一光源光谱响应值Sw1，以及在第二特征波段的第二光源光谱响应值Sw2。

在实施例三的其他一些实现方式中，可以利用实施例一和实施例二中描述的获取光源光谱响应值的方法，获取第一特征波段的第一光源光谱响应值Sw1，以及在第二特征波段的第二光源光谱响应值Sw2。此处未详细描述之处，请参见实施例一和实施例二的相关描述之处。

具体地，首先，获取多光谱图像对应的RGB图像，可以根据多光谱图像重构RGB图像(参见实施一)，也可以在拍摄多光谱图像时针对相同场景拍摄RGB图像(参见实施例二)；然后，将RGB图像转换成灰度图像；再确定灰度图像中灰度值小于阈值的目标区域；最后，计算目标区域的多通道中第一特征波段这一通道的多光谱响应值的平均值，即第一光源光谱响应值Sw1；计算目标区域的多通道中第二特征波段这一通道的多光谱响应值的平均值，即第二光源光谱响应值Sw2。

需要说明的是，一方面，由于实施例一和实施例二的方法预估的光源光谱更准确，因而基于实施例一和实施例二的相关描述得到的第一光源光谱响应值Sw1和第二光源光谱响应值Sw2更准确，从而可以提高后续活体检测结果的精度。另一方面，由于实施例一和实施例二中预估光源光谱的方法可适用于不同光源的应用场景，因此，活体检测方案在不同光源下应用时可以有更好的鲁棒性。

S340，计算Dw1/Dw2与Sw2/Sw1的乘积，并将乘积与阈值k比较，若乘积小于阈值k，则判断所述人体为活体。

至少一个像素点在第一特征波段w1的多光谱响应值为Dw1，光源光谱的第一特征波段w1的预估响应值为第一光源光谱响应值Sw1；至少一个像素点在第二特征波段w2的多光谱响应值为Dw2，光源光谱的第二特征波段w2的预估响应值为第二光源光谱响应值Sw2。

计算(Dw1/Dw2)*(Sw2/Sw1)得到乘积Rw，并将乘积Rw与阈值k进行比较，根据比较结果获得活体检测结果。

作为一实现方式，首先，计算Dw1与Sw1的比值，即计算至少一个像素在第一特征波段的反射率值，可以记为Rw1，Rw1＝Dw1/Sw1；计算Dw2与Sw2的比值，即计算至少一个像素在第二特征波段的反射率值，可以记为Rw2，Rw2＝Dw2/Sw2。然后，计算Rw1与Rw2的比值，可以记为Rw，Rw＝Rw1/Rw2＝(Dw1/Dw2)*(Sw2/Sw1)。因此，本实现方式可以称作波段比的活体检测方法。

在一些实施例中，若乘积Rw小于阈值k，则判定人体为活体；若乘积Rw等于或大于阈值k，则判断人体为假体。在其他一些实施例中，根据活体检测的实际精度要求，调整比较条件，例如，当乘积Rw等于阈值k时，对应的活体检测结果可以设置为：判定人体为活体。本申请对此不予具体限制。

在图4所示实施例的基础上，在其他实施例中，步骤S340之前还包括确定阈值k的步骤。

作为一实现方式，确定阈值k的过程包括：获取多个真实皮肤样本在第一特征波段和第二特征波段的第一样本反射率R1和第二样本反射率R2，计算多个真实皮肤样本的第一样本反射率比值，即R1/R2，确定多个真实皮肤样本中第一样本反射率比值的最大值a。此外，获取多个不同种类的假体样本在第一特征波段和第二特征波段的第三样本反射率R3和第四样本反射率R4，计算多个假体样本的第二样本反射率比值R3/R4，确定多个假体样本中第二样本反射率比值的最小值b。最后，根据最大值a和最小值b确定阈值k。

作为一非限制性示例，通过光谱仪采集M个(M为大于1的整数)不同的真实皮肤样本各自在第一特征波段的第一样本反射率R1，以及在第二特征波段的第二样本反射率R2，计算M个真实皮肤样本中每个真实皮肤样本的第一样本反射率比值，即R1/R2，找到M个真实皮肤样本的第一样本反射率比值的最大值a。此外，采用与真实皮肤样本相同的处理方法，通过光谱仪采集N个(N为大于1的整数)不同种类的假体样本各自在第一特征波段的第三样本反射率R3，以及在第二特征波段的第四样本反射率R4，计算N个假体样本中每个假体样本的第二样本反射率比值，即R3/R4，找到N个假体样本的第二样本反射率比值的最小值b。然后根据a和b确定阈值k的取值范围。例如阈值k的取值范围为：(a+b)/2>＝k>＝min(a,b)，其中min表示取最小值函数。即阈值k大于或等于a和b两个值中较小的值，阈值k小于或等于a和b两个值的均值。阈值k的具体取值可以跟实际应用的需求确定，本实施例通过简单的阈值k设计即可将较多的活体和假体进行区分。

实施例四

图5是本申请另一实施例提供的一种活体检测方法的实现流程示意图，本实施例中的活体检测方法可由电子设备执行。如图5所示，活体检测方法可以包括步骤S410至步骤S440。应理解，实施例四与实施例三中相同之处，请详见前述实施例三的描述，此处不再赘述。

S410，获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素。

S420，确定所述至少一个像素分别在第一特征波段和第二特征波段的第一多光谱响应值Dw1及第二多光谱响应值Dw2。

S430，根据所述多光谱图像分别获取第一特征波段和第二特征波段的第一光源光谱响应值Sw1以及第二光源光谱响应值Sw2。

S440，计算Dw1与Sw1的第一比值，计算Dw2与Sw2的第二比值，计算第一比值与第二比值的第三比值。

其中，计算Dw1与Sw1的第一比值，即计算至少一个像素在第一特征波段的反射率值，第一比值可以记为Rw1，Rw1＝Dw1/Sw1；计算Dw2与Sw2的第二比值，即计算至少一个像素在第二特征波段的反射率值，第二比值可以记为Rw2，Rw2＝Dw2/Sw2。然后，计算Rw1与Rw2的第三比值，第三比值可以记为Rw，Rw＝Rw1/Rw2＝(Dw1/Dw2)*(Sw2/Sw1)。

S450，将第一比值、第二比值和第三比值输入活体检测模型获得活体检测结果。

其中，活体检测模型为经训练的用于判断待测人体是否为活体的检测模型。将第一比值Rw1、第二比值Rw2和第三比值Rw输入活体检测模型，模型可以输出待测人体为活体或假体的分类结果。

在本实施例中，活体检测模型可以包括机器学习或深度学习模型。例如，支持向量机模型、神经网络模型、贝叶斯分类器或随机森林等模型。本申请对活体检测模型不予具体限制。

在一些实现方式中，活体检测模型可以包括二分类模型，二分类模型的二分类结果包括待测人体为活体和待测人体为假体。例如，将[Rw1，Rw2，Rw1/Rw2]输入活体检测模型，模型的输出为1，表示待测人体为活体；输出为0，表示待测人体为假体。

在其他一些实现方式中，活体检测模型可以包括多分类活体检测模型，在这些实现方式中，活体检测模型可以对活体和/或假体进行更细的分类。例如，对假体作进一步细分以区分假体的不同种类或类别(如不同种类或类别的假体对应不同材质的假体)。本申请对活体检测模型的分类数量不作具体限制。

需要说明的是，在利用活体检测模型之前，还需要获取经训练的活体检测模型。作为一非限制性示例，获取经训练的活体检测模型的过程包括：获取多个真实皮肤样本各自的第一样本向量及对应标签，第一样本向量包括真实皮肤样本在第一特征波段的第一样本反射率值、在第二特征波段的第二样本反射率值、第一样本反射率值与第二样本反射率值的比值这三个特征；获取多个不同种类的假体样本的第二样本向量及对应标签，第二样本向量包括假体样本在第一特征波段的第三样本反射率值、在第二特征波段的第四样本反射率值、第三样本反射率值与第四样本反射率值的比值这三个特征；利用第一样本向量及对应标签和第二样本向量及对应标签作为训练样本，对活体检测模型进行训练，得到经训练的活体检测模型。这样，经训练的活体检测模型可以实现活体和假体的分类，也就是说，经训练的活体检测模型可以用于识别待测人体是否为活体。应理解，作为一非限制性示例，获取第一样本反射率值、第二样本反射率值、第一样本反射率值与第二样本反射率值的比值、第三样本反射率值、第四样本反射率值、以及第三样本反射率值与第四样本反射率值的比值的过程，可以参见确定阈值k的相关描述。

在本实施例中，将波段比值Rw1/Rw2加入到两个特征波段的反射率特征中，组成了三维的特征组合向量，即[Rw1，Rw2，Rw1/Rw2]，增加了特征的维度。将该特征组合向量输入活体检测模型输出活体检测结果，该活体检测结果由[Rw1，Rw2，Rw1/Rw2]中的三个特征共同决定，可以获得更准确的结果。

应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本申请实施例的实施过程构成任何限定。

本申请一实施例还提供一种活体检测装置。该活体检测装置中未详细描述之处请详见前述实施例中方法的描述。

参见图6，图6是本申请实施例提供的一种活体检测装置的示意框图。所述活体检测装置包括：获取模块81，确定模块82和检测模块83。

获取模块81，用于获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素；

确定模块82，用于根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一反射率值，在第二特征波段的第二反射率值；

检测模块83，用于将所述第一反射率值，所述第二反射率值，以及所述第一反射率值与所述第二反射率值的比值输入活体检测模型获得活体检测结果。

可选的，作为一种实现方式，如图7所示，所述确定模块82包括：多光谱响应值确定子模块821，光源光谱响应值确定子模块822和反射率确定子模块823。

其中，多光谱响应值确定子模块821，用于根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一多光谱响应值Dw1以及在第二特征波段的第二多光谱响应值Dw2。

光源光谱响应值确定子模块822，用于根据所述多光谱图像获取第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2。

反射率确定子模块823，用于计算所述至少一个像素在第一特征波段的第一反射率值Dw1/Sw1，计算所述至少一个像素在第二特征波段的第二反射率值Dw2/Sw2。

可选的，作为一种实现方式，如图8所示，所述光源光谱响应值确定子模块822包括：确定孙模块8221，重构孙模块8222，转换孙模块8223和计算孙模块8224。

其中，确定孙模块8221，用于确定所述多光谱图像中每个像素的多光谱响应值。

重构孙模块8222，用于根据所述多光谱图像重构RGB图像。

转换孙模块8223，用于将所述RGB图像转换成灰度图像。

计算孙模块8224，用于确定所述灰度图像中灰度值小于阈值或灰度值小于或等于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2。

可选的，作为一种实现方式，所述灰度图像中各像素对应的灰度值根据该像素在所述RGB图像中的三通道数值计算得到。

可选的，作为一种实现方式，根据公式deta＝abs(1-G/B)+abs(1-R/B)计算各像素对应的灰度值，其中，R、G和B表示RGB图像中各像素的三通道数值，即R值、G值和B值，abs表示绝对值函数。

可选的，作为一种实现方式，在图8所示实现方式的基础上，如图9所示，所述光源光谱响应值确定子模块822，还包括：阈值确定孙模块8225。

阈值确定孙模块8225，用于根据所述灰度图像确定阈值。

可选的，作为一种实现方式，阈值确定孙模块8225具体用于：对所述灰度图像进行直方图统计，根据直方图统计结果中最小数值区间的区间参数确定阈值。

可选的，作为一种实现方式，所述根据直方图统计结果中最小数值区间的区间参数确定阈值，包括：

根据直方图统计结果中最小数值区间的区间边界数值和像素占比确定阈值。

可选的，作为一种实现方式，所述计算孙模块8224，具体用于：

计算所述目标区域在所述多光谱图像中对应的各像素的第一特征波段的多光谱响应值的平均值，得到第一光源光谱响应值Sw1；计算所述目标区域在所述多光谱图像中对应的各像素的第二特征波段的多光谱响应值的平均值，得到第二光源光谱响应值Sw2。

可选的，作为另一种实现方式，如图10所示，所述光源光谱响应值确定子模块822包括：确定孙模块8221’，匹配孙模块8222’，转换孙模块8223’和计算孙模块8224’。

其中，确定孙模块8221’，用于确定所述多光谱图像中每个像素的多光谱响应值。

匹配孙模块8222’，用于获取RGB图像，将RGB图像与多光谱图像进行匹配，获得匹配后的RGB图像。

转换孙模块8223’，用于将匹配后的RGB图像转换成灰度图像。

计算孙模块8224’，用于确定所述灰度图像中灰度值小于阈值或灰度值小于或等于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2。

可选的，作为一种实现方式，所述灰度图像中各像素对应的灰度值根据该像素在匹配后的RGB图像中的三通道数值计算得到。

可选的，作为一种实现方式，根据公式deta＝abs(1-G/B)+abs(1-R/B)计算各像素对应的灰度值，其中，R、G和B表示匹配后的RGB图像中各像素的三通道数值，即R值、G值和B值，abs表示绝对值函数。

可选的，作为一种实现方式，在图10所示实现方式的基础上，如图11所示，所述光源光谱响应值确定子模块822，还包括：阈值确定孙模块8225’。

阈值确定孙模块8225’，用于根据所述灰度图像确定阈值。

可选的，作为一种实现方式，阈值确定孙模块8225’具体用于：对所述灰度图像进行直方图统计，根据直方图统计结果中最小数值区间的区间参数确定阈值。

可选的，作为一种实现方式，所述根据直方图统计结果中最小数值区间的区间参数确定阈值，包括：

根据直方图统计结果中最小数值区间的区间边界数值和像素占比确定阈值。

可选的，作为一种实现方式，所述计算孙模块8224’，具体用于：

计算所述目标区域在所述多光谱图像中对应的各像素的第一特征波段的多光谱响应值的平均值，得到第一光源光谱响应值Sw1；计算所述目标区域在所述多光谱图像中对应的各像素的第二特征波段的多光谱响应值的平均值，得到第二光源光谱响应值Sw2。

可选的，作为一种实现方式，所述第一特征波段包括真实人体皮肤的吸收峰波段；所述第二特征波段包括真实人体皮肤的非吸收峰波段。

本申请实施例还提供了一种电子设备，如图12所示，电子设备可以包括一个或多个处理器120(图12中仅示出一个)，存储器121以及存储在存储器121中并可在一个或多个处理器120上运行的计算机程序122，例如，获取光源光谱和/或多光谱反射率图像的程序。一个或多个处理器120执行计算机程序122时可以实现光源光谱获取方法和/或多光谱反射率图像获取方法实施例中的各个步骤。或者，一个或多个处理器120执行计算机程序122时可以实现光源光谱获取装置和/或多光谱反射率图像获取装置实施例中各模块/单元的功能，此处不作限制。

本领域技术人员可以理解，图12仅仅是电子设备的示例，并不构成对电子设备的限定。电子设备可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如电子设备还可以包括输入输出设备、网络接入设备、总线等。

在一个实施例中，所称处理器120可以是中央处理单元(Central ProcessingUnit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。

在一个实施例中，存储器121可以是电子设备的内部存储单元，例如电子设备的硬盘或内存。存储器121也可以是电子设备的外部存储设备，例如电子设备上配备的插接式硬盘，智能存储卡(smart media card，SMC)，安全数字(secure digital，SD)卡，闪存卡(flash card)等。进一步地，存储器121还可以既包括电子设备的内部存储单元也包括外部存储设备。存储器121用于存储计算机程序以及电子设备所需的其他程序和数据。存储器121还可以用于暂时地存储已经输出或者将要输出的数据。

所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本申请的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

本申请实施例还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现可实现光源光谱获取方法实施例和/或多光谱反射率图像获取方法实施例中的步骤。

本申请实施例提供了一种计算机程序产品，当计算机程序产品在电子设备上运行时，使得电子设备可实现光源光谱获取方法实施例和/或多光谱反射率图像获取方法实施例中的步骤。

在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。

本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。

在本发明所提供的实施例中，应该理解到，所揭露的装置/电子设备和方法，可以通过其它的方式实现。例如，以上所描述的装置/电子设备实施例仅仅是示意性的，例如，模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。

作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。

集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，计算机程序包括计算机程序代码，计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读介质可以包括：能够携带计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(read-only memory，ROM)、随机存取存储器(random access memory，RAM)、电载波信号、电信信号以及软件分发介质等。需要说明的是，计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。

以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。

Claims (8)

1.一种活体检测方法，其特征在于，包括：

获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素；

根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一多光谱响应值Dw1以及在第二特征波段的第二多光谱响应值Dw2；

确定所述多光谱图像中每个像素的多光谱响应值；

根据所述多光谱图像重构RGB图像，或者，从获取的RGB图像中确定与所述多光谱图像相匹配的RGB图像；

将RGB图像转换成灰度图像；

确定所述灰度图像中灰度值小于或等于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2；

计算所述至少一个像素在第一特征波段的第一反射率值Dw1/Sw1，计算所述至少一个像素在第二特征波段的第二反射率值Dw2/Sw2；

将所述第一反射率值，所述第二反射率值，以及所述第一反射率值与所述第二反射率值的比值输入活体检测模型获得活体检测结果。

2.如权利要求1所述的活体检测方法，其特征在于，所述转换成灰度图像之后，还包括：

对所述灰度图像进行直方图统计，根据直方图统计结果中最小数值区间的区间参数确定阈值。

3.如权利要求2所述的活体检测方法，其特征在于，所述根据直方图统计结果中最小数值区间的区间参数确定阈值，包括：

根据直方图统计结果中最小数值区间的区间边界数值和像素占比确定阈值。

4.如权利要求2所述的活体检测方法，其特征在于，所述根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2，包括：

计算所述目标区域在所述多光谱图像中对应的各像素的第一特征波段的多光谱响应值的平均值，得到第一光源光谱响应值Sw1；计算所述目标区域在所述多光谱图像中对应的各像素的第二特征波段的多光谱响应值的平均值，得到第二光源光谱响应值Sw2。

5.如权利要求1所述的活体检测方法，其特征在于，所述第一特征波段包括真实人体皮肤的吸收峰波段；所述第二特征波段包括真实人体皮肤的非吸收峰波段。

6.一种活体检测装置，其特征在于，包括：

获取模块，用于获取包含人体皮肤的多光谱图像，所述多光谱图像包含至少一个像素；

确定模块，用于根据所述多光谱图像确定所述至少一个像素在第一特征波段的第一多光谱响应值Dw1以及在第二特征波段的第二多光谱响应值Dw2；确定所述多光谱图像中每个像素的多光谱响应值；根据所述多光谱图像重构RGB图像，或者，从获取的RGB图像中确定与所述多光谱图像相匹配的RGB图像；将RGB图像转换成灰度图像；确定所述灰度图像中灰度值小于或等于阈值的目标区域，根据所述目标区域在所述多光谱图像中对应的各像素的多光谱响应值，计算第一特征波段的第一光源光谱响应值Sw1和第二特征波段的第二光源光谱响应值Sw2；计算所述至少一个像素在第一特征波段的第一反射率值Dw1/Sw1，计算所述至少一个像素在第二特征波段的第二反射率值Dw2/Sw2；

检测模块，用于将所述第一反射率值，所述第二反射率值，以及所述第一反射率值与所述第二反射率值的比值输入活体检测模型获得活体检测结果。

7.一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现如权利要求1至5任一项所述活体检测方法。

8.一种计算机存储介质，所述计算机存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至5任一项所述活体检测方法。