发明名称 应用二维非负稀疏偏最小二乘进行人脸识别的方法
摘要 本发明涉及应用二维非负稀疏偏最小二乘进行人脸识别的方法,它包括如下步骤,首先构造人脸训练样本集的类别矩阵,再构建使投影人脸训练样本集和类别矩阵时信息量损失最少的目标函数,然后在目标函数中加入非负性约束和稀疏性约束得到收敛的非负的基矩阵,将人脸训练样本集投影在基矩阵上获得测试样本系数矩阵,对测试样本也进行前述操作得到测试样本的系数矩阵,使用最近邻策略对测试样本的系数矩阵和测试样本系数矩阵中的某个要素矩阵是同一类,则认为测试样本的系数矩阵所对应的测试样本上的人与该要素矩阵所对应的训练样本上的人为同一人。该方法识别率和鲁棒性,由于只需对基矩阵进行迭代求解,简化了运算,降低了时间复杂度,识别速度快。
申请公布号 CN103971132A 申请公布日期 2014.08.06
申请号 CN201410227645.9 申请日期 2014.05.27
申请人 重庆大学 发明人 葛永新;徐玲;杨梦宁;洪明坚;张小洪;杨丹;步文斌
分类号 G06K9/64(2006.01)I 主分类号 G06K9/64(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 王海凤;穆祥维
主权项 应用二维非负稀疏偏最小二乘进行人脸识别的方法,其特征在于具体包括如下步骤:步骤a:构造人脸训练样本集的类别矩阵Y:a1:构建人脸训练样本集:在人脸库中任意选取N张人脸图像,该N张人脸图像属于m个人,其中每个人有n张人脸图像,并记为n<sub>i</sub>,i=1,...,m,每张人脸图像为一个训练样本,则人脸训练样本集为<img file="FDA0000511930110000014.GIF" wi="760" he="81" />其中x<sub>ij</sub>表示第i个人的第j个训练样本,j=1,...,n,h表示人脸图像的纵向像素数,l表示人脸图像的横向像素数;a2:构造步骤a1中的人脸训练样本集X的类别矩阵Y,<img file="FDA0000511930110000011.GIF" wi="997" he="586" />其中,P<sub>i</sub>对应的是第i个人的n个训练样本,每个训练样本对应大小为l×l的单位矩阵I<sub>l</sub>,因此,<img file="FDA0000511930110000012.GIF" wi="536" he="150" />步骤b:构建使投影人脸训练样本集X和类别矩阵Y时信息量损失最少的目标函数如式(1):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><munder><mi>max</mi><mrow><mi>w</mi><mo>,</mo><mi>c</mi></mrow></munder><mo>&lt;</mo><msup><mi>w</mi><mi>T</mi></msup><mi>X</mi><mo>,</mo><msup><mi>c</mi><mi>T</mi></msup><mi>Y</mi><mo>></mo><mo>=</mo><msup><mi>w</mi><mi>T</mi></msup><msup><mi>XY</mi><mi>T</mi></msup><mi>c</mi></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open='{' close=''><mtable><mtr><mtd><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>=</mo><mi>I</mi></mtd></mtr><mtr><mtd><msup><mi>c</mi><mi>T</mi></msup><mi>c</mi><mo>=</mo><mi>I</mi></mtd></mtr></mtable></mfenced></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000511930110000013.GIF" wi="766" he="289" /></maths>其中w是人脸训练样本集X的投影矩阵,c为类别矩阵Y的投影矩阵,w<sup>T</sup>表示矩阵w的转置,c<sup>T</sup>表示矩阵c的转置,I表示单位矩阵;步骤c:在步骤b构建的目标函数中加入非负性约束得式(3):<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><munder><mi>max</mi><mi>w</mi></munder><msubsup><mrow><mo>|</mo><mo>|</mo><msup><mi>w</mi><mi>T</mi></msup><msup><mi>XY</mi><mi>T</mi></msup><mi>c</mi><mo>|</mo><mo>|</mo></mrow><mi>F</mi><mn>2</mn></msubsup><mo>-</mo><mi>&alpha;</mi><msubsup><mrow><mo>|</mo><mo>|</mo><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>I</mi><mo>|</mo><mo>|</mo></mrow><mi>F</mi><mn>2</mn></msubsup><mo>,</mo></mtd><mtd><mi>I</mi></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msup><mi>c</mi><mi>T</mi></msup><mi>c</mi><mo>=</mo><msup><mi>cc</mi><mi>T</mi></msup><mo>=</mo><mi>I</mi><mo>,</mo><mi>w</mi><mo>&GreaterEqual;</mo><mn>0</mn><mo>,</mo></mtd><mtd><mi>II</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000511930110000021.GIF" wi="1160" he="300" /></maths>其中α为任意正常数将式(3)中表示非负性约束的Ⅱ式代入Ⅰ式得到非负性目标函数如式(6):<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><munder><mi>max</mi><mi>w</mi></munder><msup><mi>w</mi><mi>T</mi></msup><msup><mi>XY</mi><mi>T</mi></msup><msup><mi>YX</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>&alpha;</mi><msup><mrow><mo>(</mo><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>I</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000511930110000026.GIF" wi="968" he="118" /></maths>步骤d:在步骤c得到的非负性目标函数中再加入稀疏性约束:用l<sup>1</sup>范数的约束,对w引入稀疏性约束,式(6)转化为式(7):<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mi>max</mi><msup><mi>w</mi><mi>T</mi></msup><msup><mi>XY</mi><mi>T</mi></msup><msup><mi>YX</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>&alpha;</mi><msup><mrow><mo>(</mo><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>I</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>,</mo></mtd><mtd><mi>I</mi></mtd></mtr><mtr><mtd><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><msub><mrow><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>&le;</mo><mi>&lambda;</mi><mo>,</mo></mtd><mtd><mi>II</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000511930110000022.GIF" wi="1326" he="302" /></maths>其中λ为稀疏程度,将式(7)中表示稀疏性约束的Ⅱ式加入Ⅰ式得到最优目标函数如式(8):<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><munder><mi>max</mi><mi>w</mi></munder><msup><mi>w</mi><mi>T</mi></msup><msup><mi>XY</mi><mi>T</mi></msup><msup><mi>YX</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>&alpha;</mi><msup><mrow><mo>(</mo><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>I</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mi>&beta;</mi><msub><mrow><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>|</mo></mrow><mn>1</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow><mo>:</mo></mrow>]]></math><img file="FDA0000511930110000023.GIF" wi="1296" he="119" /></maths>将式(8)中l<sup>1</sup>约束改为l<sup>2</sup>约束得到最优目标函数如式(9):<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><munder><mi>max</mi><mi>w</mi></munder><msup><mi>w</mi><mi>T</mi></msup><msup><mi>XY</mi><mi>T</mi></msup><msup><mi>YX</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>&alpha;</mi><msup><mrow><mo>(</mo><msup><mi>w</mi><mi>T</mi></msup><mi>w</mi><mo>-</mo><mi>I</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mi>&beta;</mi><msubsup><mrow><mo>|</mo><mo>|</mo><mi>w</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000511930110000024.GIF" wi="1451" he="144" /></maths>其中β为任意的正常数;将l<sup>2</sup>范数展开并且带入式稀疏性约束得式(11);f(w)=w<sup>T</sup>XY<sup>T</sup>YX<sup>T</sup>w‑α(w<sup>T</sup>w‑I)<sup>2</sup>‑βw<sup>T</sup>w   (11);对式(11)利用乘性迭代法得出迭代格式(14);<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msup><mi>w</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msup><mo>&LeftArrow;</mo><msup><mi>w</mi><mi>n</mi></msup><mfrac><mrow><mo>(</mo><msup><mi>XY</mi><mi>T</mi></msup><msup><mi>YX</mi><mi>T</mi></msup><msup><mi>w</mi><mi>n</mi></msup><mo>+</mo><msup><mrow><mn>2</mn><mi>&alpha;w</mi></mrow><mi>n</mi></msup><mo>)</mo></mrow><mrow><msup><mrow><mn>2</mn><mi>&alpha;w</mi></mrow><mi>n</mi></msup><msup><mrow><mo>(</mo><msup><mi>w</mi><mi>n</mi></msup><mo>)</mo></mrow><mi>T</mi></msup><msup><mi>w</mi><mi>n</mi></msup><mo>+</mo><msup><mi>&beta;w</mi><mi>n</mi></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000511930110000025.GIF" wi="1183" he="183" /></maths>步骤e:通过步骤d得到迭代格式(14)迭代最后得到收敛的非负的基矩阵w',将人脸训练样本集X投影在基矩阵w'上获得测试样本系数矩阵t,t=w'<sup>T</sup>X;步骤f:人脸图像识别过程:f1:采集获取人脸图像测试样本,采用与步骤a2相同的方法构造测试样本的类别矩阵;f2:采用与步骤b相同的方法构建使投影测试样本和测试样本的类别矩阵时信息量损失最少的目标函数;f3:在步骤f2构建的目标函数中加入与步骤c相同的非负性约束得到测试样本的非负性目标函数;f4:在步骤f3构建的测试样本的非负性目标函数中加入与步骤d相同的稀疏性约束得到最优目标函数,利用乘性迭代法得到测试样本的迭代格式,然后通过迭代求得测试样本的基矩阵,最后将测试样本投影在其基矩阵上获得测试样本系数矩阵;f5:使用最近邻策略对测试样本的系数矩阵与人脸训练样本集投影在基矩阵上获得测试样本系数矩阵进行分类识别,当测试样本的系数矩阵与人脸训练样本集投影在基矩阵上获得系数矩阵中的某个要素矩阵是同一类时,则认为测试样本的系数矩阵所对应的测试样本上的人与该要素矩阵对应的训练样本上的人为同一人。
地址 400044 重庆市沙坪坝区沙正街174号