一种拥挤场景下视频异常事件检测方法
【技术领域】
[0001] 本发明属于智能视频监控技术领域，涉及一种拥挤场景下视频异常事件检测方 法。
【背景技术】
[0002] 视频异常事件检测是指能自动分析视频监控场景中发生的事件，如果发生异常事 件就能立刻发出报警信号，从而提高相关部门的响应和救援效率。例如步行街上的骑车、轮 滑和车辆通行，广场上人群的恐慌、踩踏等这些异常事件都能及时的被检测到并进行报警。 它在视频监控领域有着广泛的应用前景。
[0003] 目前现有的异常事件检测方法可以大致分为两类：1)通过对目标对象进行跟踪， 分析其运动轨迹，进而对其做出判断；2)不需要对目标对象进行跟踪，通过建立正常事件 模型，对光流、纹理等特征作出分析。第一类方法，通过对目标对象进行跟踪，获取运动方向 与目标的运动速度以及根据目标对象特征与尺寸比，该类方法针对场景中只存在少数运动 目标时得到了较好的应用，但是在拥挤的环境下，由于目标之间相互重叠，很难对目标对 象进行跟踪，因此这类方法的检测性能欠佳。本发明采用第二类方法。
【发明内容】
[0004] 本发明的目的是提供一种拥挤场景下视频异常事件检测方法，以提高异常事件检 测率。
[0005] 为解决上述问题，本发明提供的技术方案如下：
[0006] 步骤⑴特征提取，具体如下：
[0007] 将视频中每一帧分割成互不重叠且大小为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，由这些原子构成该视频的原子集合。
[0008] 步骤⑵特征学习，具体如下：
[0009] 2-1.采用拉普拉斯特征映射方法把原子集合映射到低维空间中，再对其进行聚 类。首先对原子集构建图G= (V，E)，顶点集V表示各个原子，带权重的边E表示各个原子 之间的相似度，图中第i个原子与第j个原子之间边的权重按式（1)计算，1彡i彡PXf_ block，1 < j < PXf-block:
[0011] 公式⑴右边第一项中#的数学表达式为：
[0013] 余弦距离
_其中〈Xi, x.j>表示求Xi与X .j的内积，定义
，的形式就是式⑵中i = r， j = G时，％是尺度因子，其中x ^是x ^第G个近邻点，近邻的距离度量采用欧氏距离，r = i或j。
[0014] 公式⑴右边第二项中，< 7表示第i个原子与第j个原子空间上的欧氏距离；O s 为空间尺度因子。构建完图后，借助于图对原子进行谱聚类。在图论中，聚类的问题转变为 图割的问题。其原则是子图内的边权重最大化和各子图间的边权重最小化。被切断的边权 重之和最小，即使式（3)所示目标函数最小化。
[0016] 其中&由公式⑴给出；y ;和y」分别是X ;和X」映射到目标空间上的坐标向量， Y由向量yi组成，1 < i < PXf_block。目标函数最小化的问题等价于求解最优的Y :
[0017] 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)计算：
[0019] Nk表示属于第k类事件的原子个数。类中心作为码字，位置loci上所有可能的码 字组成码本。
[0020] 2-2.在loci位置块上建立码本的过程如下：
[0021] (a)建立初始码本，将训练数据进行步骤2-1的特征学习得到Num类事件，Num = NumO。按式（6)计算wMcK：i，如果wkikK：i>0,则把第k类的类中心作为码字添加到位置loci 上的码本中，同时保存第k类局部特征的训练数据，k= {1，2,...........NumO};
[0023] 其中1\1(^表示位置loci上第k类事件发生的次数，1彡loci彡f_block。
[0024] (b)输入新的训练数据〇将它与loci位置块上的码本比较，若满足特征距离相 似度4^(0<屈时，th为设定的阈值；将λ:加入到一个新的集合U中，否则把加入 到位置loci上最相似类局部特征的训练数据中，然后重新计算这类局部特征的类中心，更 新码本中码字，其中是按式（7)计算的特征距离的相似度；
[0026] 其中表示在位置loci上的码本中第k个码字，1彡k彡Num。
[0027] (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)。
[0028] 步骤（3)视频异常事件检测，具体是：
[0029] 将测试数据与训练阶段建立的loci位置块的码本进行比较，1彡loci彡f_ block，若满足爲_(0 <咖寸，则初步认为丨在l〇ci位置块上有异常事件发生，否则没 有发生；
[0030] 步骤（4)时空后处理，具体是：
[0031] 在t时刻，若loci位置块的初始被判有异常事件发生，则考虑loci位置块在t_l 时刻的8-邻域，如果在邻域中至少有两个位置存在异常，则判loci位置块上有异常事件发 生，否则没有异常事件发生。
[0032] 本发明只需对光流特征构成的原子集合进行谱聚类即可建立正常事件模型。时空 后处理考虑了判决位置的历史时刻的邻域信息，进一步提高了异常事件检测的准确性。
【附图说明】
[0033] 图1为异常事件检测模型训练过程图；
[0034] 图2为异常事件检测模型检测过程图。
【具体实施方式】
[0035] 以下结合附图和实施实例对本发明加以详细说明。
[0036] 本发明的异常检测模型训练和检测过程分别如图1、图2所示，具体步骤如下：
[0037] 步骤⑴特征提取，具体是：
[0038] 将视频中每一帧分割成互不重叠且大小为NXN的块，并连续取Μ帧，得到大小为 NX NX Μ的立体块，因此Μ帧长的视频可以由若干个立体块组成，每个立体块称为原子，取Ν =10，Μ = 3。设视频帧的分辨率大小为320X240,每一帧得到块的数量f_block = 768。t 时刻loci位置块的运动信息用一个直方图表示（，1彡loci彡f_block， loci是整数，其中h，1 < i < 4,是loci位置块根据光流方向按90度间隔量化获得的4个 方向上的光流幅值之和。当前时刻为t，结合t时刻的前1帧与后1帧的直方图信息，t时 刻loci位置上的原子可以表示为，
3取？= 10,得到的原子数量为 PXf_block = 7680,由这些原子构成该视频的原子集合。
[0039] 步骤⑵特征学习，具体如下：
[0040] 2-1.采用拉普拉斯特征映射方法把原子集合映射到低维空间中，再对其进行聚 类。首先对原子集构建图G= (V，E)，顶点集V表示各个原子，带权重的边E表示各个原子 之间的相似度，图中第i个原子与第j个原子之间边的权重按式（1)计算，1彡i彡PXf_ block，1 < j < PX f_