一种基于局部特征分析的人脸识别方法

一种基于局部特征分析的人脸识别方法
焦峰 1 山世光 1 崔国勤 1 高文 1,2 李锦涛 1
1
(中国科学院计算技术研究所 ICT-YCNC 面像识别联合实验室，北京，100080) 2 (哈尔滨工业大学计算机学院智能人机接口实验室，哈尔滨，150001)
摘要
本文在传统的弹性图匹配基础上， 提出了一种基于局部特征分析的人脸识别算法。 利用 人脸的先验结构和人脸图像的灰度分布知识， 先粗略地找出图像的特征点， 然后利用人脸弹 性图对特征点的位置进行调整， 最后在各个特征点处计算 Gabor 变化的系数。 整个人脸相应 地被表示成特征点处的 Gabor 系数集合。 我们用几种不同的距离来进行分类， 并给出测试结 果。相对于传统的弹性图匹配方法，由于特征点预先被估算出，而不是在整个图像上搜索， 所以大大减少了计算量。
图 1 自动人脸识别系统框图
对于人脸图像，系统先根据模板匹配和特征子脸相结合的方法检测出人脸, 然后利用人 脸的先验结构和灰度分布信息，粗略定位特征点。为了精确定位特征点，我们用引入一种称 为人脸弹性图的数据结构来调整特征点。特征点被精确定位后，我们计算特征点处的 Gabor 变换系数向量集合, 并以此来表示整个人脸。最后计算测试样本和训练样本的 Gabor 变换系 数向量集合距离，得到并输出结果。 人脸检测是人脸识别的前期工作，本文继承了文献[6],的工作，采用了一种模板匹配和 特征子脸相结合的两级检测模型， 该模型在各种图像条件下可以很好地工作。 由于人脸检测 不是本文的重点，而且由于篇幅限制，本文不做详细描述，感兴趣者可查阅相关文献。 以下各章将详细介绍系统以后每个步骤。
(1)
从而得到两个眼睛区域的粗略点集， 分别计算这两个点集的重心， 以该重心为中心得到一些 候选的瞳孔点，对每一个候选点 P，在边缘图像上计算下述支持函数：
Sp =
1 N
( x , y )∈ A
∑E
( x, y )
(2)
其中：E 为边缘图像，A 为由半径范围 R1 ≤ r ≤ R2 所决定的圆环区域， R1 、 R2 是根据人 瞳孔大小取的一个经验值，N 为其中的点个数。分别取两个候选点集中具有最大支持函数的 点作为粗定位阶段的虹膜中心点：
眼睛宽 0.42 0.03
鼻度 0.57 0.05
嘴宽 0.70 0.06
嘴高 0.27 0.02
表 1 面部特征结构关系统计
对整个矫正后的图像做积分变换，可以发现人脸特征点处有一定的灰度分布特点， 利用这
3
些特点和估计出的特征点位置， 我们可以找出人脸的特征点。
3.3.1 嘴部的检测
观察表明：首先根据得到的统计参数，估计出嘴部区域在垂直方向的大致范围，然后在 此范围内取窗口区域进行水平方向的积分投影。一个典型的过程如图 2、3 所示。实验表明： 除非有很大的噪声或者有浓厚的胡须， 在其他情况下都可以得到一条类似的曲线。 对该曲线 进行高斯平滑，计算曲线的极小点，具有最小值的极小点可以作为中唇点定位结果。高斯平 滑通过卷积实现：
h ( i ) = f (i ) * g (i ) =
∑
j 2
f (i ) g (i − j )
2
(8)
其中 : g ( x ) = exp( − x /( 2 s ))
曲线极小点则通过计算平滑后曲线的一阶和二阶差分得到：
h ′ ( i ) = ( h ( i + 1) − h ( i − 1)) / 2 h ′′ ( i ) = h ( i + 1) + h ( i − 1) − 2 h ( i )
其它特征点可以通过类似的方法得到[6][7]。 本文一共选取了 26 个特征点。其中眼睛处共 10 个点，每个眼睛 5 个点，分别为瞳孔及 上下左右眼角；眉心处 1 个点；瞳孔垂直往上和眉毛交点处左右各 1 个点；鼻子处 3 个点， 分别为鼻尖和两个鼻孔；嘴巴处 5 个点，分别为嘴唇中心和嘴唇上下左右端 4 个点；脸颊处 5 个点， 分别为下巴、 嘴唇水平位置与脸颊左右处的交点、 眼睛水平位置与脸颊左右的交点。 在这些特征点处人脸图像的变化比较大，基本代表了人脸图像的变化。 这些点主要分步在人脸内部，而且人脸的上半部的上半部比下半部居多， 因为根据统 计和实验结果， 在识别过程中内部的特征点比外部的特征点更重要， 上半部的点比下半部的 点更重要[2]。
ϕ j (x ) = k 2 j exp −
r
2 k2 j x 2σ 2

r exp i k j
σ2 r exp − x − 2

(12)
为了平滑地在频域中采样，我们取 5 个频率,对应 v=0,1,2,3,4； 8 个相位，对应 u=0,…,7
2. 人脸自动识别系统框架
1
模板匹配 和特征子 脸
人脸先验 结构和灰 度分布信 息
人脸弹性 图 训 练 样 本 Gabor 变 换 系 数向量集合
训练样本 图像 测试样本 图像 人脸检测 特征点粗 略定位 特征点精 确调整
计算系数向量 集合的距离, 输出结果
测 试 样 本 Gabor 变 换 系 数向量集合
PLeftIris = arg( Max S p )
p∈Rngleft
(3)
(4)
PRightIris = arg( Max S p )
p∈Rng Right
3.2 平面旋转矫正
对于人脸图像可能存在的图像平面内的旋转， 可以通过定位出的两个虹膜中心点的位置 用仿射变换进行矫正。 这里对面部图像的旋转以人脸的中心为原点， 而不是图像的中心为原 点， 这样可以保证当人脸落在图像的边缘部分时不致因为旋转而使人脸部分丢失。 人脸的中 心是利用两个虹膜中心点和人脸结构的统计知识估计得到的。
( ) (
)
(13)
通过对不同方向和不同频率的 Gabor 分解变换， 对于图像上的每一个点可以得到一组系数向 量，包含 40 个变换系数（5 个频率， 8 个相位） ，它可以写成 J j = a j exp iφ j ， 1 ≤ j ≤ 40 我们引入两种相似度，对于系数向量 J和J ,如果不考虑其相位差别，则其相似度公式为
得到极小点集合为：
(9)
S = {i | h ′( i + 1) h ′(i + 1) <= 0 & & h ′′( i ) >= 0}
从而可以得到中唇点：
(10)
m = arg min { f ( i ), ∀ i ∈ S }
i
灰度均值 位置坐标
(11)
图2
图3
图 2-3 嘴部积分投影及定位结果 2 积分区域及结果 3 积分投影曲线
2
实验表明通过上述五条规则的限制， 总可以找到一对满足条件的连通区， 而这对连通区恰好 就是由两个眼睛所形成的区域。在该算法检测得到的区域 L、R 区域的二值图像 B 和二值化 的边缘图像 E 中，利用虹膜区域的谷特征和眼睛区域的高频特性，分别形成两个眼睛的区 域点集描述为：
Rng = ( x, y ) | B( x, y ) = 0, ∑ E (i , j ) > θ , Ω为(x, y)的8邻域 ( i , j )∈Ω
r kj =
( )= (
k jx k jy
cosφu kv sinφu
)
kv = 2
→
−
v+2 π 2
φ
u
= u
π
8
对一个灰度分布为 L x 的图，在给定一点 x = ( x, y ) ，Gabor 函数的变换为
→
r r r r r J j (x ) = ∫ L x ' ϕ j x − x ' d 2 x '
关键字
人脸识别, 特征定位, Gabor 变换系数向量集合, 人脸弹性图
1 引言
人脸识别是用计算机分析人脸图像，从人脸图像中提取有效信息，用于识别身份的技 术,因其在公安（罪犯识别等） 、安全验证系统、信用卡验证、医学、档案管理、视频会议、 人机交互系统等方面的巨大应用前景而越来越成为当前模式识别和人工智能领域的一个研 究热点[1][2]。人脸有复杂的三维表面结构，因而人脸知识的表达相对简单物体来说较为复 杂；同时面部肌肉的运动使得人脸成为一种非刚体目标，与刚性目标识别相比，识别更困 难。此外，所有人脸都具有相似的结构特征，而同一人的不同面部图像由于采集条件的变 化会发生较大的变化（对计算机而言） ，所以人脸识别算法必须挖掘不同类别间微妙而可靠 的差别。人脸识别的一些算法可以推广到一般的三维非刚体目标的识别上去，促进其它模 式识别问题的解决。 人脸识别作为图像匹配的一个子集， 和一般的图像匹配有一些区别。 一般的图像匹配是 从颜色、形状、纹理中提取特征，而人脸匹配是根据人脸的特征来识别。在人脸识别中，如 何描述人脸，即如何提取特征是最重要的一个环节。在现有的人脸识别算法中，比较流行的 有模板匹配法[3]，特征子脸法(即 Eigenface)[4], 基于弹性图匹配方法[5]。 本文提出了一种新的识别算法， 利用人脸的先验结构知识和人脸图像的灰度分布知识， 先粗略找出人脸的特征点， 然后利用人脸弹性图来对其进行调整， 最后在各个特征点处计算 Gabor 变化系数集合, 并以此来表示人脸的特征。相对于传统的弹性图匹配方法，特征点预 先被估算出来，而不是在整个图中搜索，因此大大地降低了计算量。 本文组织如下：第一章为引言，第二章介绍本文自动人脸识别系统框图，第三章介绍人 脸的先验结构知识和人脸图像的灰度分布知识来粗略定位特征点， 第四章介绍利用人脸弹性 图对特征点精确定位，第五章介绍识别过程和实验结果，第六章是结论和展望。
3 特征点定位
人脸图像在灰度空间总是满足一定的分步，我们可以利用这一点来指导进行面部特征点 的粗略定位，以此作为后续特征点精确定位的基础。