[发明专利]基于混合误差编码的人脸识别方法

摘要

一种基于混合误差编码的人脸识别方法，为了排除遮挡和光照变化所形成的局部显著特征的影响，利用遮挡和光照变化的空间结构，即连续性和局部方向性，构造结构化误差，以选取待识别图像的有效特征；为了弥补遮挡和光照变化所造成的特征损失，由稀疏表示理论和图像在高维特征空间中的表示，构造判别误差，以有效地识别人脸图像。结构化误差和判别误差既相互独立又相互影响，本发明将基于这两种误差编码的人脸识别方法称为混合误差编码方法。实验表明，混合误差编码方法，可以有效地识别现实中的人脸图像，为面向现实的人脸识别系统提供了有力的支撑。本发明识别率高、可行性好，特别地，对于有遮挡或光照变化的人脸图像，具有良好的识别性能。

主权项

一种基于混合误差编码的人脸识别方法，其特征在于：所述人脸识别方法包括以下步骤：步骤1将实数域上的m×n维的人脸图像y作为待识别人脸图像，将来自于K个人的N个不含遮挡的且在正常光照条件下采集到的人脸图像集A={a11,...,aN11,a12,...,aN22...a1K,...,aNKK}]]>作为训练样本，这里，为实数域上的m×n阶矩阵，表示第k个人的第i个训练样本，k＝1，2…K，Nk为第k个人的人脸图像的样本数；步骤2将y和A分别变换到梯度方向域中：φ(y)=arctan(y*hdTy*hd),]]>Φ(A)={φ(a11),...,φ(aN11),φ(a12),...,φ(aN22)...φ(a1K),...,φ(aNKK)},]]>φ(aik)=arctan(aik*hdTaik*hd),]]>这里，*表示二维卷积，arctan(·)为反正切函数，hd为差分滤波器的系数矩阵，符号T表示矩阵的转置；步骤3将待识别图像y和训练集A变换到高维特征空间中：其中，表示待识别图像y在高维特征空间中的表示，表示训练集A在高维特征空间中的表示，V(·)为拉伸变换，即：将m×n阶矩阵拉伸为D＝m×n维列向量，或者将所有的m×n阶矩阵的集合A中的元素依次拉伸为D＝m×n维列向量，然后再将拉伸后的列向量组装成矩阵，即：V(cos(Φ(A)))=[V(b11),...,V(bN11),V(b12),...,V(bN22)...V(b1K),...,V(bNKK)]V(sin(Φ(A)))=[V(c11),...,V(cN11),V(c12),...,V(cN22)...V(c1K),...,V(cNKK)],]]>其中，bik=cos(φ(aik)),cik=sin(φ(aik));]]>步骤4初始化待识别人脸图像y相对于训练集A的编码系数为这里，1N×1表示N维列向量，其全部元素为1；步骤5初始化参数：t＝1，ζ*＝+∞，这里，t表示迭代次数，ζ*表示历次迭代过程中产生的判别误差的最小质量；步骤6在像素空间中，按照如下公式，计算待识别图像y的重构图像y‾(t)=Mm×n(V(A)×x(t)),]]>其中，Mm×n(·)是V(·)的逆运算，即把D＝m×n维列向量还原为m×n阶矩阵，x(t)为实数域上的N维列向量，表示在第t次迭代时所产生的待识别人脸图像y相对于训练集A的编码系数；V(A)为集合A的拉伸变换：V(A)=[V(a11),...,V(aN11),V(a12),...,V(aN22)...V(a1K),...,V(aNKK)];]]>步骤7将变换到梯度方向域中：步骤8计算待识别图像y与其重构图像的梯度方向之差△θ(t)：Δθ(t)=φ(y)-φ(y‾(t));]]>步骤9将△θ(t)在cos域中做如下变换：步骤10对进行均值滤波：这里，hm为均值滤波器的系数；步骤11将做如下变换，得到结构化误差：这里，D＝m×n为的维数，|·|表示取绝对值，||·||1表示l1范数；步骤12对结构化误差进行阈值聚类，得到待识别图像y的有效特征支撑的初步估计sτ(t)：siτ(t)=κτ(e~i(t))=1e~i(t)<τ0e~i(t)≥τ,]]>其中，κτ(·)表示阈值滤波算子，sτ(t)和的下标i∈{1，2，…，m×n}表示sτ(t)和中元素的索引，τ为给定的阈值；步骤13通过如下公式，进一步估计待识别图像y的有效特征支撑s(t)：其中，s表示优化式中的与s(t)相对应的有效特征支撑的变元，s和sτ(t)的维数相同，下标i，j∈{1，2，…，m×n}表示s和sτ(t)中元素的索引，表示i的邻域节点，λs为光滑参数，λu为数据参数；步骤14由待识别图像y的有效特征支撑s(t)，计算y在高维特征空间中的表示相对于训练集A在高维特征空间中的表示的稀疏编码系数x(t)，以及与x(t)相关联的判别误差其中，s.t.表示约束条件，||·||2表示l2范数，表示点乘，x表示优化式中的与x(t)相对应的稀疏编码系数的变元，表示优化式中的与相对应的判别误差的变元；步骤15度量判别误差的质量ζ(t)：如果ζ(t)＜ζ*，则令ζ*＝ζ(t)，x*＝x(t)，这里，ζ*，x*，分别记录ζ(t)，x(t)，的最优表示；步骤16令t＝t+1，迭代上述步骤6～步骤15，直至给定的最大迭代次数；步骤17在高维特征空间中，按如下公式，计算待识别图像y在训练集中的每一个人的样本空间中的投影残差rk：其中，k∈{1，2，…，K}，δk(x*)表示将x*中的所有除了对应于第k个人的编码系数外的元素全部置为0；步骤18按如下公式对待识别图像y进行分类：ID(y)=argminkrk.]]>