[发明专利]一种低秩分块稀疏表示的人脸识别方法

摘要

发明公开了一种低秩分块稀疏表示的人脸识别方法，采用低秩矩阵分解并引入参考项，采用DCT算法，实现了图像的归一化，有效地解决了人脸图像中光照不均匀的问题。在分类阶段，利用聚类的思想，行之有效地提高了识别速度。将该算法用于标准人脸数据库上进行多次实验，实验结果表明：与现有人脸识别算法相比，本文算法的识别准确率和计算效率都得到了一致提高。在人脸图像出现遮挡、伪装、光照变化等复杂情况下以提高人脸识别的精度和稳健性。

主权项

1.一种低秩分块稀疏表示的人脸识别方法，其特征在于，具体步骤如下：步骤一、对于人脸图像数据库中的每一个主体，随机选择其一部分作为训练图像，另一部分作为测试图像，将所有主体的训练图像和测试图像分别整合构成初始的训练数据矩阵和测试矩阵；步骤二、将训练数据矩阵D分解成A+E，其中A代表低秩分解矩阵，E代表分解后的稀疏误差，通过最小化低秩分解矩阵A的秩，同时减小零范数||E||₀的值来达到训练数据矩阵D的最佳低秩逼近，低秩分解矩阵分解的公式：minA,E||A||*+λ||E||1s.t.D=A+E---(1)]]>公式（1）中，核范数||A||_*是低秩分解矩阵A秩的近似值，一范数||E||₁是零范数||E||₀的替代值，λ为参数；步骤三、引入参考项，根据低秩分解矩阵分解的公式（1），建立目标函数：minA,EΣi=1c{||Ai||*+λ1||Ei||1}+λ2Ψ(A1,A2,....,Ac)s.t.Di=Ai+Ei---(2)]]>公式（2）中，i＝1,2,...,c，c为训练数据矩阵中类的个数，D_i为第i个训练数据矩阵，A_i为第i个低秩分解矩阵，E_i为第i个稀疏误差矩阵，Ψ(A₁,A₂,...,A_c)为提高低秩分解矩阵区分能力的参考项，参数λ₁是正的权重系数，参数λ₂为常数且λ₂≥0；步骤四、对目标函数进行化简，将目标函数最终整理为：minAi,Ei||Ai||*+λ1||Ei||1+λ2||Ai-mi||F2s.t.Di=Ai+Ei---(10)]]>其中d为人脸图像化成列向量后的行数，m_i代表第i个主体的均值向量，n_i代表第i个主体中样本的个数，||·||_F代表是是f范数；步骤五、利用增广拉格朗日乘法对公式（10）进行处理：L(Ai,Ei,Yi,μ,λ2)=||Ai||*+λ1||Ei||1+λ2||Ai′||F2+<Yi,Di-Ai-Ei>+u2||Di-Ai-Ei||F2---(12)]]>其中，A_i'＝A_i-M_i，Y_i代表第i个主体所对应的拉格朗日乘数，目的是通过迭代求取函数的极值，μ代表一个正数参数；步骤六、对于公式（12）通过迭代求解A，对低秩分解矩阵A进行更新：Aik+1=argminAiL(Ai,Eik,Yikμk,λ2)=argminAi||Ai||*+λ2||Ai′||F2+<Yik,Di-Ai-Eik>+μk2||Di-Ai-Eik||F2=argminAi||Ai||*+(λ2+μk2)<Ai,Ai′>-μk<Di-Eik+1μkYik,Ai′>=argminAiϵ||Ai||*+12||Xa-Ai||F2]]>Aik+1=USϵVTwhere(U,S,VT)=SVD(Xa)]]>式中ε＝(2λ₂+μ^k)^-1，k代表迭代次数，代表第i个稀疏误差矩阵经过k次迭代后的值，Y_i^k代表第i个主体所对应的拉格朗日乘数经过k次迭代后的值，μ^k代表正数参数经过k次迭代后的值，通过调节Y_i^k的值控制迭代过程，SVD代表矩阵奇异值分解，其中U、V为酉矩阵，S为对角矩阵；步骤七、对误差矩阵E进行更新：Eik+1=argminEiL(Aik+1,Ei,Yik,μk,λ2)=argminEi||Ei||1+<Yik,Aik+1+Ei-Di>+μk2||Aik+1+Ei-Di||F2=argminEiϵ′||Ei||1+12||Xe-Ei||F2]]>式中ϵ′=1μk,Xe=Di-Aik+1+1μkYik;]]>步骤八、对经过低秩分解矩阵分解的A应用离散余弦变换，实现光照归一化处理；步骤九、对训练图像进行重叠分块；步骤十、对照步骤二至步骤九训练图像的分块方式，将测试图片也相应地进行分块处理；步骤十一、利用基于稀疏表示的人脸识别算法，结合l₁范数最小化，求解模块的稀疏系数：xp=argminx{||Apx-yp||22+λ||x||1}for1≤p≤12]]>式中，x为稀疏系数，x_p代表第p个模块相应的稀疏系数，y_p代表测试图片第p个模块所组成的列向量；步骤十二、利用不同模块的稀疏系数，将不同模块同一类稀疏系数合并在一起，再对同一类稀疏系数从大到小进行排序，对于一个有c个主体的训练数据矩阵，每个类有N个训练样本的训练集，它的测试图片第p个模块所对应的稀疏系数为：x_p∈R^cN×1，1≤p≤12v_pj＝δ_j(x_p)₁v_pj代表测试图片第p个模块第j个主体所对应的同一类的稀疏系数绝对值之和，δ_j(x_p)代表在x_p中与第j个主体有关联的特征函数，n为训练样本的总数，R^cN×1是一个cN维的列向量，1≤p≤12；因此：v_p＝[v_p1,v_p2,...,v_pc]^T将各个模块的稀疏系数联系起来有：f=112Σp=112vp]]>f∈R^c×1，R^c×1是一个c维的列向量，通过求解f中最大值的坐标确定测试图片属于哪一个主体，从而实现正确分类与识别。