发明名称 一种基于加速度计的人体行为识别方法
摘要 本发明公开一种基于加速度计的人体行为识别分类方法,包括如下步骤:1)收集人体行为样本作为训练集;2)寻找对该训练集识别分类最优的投影矩阵U;3)对无标注数据进行投影;4)对投影后的数据采用最小距离分类器分类,获得识别结果。本发明对标注数据形成的近邻块做局部近似线性的假设,并使块上不同类别之间样本距离足够大,相同类别样本位置顺序信息通过类sigmoid函数惩罚因子尽可能的保留,最后在所有块上目标函数的基础上建立全局目标函数。利用本发明提出的方法能够合适的保留高维空间中样本之间距离的信息,减少识别模型对人工标注样本的依赖,其识别效果优于有代表性的基于线性判别分析的人体行为识别方法。
申请公布号 CN103500342B 申请公布日期 2017.01.04
申请号 CN201310428353.7 申请日期 2013.09.18
申请人 华南理工大学 发明人 陶大鹏;金连文;黎小凤
分类号 G06K9/62(2006.01)I 主分类号 G06K9/62(2006.01)I
代理机构 广州市华学知识产权代理有限公司 44245 代理人 蔡茂略
主权项 一种基于加速度计的人体行为识别方法,包括以下步骤:1)收集人体行为样本作为训练集X,即:X=[x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>]∈R<sup>D×N</sup>,样本维数为D,样本个数为N,每个样本有相应的类别标志C<sub>i</sub>∈Z<sup>n</sup>;2)寻找基于加速度计的人体行为识别分类最优的投影矩阵U;3)通过对无标注数据X<sub>u</sub>进行投影,即:Y<sub>u</sub>=U<sup>T</sup>X<sub>u</sub>;4)对Y<sub>u</sub>采用最小距离分类器分类,以获得人体行为识别的结果;其特征在于,所述寻找基于加速度计的人体行为识别分类最优的投影矩阵的方法包括以下步骤:步骤1:建立局部优化目标函数;步骤2:建立全局优化目标函数;步骤3:利用拉格朗日乘数法,投影矩阵U由式子XLX<sup>T</sup>的前d个最小特征值对应的d个特征向量得到;所述步骤1中,建立局部优化目标函数的方法为:对每一个已标注的样本x<sub>i</sub>,找到同类样本<img file="FDA0001022868550000011.GIF" wi="186" he="71" />的k<sub>1</sub>近邻和不同类别样本<img file="FDA0001022868550000012.GIF" wi="202" he="63" />的k<sub>2</sub>近邻来形成一个局部块,即<img file="FDA0001022868550000013.GIF" wi="922" he="94" />另外,定义R<sub>ij</sub>为第j个样本相对第i个样本里的位置顺序;当<img file="FDA0001022868550000014.GIF" wi="893" he="103" />时,满足不同类别样本距离足够大,同时,相同类样本位置顺序信息尽量保留;对类间样本距离建立(1)式:<maths num="0001"><math><![CDATA[<mrow><mi>M</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mn>2</mn></msub></munderover><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><msub><mi>i</mi><mi>p</mi></msub></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000015.GIF" wi="1230" he="158" /></maths>对类内样本的位置顺序信息建立(2)式:<maths num="0002"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mn>1</mn></msub></munderover><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><msup><mi>i</mi><mi>j</mi></msup></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msub><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>j</mi></msub><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000021.GIF" wi="1302" he="158" /></maths>式中,(w<sub>i</sub>)<sub>j</sub>表示区别小距离和大距离之间的惩罚因子,当原始空间的距离小时,赋予低维子空间较大的权重,当距离大时,赋予较小的权重;所述惩罚因子定义为类‑sigmoid函数,所述惩罚因子的表达式如下:<img file="FDA0001022868550000022.GIF" wi="1493" he="222" />式中,f(u|μ,σ)是均值为μ,标准差为σ的高斯概率密度函数,均值μ和标准差σ可以分别由以下两条公式估计得到:<maths num="0003"><math><![CDATA[<mrow><mi>&mu;</mi><mo>=</mo><munder><mrow><mi>m</mi><mi>e</mi><mi>a</mi><mi>n</mi></mrow><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&lt;</mo><mi>j</mi><mo>&le;</mo><mi>N</mi></mrow></munder><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&lambda;</mi><mo>)</mo></mrow><munder><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&lt;</mo><mi>j</mi><mo>&le;</mo><mi>N</mi></mrow></munder><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000023.GIF" wi="1508" he="119" /></maths><maths num="0004"><math><![CDATA[<mrow><mi>&sigma;</mi><mo>=</mo><mn>2</mn><mi>&lambda;</mi><munder><mrow><mi>s</mi><mi>t</mi><mi>d</mi></mrow><mrow><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>&lt;</mo><mi>j</mi><mo>&le;</mo><mi>N</mi></mrow></munder><mrow><mo>(</mo><msub><mi>d</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000024.GIF" wi="1150" he="111" /></maths>式中,d<sub>ij</sub>是高维空间样本之间的距离,参数λ∈[0,1],并且,参数λ通过交叉验证得到;由于局部块X<sub>i</sub>是近似线性的,由公式(1)(2)和一个权衡参数γ,得到局部优化的目标函数(6):<maths num="0005"><math><![CDATA[<mrow><mi>arg</mi><munder><mi>min</mi><msub><mi>y</mi><mi>i</mi></msub></munder><mrow><mo>(</mo><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mn>1</mn></msub></munderover><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><msup><mi>i</mi><mi>j</mi></msup></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msub><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>j</mi></msub><mo>-</mo><mi>&gamma;</mi><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mn>2</mn></msub></munderover><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><msub><mi>i</mi><mi>p</mi></msub></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000025.GIF" wi="1534" he="181" /></maths>式中,γ∈[0,1]是一个用于整合类内样本和类间样本贡献值的权衡系数;公式(6)进一步化为以下形式:<maths num="0006"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><munder><mrow><mi>arg</mi><mi>min</mi></mrow><msub><mi>y</mi><mi>i</mi></msub></munder><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mn>1</mn></msub></munderover><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><msub><mi>i</mi><mi>j</mi></msub></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><msub><mrow><mo>(</mo><msub><mi>w</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>j</mi></msub><mo>-</mo><mi>&gamma;</mi><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>k</mi><mn>2</mn></msub></munderover><mo>|</mo><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><msub><mi>i</mi><mi>j</mi></msub></msub><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><msub><mi>Y</mi><mi>i</mi></msub></munder><mi>t</mi><mi>r</mi><mrow><mo>(</mo><mrow><msub><mi>Y</mi><mi>i</mi></msub><msub><mi>L</mi><mi>i</mi></msub><msubsup><mi>Y</mi><mi>i</mi><mi>T</mi></msubsup></mrow><mo>)</mo></mrow><mo>,</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000026.GIF" wi="1550" he="310" /></maths>式中,tr(·)是迹算子,<img file="FDA0001022868550000027.GIF" wi="883" he="191" /><img file="FDA0001022868550000031.GIF" wi="1852" he="211" /><img file="FDA0001022868550000032.GIF" wi="878" he="239" />所述步骤2中,建立全局优化目标函数的方法为:通过样本选择矩阵,低维空间表达Y<sub>i</sub>的坐标是从全局坐标Y=U<sup>T</sup>X=[y<sub>1</sub>,y<sub>2</sub>,…y<sub>N</sub>]∈R<sup>d×N</sup>中选择出来的,即:Y<sub>i</sub>=YS<sub>i</sub>,                         (8)式中,S<sub>i</sub>∈R<sup>N×(K+1)</sup>是选择矩阵,令F<sub>i</sub>={i,i<sub>1</sub>,…i<sub>K</sub>}为指示集合,则选择矩阵的定义如下:<img file="FDA0001022868550000033.GIF" wi="1294" he="182" />根据公式(9),公式(7)写为:<maths num="0007"><math><![CDATA[<mrow><munder><mi>argmin</mi><mi>Y</mi></munder><mi>t</mi><mi>r</mi><mrow><mo>(</mo><msub><mi>YS</mi><mi>i</mi></msub><msub><mi>L</mi><mi>i</mi></msub><msubsup><mi>S</mi><mi>i</mi><mi>T</mi></msubsup><msup><mi>Y</mi><mi>T</mi></msup><mo>)</mo></mrow><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000034.GIF" wi="1259" he="111" /></maths>通过对公式(10)的局部优化求和,得到整体调整公式(11):<maths num="0008"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>Y</mi></munder><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>t</mi><mi>r</mi><mrow><mo>(</mo><mrow><msub><mi>YS</mi><mi>i</mi></msub><msub><mi>L</mi><mi>i</mi></msub><msubsup><mi>S</mi><mi>i</mi><mi>T</mi></msubsup><msup><mi>Y</mi><mi>T</mi></msup></mrow><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munder><mrow><mi>arg</mi><mi>min</mi></mrow><mi>Y</mi></munder><mi>t</mi><mi>r</mi><mrow><mo>(</mo><mrow><msup><mi>YLY</mi><mi>T</mi></msup></mrow><mo>)</mo></mrow><mo>,</mo></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000035.GIF" wi="1326" he="277" /></maths>令Y=U<sup>T</sup>X,而U<sup>T</sup>U=I<sub>d</sub>,这里I<sub>d</sub>是d×d维的单位矩阵,因此公式(9)写成:<maths num="0009"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><munder><mi>argmin</mi><mi>Y</mi></munder><mi>t</mi><mi>r</mi><mrow><mo>(</mo><msup><mi>U</mi><mi>T</mi></msup><msup><mi>XLX</mi><mi>T</mi></msup><mi>U</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo></mrow></mtd><mtd><mrow><msup><mi>U</mi><mi>T</mi></msup><mi>U</mi><mo>=</mo><msub><mi>I</mi><mi>d</mi></msub><mo>,</mo></mrow></mtd></mtr></mtable></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001022868550000036.GIF" wi="1310" he="213" /></maths>式中,U为投影矩阵;X为训练集;<img file="FDA0001022868550000037.GIF" wi="310" he="151" />为全局校准矩阵;s.t.为subject to,即满足约束;I<sub>d</sub>为d×d维的单位矩阵。
地址 510640 广东省广州市天河区五山路381号