[发明专利]一种基于稀疏低秩的人体行为识别方法

摘要

本发明属于数字图像处理技术领域，涉及计算机视觉、模式识别、机器学习和数据挖掘等相关理论知识。本发明首先利用光流直方图对相邻两帧图片提取光流特征，并对单帧图像提取梯度直方图信息，获得监测场景内的运动特征信息，采用低维空间的特征信息并按照[动作1|动作2|动作3|……]的方式排列；然后用K均值的方式聚类，得到聚类中心后，将聚类中心作为过完备字典，求解测试样本在过完备字典下的稀疏低秩表达，得到表达矩阵；最后根据表达矩阵中的最大值求解出该测试样本所属的行为类别。本发明采用基于低秩稀疏的人体动作识别，采用交叉验证的方法，识别率92.3～98.79％，误识率1.21～7.6％。本发明具有具有低秩特点，且识别率上达到92.3～98.79％，误识率1.21～7.6％。

主权项

一种基于稀疏低秩的人体行为识别方法，具体步骤如下：步骤1、视频图像灰度转换并提取场景特征步骤1‑1：首先接收监测设备采集到的视频流，对得到的单帧图像进行灰度化处理。步骤1‑2：将每帧图片归一化成96*48的大小；每个cell的大小为8*8的像素点，每个block中包含2*2个cell即16*16的像素点；相邻block有一半的区域是重叠的；梯度方向分为9个方向，得到的一张图片的梯度直方图HOG特征向量为1980维；对连续两帧灰度图像在3*3的网格中提取光流直方图特征HOF，即每3*3的网格求出一个光流特征向量，根据视频帧的不同大小得到N个光流特征向量，将光流范围0～2π等分成32个子空间，构建32维的直方图B，在统计直方图时，根据角度值对应位置投票，投票所得结果作为权值，因此得权值向量h＝[h_k],k＝{1,2,3…,K},其中k为每个网格的索引值，K是网格总数，h_k为第k个网格中32维直方图B对应的权值；步骤2、特征聚类设光流特征向量x＝[x_i],i＝{1,2,3…,m},对于每一个特征向量，动作类型计算公式：$C^{\left(j\right)}=\underset{j}{\mathrm{argmin}}\left|\right|x^{\left(i\right)}-{\mathrm{\μ}}_j|{|}^2,j=\{1,2,3...,K\}$对于每类动作C^(j)j,计算该类质心这里K表示聚类中心，取值500～2000；步骤3、低秩稀疏表达步骤3‑1：将聚类中心构成的字典D＝[d_k],k＝{1,2,3…,n},则X为一个d_s×n的矩阵，d_s为特征维度，识别样本集X＝[x_k],k＝{1,2,3…,n}，X中的任意一列x_i都可以用D进行线性表达，将这种关系式写为：X＝DZ+E,其中Z的每一列表示X的对应列用D进行线性表达的系数，矩阵Z称之为表达子，E为由噪声引起的误差；步骤3‑2：据此建立目标方程：argmin_Z,E||Z||_*+β||Z||₁+γ||E||₁ s.t.X＝DZ+E式1其中β和γ分别表示各项的关注因子，它们的值越大，表示相关项越受关注，这里β＝2,γ＝2，公式中分别用核范数和L1范数求解低秩和稀疏问题；步骤3‑3：求解步骤3‑2中的式1，这里用到了增广拉格朗日乘数法，得到目标函数后分别对各个参数进行迭代优化，则参数的迭代过程为：$Z_{k+1}=J_{{\left({\mathrm{\η\μ}}_k\right)}^{-1}}(Z_k+\frac{Q(D,X,Z_k,Y_{1,k},Y_{2,k},{\mathrm{\μ}}_k)}{\mathrm{\η}});$其中$W_{k+1}={\mathrm{\Θ}}_{{\mathrm{\β\μ}}_k^{-1}}(Z_{k+1}+\frac{Y_{2,k}}{{\mathrm{\μ}}_k}),$$E_{k+1}={\mathrm{\Θ}}_{{\mathrm{\γ\μ}}_k^{-1}}(X-{\mathrm{DZ}}_{k+1}+\frac{Y_{1,k}}{{\mathrm{\μ}}_k}),$Y_1,k+1＝Y_1,k+μ_k(X‑DZ_k+1‑E_k+1),Y_2,k+1＝Y_2,k+μ_k(Z_k+1‑W_K+1),μ_k+1＝ρμ_k,其中函数J和Θ分别表示奇异值压缩算子和压缩算子，ρ＝1.1,各参量的初始值为Z₀＝W₀＝E₀＝Y_1,0＝Y_2,0＝0；步骤4、行为识别步骤4‑1、通过步骤3求得的稀疏低秩矩阵Z，其中Z的每一列就是待识别样本在D中的线性表达，现在只需要找出表达最大的那个系数即可知道待识别样本所属的类别，采用最大池化法，即找出矩阵Z中每一行的最大值；令y_i＝max(|z_i1|,|z_i2|,…|z|_in),i＝1,2,…l；步骤4‑2、得到矩阵Z中每一行的最大值y_i，取最大值中的前十个，分别找出这十个值所对应的相应类别，并计算该类别在矩阵Z中对应的系数，并把相同类别的系数叠加起来，得到的最大值所对应的类别即为最终的识别结果，相应的计算公式如下：y_i＝max(|z_i1|,|z_i2|，...|zin|),i＝1,2，...T上式表示取Z中行向量中最大值的前T个；分别找出对应的P_j＝C_i(y_i)，i＝1，2，…T，其中C_i(y_i)为步骤2的聚类结果；其中δ(·)是Dirac Delta函数，PreC即为最终的识别结果，T取5～10。