发明名称 一种人脸图像识别方法
摘要 本发明属于图像处理和模式识别技术领域,具体涉及一种人脸图像识别方法,该方法主要包括如下步骤,步骤a:构建人脸样本图像数据库;步骤b:构建人脸样本图像的训练样本矩阵;步骤c:训练样本矩阵<img file="DDA00002395550900011.GIF" wi="42" he="55" />的近似分解,在该步骤中,在基矩阵W中增加常数矩阵C,并在损失函数上,将系数矩阵H不同列之间的方差作为罚项;步骤d:人脸图像识别过程。由于在基矩阵W中增加光滑常数矩阵C,从而增强了基矩阵的光滑性,消弱了噪声点的影响,使得迭代过程更加快速,大大减少迭代次数,另外将系数矩阵H不同列之间的方差作为罚项,增大系数矩阵H不同列之间的区分度,更好地区分不同人脸图像,提高了人脸识别的准确率。
申请公布号 CN102930258A 申请公布日期 2013.02.13
申请号 CN201210457794.5 申请日期 2012.11.13
申请人 重庆大学 发明人 葛永新;杨梦宁;马鹏;张小洪;洪明坚;徐玲;胡海波
分类号 G06K9/00(2006.01)I;G06K9/66(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 重庆博凯知识产权代理有限公司 50212 代理人 张先芸;梁展湖
主权项 1.一种人脸图像识别方法,其特征在于:具体包括如下步骤:步骤a:构建人脸样本图像数据库:在人脸库中任意选取n张人脸图像,对选取的人脸图像归一化为a×b像素的人脸样本图像;步骤b:构建人脸样本图像的训练样本矩阵:b1:把人脸样本图像转化为高维向量,将步骤a构建的人脸样本图像数据库中的n张人脸样本图像分别转化为m维向量,且m=a×b;b2:通过b1把人脸样本图像转化为高维向量后得到人脸样本图像的训练样本矩阵V,所述训练样本矩阵V大小为m×n,记为<img file="FDA00002395550600011.GIF" wi="365" he="47" />其中<img file="FDA00002395550600012.GIF" wi="28" he="47" />代表训练样本矩阵V中第i列的列向量,且<img file="FDA00002395550600013.GIF" wi="149" he="48" />i=1,2,…,n,R<sup>m</sup>表示m维欧式空间,该训练样本矩阵V的每一列代表一张人脸样本图像;步骤c:训练样本矩阵V的近似分解:c1:输入特征维数r和常数矩阵C,其中<maths num="0001"><![CDATA[<math><mrow><mi>C</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mn>10</mn><mi>m</mi></mrow></mfrac><msub><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mrow><mi>m</mi><mo>&times;</mo><mi>r</mi></mrow></msub><mo>,</mo></mrow></math>]]></maths>且r满足(m+n)r&lt;mn,并由计算机随机生成非负的基矩阵W和系数矩阵H,其中,基矩阵W的大小为m×r,系数矩阵H的大小为r×n,且基矩阵W满足<img file="FDA00002395550600015.GIF" wi="206" he="111" />其中w<sub>ik</sub>为基矩阵第i行,第k列的元素;c2:定义损失函数f,如式(1):<maths num="0002"><![CDATA[<math><mrow><mi>f</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msub><mi>v</mi><mi>ij</mi></msub><mo>-</mo><msub><mi>u</mi><mi>ij</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mi>&alpha;</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msup><mrow><mo>(</mo><msub><mi>h</mi><mi>kj</mi></msub><mo>-</mo><msub><mi>l</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mi>&beta;</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>r</mi></munderover><msubsup><mi>h</mi><mi>kj</mi><mn>2</mn></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>式(1)中l<sub>k</sub>表示系数矩阵H第k行的均值,且<img file="FDA00002395550600017.GIF" wi="259" he="118" />α和β为常数,且α,β∈[10<sup>-3</sup>,10<sup>-6</sup>],v<sub>ij</sub>表示训练样本矩阵V的第i列,第j行的元素;U表示中间矩阵,且U=(W+C)H,u<sub>ij</sub>表示矩阵U中第i行,第j列的元素;c3:采用乘性迭代法求解式(1),具体迭代方式如下:c31:固定系数矩阵H不变,以<img file="FDA00002395550600021.GIF" wi="398" he="130" />更新元素w<sub>ik</sub>;c32:在更新完w<sub>ik</sub>的基础上,固定基矩阵W不变,以<maths num="0003"><![CDATA[<math><mrow><msub><mi>h</mi><mi>kj</mi></msub><mfrac><mrow><msub><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>W</mi><mo>+</mo><mi>C</mi><mo>)</mo></mrow><mi>T</mi></msup><mi>V</mi><mo>)</mo></mrow><mi>kj</mi></msub><mo>+</mo><mi>&alpha;</mi><msub><mi>h</mi><mi>kj</mi></msub></mrow><mrow><msub><mrow><mo>(</mo><msup><mrow><mo>(</mo><mi>W</mi><mo>+</mo><mi>C</mi><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><mi>W</mi><mo>+</mo><mi>C</mi><mo>)</mo></mrow><mi>H</mi><mo>)</mo></mrow><mi>kj</mi></msub><mo>+</mo><mi>&beta;</mi><msub><mi>h</mi><mi>kj</mi></msub><mo>+</mo><mfrac><mi>&alpha;</mi><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>h</mi><mi>kj</mi></msub></mrow></mfrac></mrow></math>]]></maths>更新元素h<sub>kj</sub>;c33:按所述式(1)计算更新完各个元素w<sub>ik</sub>和各个元素h<sub>kj</sub>后所对应的损失函数f;重复执行步骤c31~c33所述的迭代过程,直至满足<img file="FDA00002395550600023.GIF" wi="257" he="118" />的条件时停止迭代;其中,f<sub>x</sub>和f<sub>x-1</sub>分别表示第x次迭代计算的损失函数值和第x-1次迭代计算的损失函数值,x≥2;记录停止迭代时更新产生的新的基矩阵W′和新的系数矩阵H',由所述新的基矩阵W'和新的系数矩阵H'近似地分解训练样本矩阵V,即如式(2)所示:V≈W′H′    (2)记<img file="FDA00002395550600024.GIF" wi="465" he="59" />其中<img file="FDA00002395550600025.GIF" wi="33" he="60" />代表H’中第i列的列向量,且<img file="FDA00002395550600026.GIF" wi="146" he="60" />i=1,2,…,n,R<sup>r</sup>代表r维欧式空间,故此,训练样本矩阵V的第i列的列向量对应于系数矩阵H'的第i列列向量<img file="FDA00002395550600027.GIF" wi="58" he="61" />步骤d:人脸图像识别过程:d1:采集获取人脸图像测试样本,对人脸图像测试样本归一化为a×b像素的测试样本图像;d2:把测试样本图像转化为高维向量,将步骤d1获得的测试样本图像转化为m维向量,m=a×b,得到测试样本图像的高维向量记为<img file="FDA00002395550600028.GIF" wi="60" he="48" />且<img file="FDA00002395550600029.GIF" wi="150" he="48" />t=1,2,…,n;d3:将步骤d2得到的测试样本图像的高维向量<img file="FDA000023955506000210.GIF" wi="29" he="47" />投影到步骤c33中新的基矩阵W’上,得到该测试样本图像的系数向量h<sub>t</sub>,其中<img file="FDA000023955506000211.GIF" wi="191" he="54" />且h<sub>t</sub>∈R<sup>r</sup>,R<sup>r</sup>表示m维欧式空间;d4:分别计算测试样本图像的系数向量h<sub>t</sub>和系数矩阵H'各列的列向量<img file="FDA000023955506000212.GIF" wi="37" he="61" />之间的欧式距离,i=1,2,…,n,找出其中的欧式距离最小值对应的列向量<img file="FDA00002395550600031.GIF" wi="37" he="60" />作为测试样本图像的系数向量h<sub>t</sub>的识别匹配向量,进而判定测试样本图像与所述识别匹配向量对应的人脸样本图像来自于同一人。
地址 400044 重庆市沙坪坝区沙正街174号