CN102142083A - 一种基于lda子空间学习的人脸识别方法 - Google Patents

Abstract

本发明公开了一种基于LDA子空间学习的人脸识别方法，其特征在于：其包括以下步骤：获取人脸图像，并进行预处理；计算所述人脸图像的梯度多尺度局部相位量化GMLPQ特征集；应用Adaboost选择器筛选其中具有鉴别能力的特征组成候选特征子集；应用LDA子空间分析器，得到人脸特征模板；将所述人脸特征模板与预建人脸特征模板库进行匹配，获取识别人身份信息。本方法所述人脸识别技术比现有其它人脸识别技术具有更强的环境适应性，在图像模糊（失焦、运动等）、低分辩率、各种光照条件（红外、可见光）下具有较好的识别率和误识率，并且计算速度快，特别适合于嵌入式产品，可大规模推广应用。

Description

一种基于LDA子空间学习的人脸识别方法

技术领域

本发明属于人脸识别技术，尤其是基于LDA子空间识别方法的人脸识别技术。

背景技术

人脸识别技术是当前大力发展的生物识别技术之一。人脸识别系统主要包括数据采集子系统、人脸检测子系统和人脸识别子系统。人脸特征提取是人脸识别子系统最为关键的技术，好的人脸特征提取技术将使提取的人脸特征值更小、辨别性能更好，可以提高识别率和降低误识率。目前已存在的人脸特征提取方法主要有：基于几何特征方法、基于子空间分析方法、基于小波理论方法、基于神经网络方法、基于隐马尔可夫模型方法、基于支持向量机方法和基于三维模型方法。基于几何特征方法的思想是提取人脸面部具有代表性的部位(例如眉毛、眼睛、鼻子、嘴巴等)的相对位置和相对大小作为特征，再辅以人脸轮廓的形状信息作为特征，该方法容易受光照、表情、遮挡等因素的影响，稳定性不高。基于小波理论的人脸识别方法主要思想是人脸图像在经过小波变换后得到的低频图像可用于表示人脸。人工神经网络ANN是由大量简单的处理单元相互连接而成的网络系统，在自学习、自组织、联想及容错方面具有较强的能力，在学习过程中提取得到的特征可作为人脸特征进行识别。基于子空间分析的方法是目前比较流行的人脸识别方法，基本思想是把高维空间中松散分布的人 脸图像通过线性或非线性变换投影到一个低维的子空间中，使得人脸图像的分布在低维子空间中更加紧凑，更有利于分类，并且从高维计算变为低维计算可解决“维数灾难”问题。线性子空间方法有：主元分析PCA、奇异值分解SVD、线性判别分析LDA、独立主元分析ICA和非负矩阵因子NMF等，非线性子空间方法有：核主成份分析、核Fishe判别分析、流形学习方法等。子空间分析的方法具有计算代价小、描述能力强、分类性强等特点，基于可分性准则的线性判别分析方法是当前人脸识别的主流方法之一。

发明内容

本发明提供了一种改善上述问题的方案，提供一种性能更好、鲁棒性更强的应用于人脸识别方法。

本发明的技术方案是提供一种基于LDA子空间学习的人脸识别方法，其特征在于：其包括以下步骤：

1)获取人脸图像，并对所述人脸图像进行转正、滤波、规定化分辩率等预处理；

2)计算1)中所述人脸图像的梯度多尺度局部相位量化GMLPQ特征集；

3)应用Adaboost选择器筛选2)中所述梯度多尺度局部相位量化GMLPQ特征集中具有鉴别能力的特征组成候选特征子集；

4)应用LDA子空间分析器对3)中所述候选特征子集进行分析，得到一个低维特征向量作为人脸特征模板；

5)将4)中所述人脸特征模板与预建人脸特征模板库进行匹配， 获取识别人身份信息。

优选的，步骤4)中所述LDA子空间分析器基于改进的目标函数训练而成，所述改进的目标函数如下：

$J=w^T{S}_{w}w-\mathrm{\λ}{w}^T{S}_{b}w$

$+\mathrm{\η}(\underset{i=1}{\overset{C}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{C}{\mathrm{\Σ}}}{d}_{\mathrm{ij}}^2-C(C-1){\stackrel{\‾}{d}}^2)$

λ，η：权重系数

优选的，所述改进的目标函数基于梯度下降法或梯度共轭法进行优化求解。

传统的子空间学习方法(如LDA)中，通常考虑寻找一个子空间，使得类内样本离散度小，类间样本离散度大。在解决实际问题过程中，为了满足优化目标函数，经常会有两个类别的样本聚集在一起的情况发生，而这种现象对于分类来说是极为不利的。本发明的创新点就是提出了一种解决方法，使得在求解得到的子空间中，在除了保持原来的优化目标函数外，尽可能避免这种不同类别的样本聚集到一起的情况，从而提高分类正确率。本方法所述人脸识别技术比现有其它人脸识别技术具有更强的环境适应性，在图像模糊(失焦、运动等)、低分辩率、各种光照条件(红外、可见光)下具有较好的识别率和误识率，并且计算速度快，特别适合于嵌入式产品，可大规模推广应用。

附图说明

图1是本发明的算法原理框图。

具体实施方式

下面对本发明的具体实施方式作进一步详细的描述。

如图1所示，本发明所述应用于人脸识别的改进目标函数的LDA子空间学习方法步骤是：人脸图像获取、提取GMLPQ特征集、Adaboost选择器、LDA子空间分析器，最后进行人脸特征比对。

以下结合附图1所示的算法原理图，详细说明本方法的具体实施方式。

①获取人脸图像，并进行转正、滤波、规定化分辩率等预处理。

②计算①中所述人脸图像的梯度多尺度局部相位量化(GMLPQ)特征集。

GMLPQ特征提取原理：

GMLPQ特征就是基于梯度图像提取MLPQ特征，梯度图像包括水平梯度图像和垂直梯度图像。

MLPQ特征提取原理：

MLPQ特征就是多尺度多频域的LPQ融合特征。

对于图像来说，假设图像受到一定噪声的模糊影响，在频域上可以表示为：

G＝F·H

其中F是原来图像的傅立叶变换，G是模糊后的图像的傅立叶变换，H是模糊函数(点扩散函数的傅立叶变换)。他们的幅度和相位分量分别满足：

|G|＝|F|·|H|

∠G＝∠F+∠H

假设点扩散函数在空域中是中心对称的，则H是处在实数域的，也就是在频域中H的相位为0或者pi。又假设在低频区域，H的值为正，所以F和G的具有共同的相位。LPQ就是基于这个原理提出的，它对图像模糊具有比较好的鲁棒性。MLPQ是对LPQ的一个扩展，采用多个不同的尺度来提取不同频域下的LPQ特征，从而能够更好地刻画人脸。

LPQ提取流程：

(1)将图像分成多个重叠的一定大小的小块。

(2)在每个小块中分别进行短时傅立叶变换。

提取某些特定频率下的相位信息，对其进行编码，得到LPQ特征。原始的LPQ特征提取方法计算量较大，中间存在着比较多的重复计算，难以达到特征的实时提取。本发明提出一种基于积分图的快速LPQ特征提取方法。

LPQ特征编码：

(1)每个相位象限用2个bit来进行编码，例如00表示第一象限，01表示第二象限，10表示第三象限，11表示第四象限。将4个不同频域的相位编码连接在一起，得到8个bit的二进制串，然后转换成十进制的0-255的编码。

(2)原始的LPQ编码维数较高，且易受到噪声的影响。采用典型模式的方法，从样本中统计出出现频率最高的K个LPQ码，将 剩下的码全部合并到一个码中，从而减少了LPQ特征码维度，加快了计算速度，同时保证了LPQ码的鲁棒性。

③应用Adaboost选择器筛选②中所述梯度多尺度局部相位量化(GMLPQ)特征集中具有鉴别能力的特征组成候选特征子集。

本方法所述算法采用Adaboost训练特征选择器。Adaboost是一种迭代算法，其核心思想是针对同一个训练集训练不同的分类器(弱分类器)，然后把这些弱分类器集合起来，构成一个更强的最终分类器(强分类器)。其算法本身是通过改变数据分布来实现的，它根据每次训练集之中每个样本的分类是否正确，以及上次的总体分类的准确率，来确定每个样本的权值。将修改过权值的新数据集送给下层分类器进行训练，最后将每次训练得到的分类器最后融合起来，作为最后的决策分类器。使用adaboost分类器可以排除一些不必要的训练数据特征，将注意力放在关键的训练数据上面。

Adaboost训练过程如下：

1.先通过对N个训练样本的学习得到第一个弱分类器；

2.将分错的样本和其他的新数据一起构成一个新的N个的训练样本，通过对这个样本的学习得到第二个弱分类器；

3.将1和2都分错了的样本加上其他的新样本构成另一个新的N个的训练样本，通过对这个样本的学习得到第三个弱分类器；

4.最终经过提升得到强分类器。

④应用LDA子空间分析器对③中所述候选特征子集进行分析，得到一个低维特征向量作为人脸特征模板。

所述LDA子空间由多个子空间线性辨别分析器组合而成，包括：原始LDA、增强LDA(E-FLDA)、直接LDA(D-LDA)、零空间LDA(N-LDA)、边缘LDA(MFA)。在每种方法的过程中，都随机选择训练样本，使得不同方法具有更多的互补性，提高模型泛化能力。

LDA子空间分析方法简述。

LDA子空间分析方法就是线性鉴别分析方法，目标是从高维特征空间里提取出具有鉴别能力的低维特征，这些特征有助于将属于同一个类的样本更加聚集在一起，属于不同类的样本更加分开。数学描述方式如LDA目标函数定义，就是找到一个投影矩阵W，使得类内散度矩阵Sw与类间散度矩阵Sb的比值最大化。

类内散度矩阵Sw定义：

$S_w=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1}{\overset{N_k}{\mathrm{\Σ}}}(X_{\mathrm{ij}}-m_i){(X_{\mathrm{ij}}-m_i)}^T$

类间散度矩阵Sb定义：

$S_b=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{N}{\mathrm{\Σ}}}(m_i-m_j){(m_i-m_j)}^T$

LDA目标函数定义：

$J=\frac{w^T{S}_{b}w}{w^T{S}_{w}w}$

基于LDA子空间的人脸识别方法是一个多类学习问题，一个人就是一个类，同一个人的人脸样本属于同一个类，不同人的人脸样本属 于不同类。LDA子空间人脸识别方法就是学习得到一个低维子空间，使得人脸特征点从高维特征空间投影到该低维子空间中，属于同一个人的人脸特征点更加聚集，属于不同人的人脸特征点更加分开。对于多类学习问题，一个存在的问题就是投影到低维子空间后类间距离并不均匀，并且某些类可能存在交叠，这样使得人脸识别中不同人特征比对分数出现异常抖动现象，增大了误识率。传统的LDA目标函数学习仅仅解决了两个类间的分离，而解决不了多类的类间距保持均匀问题。本发明改进了传统LDA目标函数，并提出基于梯度下降法或共轭梯度法的优化类间距的学习方法，使得学习得到的低维子空间不仅能够很好的分开不同类点，而且使得不同类中心之间距离更加均匀，有效降低了人脸误识率。

本发明改进LDA目标函数的方法，就是在保证类内样本更加聚集和类间样本更加分开的同时尽可能的使不同类中心之间距离更加均匀，用数学描述方式就是，在保证类内离散度最小和类间离散度最大的同时尽可能的使不同类均值向量之间的距离方差更加均匀，从而避免某几类样本聚集在一起的现象，有效降低人脸误识率。

假设dij为第i类中心与第j类中心的距离，则所有类间距离平均值为：

$\stackrel{\‾}{d}=\underset{i=1}{\overset{C}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{C}{\mathrm{\Σ}}}\frac{1}{C(C-1)}{d}_{\mathrm{ij}}$

类间距离分布目标函数：

$J^{\′}={\mathrm{\Σ}}_{i=1}^C{\mathrm{\Σ}}_{j=i+1}^C{(d_{\mathrm{ij}}-\stackrel{\‾}{d})}^2={\mathrm{\Σ}}_{i=1}^C{\mathrm{\Σ}}_{j=i+1}^C(d_{\mathrm{ij}}^2+{\stackrel{\‾}{d}}^2-2\stackrel{\‾}{d}{d}_{\mathrm{ij}})=$

${\mathrm{\Σ}}_{i=1}^C{\mathrm{\Σ}}_{j=i+1}^C{d}_{\mathrm{ij}}^2+C(C-1){\stackrel{\‾}{d}}^2-2C(C-1){\stackrel{\‾}{d}}^2={\mathrm{\Σ}}_{i=1}^C{\mathrm{\Σ}}_{j=i+1}^C{d}_{\mathrm{ij}}^2-$

$C(C-1){\stackrel{\‾}{d}}^2$

本发明改进的LDA目标函数：

$J=w^T{S}_{w}w-\mathrm{\λ}{w}^T{S}_{b}w$

$+\mathrm{\η}(\underset{i=1}{\overset{C}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{C}{\mathrm{\Σ}}}{d}_{\mathrm{ij}}^2-C(C-1){\stackrel{\‾}{d}}^2)$

λ，η：权重系数

应用梯度下降法或共轭梯度法等方法优化改进的LDA目标函数，对各个目标函数J求偏导。改变投影矩阵w，使得目标函数J向着梯度降低最快的方向走，在确定全局最佳的目标函数时获取最佳的投影矩阵w。

$\frac{\∂J}{\∂w}=S_{w}w-\mathrm{\λ}{S}_{b}w$

$+\mathrm{\η}(\underset{i=1}{\overset{C}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{C}{\mathrm{\Σ}}}{d}_{\mathrm{ij}}\frac{\∂d_{\mathrm{ij}}}{\∂w}-\stackrel{\‾}{d}\frac{\∂\stackrel{\‾}{d}}{\∂w})$

$\frac{\∂d_{\mathrm{ij}}}{\∂w}=\frac{1}{2\sqrt{{(m_i-m_j)}^{T}w{w}^T(m_i-m_j)}}(m_i-m_j){(m_i-m_j)}^{T}w$

$\frac{\∂\stackrel{\‾}{d}}{\∂w}=\underset{i=1}{\overset{C}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{C}{\mathrm{\Σ}}}\frac{\∂d_{\mathrm{ij}}}{\∂w}$

LDA传统目标函数与优化目标函数投影效果相比较，LDA传统目标函数学习到的低维子空间为水平超平面，本发明LDA改进目标函数学习到的低维子空间为垂直超平面。在水平超平面的投影会出现蓝色点与红色点交叠的情况，而在垂直超平面的投影不会出现不同类点交叠 的情况，说明本发明改进的LDA目标函数使得学习到的低维子空间中的不同类点之间距离更加均匀，抑制了不同类点出现交叠的情况，降低了误识率。

⑤将④中所述人脸特征模板与预建人脸特征模板库进行匹配，获取识别人身份信息。

以上详细说明了本发明所述应用于人脸识别的改进目标函数的LDA子空间学习方法的执行过程及其实现细节。

以上实施例仅为本发明其中的一种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (3)

1.一种基于LDA子空间学习的人脸识别方法，其特征在于：其包括以下步骤：

1)获取人脸图像，并对所述人脸图像进行转正、滤波、规定化分辩率等预处理；

2)计算1)中所述人脸图像的梯度多尺度局部相位量化GMLPQ特征集；

3)应用Adaboost选择器筛选2)中所述梯度多尺度局部相位量化GMLPQ特征集中具有鉴别能力的特征组成候选特征子集；

4)应用LDA子空间分析器对3)中所述候选特征子集进行分析，得到一个低维特征向量作为人脸特征模板；

5)将4)中所述人脸特征模板与预建人脸特征模板库进行匹配，获取识别人身份信息。

2.根据权利要求1所述一种基于LDA子空间学习的人脸识别方法，其特征在于：步骤4)中所述LDA子空间分析器基于改进的目标函数训练而成，所述改进的目标函数如下：

$J=w^T{S}_{w}w-\mathrm{\λ}{w}^T{S}_{b}w$

$+\mathrm{\η}(\underset{i=1}{\overset{C}{\mathrm{\Σ}}}\underset{j=i+1}{\overset{C}{\mathrm{\Σ}}}{d}_{\mathrm{ij}}^2-C(C-1){\stackrel{\‾}{d}}^2)$

λ，η：权重系数

3.根据权利要求2所述一种基于LDA子空间学习的人脸识别方法，其特征在于：所述改进的目标函数基于梯度下降法或梯度共轭法进行优化求解。