发明名称 一种基于图嵌入低秩稀疏表示恢复的稀疏表示人脸识别方法
摘要 本发明公开了一种图嵌入低秩稀疏表示恢复稀疏表示人脸识别方法,属于计算机视觉和模式识别技术领域。本发明包括以下步骤:首先,提出一种图嵌入低秩稀疏表示恢复方法,能够从训练样本数据矩阵恢复出判别力强的干净训练样本数据矩阵,同时得到训练样本数据误差矩阵;然后,以干净训练样本数据矩阵为字典,以训练样本数据误差矩阵为误差字典,采用<img file="dest_path_image002.GIF" wi="10" he="17" />范数最优化技术求解待识别人脸数据的稀疏表示系数;更进一步,利用待识别人脸数据的稀疏表示系数,对待识别人脸数据进行类关联重构;最后,基于待识别人脸数据的类关联重构误差,完成待识别人脸图像的识别。本发明能够解决训练样本图像和待识别图像都受噪声污染或局部被遮挡情况下的人脸识别问题。
申请公布号 CN104318261B 申请公布日期 2016.04.27
申请号 CN201410607957.2 申请日期 2014.11.03
申请人 河南大学 发明人 杜海顺;王俊;张延宇;杜晓玉;胡青璞;蒋曼曼
分类号 G06K9/64(2006.01)I;G06K9/00(2006.01)I 主分类号 G06K9/64(2006.01)I
代理机构 郑州联科专利事务所(普通合伙) 41104 代理人 刘建芳;崔卫琴
主权项 一种基于图嵌入低秩稀疏表示恢复的稀疏表示人脸识别方法,其特征在于:包括以下步骤:S01:假设有K类训练样本图像,每类有n张训练样本图像,共N=K×n张训练样本图像,设每张训练样本图像分辨率为r×c,将每张训练样本图像转换为M=r×c维向量,则训练样本数据矩阵记为<img file="FDA0000904911070000011.GIF" wi="574" he="79" />S02:利用训练样本数据矩阵X及其对应的类别标签,构建一个包含N个节点的有监督无向近邻图G;其中,图G的节点连接关系及其权重为:如果数据x<sub>i</sub>与数据x<sub>j</sub>属于同一类,且x<sub>i</sub>为x<sub>j</sub>的k近邻或x<sub>j</sub>为x<sub>i</sub>的k近邻,k为正整数,取值范围3≤k≤n‑1,则图G的节点i和节点j连接,且其连接权重W<sub>ij</sub>为1;否则,图G的节点i和节点j不连接,其权重W<sub>ij</sub>为0,定义数据x<sub>i</sub>的度为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>d</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>W</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>;</mo></mrow>]]></math><img file="FDA0000904911070000012.GIF" wi="614" he="102" /></maths>定义数据度矩阵D为:D=diag(d<sub>1</sub>,d<sub>2</sub>,…,d<sub>N</sub>),其中diag(d<sub>1</sub>,d<sub>2</sub>,…,d<sub>N</sub>)表示以d<sub>1</sub>,d<sub>2</sub>,…,d<sub>N</sub>为对角元素的对角矩阵;S03:根据步骤S02求得的权重矩阵W和数据度矩阵D,定义拉普拉斯矩阵<img file="FDA0000904911070000015.GIF" wi="70" he="55" /><img file="FDA0000904911070000016.GIF" wi="262" he="70" />假设Z=[z<sub>1</sub>,z<sub>2</sub>,…,z<sub>N</sub>]为X的表示系数矩阵,为了使训练样本数据的表示系数能够保持数据空间的局部几何结构并具有强的判别力,基于拉普拉斯矩阵L,定义一个有监督的图嵌入正则项:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mn>1</mn><mn>2</mn></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>|</mo><mo>|</mo><msub><mi>z</mi><mi>i</mi></msub><mo>-</mo><msub><mi>z</mi><mi>j</mi></msub><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><msub><mi>W</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mi>T</mi><mi>r</mi><mrow><mo>(</mo><msup><mi>ZLZ</mi><mi>T</mi></msup><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000904911070000013.GIF" wi="654" he="143" /></maths>式中,Tr(·)表示矩阵的迹;S04:令E为训练样本数据误差矩阵,α>0、β>0、γ>0为正则化常数,定义图嵌入低秩稀疏表示恢复方法目标函数为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>Z</mi><mo>,</mo><mi>E</mi></mrow></munder><mo>|</mo><mo>|</mo><mi>Z</mi><mo>|</mo><msub><mo>|</mo><mo>*</mo></msub><mo>+</mo><mi>&alpha;</mi><mo>|</mo><mo>|</mo><mi>E</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>+</mo><mi>&beta;</mi><mo>|</mo><mo>|</mo><mi>Z</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>+</mo><mfrac><mi>&gamma;</mi><mn>2</mn></mfrac><mi>T</mi><mi>r</mi><mrow><mo>(</mo><msup><mi>ZLZ</mi><mi>T</mi></msup><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><mi>X</mi><mo>=</mo><mi>X</mi><mi>Z</mi><mo>+</mo><mi>E</mi><mo>;</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000904911070000014.GIF" wi="1244" he="110" /></maths>式中||·||<sub>*</sub>表示矩阵的核范数,即矩阵奇异值之和,||·||<sub>1</sub>表示矩阵的l<sub>1</sub>范数,Tr(·)表示矩阵的迹;根据图嵌入低秩稀疏表示恢复方法目标函数求得训练样本数据矩阵的表示系数矩阵Z和训练样本数据误差矩阵E;S05:根据步骤S04求得的训练样本数据矩阵X的表示系数矩阵Z,采用下式恢复干净训练样本数据矩阵C:C=XZ;S06:对于任一待识别人脸图像,将其转化为M维向量,记为待识别人脸图像数据y,<img file="FDA0000904911070000021.GIF" wi="187" he="71" />S07:以C为字典、E为误差字典,采用l<sub>1</sub>范数最优化技术,通过下式求解待识别人脸图像数据y的稀疏表示系数:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>min</mi><msup><mrow><mo>&lsqb;</mo><msup><mi>&alpha;</mi><mi>T</mi></msup><mo>,</mo><msup><mi>&beta;</mi><mi>T</mi></msup><mo>&rsqb;</mo></mrow><mi>T</mi></msup></msub><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><mi>C</mi><mi>&alpha;</mi><mo>-</mo><mi>E</mi><mi>&beta;</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><msup><mrow><mo>&lsqb;</mo><msup><mi>&alpha;</mi><mi>T</mi></msup><mo>,</mo><msup><mi>&beta;</mi><mi>T</mi></msup><mo>&rsqb;</mo></mrow><mi>T</mi></msup><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>;</mo></mrow>]]></math><img file="FDA0000904911070000022.GIF" wi="869" he="134" /></maths>式中,正则化常数λ=0.05;S08:计算各类训练样本对待识别人脸图像数据y的类重构误差:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><mi>y</mi><mo>-</mo><msub><mi>C&delta;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>&alpha;</mi><mo>)</mo></mrow><mo>-</mo><mi>E</mi><mi>&beta;</mi><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>K</mi><mo>;</mo></mrow>]]></math><img file="FDA0000904911070000023.GIF" wi="813" he="87" /></maths>式中,δ<sub>i</sub>(α)是仅与第i类对应系数保留,与其它类对应系数置为0的系数向量;S09:根据类重构误差e(i),i=1,2,…,K,计算待识别人脸图像数据y的类别标签:Class(y)=arg min<sub>i</sub>e(i);式中,Class(y)表示待识别人脸图像数据y的类别标签,argmin<sub>i</sub>e(i)表示值最小的e(i)对应的i;S10:输出待识别人脸图像的识别结果。
地址 475001 河南省开封市明伦街85号