一种拥挤场景下视频异常事件检测方法技术

本发明专利技术涉及一种拥挤场景下视频异常事件检测方法。具体步骤是：训练阶段，首先对视频分块并提取块的光流信息表示块的局部特征，然后以块的局部特征作为原子构建图，采用拉普拉斯特征映射降维并采用自适应聚类方法对局部特征分类，获得类中心作为码字；最后由码字组成码本。测试阶段，先对视频分块提取局部特征，然后计算与码本的特征距离相似度，并结合前一时刻的8-邻域信息检测事件。本发明专利技术只需对光流特征构成的原子集合进行聚类即可建立正常事件模型，历史时刻的邻域信息进一步提高了异常事件检测的准确性。

【技术实现步骤摘要】

本专利技术属于智能视频监控
，涉及一种拥挤场景下视频异常事件检测方 法。
技术介绍
视频异常事件检测是指能自动分析视频监控场景中发生的事件，如果发生异常事 件就能立刻发出报警信号，从而提高相关部门的响应和救援效率。例如步行街上的骑车、轮 滑和车辆通行，广场上人群的恐慌、踩踏等这些异常事件都能及时的被检测到并进行报警。 它在视频监控领域有着广泛的应用前景。 目前现有的异常事件检测方法可以大致分为两类：1)通过对目标对象进行跟踪， 分析其运动轨迹，进而对其做出判断；2)不需要对目标对象进行跟踪，通过建立正常事件 模型，对光流、纹理等特征作出分析。第一类方法，通过对目标对象进行跟踪，获取运动方向 与目标的运动速度以及根据目标对象特征与尺寸比，该类方法针对场景中只存在少数运动 目标时得到了较好的应用，但是在拥挤的环境下，由于目标之间相互重叠，很难对目标对 象进行跟踪，因此这类方法的检测性能欠佳。本专利技术采用第二类方法。
技术实现思路
本专利技术的目的是提供，以提高异常事件检 测率。 为解决上述问题，本专利技术提供的技术方案如下： 步骤⑴特征提取，具体如下： 将视频中每一帧分割成互不重叠且大小为NXN的块，并连续取Μ帧，得到大小 为ΝΧΝΧΜ的立体块，因此Μ帧长的视频由若干个立体块组成，每个立体块称为原子。设 视频帧的分辨率大小为W X Η，每一帧得到块的数量I其中 L·」表示向下取整。t时刻loci位置块的运动信息用一个直方图表示1 < loci < f_block，loci是整数，其中h，1 < i < 4,是loci位置块根据光流方向按90度 间隔量化获得的4个方向上的光流幅值之和。当前时刻为t，结合t时刻的前（M-1) /2帧与 后（M-1)/2帧的直方图信息，t时刻loci位置上的原子表示为.Μ取奇数，1 < loci < f_block。对于一段视频，以Μ帧为单位分成P小段，得到的原子数量 为PXf_block，由这些原子构成该视频的原子集合。 步骤⑵特征学习，具体如下： 2-1.采用拉普拉斯特征映射方法把原子集合映射到低维空间中，再对其进行聚 类。首先对原子集构建图G= (V，E)，顶点集V表示各个原子，带权重的边E表示各个原子 之间的相似度，图中第i个原子与第j个原子之间边的权重按式（1)计算，1彡i彡PXf_ block，1 < j < PXf-block: 公式⑴右边第一项中#的数学表达式为： 余弦距离_其中〈Xi, x.j>表示求Xi与X .j的内积，定义，的形式就是式⑵中i = r， j = G时，％是尺度因子，其中x ^是x ^第G个近邻点，近邻的距离度量采用欧氏距离，r = i或j。 公式⑴右边第二项中，< 7表示第i个原子与第j个原子空间上的欧氏距离；O s 为空间尺度因子。构建完图后，借助于图对原子进行谱聚类。在图论中，聚类的问题转变为 图割的问题。其原则是子图内的边权重最大化和各子图间的边权重最小化。被切断的边权 重之和最小，即使式（3)所示目标函数最小化。 其中&由公式⑴给出；y ;和y」分别是X ;和X」映射到目标空间上的坐标向量， Y由向量yi组成，1 < i < PXf_block。目标函数最小化的问题等价于求解最优的Y : Yopt= argmin (YLY) s. t. Y TDY = 1 (4)公式⑷中，拉普拉斯矩阵 L = D_W。D 是对角阵，其对角线上元素值dii= 由w ^构成。计算L相对于D的广义特征值 和特征向量，选取1个最小非零特征值以及对应的特征向量。将求得的1个特征向量并 成一个（PXf_block) XI的特征向量空间，其中每一行代表原子在1维空间的坐标。最 后根据每个原子的1维空间坐标，采用Lihi Zelnik-Manor的自适应聚类方法（Lihi Zelnik-Manor, Self-Tuning Spectral Clustering. In Proceedings of the 18th Annual Conference on Neural Information Processing Systems 2004)在 1 维空间对原子进行 聚类，自适应确定聚类的类数。聚类后得到NumO个类，每一类的类中心表示一类事件的局 部特征，类中心按式（5)计算： Nk表示属于第k类事件的原子个数。类中心作为码字，位置loci上所有可能的码 字组成码本。 2-2.在loci位置块上建立码本的过程如下： (a)建立初始码本，将训练数据进行步骤2-1的特征学习得到Num类事件，Num = NumO。按式（6)计算wMcK：i，如果wkikK：i>0,则把第k类的类中心作为码字添加到位置loci 上的码本中，同时保存第k类局部特征的训练数据，k= {1，2,...........NumO}; 其中1\1(^表示位置loci上第k类事件发生的次数，1彡loci彡f_block。 (b)输入新的训练数据〇将它与loci位置块上的码本比较，若满足特征距离相 似度4^(0<屈时，th为设定的阈值；将λ:加入到一个新的集合U中，否则把加入 到位置loci上最相似类局部特征的训练数据中，然后重新计算这类局部特征的类中心，更 新码本中码字，其中是按式（7)计算的特征距离的相似度； 其中表示在位置loci上的码本中第k个码字，1彡k彡Num。 (c)当集合U中的数据数量未达到Q时，返回步骤（b);当集合U中的数据数量达 到Q时，重新对集合U进行步骤2-1操作，聚类成Numl类，更新NumO为NumO = Numl，如果 wk,lcicl>〇,贝1J就把第k类事件类中心作为码字添加到位置loci上的码本中，同时保存第k类 局部特征的训练数据，k = {1，2,...........NumO};清空集合U，类目总数Num更新为Num =Num+Numl，判断所有输入的训练数据是否训练完毕，若未完毕则返回步骤（b)。 步骤（3)视频异常事件检测，具体是： 将测试数据与训练阶段建立的loci位置块的码本进行比较，1彡loci彡f_ block，若满足爲_(0 <咖寸，则初步认为丨在l〇ci位置块上有异常事件发生，否则没 有发生； 步骤（4)时空后处理，具体是： 在t时刻，若loci位置块的初始被判有异常事件发生，则考虑loci位置块在t_l 时刻的8-邻域，如果在邻域中至少有两个位置存在异常，则判loci位置块上有异常事件发 生，否则没有异常事件发生。 本专利技术只需对光流特征构成的原子集合进行谱聚类即可建立正常事件模型。时空 后处理考虑了判决位置的历史时刻的邻域信息，进一步提高了异常事件检测的准确性。【附图说明】 图1为异常事件检测模型训练过程图； 图2为异常事件检测模型检测过程图。【具体实施方式】 以下结合附图和实施实例对本专利技术加以详细说明。 本专利技术的异常检测模型训练和检测过程分别如图1、图2所示，具体步骤如下： 步骤⑴特征提取，具体是： 将视频中每一帧分割成互不重叠且大小为NXN的块，并连续取Μ帧，得到大小为 NX NX Μ的立体块，因此Μ帧长的视频本文档来自技高网...

【技术保护点】
一种拥挤场景下视频异常事件检测方法，其特征在于包括如下步骤：步骤(1)特征提取，具体如下：将视频中每一帧分割成互不重叠且大小为N×N的块，并连续取M帧，得到大小为N×N×M的立体块，因此M帧长的视频由若干个立体块组成，每个立体块称为原子；设视频帧的分辨率大小为W×H，每一帧得到块的数量其中表示向下取整；t时刻loci位置块的运动信息用一个直方图表示1≤loci≤f_block，loci是整数，其中hi，1≤i≤4，是loci位置块根据光流方向按90度间隔量化获得的4个方向上的光流幅值之和；当前时刻为t，结合t时刻的前(M‑1)/2帧与后(M‑1)/2帧的直方图信息，t时刻loci位置上的原子表示为xloci=[hlocit-(M-1)/2...hlocit+(M-1)/2],]]>M取奇数，1≤loci≤f_block；对于一段视频，以M帧为单位分成P小段，得到的原子数量为P×f_block，由这些原子构成该视频的原子集合；步骤(2)特征学习，具体如下：2‑1.采用拉普拉斯特征映射方法把原子集合映射到低维空间中，再对其进行聚类；首先对原子集构建图G＝(V,E)，顶点集V表示各个原子，带权重的边E表示各个原子之间的相似度，图中第i个原子与第j个原子之间边的权重按式(1)计算，1≤i≤P×f_block，1≤j≤P×f_block:wij=exp(-(1-dfi,j)σiσj)exp(-dsi.jσs)---(1)]]>公式(1)右边第一项中的数学表达式为：dfi,j=αi,j×df(xi,xj)---(2)]]>余弦距离其中<xi,xj>表示求xi与xj的内积，定义df(xi,xj)=Σk=04M-1min[xi(k),xj(k)]max[Σk=04M-1xi(k),Σk=04M-1xj(k)];σr=(1-dfr,G),]]>σr是尺度因子，其中xG是xr第G个近邻点，近邻的距离度量采用欧氏距离，r＝i或j；公式(1)右边第二项中，表示第i个原子与第j个原子空间上的欧氏距离；σs为空间尺度因子；构建完图后，借助于图对原子进行谱聚类；在图论中，聚类的问题转变为图割的问题；其原则是子图内的边权重最大化和各子图间的边权重最小化；被切断的边权重之和最小，即使式(3)所示目标函数最小化；E(r)=Σijwij||yi-yj||2---(3)]]>其中wij由公式(1)给出；yi和yj分别是xi和xj映射到目标空间上的坐标向量，Y由向量yi组成，1≤i≤P×f_block；目标函数最小化的问题等价于求解最优的Y：Yopt＝argmin(YLY)s.t.YTDY＝1   (4)公式(4)中，拉普拉斯矩阵L＝D‑W；D是对角阵，其对角线上元素值dii＝∑jwij；W由wij构成；计算L相对于D的广义特征值和特征向量，选取l个最小非零特征值以及对应的特征向量；将求得的l个特征向量并成一个(P×f_block)×l的特征向量空间，其中每一行代表原子在l维空间的坐标；最后根据每个原子的l维空间坐标，采用Lihi Zelnik‑Manor的自适应聚类方法在l维空间对原子进行聚类，聚类后得到Num0个类，每一类的类中心表示一类事件的局部特征，类中心按式(5)计算:μk=1NkΣi=1Nkxi---(5)]]>Nk表示属于第k类事件的原子个数；类中心作为码字，位置loci上所有可能的码字组成码本；2‑2.在loci位置块上建立码本的过程如下：(a)建立初始码本，将训练数据进行步骤2‑1的特征学习得到Num类事件，Num＝Num0；按式(6)计算wk,loci，如果wk,loci>0，则把第k类的类中心作为码字添加到位置loci上的码本中，同时保存第k类局部特征的训练数据，k＝{1,2,...........Num0}；wk,loci=nk,lociΣknk,loci,k={1,2,.........··Num0}---(6)]]>其中nk,loci表示位置loci上第k类事件发生的次数，1≤loci≤f_block；(b)输入新的训练数据将它与loci位置块上的码本比较，若满足特征距离相似度时，th为设定的阈值；将加入到一个新的集合U中，否则把加入到位置loci上最相似类局部特征的训练数据中，然后重新计算这类局部特征的类中心，更新码本中码字,其中是按式(7)计算的特征距离的相似度；dmax(xlocinew)=maxk(df(xlocinew,clocik))---(7)]]>其中表示在位置...

【技术特征摘要】

【专利技术属性】
技术研发人员：陈华华，周灵娟，郭春生，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33