发明名称 基于监督等度规投影的人脸识别方法
摘要 本发明提供的是一种基于监督等度规投影的人脸识别方法。包括人脸样本训练过程和人脸样本测试过程。人脸样本训练过程有包括首先对人脸训练图像进行预处理、采用Gabor小波对图像进行滤波、提出新的距离公式计算训练样本的邻接矩阵、由训练样本邻接矩阵DG计算训练样本间的最短路径距离矩阵D、求取描述人脸训练样本数据的低维投影矩阵、通过投影转换矩阵A计算训练样本在低维空间的投影等步骤;人脸样本测试过程又包括对人脸测试图像进行预处理、采用Gabor小波对图像进行滤波、计算测试样本在低维空间的投影和采用最近邻算法判断测试样本的类别等过程。本发明具有对样本数据的结构描述更强,可以消除高阶冗余、计算代价小,更适合于模式分类任务等特点。
申请公布号 CN101673348B 申请公布日期 2012.05.09
申请号 CN200910073077.0 申请日期 2009.10.20
申请人 哈尔滨工程大学 发明人 张汝波;王庆军;刘冠群;徐东;杨歌;史长亭;刘佰龙;张子迎
分类号 G06K9/66(2006.01)I 主分类号 G06K9/66(2006.01)I
代理机构 代理人
主权项 1.基于监督等度规投影的人脸识别方法,其特征是:(1)人脸样本训练过程①首先对人脸训练图像进行预处理,获取高维空间中原始训练样本矩阵X=[x<sub>1</sub>,x<sub>2</sub>…x<sub>N</sub>],N为训练样本个数;②采用Gabor小波进行滤波得到训练样本的小波特征;所述Gabor小波对图像进行滤波,得到新的图像特征,处理过程是,将原始训练图像I(x,y)与Gabor小波滤波器进行卷积计算,二维Gabor小波滤波器定义为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>=</mo><mfrac><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><msup><mi>&sigma;</mi><mn>2</mn></msup></mfrac><msup><mi>e</mi><mrow><mo>(</mo><mo>-</mo><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>k</mi><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><msup><mrow><mo>|</mo><mo>|</mo><mi>z</mi><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>/</mo><msup><mrow><mn>2</mn><mi>&sigma;</mi></mrow><mn>2</mn></msup><mo>)</mo></mrow></msup><mrow><mo>(</mo><msup><mi>e</mi><mrow><msub><mi>ik</mi><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mi>z</mi></mrow></msup><mo>-</mo><msup><mi>e</mi><mrow><mo>-</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>/</mo><mn>2</mn></mrow></msup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>k</mi><mrow><mi>u</mi><mo>,</mo><mi>v</mi></mrow></msub><mo>=</mo><msub><mi>k</mi><mi>v</mi></msub><msup><mi>e</mi><msub><mi>i&phi;</mi><mi>u</mi></msub></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,u,v表示Gabor核的方向和尺度,z=(x,y)为坐标值,||·||表示向量模;k<sub>v</sub>=k<sub>max</sub>/f<sub>v</sub>,φ<sub>u</sub>=πu/8;k<sub>max</sub>为最高采样频率,f<sub>v</sub>为频率域中的采样步长;得到新的图像特征是将图像进行Gabor小波滤波,若x<sub>i</sub>′为滤波后的新特征,将x<sub>i</sub>′堆叠为一M维长向量,即构成一个原始训练样本,所有训练图像构成一个原始训练样本矩阵X=[x<sub>1</sub>,x<sub>2</sub>…x<sub>N</sub>]∈R<sup>M×N</sup>,N为训练样本个数,矩阵X的每一列代表了一个原始训练样本;③引入训练样本的类别信息,采用距离公式计算训练样本的邻接矩阵,得到训练样本的权值距离,所述距离公式为:<maths num="0003"><![CDATA[<math><mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msqrt><mn>1</mn><mo>-</mo><msup><mi>e</mi><mfrac><mrow><msup><mrow><mo>-</mo><mi>d</mi></mrow><mn>2</mn></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mi>&beta;</mi></mfrac></msup></msqrt></mtd><mtd><msub><mi>C</mi><mi>i</mi></msub><mo>=</mo><msub><mi>C</mi><mi>j</mi></msub></mtd></mtr><mtr><mtd><msqrt><msup><mi>e</mi><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow></mrow><mi>&beta;</mi></mfrac></msup></msqrt></mtd><mtd><msub><mi>C</mi><mi>i</mi></msub><mo>&NotEqual;</mo><msub><mi>C</mi><mi>j</mi></msub></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>在得到样本点权值距离矩阵后,应用k近邻或ε邻域算法得到训练样本的邻接矩阵D<sub>G</sub>,其中d(x<sub>i</sub>,x<sub>j</sub>)为样本间的欧式距离,β为经验值,C<sub>i</sub>为样本x<sub>i</sub>所属的类别;④由训练样本邻接矩阵D<sub>G</sub>计算训练样本间的最短路径距离矩阵D,用两点间最短路径近似两点间测地线距离;⑤计算投影转换矩阵,通过保持映射前后人脸训练样本测地线距离给出目标函数,通过求解目标函数求取描述人脸训练样本数据的低维投影矩阵;所述计算投影转换矩阵的具体过程为:函数f为高维空间到低维空间的投影函数,y<sub>i</sub>为x<sub>i</sub>的低维投影,y<sub>i</sub>=f(x<sub>i</sub>),目标函数为保持映射前后训练样本的测地线距离、即最短路径,即:<maths num="0004"><![CDATA[<math><mrow><msup><mi>f</mi><mi>opt</mi></msup><mo>=</mo><mi>arg</mi><mi>min</mi><munder><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></munder><msup><mrow><mo>(</mo><msub><mi>D</mi><mi>ij</mi></msub><mo>-</mo><mi>d</mi><mrow><mo>(</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>,</mo><mi>f</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>其中d(f(x<sub>i</sub>),f(x<sub>j</sub>))为对应低维空间的欧式距离;令S<sub>ij=</sub>D<sub>ij</sub><sup>2</sup>,采用与原始ISOMAP算法相同的处理方式,用中心化矩阵<img file="FSB00000710198200022.GIF" wi="277" he="120" />中心化后得τ(D)=-HSH/2;D<sub>Y</sub>表示降维后子空间中的欧式距离矩阵,τ(D<sub>Y</sub>)为对应的内积矩阵,则上述的目标函数等价于:<img file="FSB00000710198200023.GIF" wi="552" he="79" />其中,<img file="FSB00000710198200024.GIF" wi="471" he="94" />考虑映射函数f(x)=a<sup>T</sup>x,则τ(D<sub>Y</sub>)=Y<sup>T</sup>Y=X<sup>T</sup>aaT<sup>X</sup>;则<img file="FSB00000710198200025.GIF" wi="503" he="68" />式的目标函数转化为:X[τ(D<sub>G</sub>)]X<sup>T</sup>a=λXX<sup>T</sup>a对公式X[τ(D<sub>G</sub>)]X<sup>T</sup>a=λXX<sup>T</sup>a求解得到的前d个最大特征值对应的特征向量构成了投影转换矩阵A=[a<sub>1</sub>,a<sub>2</sub>…a<sub>d</sub>];⑥通过投影转换矩阵A计算训练样本在低维空间的投影,令Y=[y<sub>1</sub>,y<sub>2</sub>…y<sub>N</sub>]为训练样本的低维投影,则Y=A<sup>T</sup>X;(2)人脸样本测试过程①对人脸测试图像进行预处理,获取高维空间中原始测试样本x;②采用Gabor小波进行滤波得到测试样本的小波特征;③计算测试样本在低维空间的投影,通过投影转换矩阵A,测试样本x的低维投影为y=A<sup>T</sup>x;④采用最近邻算法判断测试样本的类别,在低维空间中对y与Y采用最近邻算法,则与y最近的训练样本所对应的类别即为测试样本类别。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号1号楼哈尔滨工程大学科技处知识产权办公室