视频分析在交通控制中的应用

1、承 诺 书我们仔细阅读了第四届文鼎创杯华中地区大学生数学建模邀请赛的选手须知。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们的竞赛编号为： 11784001 我们的选择题号为： B 参赛队员（打印并签名）：队员1： 队员2： 队员3

2、： （以下内容参赛队伍不需要填写） 评阅编号： 武汉工业与应用数学学会第四届文鼎创杯华中地区大学生数学建模邀请赛竞赛组委会 题目： 视频分析在交通控制中的应用 摘要传统的视频监控由人工进行视频监测发现安全隐患或异常状态，或者用于事后分析，这种应用具有其固有的缺点，难以实现实时的安全监控和检测管理。本文通过对提取视频模型的建立，计算车辆的模型来进行计算机监控来解决人工对视频监控的缺点。对问题进行分析，主要分为如下步骤：对于问题一，主要是要建立怎么提取视频的数学模型，我们通过分析对场景中的每个像素点建立高斯模型：，对新来的帧进行高斯模型的拟合来提取背景图像。我们为了提高视频的清晰度和正确性，在高斯

3、模型的基础上还建立了混合高斯模型：，我们通过单高斯背景建模算法和混合高斯模型的k-mean聚类算法以及利用matlab程序对视频图像进行处理分析，得出清晰的图像。对于问题二，根据对视频的观察（1）一段时间车流的密度不大，车辆之间相互影响很小，其它外界干扰因素几乎不存在，即车流是随机的；（2）另一段时间车辆比较拥挤，自由行驶不多。这样我们根据（1）想到了泊松分布，我们根据泊松分布的公式建立泊松模型：；根据（2）我们想到了二项分布，所以我们根据二项分布公式建立二项分布模型：，然后我们根据上述的模型建立模型三（车流量的模型）：从而我们通过题目提供的视频进行计算，得出该视频的车流量是3辆/秒。关键词：

4、高斯模型、混合高斯模型、泊松分布、二项分布、k-mean聚类算法一、问题重述1.1基本情况传统的视频监控由人工进行视频监测发现安全隐患或异常状态，或者用于事后分析，这种应用具有其固有的缺点，难以实现实时的安全监控和检测管理。视频分析是指将先进的信息技术、数据通讯传输技术、电子控制技术、计算机处理技术等应用于工业生产，从而提高生产效率。在交通管理中，带有智能分析功能的监控系统可以通过区分监控对象的外形、动作等特征，做到主动收集、分析数据，并根据预设条件执行报警、记录、分析等动作。对于视频分析而言，其关键技术是建立合理的数学模型。视频附件土豆网：1.2需要解决的问题问题一：根据提供的视频，建立提取

5、附录提供的视频图像的数学模型。问题二：结合已提取的背景，建立计算车流量数学模型。2、 问题分析对于问题一，主要是建立提取视频图像的模型，我们想到了高斯模型，所以我们对场景中的每个像素点建立高斯模型，对新来的帧进行高斯模型的拟合来提取背景图像。从而建立模型：但是我们又想到高斯模型的单一性，对视频图像的提取误差会大一点，那么我们就还建立了混合高斯模型：，然后我们考虑到解法，运用了k-mean聚类算法，从而得出视频图像。对于问题二，我们刚开始想到泊松分布，但是后来再次仔细看了视频，发现有些时间不符合泊松分布而是符合二项分布，那么我们就建立泊松模型和二项分布模型，然后综合两个模型建立车流量模型：通过题

6、目提供的视频进行计算，得出该视频的车流量是3辆/秒。3、 符号说明符号表示均值方差t每个计数间隔持续的时间视频时间i第i帧I一个图像序列高斯模型第i个高斯模型的数学期望相应高斯函数的权重系数第i个高斯模型的协方差矩阵k混合高斯所包括高斯函数的个数权重分布在技数间隔t内到达x辆车的概率在计数间隔t内到达的k辆车的概率单位时间的车辆平均到达率（辆/秒） n正整数Z车流量某一个像素点的H、S、V分量新一帧相应像素点的值g观测数据的分组数m泊松分布的参数4、 模型假设 假设光线等因素对视频图像的背景影响不大。 背景颜色与车的颜色差别不太大。 外界环境的关照变化不大，外界景观和花草数目摆动不大。5、 模

7、型建立5.1模型一我们对场景中的每个像素点建立高斯模型，对新来的帧进行高斯模型的拟合来提取背景图像。对每一个像素利用高斯模板建模，每一像素点都认为服从均值和标准方差的分布，且每一个点的高斯分布是独立的。假设任意时间t，像素点，I表示这一个图像序列，i表示第i帧。我们建立模型为：5.2模型二视频图像序可以表达为三维空间其中t表示时间序列，表示图像位置点。对于一特定点，其时间序为。对任意点的时间序列，可看作K个高斯模型的叠加逼近如图一，对于当前点的概率可表示为：图一其中为高斯模型，表示 第i个高斯模型的数学期望；表示相应高斯函数的权重系数，表示第i个高斯模型的协方差矩阵，k为混合高斯所包括高斯函数

8、的个数，其具体形式为：其中k的取值由实际情况确定，通常k=3，为了提高运算速度，降低计算复杂度，假定协方差矩阵为对角阵，且。进一步确定混合高斯模型背景成分将进行排序生成权重分布。若权重分布的比值系数越大，则其对应的高斯所占混合高斯分布的比例就越大。其中T为阀值，当时，其对应的高斯成分为背景，否则当前含有前景部分。5.3模型三5.3.1泊松分布模型：在技数间隔t内到达x辆车的概率；：单位时间的车辆平均到达率（辆/秒）； t:每个计数间隔持续的时间（秒）。若令为计数间隔t内平均到达的车辆数，则m又称为泊松分布的参数。公式可以写成：g：观测数据的分组数5.3.2二项分布模型 （k=0,1,2，n）：

9、在计数间隔t内到达的k辆车的概率；：单位时间的车辆平均到达率（辆/秒）； t:每个计数间隔持续的时间（秒）。 n：正整数我们计，则 （k=0,1,2，n）5.3.2车流量模型 我们根据上述的两个模型可以得出车流量模型：车流量6、 模型求解6.1模型一和模型二的求解：设表示当前模板上某一个像素点的H、S、V分量。表示新一帧相应像素点的值，ST表示色彩的阀值。我们可以得出背景像素的四种见表，如下表：情况环境条件分类检验判断条件1234 表一我们在对上述算法模型进一步分析在高斯模型求解的基础上我们为了降低混合高斯参数求解的计算方式的复杂性，我们运用了k-means聚类算法实现高斯模型。下图是k-me

10、ans聚类模型的流程图（图二） 图二我们从如下图一个聚类样本的分布随机提取五个数据 图三那么随机的五个数据如下图 图四那么我们把这五个随机数据称作一个集合，我们计算数据集合的平均值，把该平均值作为相应的中心点，进行中心点的两次分配：如下图（图五和图六）图五图六最终我们可以得到聚类的结果如下图：图七图八我们编写程序（附录一：程序）我们可以得到视频图像的处理分析如下（图九）图九从而得出本文的算法能够满足光线等环境变化得出清晰的画面。6.2模型三的求解：6.2.1泊松分布若令为计数间隔t内平均到达的车辆数，则m又称为泊松分布的参数。公式可以写成：g：观测数据的分组数6.2.2二项分布 （k=0,1,

11、2，n）若，则 （k=0,1,2，n）6.2.3车流量根据视频我们可以得出一些数据如下（表二）时间1s5s10s15s20s100s200s300s车辆数4527615910 表二根据上面数据计算得：Z=3辆/秒7、 模型的及改进7.1模型的评价：优点： 建立的模型有成熟的理论基础，又有相应的软件支持，可信度较高。 给出的算法合理，具有实用性。 计算方便简单。缺点： 模型虽然综合考虑到很多因素，但为了建立模型，理想化了很多影响因素，具有一定的局限性，得到的结果可能与实际有一定的出入。 单一算法误差会比较大一些。7.2 模型的改进方向争取在算法上看看能不能找到更精确的算法，来减少局限的误差；或者

12、用多种算法 进行综合。8、 参考文献1、姜启源 数学模型(第三版) 高等教育出版社 2003年第3版2、赵静.但琦 数学建模与数学实验.高等教育出版社 2003年第2版 3、吴清等 MATLAB7简明教程.清华大学出版社 20064、陈超 基于视频的车辆检测技术研究 武汉理工大学硕士学士论文2010年5、王芯 基于视频的运动物体检测算法研究 吉林大学硕士学士论文2009年9、 附录9.1附录一（程序） include <stdio.h> #include <cv.h> #include <cxcore.h> #include <highgui.h>

13、; #include <cvaux.h>int main( int argc, char* argv ) IplImage* pFrame = NULL; IplImage* pFrImg = NULL; IplImage* pBkImg = NULL; CvCapture* pCapture = NULL; int nFrmNum = 0; cvNamedWindow("video", 1); cvNamedWindow("background",1); cvNamedWindow("foreground",1); cv

14、MoveWindow("video", 30, 0); cvMoveWindow("background", 360, 0); cvMoveWindow("foreground", 690, 0); if( argc > 2 ) fprintf(stderr, "Usage: bkgrd video_file_namen"); return -1; if(argc = 2) if( !(pCapture = cvCaptureFromFile(argv1) fprintf(stderr, "Can

15、not open video file %sn", argv1); return -2; if (argc = 1) if( !(pCapture = cvCaptureFromCAM(-1) fprintf(stderr, "Can not open camera.n"); return -2; CvGaussBGModel* bg_model=NULL; while(pFrame = cvQueryFrame( pCapture ) nFrmNum+; if(nFrmNum = 1) pBkImg = cvCreateImage(cvSize(pFrame-&

16、gt;width, pFrame->height), IPL_DEPTH_8U,3); pFrImg = cvCreateImage(cvSize(pFrame->width, pFrame->height), IPL_DEPTH_8U,1); bg_model = (CvGaussBGModel*)cvCreateGaussianBGModel(pFrame, 0); else cvUpdateBGStatModel(pFrame, (CvBGStatModel *)bg_model ); cvCopy(bg_model->foreground,pFrImg,0); cvCopy(bg_model->background,pBkImg,0); pBkImg->origin=1; pFrImg->origin=1; cvShowImage("video", pFrame); cvShowImage("background", pBkImg); cvShowImage("foreground", pFrImg); if( cvWaitKey(2) >= 0 ) break; cvRele