发明名称 基于步态和二维判别分析的煤矿井下人员身份鉴别方法
摘要 基于步态和二维判别分析的煤矿井下人员身份识别方法,先进行煤矿井下人员的步态图像预处理,然后利用二维判别分析算法对步态能量图进行维数约简,提取所有训练步态能量图的识别特征,并将该识别特征输入模板数据库,每个步态能量图的识别特征与预先存入模板数据库中的煤矿井下人员的身份信息相对应,最后使用最近邻分类器进行煤矿井下人员身份识别,本发明采用了二维判别分析的步态能量图维数约简方法,不仅能够揭示步态图像数据集的内在分布和几何结构,而且能够得到较多的辨别信息,同时它对于步态图像拍摄时的光线变化、步态图像之间的位移旋转等问题具有较强的容错能力,识别速度快,识别效果稳定,实用性强。
申请公布号 CN103593651A 申请公布日期 2014.02.19
申请号 CN201310521289.7 申请日期 2013.10.28
申请人 西京学院 发明人 袁方;王旭启;张善文
分类号 G06K9/00(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 西安智大知识产权代理事务所 61215 代理人 贺建斌
主权项 1.基于步态和二维判别分析的煤矿井下人员身份识别方法,其特征在于,包括以下步骤:第一步,煤矿井下人员的步态图像预处理:由摄像机拍摄到的步态图像直接输入到计算机,拍摄没有步态的图像作为背景模型,将拍摄到的步态图像和背景模型相减后得到消除背景后的步态图像序列;然后使用一个高和宽之比为120:80的矩形框来框住步态轮廓,矩形的高为人体步态图像轮廓的高度,矩形水平方向的中心为人体步态图像轮廓质心的水平坐标,将该矩形框截取出来,并按1:1的缩放比例归一化为120×80大小,再进行二值化处理,得到二值化的步态图像序列B<sub>t</sub>(x,y),确定一个人行走时的连续2个最低高度之间的时间差作为一个步态周期,利用下式(1)将一个周期的步态图像经过平均的方法合成为一幅步态能量图,<maths num="0001"><![CDATA[<math><mrow><mi>G</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>B</mi><mi>t</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,N为完整步态周期序列的长度,t为时间,x、y为二维图像平面坐标;第二步,利用二维判别分析算法对步态能量图进行维数约简,提取所有训练步态能量图的识别特征,并将该识别特征输入模板数据库,每个步态能量图的识别特征与预先存入模板数据库中的煤矿井下人员的身份信息相对应,特征提取按以下步骤进行:设K类n幅步态能量图为{G<sub>1</sub>,G<sub>2</sub>,...,G<sub>n</sub>},类别编号用C<sub>1</sub>,C<sub>2</sub>,...,C<sub>K</sub>表示,并设C<sub>i</sub>类有n<sub>i</sub>个幅图像,计算所有步态能量图和C<sub>i</sub>类步态能量图的平均值分别为<maths num="0002"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mover><mi>G</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>G</mi><mi>i</mi></msub></mtd></mtr><mtr><mtd><mover><msub><msub><mi>G</mi><mi>C</mi></msub><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><msub><mi>n</mi><mi>i</mi></msub></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>i</mi></msub></munderover><msub><mi>G</mi><mi>j</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>计算类间散度矩阵S<sub>B</sub>和类内散度矩阵S<sub>W</sub>分别为<maths num="0003"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>S</mi><mi>B</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><mi>P</mi><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><mover><msub><msub><mi>G</mi><mi>C</mi></msub><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>-</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><msup><mrow><mo>(</mo><mover><msub><msub><mi>G</mi><mi>C</mi></msub><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>-</mo><mover><mi>G</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mi>T</mi></msup></mtd></mtr><mtr><mtd><msub><mi>S</mi><mi>W</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>n</mi><mi>i</mi></msub></munderover><mi>P</mi><mrow><mo>(</mo><msub><mi>C</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><mrow><mrow><mo>(</mo><msub><mi>G</mi><mi>j</mi></msub><mo>-</mo><mover><msub><msub><mi>G</mi><mi>C</mi></msub><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>G</mi><mi>j</mi></msub><mo>-</mo><mover><msub><msub><mi>G</mi><mi>C</mi></msub><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow><mi>T</mi></msup><mo>/</mo><msub><mi>n</mi><mi>i</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,P(C<sub>i</sub>)为模式C<sub>i</sub>的先验概率,取<img file="FDA0000403423120000023.GIF" wi="288" he="170" />由S<sub>B</sub>和S<sub>W</sub>建立目标优化函数J(A)=A<sup>T</sup>(S<sub>B</sub>-S<sub>W</sub>)A     (4)式中,A为待求的映射矩阵,T为矩阵转置运算,对式(4)中的S<sub>B</sub>-S<sub>W</sub>进行特征值分解,计算(S<sub>B</sub>-S<sub>W</sub>)a=λa的d个最大的特征值λ<sub>0</sub>,λ<sub>1</sub>,L,λ<sub>d-1</sub>(λ<sub>0</sub>≥λ<sub>1</sub>≥L≥λ<sub>d-1</sub>)对应的约简特征向量a<sub>1</sub>,a<sub>2</sub>,...,a<sub>d</sub>,其中λ为特征值,d为低维特征向量的维数,由下式(5)将[a<sub>1</sub>,a<sub>2</sub>,...,a<sub>d</sub>]转换为一个正交矩阵P=[p<sub>1</sub>,p<sub>2</sub>,...,p<sub>d</sub>],<maths num="0004"><![CDATA[<math><mrow><msub><mi>p</mi><mi>k</mi></msub><mo>=</mo><msub><mi>a</mi><mi>k</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></munderover><mfrac><mrow><msubsup><mi>p</mi><mi>i</mi><mi>T</mi></msubsup><msub><mi>a</mi><mi>k</mi></msub></mrow><mrow><msubsup><mi>p</mi><mi>i</mi><mi>T</mi></msubsup><msub><mi>p</mi><mi>i</mi></msub></mrow></mfrac><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>由下式(6)将所有的步态能量图{G<sub>1</sub>,G<sub>2</sub>,...,G<sub>n</sub>}分别映射到低维子空间,G<sub>i</sub>→Y<sub>i</sub>=P<sup>T</sup>G<sub>i</sub>     (6)由{Y<sub>1</sub>,Y<sub>2</sub>,...,Y<sub>n</sub>}构建煤矿井下人员身份识别的特征模板,由式(6)将待测试的任意一幅步态能量图G<sub>new</sub>映射为低维特征图像,即G<sub>new</sub>→Y<sub>new</sub>=P<sup>T</sup>G<sub>new</sub>;第三步,使用最近邻分类器进行煤矿井下人员身份识别:利用欧氏距离度量作为任意两个矩阵A、B之间的距离,即d(A,B)=||A-B||,由第二步得到所有训练步态能量图{G<sub>1</sub>,G<sub>2</sub>,...,G<sub>n</sub>}的低维映射矩阵分别为{Y<sub>1</sub>,Y<sub>2</sub>,...,Y<sub>n</sub>},对于任意一个待测试的步态能量图G<sub>new</sub>的低维映射矩阵Y<sub>new</sub>,分别计算Y<sub>new</sub>与{Y<sub>1</sub>,Y<sub>2</sub>,...,Y<sub>n</sub>}之间的距离,若<img file="FDA0000403423120000031.GIF" wi="602" he="98" />则将测试步态能量图G<sub>new</sub>判为第k类,由此确定待鉴别煤矿井下人员的身份。
地址 710123 陕西省西安市长安区西京路1号