《一种基于稀疏表示与多特征点的三维人脸识别方法.pdf》由会员分享,可在线阅读,更多相关《一种基于稀疏表示与多特征点的三维人脸识别方法.pdf(12页完整版)》请在专利查询网上搜索。
1、(10)申请公布号 CN 103984920 A (43)申请公布日 2014.08.13 CN 103984920 A (21)申请号 201410172494.1 (22)申请日 2014.04.25 G06K 9/00(2006.01) G06K 9/46(2006.01) (71)申请人 同济大学 地址 200092 上海市杨浦区四平路 1239 号 (72)发明人 张林 丁志轩 李力达 李宏宇 沈莹 (74)专利代理机构 上海智信专利代理有限公司 31002 代理人 吴林松 (54) 发明名称 一种基于稀疏表示与多特征点的三维人脸识 别方法 (57) 摘要 本发明公开了一种基于稀疏表。
2、示与多特征点 的三维人脸识别方法, 步骤 ( 一 ) 在尺度空间中 提取三维人脸扫描数据中的特征点并对其进行描 述 ; 步骤(二)根据步骤(一)中确定的待测人脸 的三维数据特征构建三维人脸数据库 ; 步骤(三) 运用步骤 ( 一 ) 提取得到待测人脸的特征描述算 子集合, 在步骤 ( 二 ) 中所确定的三维人脸数据 库比较以确定待测人脸身份的方法。本发明所示 的三维人脸识别方法选取固定列数的特征描述算 子构建子字典, 并进行稀疏滤波, 从而提高识别效 率与精确度, 可用于对身份识别有严格要求的场 合, 解决了三维人脸扫描数据中存在的人脸部分 缺失、 表情变化大、 面部被遮挡等问题。 (51)I。
3、nt.Cl. 权利要求书 2 页 说明书 8 页 附图 1 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书2页 说明书8页 附图1页 (10)申请公布号 CN 103984920 A CN 103984920 A 1/2 页 2 1. 一种基于稀疏表示与多特征点的三维人脸识别方法, 其特征在于 : ( 一 ) 确定待测人脸的三维数据特征 : (1) 定义人脸网状结构 M 包含顶点数据 V 与边缘数据 E, V 为各顶点在三维空间中 的具体坐标, E 为包含所有两两相连的顶点对, 对人脸网状结构 M 构建平滑尺度空间 : 在建 立 k 个尺度空间的过程中, 选取 k 。
4、个不同大小的 值, 对顶点数据 V 进行平滑操作, 得到 k 组顶点数据集, 所述 k 个不同大小的 值通过得到, 其中, 表示所 选取的第 s 个 的估计值, 0表示用于创建 k 个尺度空间的基准值, 表示边缘数据的平 均值 ; (2) 提取特征点 : (2-1) 分别计算输入的 k 组顶点数据集中各顶点的平均曲率其 中,与分别代表该第 i 个顶点在平滑尺度 s 中的最大和最小曲率 ; (2-2) 分别计算输入的人脸网状结构中所有顶点在相邻平滑尺度下的曲率差 (2-3) 选取满足顶点在相邻平滑尺度下的曲率差的值同时大于或者同时小于其所有邻 居顶点的曲率差的值的顶点为特征点 ; (3) 创建特。
5、征点描述算子, 用于描述特征点周围的信息 : (3-1) 选取每个特征点的邻居顶点, 对于每一个特征点 P, 选取以该点为球心、 r 为半径 的球形区域内的所有顶点为其邻居顶点, 并确定该点的正则方向 ; (3-2) 将所有的邻居顶点的法线投影到过 P 点的切平面上, 同时将投影后的法线向量 归纳到含 h 项的带权重的直方图中 ; (3-3) 划分若干个子区域, 确定每个子区域的形状索引信息直方图与倾斜角信息直方 图 : 将球形区域划分为 q 个子区域, 每个子区域使用两个直方图 ps和 p用于构造描述算 子, 其中, ps的值通过公式 Si 2/tan-1(ci, 1+ci, 2)/(ci,。
6、 1-ci, 2) 得到, ci, 1与 ci, 2分别表示 最大、 最小曲率, p定义为投影后的法线向量与正则方向间的倾斜角 ; (3-4)将所有子区域的两个直方图值连接在一起得到特征点的描述算子f : 将q个区域 的直方图连接成向量, 得到 f ps, 1p, 1ps, qp, qT作为特征点的描述算子 ; 其中, r,h,q 根据实际构建特征点的描述算子维度的大小来调整选取 ; (二)根据步骤(1)中确定的待测人脸的三维数据特征构建三维人脸数据库D : 对于每 一个对象, 建立其对应的特征池 Si F1, F2, , Fj 所述特征池中包含所有该对象历 次注册后经计算而得的特征点的描述算。
7、子集, 则三维人脸数据库 D 中, 将所有对象的特征 池根据其对应的身份序号依次排序而构成, 即 D S1, S2, Sc, 其中 C 为三维人脸数据 库 D 中注册成功的对象个数 ; ( 三 ) 比较待测人脸与数据库中数据信息以确定待测人脸的身份 : (I) 多特征点稀疏表示 : 将待识别人脸的特征点的描述算子集 Fprobe, 使用多特征点稀 疏表示方法求解方程式 Fprobe DX 中与 D 对应的系数矩阵 X, 由于 D 的列数太大, 对于每一 个 f 通过提取 D 中与 f 最相似的前 L 个特征描述算子构建 D 的子字典 A, 矩阵 X 通过 n 次并 权 利 要 求 书 CN 1。
8、03984920 A 2 2/2 页 3 行计算满足 f Ax 来确定 ; ( ) 稀疏滤波, 筛选待测人脸的特征点的描述算子集对应的系数向量 : 若待检测人脸 中的一个特征描述算子 f 与子字典 A 中的一个特征描述算子对应的系数 x 的各个数值分布 的非常离散, 则此特征描述算子 f 视为噪音向量, 将其排除出进一步的人脸身份识别参考 范围 ; ( ) 误差重构, 确定待识别人脸身份 : 分别计算筛选得到的待测人脸特征描述算子集 与子字典中各类的关联程度, 选取最大的关联度的对象作为待测人脸的身份, 关联度通过 确定重构误差而得到, 重构误差越小, 关联度越高。 2. 根据权利要求 1 所。
9、述的三维人脸识别方法, 其特征如下 : 所述步骤 ( 一 ) 中, 顶点数 据 V 进行平滑操作如下 : (1-1)确定所需创建的尺度空间的个数k, 并选取用于创建k个尺度空间的基准值0; (1-2) 分别确定 k 个尺度空间对应的 估计值, 第 s 个尺度空间对应的 的估计值 表示 E 的平均值 ; (1-3)初始化临时变量c0, 依次确定k组平滑后的顶点数据 : 在创建第s组顶点数 据的过程中, 使用第s-1组的顶点数据Vs-1作为第s组顶点数据Vs的初始值, 当确定第一组 顶点数据 V1时使用输入的顶点数据作为初始值, 并执行以下循环, 迭代更新 Vs: 若 则更新 c以及 Vs中所有顶。
10、点的坐标, 令 c c+1, 令表示 更新后的第 s 组顶点数据中的第 i 个顶点的坐标, Ni表示所有与第 i 个顶点相邻的顶点构 成的集合, vi表示第 s 组顶点数据中的第 i 个顶点的当前的坐标 ; 若则得到第 s 组 顶点数据 Vs, 退出更新 Vs的循环, 利用 Vs继续确定 Vs+1。 3. 根据权利要求 1 所述的三维人脸识别方法, 其特征在于 : 所述步骤 (3-1) 中, 选取直 方图中的最高峰值以及大于最高峰值的 t的峰值作为正则方向 ; 所述步骤 (3-2) 中, 直方图中各项权重为到点 P 的测地距离的高斯权重。 4. 根据权利要求 3 所述的三维人脸识别方法, 其特。
11、征在于 : 所述步骤 (3-4) 中, r 取值 为 9s, h 取值为 360, t 取值为 80, q 取值为 9。 5. 根据权利要求 1 所述的三维人脸识别方法, 其特征在于 : 所述步骤 (I) 中, 通过确定 f 与 D 中的每一列的内积而确定相似度。 6. 根据权利要求 5 所述的三维人脸识别方法, 其特征在于 : 所述步骤 (I) 中, L 根据实 际运算速度在 200 到 1000 内选取。 7. 根据权利要求 1 所述的三维人脸识别方法, 其特征在于 : 所述步骤 ( ) 中, 系 数向量 x 的稀疏度通过公式获得, 其中, 表示 x 中最大的前 u的元素之和, l 表示子。
12、字典 A 中的对象的个数。 8. 根 据 权利要求 1 所述的三维人脸识别方法, 其特征在于 : 重构误差经 由 确定, 其中, ni为筛选后的系数向量的个数, 意为选取 中只与类别 j 相关的系数。 权 利 要 求 书 CN 103984920 A 3 1/8 页 4 一种基于稀疏表示与多特征点的三维人脸识别方法 技术领域 0001 本发明涉及一种基于稀疏表示与多特征点的三维人脸识别方法, 属于模式识别领 域, 可应用于对身份识别要严格要求的场合。 背景技术 0002 人脸识别 (Face Recognition), 作为计算机生物特征 (Biometric) 识别技术中 的典范一直受到社会。
13、各界及各个科学领域的广泛关注, 其已经不再仅仅作为计算机视觉研 究性的课题。随着社会硬件技术和视频监控技术的发展, 许多监控及鉴别领域需要精准及 快速的身份确认技术 : 在国外人脸识别技术应用于国家重要部门及军事区国防区等安防部 门, 在国内这一项技术也被广泛地应用于公安, 交通, 金融以及物业管理等领域, 随着近年 来新媒体技术带动的人机交互热潮以及简易计算机摄像头在社会民众间的逐渐普及, 人们 逐渐将人脸识别看成是一种人机交互的方式, 将其运用于网络或新媒体交互产品也是今后 一个发展新媒体技术的切入点。 0003 目前已经有能够获取目标三维信息的深度摄像机, 随着各种三维摄像机和打印机 在。
14、工程界的热议和发展, 带有深度信息的三维摄像机已经成为了未来摄像相机硬件发展的 趋势, 现在生物特征界已经有一批学者在进行三维人脸识别研究。和通过普通摄像头获取 人脸数据不同, 三维人脸数据获取的是目标信息的三维形状, 而不再是目标信息的色彩图, 所以和普通的二维人脸识别相比, 三维人脸识别相对于光照以及遮挡等问题并不敏感, 在 不同的光照和遮挡下也做到较高的识别率, 故应用范围及效果更佳。 发明内容 0004 本发明的目的在于提供一种基于稀疏表示与多特征点的三维人脸识别方法, 其针 对三维人脸扫描数据中存在的人脸部分缺失、 表情变化大、 面部被遮挡等问题, 通过三维人 脸数据上局部特征描述算。
15、子的确定以及多特征点稀疏表示 (3DMKDSRC), 从而建立了一种精 确、 快速的三维人脸识别方法。 0005 为达到上述目的, 本发明的解决方案是 : 0006 一种基于稀疏表示与多特征点的三维人脸识别方法 : 0007 ( 一 ) 确定待测人脸的三维数据特征 : 0008 (1) 定义人脸网状结构 M 包含顶点数据 V 与边缘数据 E, V 为各顶点在三维空 间中的具体坐标, E 为包含所有两两相连的顶点对, 对人脸网状结构 M 构建平滑尺度空间 : 在建立 k 个尺度空间的过程中, 选取 k 个不同大小的 值, 对顶点数据 V 进行平滑操作, 得到 k 组顶点数据集, 所述 k 个不同。
16、大小的 值通过得到, 其中,表 示所选取的第 s 个 的估计值, 0表示用于创建 k 个尺度空间的基准值, 表示边缘数据 的平均值 ; 0009 (2) 提取特征点 : 说 明 书 CN 103984920 A 4 2/8 页 5 0010 (2-1) 分 别 确 定 输 入 的 k 组 顶 点 数 据 集 中 各 顶 点 的 平 均 曲 率 其中,与分别代表该第 i 个顶点在平滑尺度 s 中的最大和 最小曲率 ; 0011 (2-2) 分别确定输入的人脸网状结构中所有顶点在相邻平滑尺度下的曲率差 0012 (2-3) 选取满足顶点在相邻平滑尺度下的曲率差的值同时大于或者同时小于其所 有邻居顶。
17、点的曲率差的值的顶点为特征点 ; 0013 (3) 创建特征点描述算子, 用于描述特征点周围的信息 : 0014 (3-1) 选取每个特征点的邻居顶点, 对于每一个特征点 P, 选取以该点为球心、 r 为 半径的球形区域内的所有顶点为其邻居顶点, 并确定该点的正则方向 ; 0015 (3-2) 将所有的邻居顶点的法线投影到过 P 点的切平面上, 同时将投影后的法线 向量归纳到含 h 项的带权重的直方图中 ; 0016 (3-3) 划分若干个子区域, 确定每个子区域的形状索引信息直方图与倾斜角信息 直方图 : 将球形区域划分为 q 个子区域, 每个子区域使用两个直方图 ps和 p用于构造描述 算。
18、子, 其中, ps的值通过公式 Si 2/tan-1(ci, +ci, )/(ci, -ci, ) 得到, ci, 与 ci, 分 别表示最大、 最小曲率, p定义为投影后的法线向量与正则方向间的倾斜角 ; 0017 (3-4)将所有子区域的两个直方图值连接在一起得到特征点的描述算子f : 将q个 区域的直方图连接成向量, 得到 f ps, 1p, 1ps, qp, qT, 作为特征点的描述算子 ; 0018 其中, r,h,q 根据实际构建特征点的描述算子维度的大小来调整选取 ; 0019 (二)根据步骤(1)中确定的待测人脸的三维数据特征构建三维人脸数据库D : 对 于每一个对象, 建立其。
19、对应的特征池 Si F1, F2,Fj, 所述特征池中包含所有该对象历 次注册后经计算而得的特征点的描述算子集, 则三维人脸数据库 D 中, 将所有对象的特征 池根据其对应的身份序号依次排序而构成, 即 D S1, S2, Sc, 其中 C 为三维人脸数据 库 D 中注册成功的对象个数 ; 0020 ( 三 ) 比较待测人脸与数据库中数据信息以确定待测人脸的身份 : 0021 (I) 多特征点稀疏表示 : 将待识别人脸的特征点的描述算子集 Fprobe, 使用多特征 点稀疏表示方法求解方程式 Fprobe DX 中与 D 对应的系数矩阵 X, 由于 D 的列数太大, 对于 每一个 f 我们通过。
20、提取 D 中与 f 最相似的前 L 个特征描述算子构建 D 的子字典 A, 矩阵 X 通 过 n 次并行计算满足 f Ax 来确定 ; 0022 ( ) 稀疏滤波, 筛选待测人脸的特征点的描述算子集对应的系数向量 : 若待检测 人脸中的一个特征描述算子 f 与子字典 A 中的一个特征描述算子对应的系数 x 的各个数值 分布的非常离散, 则此特征描述算子 f 视为噪音向量, 将其排除出进一步的人脸身份识别 参考范围 ; 0023 ( ) 误差重构, 确定待识别人脸身份 : 分别确定筛选得到的待测人脸特征描述算 子集与子字典中各类的关联程度, 选取最大的关联度的对象作为待测人脸的身份, 关联度 通。
21、过确定重构误差而得到, 重构误差越小, 关联度越高。 0024 所述步骤 ( 一 ) 中, 顶点数据 V 进行平滑操作如下 : 说 明 书 CN 103984920 A 5 3/8 页 6 0025 (1-1)确定所需创建的尺度空间的个数k, 并选取用于创建k个尺度空间的基准值 0; 0026 (1-2)分别确定k个尺度空间对应的估计值, 第s个尺度空间对应的的估计 值表示 E 的平均值 ; 0027 (1-3) 初始化临时变量 0 0, 依次确定 k 组平滑后的顶点数据 : 在创建第 s 组 顶点数据的过程中, 使用第 s-1 组的顶点数据 Vs-1作为第 s 组顶点数据 Vs的初始值, 当。
22、确 定第一组顶点数据v1时使用输入的顶点数据作为初始值, 并执行以下循环, 迭代更新Vs: 若 则更新c以及Vs中所有顶点的坐标, 令cc+1, 令 表示更新后的第 s 组顶点数据中的第 i 个顶点的坐标, Ni表示所有与第 i 个顶点相邻的 顶点构成的集合, Vi表示第 s 组顶点数据中的第 i 个顶点的当前的坐标 ; 若则得到 第 s 组顶点数据 Vs, 退出更新 Vs的循环, 利用 Vs继续确定 Vs+1。 0028 所述步骤 (3-1) 中, 选取直方图中的最高峰值以及大于最高峰值的 t的峰值作 为正则方向 ; 0029 所述步骤 (3-2) 中, 直方图中各项权重为到点 P 的测地距。
23、离的高斯权重。 0030 所述步骤 (3-4) 中, r 取值为 9s, h 取值为 360, t 取值为 80, q 取值为 9。 0031 所述步骤 (I) 中, 通过确定 f 与 D 中的每一列的内积而确定相似度。 0032 所述步骤 (I) 中, L 根据实际运算速度需要在 200 到 1000 内选取。 0033 所述步骤 ( ) 中, 系数向量 x 的稀疏度通过公式 获得, 其中,表示 x 中最大的前 u的元素之和, l 表示子字典 A 中的对象的个数。 0034 重构误差的经由公式确定, 其中, ni为筛 选后的系数向量的个数,意为选取 中只与类别 j 相关的系数。 0035 由。
24、于采用上述方案, 本发明的有益效果是 : 0036 步骤 ( 一 ) 在尺度空间中提取三维人脸扫描数据中的特征点并对其进行描述, 由 于所选取的特征描述算子通过局部区域内的邻居顶点信息得到, 因此即使产生人脸表情变 化后, 特征描述算子中包含的特征信息可以被很好地保留下来。 得益于此, 在三维人脸识别 过程中, 可以排除表情的干扰, 使用提取的多个计算特征描述算子进行相似度比较。 通过本 步骤, 对于输入的三维人脸数据, 可以提取出多个特征点, 从而构成该人脸的特征描述算子 集合。 0037 步骤 ( 二 ) 根据步骤 ( 一 ) 中确定的待测人脸的三维数据特征构建三维人脸数据 库。为了保证该。
25、数据库具有支持多人多次注册的性质, 步骤 ( 二 ) 对于每一次注册的三维 人脸数据标记上所属的用户编号与注册次数, 进而可以满足三维人脸数据的注册、 识别、 验 证等基本需求。 0038 步骤 ( 三 ) 运用步骤 ( 一 ) 提取得到待测人脸的特征描述算子集合, 在步骤 ( 二 ) 中所确定的三维人脸数据库比较以确定待测人脸身份的方法。首先, 通过多特征点稀疏表 说 明 书 CN 103984920 A 6 4/8 页 7 示的方法, 为待测人脸的每一个特征描述算子选取数据库中最接近的若干特征描述算子作 为子字典, 从而大幅度减少后续步骤中的运算量。其次, 进行稀疏滤波, 筛选待测人脸的特。
26、 征点的描述算子集对应的系数向量。 根据稀疏表示的理论, 如果一个系数矩阵具有区分性, 则其系数向量中的非零元素, 应集中出现在对应某一个对象所选取的特征描述算子在子字 典中的位置上。反之, 若系数向量的非零元素均匀分布, 则该系数向量缺少区分性, 并不将 其作为人脸进一步身份识别的参考, 从而提高三维人脸的识别效果。重构误差的大小反映 了待测样本与已知标记样本集合之间的差距, 重构误差越小, 说明待测样本与标记样本集 合越接近。因此, 运用重构误差描述待测人脸与数据库的关联程度, 选取关联程度最大, 即 重构误差最小的对象作为待测人脸的身份。 0039 通过执行以上三个步骤并根据实际需求稍作。
27、调整, 即可实现三维人脸扫描数据的 注册、 识别和验证等基本功能, 同时, 也解决了三维人脸扫描数据中存在的人脸部分缺失、 表情变化大、 面部被遮挡等问题。 附图说明 0040 图 1 为本发明基于稀疏表示与局部多特征点的三维人脸识别方法的工作流程图。 具体实施方式 0041 以下结合附图所示实施例对本发明作进一步的说明。 0042 针对三维人脸扫描数据中存在的人脸部分缺失、 表情变化大、 面部被遮挡等问题, 本发明提供一种基于稀疏表示与局部多特征点的三维人脸识别方法, 具体流程如图 1 所 示 : 0043 ( 一 ) 确定待测人脸的三维数据特征 : 0044 (1) 平滑处理。对于每一张输。
28、入的人脸数据, 我们以网状结构 M 表示, 其中应该包 含的是顶点数据 V, 边缘数据 E, V 包含各顶点在三维空间中的具体坐标, E 包含所有两 两相连的顶点对。在此之上, 我们对 M 构建平滑尺度空间 : 在建立 k 个尺度空间的过程中, 选取 k 个不同大小的 值, 随后对顶点数据 V 进行一系列类似二维图像卷积的平滑操 作, 得到 k 组顶点数据集。具体过程如下 : 0045 (1-1)确定所需创建的尺度空间的个数k, 并选取用于创建k个尺度空间的基准值 0, 其特征在于 : k 可取为 5, 0可取为 5 ; 0046 (1-2)分别确定k个尺度空间对应的估计值, 第s个尺度空间对。
29、应的的估计 值表示 E 的平均值 ; 0047 (1-3)初始化临时变量00, 依次确定k组平滑后的顶点数据 : 在创建第s组顶 点数据的过程中, 使用第s-1组的顶点数据Vs-1作为第s组顶点数据Vs的初始值(当确定第 一组顶点数据 V1时可使用输入的顶点数据作为初始值 ), 并执行以下循环, 迭代更新 Vs: 若 则更新c以及Vs中所有顶点的坐标, 令cc+1令 表示更新后的第 s 组顶点数据中的第 i 个顶点的坐标, Ni表示所有与第 i 个顶点相邻的 说 明 书 CN 103984920 A 7 5/8 页 8 顶点构成的集合, Vi表示第 s 组顶点数据中的第 i 个顶点的当前的坐标。
30、 ; 反之, 若 则得到第 s 组顶点数据 Vs, 退出更新 Vs的循环, 利用 Vs继续确定 Vs+1; 0048 (2) 提取特征点 : 0049 (2-1) 分别确定平滑后的 k 组顶点数据集中各顶点的平均曲率即 k 组平滑 后的顶点数据集中的每个顶点的曲率, 第 s 组顶点数据集中的第 i 个顶点的平均曲率 其中与分别代表该顶点在第 s 个尺度空间中的最大和最小 曲率 ; 0050 (2-2) 分别确定输入的人脸网状结构中所有顶点在相邻平滑尺度下的曲率差 0051 (2-3) 选取满足顶点在相邻平滑尺度下的曲率差的值同时大于或者同时小于其所 有邻居顶点的曲率差的顶点为特征点 ; 005。
31、2 (3) 创建特征点描述算子, 用于描述特征点周围的信息 : 0053 (3-1) 选取每个特征点的邻居顶点, 对于每一个特征点 P, 选取以该点为球心、 r 为 半径的球形区域内的所有顶点为其邻居顶点, 并确定该点的正则方向 ; 本实施例中, 选取直 方图中的最高峰值以及大于最高峰值的 t的峰值作为正则方向。 0054 (3-2) 将所有的邻居顶点的法线投影到过 P 点的切平面上, 同时将投影后的法线 向量归纳到含h项的带权重的直方图中 ; 本实施例中, 直方图中的各项权重为到点P的测地 距离的高斯权重。 0055 (3-3) 划分若干个子区域, 确定每个子区域的形状索引信息直方图与倾斜角。
32、信息 直方图 : 将球形区域划分为 q 个子区域, 每个子区域使用两个直方图 ps和 p用于构造描述 算子, 其中, ps的值通过公式 Si 2/tan-1(ci, +ci, 2)/(ci, 1-ci, 2) 得到, ci, 1与 ci, 2分别 表示最大、 最小曲率, p定义为投影后的法线向量与正则方向间的倾斜角 ; 0056 (3-4)将所有子区域的两个直方图值连接在一起得到特征点的描述算子f : 将q个 区域的直方图连接成向量, 得到作为特征点的描述算子 ; 0057 r,h,t,q 可根据实际构建特征描述算子维度的大小来调整选取, 本实施例中, r 可 取为 9s, h 可取为 360。
33、, t 可取为 80, q 可取为 9。 0058 (二)根据步骤(1)中确定的待测人脸的三维数据特征构建三维人脸数据库D。 对 于任意一个给定的三维人脸数据, 若其满足步骤 A 中所阐述的条件, 则可以通过如上步骤, 可得出一组极值特征描述算子 F f1,f2,fnT用于表示该张人脸。 0059 对于构造的三维人脸数据库, 其应满足支持多人注册, 同时也需要支持每人多次 注册的需求。为此, 我们建议为每一个对象, 即每一个人, 建立其对应的特征池 S, 该特征池 中包含所有该对象历次注册后, 经计算机计算而得的极值特征算子集。 比如, 对象Hi曾在数 据库中有 j 次注册记录, 则其特征池 。
34、Si F1, F2, ., Fj。若三维人脸数据库 D 中有 c 个 注册成功的对象, 我们建议将所有对象的特征池根据其对应的身份序号依次排序而构成, 即 D S1, , ., Sc。 0060 ( 三 ) 比较待测人脸与数据库中数据信息以确定待测人脸的身份 : 说 明 书 CN 103984920 A 8 6/8 页 9 0061 ( ) 多特征点稀疏表示, 为了比较待测身份的人脸与三维数据人脸数据库, 我们 将该人脸的极值特征描述算子集Fprobe, 通过使用多特征点稀疏表示方法, 求解方程式probe DX中与D对应的系数矩阵X, 并根据稀疏表示理论进一步确定待测人脸的身份。 其中, 求。
35、 解方程的过程为 : 对于 Fprobe中的每一个特征描述算子 f, 分别求解范数 1 最小化的优化解, 即满足 f Dx。 0062 由于 D 的列数太大, 通常情况下超过了百万, 为了优化计算效率, 对于每一个 f 我 们通过提取 D 中与 f 最相似的前 L 个特征描述算子构建 D 的子字典 A。其中, 确定相似度的 方法可通过确定 f 与 D 中的每一列的内积而确定。由于 A 的列数是固定值, 矩阵 X 通过数 次并行计算满足 f Ax 而完成, 使得整个计算过程更加有效。在实际操 作中, L 可根据实际运算速度需要在 200 到 1000 内选取, 本实施例中 L 取值为 400。 。
36、0063 ()稀疏滤波, 筛选待测人脸的特征描述算子集对应的系数向量。 为了提高识别 的准确性, 我们在进一步运算前, 根据待测人脸的特征描述算子集对应的系数向量, 保留其 中具有区分性的系数向量, 排除缺乏区分性的向量。 0064 根据稀疏表示的理论, 如果一个系数矩阵具有区分性, 则其在本步骤中得到的系 数向量 x 中的非零元素, 应集中出现在某一个对象在子字典 A 中的所有选取的特征描述算 子在 x 中所对应的位置上 ; 反之, 若系数向量 x 的非零元素均匀分布, 则该系数向量缺少区 分性, 则不将其作为人脸进一步身份识别的参考。本实施例中, 系数向量 x 的稀疏度可以通 过计算公式获。
37、得, 其中,表示 x 中 最大的前 u的元素之和, l 表示子字典 A 中的对象的个数。在实际操作中, u 可取 5。 0065 ()误差重构, 确定待测人脸身份。 通过上一步骤筛选得到的待测人脸特征描述 算子集, 分别计算其与子字典中各类的关联程度, 选取最大的关联度的对象作为待测人脸 的身份。关联度通过计算重构误差而得到, 重构误差越小, 关联度越高。重构误差的计算公 式为其中, ni为筛选后的系数向量的个数, 意为选取中只与类别 j 相关的系数, 即保留系数向量 x 中, 对应于子字典 A 中所选取的、 属于类别 j 的极值特征描述算子, 在系数向量 x 中对应位置的系数值, 将 x 中。
38、其余位置的元 素的值都设为 0。 0066 以下结合具体实验对本发明的有益效果进行说明 : 0067 实验一 : Bosphorus 基准数据库中包含 105 人的 4666 张人脸深度扫描数据。这些 数据中包含面部表情变化、 头部旋转以及遮挡问题。 0068 在实验中, 对于数据库中的每一人, 各选取三张正脸的扫描数据作为已知样本集 合, 其余的正脸扫描数据用于测试。实验采用一阶识别率作为识别的评判标准, 表 1 总结了 包括 meshSIFT 等方法在内的其他方法的实验结果, 其中, Dibeklioglu 等提出的方法以及 Alyuz 等提出的方法针对表情变化的问题, 在面部的部分区域 。
39、( 主要是鼻子区域 ) 中, 具有 很好的鲁棒性。 此外, 表1中所列的迭代最近点算法(ICP)以及基于主成分分析(PCA)处理 深度图像的方法是在 Bosphorus2.0 基准数据库中验证的, 该数据库只包含 47 个人的 2491 张面部扫描图像, 且所有采样都是正脸。 0069 说 明 书 CN 103984920 A 9 7/8 页 10 0070 0071 表 1 0072 从表 1 中可以看出, 在仅包含正脸和使用各类人脸数据情形下, 本方法的一阶识 别率都高于相比较的方法。 由此可以说明, 应用本方法识别三维人脸, 可以克服个体采样中 的表情变化、 头部旋转、 遮挡等问题, 同。
40、时对比现有方法, 取得了更高的准确率。 0073 实验二 : GavabDB 基准数据库是包含表情最丰富、 噪音最多的三维人脸数据库。通 过激光扫描的方式采集了 61 人的三维人脸数据。对于每个人, 分别采集不同的头部旋转及 不同面部表情在内的 9 张人脸扫描数据。 0074 在实验中, 去除了两类极大的头部旋转情形下的三维人脸扫描数据 (90 ), 并 针对每个人选取三张正常表情的三维人脸数据用于构造已知样本集合。 实验采用两种情形 分别作测试。一种为选取所有剩余的样本用于测试, 另一种为选取剩余的自然表情的样本 作测试。实验同样选取一阶识别率作为评判标准。表 2 总结了本方法以及其他一些具。
41、有代 表性的算法的实验结果。 0075 说 明 书 CN 103984920 A 10 8/8 页 11 0076 表 2 0077 从表 2 中可以看出, 本方法相比其余方法所产生的实验结果, 具有明显的优势。特 别是在识别仅包含自然表情的人脸数据中, 其一阶准确率达到了 100。 0078 实验三 : FRGC2.0 基准数据库包含 4007 张, 在受控制的光照条件下, 采集得到的三 维深度扫描数据。扫描对象共有 466 人。 0079 在实验中, 分别从每人所采集的扫描数据中随机抽取三张, 组成已知样本集合。 对 于采集数量不满三张的人, 则选取其全部的扫描数据, 一并归入已知样本集合。
42、中。 在构建完 已知样本集合之后, 选取所有剩余的人脸作为测试集合。表 3 总结了本方法与 meshSIFT 方 法在上述实验条件下取得的实验结果。 0080 0081 表 3 0082 从表 3 中可以看出, 本方法在同等实验条件下, 所取得的一阶识别率高于 meshSIFT方法。 需要注意的是, 实验是在未执行任何预处理的条件下展开的。 实际上, 目前 的一些较好的方法在 FRGC2.0 上验证前, 通常需要执行复杂的数据预处理过程, 比如填补 三维人脸数据中的孔洞, 以改善数据的质量, 因此会取得更好的结果。 在不执行预处理的条 件下, 本方法在 FRGC2.0 基准数据库中取得识别率仍能接近 90。 0083 上述的对实施例的描述是为便于该技术领域的普通技术人员能理解和使用本发 明。熟悉本领域技术的人员显然可以容易地对这些实施例做出各种修改, 并把在此说明的 一般原理应用到其他实施例中而不必经过创造性的劳动。 因此, 本发明不限于上述实施例, 本领域技术人员根据本发明的揭示, 不脱离本发明范畴所做出的改进和修改都应该在本发 明的保护范围之内。 说 明 书 CN 103984920 A 11 1/1 页 12 图 1 说 明 书 附 图 CN 103984920 A 12 。