本发明涉及图像处理技术领域，特别是一种监控视频中基于运动方向的异常行为检测方法。
背景技术：
近年来，随着各种安全突发事件的增多，大众安全意识的提升，人工智能理念的渗透，智能监控越来越受到人们的关注。传统的监控系统主要通过人工监控的方式实现对公共场合的安全管理，缺乏实时性和主动性。在很多情况下，视频监控由于无人管理只是起到了视频备份的作用没有做到监管的职责。此外，随着监控摄像头的普及和广泛布设，传统的人工监控方式已经不能满足现代监控的需求。
针对这一问题，大众都致力于寻求解决方案来代替人工操作。目前，随着视频监控技术以及信息科学的不断发展，视频监控、人机交互、视频搜索等领域有了长足的发展，自动化监控逐渐成为一个具有广泛应用前景的研究课题。异常行为检测是自动监控的重要内容，相较于一般地人体行为识别集中在人的常规动作的识别上，异常行为通常具有高突发性，并且持续时间较短，很难获取行为特征的特点。
近年来，对于异常行为的检测，研究者们提出了不同的方法，主要分为基于人体运动特征的检测方法和基于前景目标运动轨迹的检测方法。基于人体运动特征的检测方法通过提取前景目标的运动特征达到高准确度，但是目标的运动特征不易提取且计算复杂；基于前景目标运动轨迹的检测方法计算量小，但是精度低，一般用于特定的场景。
总的来说，对于异常行为的检测，现在并特别没有高效的方法，还需要不断的探索。
技术实现要素：
本发明所要解决的技术问题是克服现有技术的不足而提供一种监控视频中基于运动方向的异常行为检测方法，以有效的提高对视频中异常行为的检测效果。
本发明为解决上述技术问题采用以下技术方案：
根据本发明提出的一种监控视频中基于运动方向的异常行为检测方法，包括以下步骤：
步骤1、输入3个原始视频帧和1个背景视频帧，所述原始视频帧是用户输入的监控视频的视频帧，所述背景视频帧是用户输入的监控视频的背景的视频帧，每个原始视频帧和背景视频帧均由k个像素点构成；
步骤2、比较背景视频帧与原始视频帧，具体步骤如下：
步骤21、取原始视频帧的第一帧、第二帧和第三帧中每个像素点的灰度值，原始视频帧的第一帧在位置(xi,yi)处的像素点的灰度值为f(xi,yi,t-1)，原始视频帧的第二帧在位置(xi,yi)处的像素点的灰度值为f(xi,yi,t)，原始视频帧的第三帧在位置(xi,yi)处的像素点的灰度值为f(xi,yi,t+1)；再取背景视频帧像素点的灰度值，背景视频帧在位置(xi,yi)处的像素点的运动矢量为f0(xi,yi,t)，i依次取1到k中的整数值；然后分别计算背景视频帧与原始视频帧第一帧、第二帧和第三帧在位置(xi,yi)处的像素点的运动矢量的差值d1(xi,yi,t)、d2(xi,yi,t)和d3(xi,yi,t)，采用以下计算公式得到三组差值：
d1(xi,yi,t)＝|f(xi,yi,t-1)-f0(xi,yi,t)|，
d2(xi,yi,t)＝|f(xi,yi,t)-f0(xi,yi,t)，
d3(xi,yi,t)＝|f(xi,yi,t+1)-f0(xi,yi,t)；
接着计算上述三组差值的平均值d(xi,yi,t)：
d(xi,yi,t)＝(d1(xi,yi,t)+d2(xi,yi,t)+d3(xi,yi,t))/3；
步骤22、设置背景区域边界阈值为m1，当d(xi,yi,t)≥m1，设置b(xi,yi,t)＝1，当d(xi,yi,t)＜m1，设置b(xi,yi,t)＝0，设定原始视频帧的第一帧中满足b(xi,yi,t)＝1的所有像素点组成的区域为初步视频目标，b(xi,yi,t)为在位置(xi,yi)处的像素点的二值化表示；
步骤3、输入初步视频目标的所有像素点，根据光流法得到光流约束方程，对光流约束方程做泰勒展开，得到
计算初步视频目标的每个像素点的速度t(xi,yi)：
步骤4、设置噪声区域边界阈值为m2，当初步视频目标中存在位置(xi,yi)处的像素点的速度t(xi,yi)≤m2，设定满足t(xi,yi)≤m2的所有像素点组成的区域为视频目标，转入步骤5，当初步视频目标的所有像素点的速度均满足t(xi,yi)＞m2，设定原始视频帧中不存在视频目标，转入步骤1；
步骤5、根据步骤3所述的光流约束方程
转换u(xi,yi)和v(xi,yi)为极坐标表示形式，转换公式为
计算视频目标的每个像素点的运动方向r(xi,yi)：
步骤6、查询视频目标的所有像素点，当视频目标中位置为(xi,yi)处的像素点存在二值化表示b(xi,yi,t)＝0的毗邻像素点，设定该像素点为视频目标的轮廓，设置该像素点的权值a(xi,yi,t)＝2，当视频目标中位置为(xi,yi)处的像素点的毗邻像素点的二值化表示均为b(xi,yi,t)＝1，设置该像素点的权值a(xi,yi,t)＝1，所述的毗邻像素点是位置在(xi-1,yi-1)、(xi-1,yi)、(xi-1,yi+1)、(xi,yi-1)，(xi,yi+1)、(xi+1,yi-1)、(xi+1,yi)、(xi+1,yi+1)处的像素点；
步骤7、统计视频目标的运动方向，具体步骤如下：
步骤71、将方向范围360°均分为c个区间进行运动方向统计，每个区间大小为
步骤72、采用方向直方图模型h统计视频目标的像素点的运动方向，所述的h＝{hj},j＝1,2,…,c，hj表示第j区间的方向直方图统计；接着对每一区间进行归一化处理，定义归一化后的方向直方图模型为hnorm，hnorm＝{hnorm(j)}，
步骤73、计算归一化方向直方图模型的运动方向熵e：
步骤8、设置异常行为边界阈值为m3，当e≥m3，输出该视频目标的行为为异常行为，当e＜m3，设定该视频目标行为正常，转入步骤1。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤1中，输入的监控视频是在摄像机没有抖动的情况下获得。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤1中，k取720×560。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤1中，相邻两个原始视频帧的时间间隔为0.5秒。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤22中，m1取10，表示背景区域边界阈值。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，步骤3中所述的光流法指设经过时间dt，在位置(xi,yi)处的像素点移动至位置(xi+dx,yi+dy)处，由此存在光流约束方程f(xi,yi,t)＝f(xi+dx,yi+dy,t+dt)的光流处理办法，dx表示dt时间内位置(xi,yi)处的像素点在水平方向移动的距离，dy表示dt时间内位置(xi,yi)处的像素点在垂直方向移动的距离。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤4中，m2取0.5，表示噪声区域边界阈值。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤71中，c取9。
作为本发明所述的一种监控视频中基于运动方向的异常行为检测方法进一步优化方案，所述步骤8中，m3的取值在0.5～0.6之间，表示异常行为边界阈值。
本发明采用以上技术方案与现有技术相比，具有以下技术效果：
本发明使用背景相减法处理原始视频帧并获取初步视频目标，并对获取的视频目标进行优化，通过对视频目标运动方向的统计分析达到检测异常行为的效果。能够检测出大部分异常情况，具有良好的实时性和有效性，具体来说：
(1)监控视频中存在大量静止画面，不需要进行行为检测，通过背景减法快速筛选视频帧，只对存在灰度值变化的视频帧进行检测，提高算法的检测速度；
(2)监控视频中受光照和天气变化的影响，只对运动的视频目标进行检测，减少光照和天气对检测效果的影响；
(3)视频前景目标根据运动的剧烈程度相异，对于视频帧像素点速度大于阈值的区域不予考虑，大大的减少了运算量，提供了算法的时效性；
(4)使用方向直方图对运动方向进行划分，可以使行人细微的肢体动作被忽略而不影响检测效果，大大降低了计算的繁琐度。
附图说明
图1是监控视频中基于运动方向的异常行为检测方法流程。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明：
在具体实施中，图1是监控视频中基于运动方向的异常行为检测方法流程。首先用户输入3个原始视频帧和1个背景视频帧，这里的原始视频帧取自监控视频的视频帧，时间间隔为0.5秒，背景视频帧是监控视频的背景的视频帧，每个视频帧由k个像素点构成。
背景相减法通过比较背景图像和视频图像的差别可以快速并且较精确的获取覆盖视频背景的区域，所以选取原始视频帧的第一帧、第二帧和第三帧中每个像素点的运动矢量。原始视频帧第一帧在位置(xi,yi)处的像素点的运动矢量为f(xi,yi,t-1)，原始视频帧第二帧在位置(xi,yi)处的像素点的运动矢量为f(xi,yi,t)，原始视频帧第三帧在位置(xi,yi)处的像素点的运动矢量为f(xi,yi,t+1)；再取背景视频帧像素点的运动矢量，背景视频帧在位置(xi,yi)处的像素点的运动矢量为f0(xi,yi,t)，所述i依次取1到k中的值，这里的f(xi,yi,t)和f0(xi,yi,t)均表示时间t时位置(xi,yi)处的像素点的灰度值，f(xi,yi,t-1)表示时间t-1时在位置(xi,yi)处的像素点的灰度值，f(xi,yi,t+1)表示时间t+1时位置(xi,yi)处像素点的灰度值。再计算背景视频帧与上述三个视频帧在位置(xi,yi)处的像素点的运动矢量的差值d1(xi,yi,t)、d2(xi,yi,t)和d3(xi,yi,t)，采用以下计算公式得到三组差值：
d1(xi,yi,t)＝|f(xi,yi,t-1)-f0(xi,yi,t)|，
d2(xi,yi,t)＝|f(xi,yi,t)-f0(xi,yi,t)，
d3(xi,yi,t)＝|f(xi,yi,t+1)-f0(xi,yi,t)；
通过多帧对比可以消除单帧对比带来的误差，所以计算背景视频帧与上述三个视频帧在位置(xi,yi)处的像素点的运动矢量的差值d(xi,yi,t)，采用以下公式计算d(xi,yi,t)：d(xi,yi,t)＝(d1(xi,yi,t)+d2(xi,yi,t)+d3(xi,yi,t))/3。
为了对消除视频背景的影响，所有像素点设置二值化表示b(xi,yi,t)标识位置(xi,yi)处的像素点是否属于背景，那么设置背景区域边界阈值为m1，m1按照经验值取0.5。如果d(xi,yi,t)＜m1，设置b(xi,yi,t)＝0，设定满足b(x,y,t)＝0的所有像素点组成的区域为背景区域；如果d(xi,yi,t)≥m1，设置b(xi,yi,t)＝1，设定满足b(x,y,t)＝1的所有像素点组成的区域为初步视频目标，只对初步视频目标作进一步处理。
假设摄像头不抖动并且视频帧的背景不变，初步视频目标与背景不存在相对运动，那么初步视频目标运动产生速度。对于上一步获取的初步视频目标，可以计算像素点的速度判断初步视频目标是运动物体还是由于光照、天气变化影响误判的背景区域。首先根据光流法得到光流约束方程f(xi,yi,t)＝f(xi+dx,yi+dy,t+dt)，对光流约束方程做泰勒展开，得到
为了有效的消除噪声对检测效果影响，那么设置噪声区域边界阈值为m2，如果初步视频目标的所有像素点的速度均满足t(xi,yi)＞m2，设定原始视频帧中不存在视频目标，回到开始重新进行检测。如果初步视频目标中存在位置(xi,yi)处的像素点速度t(xi,yi)≤m2，设定满足t(xi,yi)≤m2的所有像素点组成的区域为视频目标，继续下面的操作。
对于视频目标区域的像素点，求解光流约束方程
快速查询视频目标的所有像素点，如果视频目标中位置为(xi,yi)处的像素点存在二值化表示的毗邻像素点，说明该像素点为视频目标的轮廓。异常行为的发生通常在人体轮廓更加剧烈，比如投掷、打斗等，所以设置该像素点的权值a(xi,yi,t)＝2。如果视频目标中位置为(xi,yi)处的像素点的毗邻像素点的二值化表示均为b(xi,yi,t)＝1，设置该像素点的权值a(xi,yi,t)＝1。这里的毗邻像素点是位置(xi-1,yi-1)、(xi-1,yi)、(xi-1,yi+1)、(xi,yi-1)，(xi,yi+1)、(xi+1,yi-1)、(xi+1,yi)、(xi+1,yi+1)处的像素点。
本发明主要通过视频目标的运动方向混乱程度判断异常行为，运动方向越混乱说明越可能发生异常行为，那么需要对视频目标的运动方向进行统计分析。将方向范围360°均分为9个区间，每个区间大小为40°，所得的视频目标的每个像素点的运动方向按照权值a(xi,yi,t)的大小分别归类到相应区间，如果a(xi,yi,t)＝1，那么运动方向归类到相应区间1次，如果a(xi,yi,t)＝2，运动方向归类到相应区间2次。
使用方向直方图模型可以对方向进行简明有效的分析，定义方向直方图模型为h，h＝{hj},j＝1,2,…,9，hj表示第j区间的方向直方图统计，对每一区间进行归一化处理，定义归一化后的方向直方图模型为hnorm，hnorm＝{hnorm(j)}，
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替代，都应当视为属于本发明的保护范围。