基于计算机视觉的汽车流量检测统计
- 格式:pdf
- 大小:321.33 KB
- 文档页数:4
第34卷 第5期2006年 5月 华 中 科 技 大 学 学 报(自然科学版)J.HuazhongUniv.ofSci.&Tech.(NatureScienceEdition)Vol.34No.5 May 2006
收稿日期:2005205216.作者简介:陈振学(19772),男,博士研究生;武汉,华中科技大学图像识别与人工智能研究所(430074).E2mail:chenzhenxue@163.com基金项目:总装备部基金资助项目(51401020201JW0521).基于计算机视觉的汽车流量检测统计陈振学1 汪国有1 刘成云2(1华中科技大学图像识别与人工智能研究所,湖北武汉430074;2武汉科技大学理学院,湖北武汉430081)摘要:利用计算机视觉和图像处理的方法对交通路口通行车辆进行流量统计.首先,由摄像机视场的标定原理和标定公式,得到了摄像机成像平面内的像素距离与实际三维空间距离的对应关系;然后,采用基于光流场的车辆检测方法,对所有经过设置的虚拟线圈的车辆进行实时检测和信息提取;最后,对提取的车辆信息进行统计和计算,得到实时路况信息.实验结果表明,该算法具有较高的检测概率,能够满足系统的精度和实时处理要求.关 键 词:交通流量;计算机视觉;光流场;流量统计中图分类号:TP391 文献标识码:A 文章编号:167124512(2006)0520046204StatisticsofvehicleflowsbasedoncomputervisionChenZhenxue WangGuoyou LiuChenyunAbstract:Onthebasisofcomputervisionandimageprocessing,amethodwasproposedtoimprovethestatisticalresultsofvehicleflowdetectionattrafficcrossing.Accordingtocamerascenecalibra2tiontheoryandexpressions,thecorrespondingrelationbetweenthepixelsdistanceonthecameraima2gingplaneandthefactualdistanceof3Dspacewasobtained.Thevehicledetectionmethodbasedonopticalflowwasutilizedtoextracttheallvehiclespassingthesetdummyloop.Thereal2timeroutein2formationwouldbeobtainedwithcomputingtheextractedvehicles.Theexperimentalresultsshowedthatthismethodhashigherdetectionprobabilityandsatisfiedwithsystemaccuracyandrobustnessre2quire.Keywords:trafficflow;computervision;opticalflow;flowstatisticsChenZhenxue DoctoralCandidate;InstituteofPatternRecognition&AI,HuazhongUniv.ofSci.&Tech.,Wuhan430074,China. 交通流量信息是交通控制中的重要信息.要实现交通流数据的提取,先要进行车辆检测.目前国内常见的交通流检测的方法有超声波检测、红外检测、环形感应圈检测和计算机视觉检测[1~3],超声波检测精度不高,容易受车辆遮挡和行人的影响,检测的距离短(一般不超过12m);红外检测受车辆本身热源的影响,抗噪声的能力不强,检测精度也不高;环形感应圈检测精度高,但要求设置于路面土木结构中,对路面有损坏,施工和安装不便,而且安装的数量多.现在,有学者将视频检测应用到交通流量的检测和交通控制上,获得了很大的成功.它是一种基于视频图像的检测技术,也称为计算机视觉或机器视觉检测,是一种结合数字视频图像和人工模式识别的技术.本文借鉴国内外相关技术研究成果,提出了一套视频检测系统的设计方案.1 摄像机视场距离的标定算法摄像机视场距离的标定涉及成像几何的理论,将客观世界的三维场景投射到摄像机的二维平面,需要建立摄像机的模型来描述.通过确定摄像机的物理参数和方位参数,找到两者的对应关系,才能把像平面的度量确定下来.成像变换涉及不同坐标系之间的变换.在摄像机的成像系统中,涉及到的坐标系有4个(如图1所示):
图1 成像系统坐标系间的关系a.世界坐标系,也称真实世界或绝对坐标系XYZ;b.摄像机坐标系,以摄像机为中心制定的坐标系xyz;c.像平面坐标系,在摄像机内所形成的像平面坐标系x′y′;d.计算机图像坐标系,计算机内部数字图像所用的坐标系MN,以像素为单位.根据以上几个坐标系不同的转换关系,就可以得到本系统所需要的摄像机成像模型,换算出二维图像到三维场景的对应关系.摄像机成像变换的过程由三部分组成,分别对应坐标系的相互转换(见图2).图2 成像变换的过程从世界坐标(X,Y,Z)到摄像机坐标(x,y,z)的变换中,若仅考虑刚体目标的情况,则可表示为 [x,y,z,1]T=R[T・[X,Y,Z,1]T](1)式中:T是平移矩阵;R是旋转矩阵,表示像平面在世界坐标系里的旋转变换.旋转可以借助刚体转动的3个欧拉角来表示(如图3),利用欧拉角可将旋转矩阵R表示成θ,φ和
图3 欧拉角示意图R=r1r2r30r4r5r60r7r8r900001,式中:r1=cos
Ty)+r3(Z+Tz)r4(X+Tx)+r5(Y+Ty)+r6(Z+Tz)r7(X+Tx)+r8(Y+Ty)+r9(Z+Tz)1. (2)根据透镜成像原理,从摄像机的三维空间坐标(x,y,z)到像平面坐标(x′,y′)的变换为x′=fx/z; y′=fy/z (f>0).(3)计算机中使用的图像坐标单位是存储器中离散像素的个数,所以对实际像平面的坐标还需取整数转换才能映射到计算机的成像平面,从像平面坐标(x′,y′)到计算机图像坐标(M,N)的变换为M=Om-x′/Sx; N=On-y′/Sy,(4)式中:Om和On分别为像平面的原点在计算机图像平面上所映射的点像素所在的行数和列数;Sx和Sy分别为在x和y方向上的尺度因子.将式(2)~(4)映射方程结合起来,就可推导出计算机图像坐标(M,N)与实际的三维空间坐标(X,Y,Z)的变换方程 M=Om-(f/Sx)・r1(X+Tx)+r2(Y+Ty)+r3(Z+Tz)r7(X+Tx)+r8(Y+Ty)+r9(Z+Tz), N=On-(f/Sy)・r4(X+Tx)+r5(Y+Ty)+r6(Z+Tz)r7(X+Tx)+r8(Y+Ty)+r9(Z+Tz).(5) 根据建立的摄像机模型(5)可知,最终要确定的参数有R,T,f,Sx,Sy,Om,On,其中旋转矩阵R有9个未知元素,但实际上只有3个自由度(θ,φ,<);平移矩阵T有3个未知元素,记为(Tx,Ty,・74・第5期 陈振学等:基于计算机视觉的汽车流量检测统计 h),h是摄像机的高度.R和T是摄像机的方位参数;f,Sx,Sy,Om,On是摄像机的物理参数.由此可知需确定的参数为11个,把它们记成向量的形式如下P=[θ,φ,<,h,f,Sx,Sy,Om,On,Tx,Ty]T. 在实际标定时,取6个样本点的坐标作为式(5)的输入,就可以求解摄像机的标定参数[4].注意到需要求解的方程是非线性的,直接求解将很难得出结果,因此必须采用迭代的数值计算方法,最终得出近似值.利用迭代法,只要选择合适的样本点和解的初值,在给定的误差范围内,就可以得到满足精度要求的参数值.2 通行车辆的检测与识别汽车流量的统计首先要对摄像机视场内的通行车辆进行检测和识别.由于图像光流的计算不需要在图像序列中建立特征之间的对应关系,属于较高层次的计算机视觉表达,比背景差影法、边缘检测法具有更好的鲁棒性.因此,本文采用光流场的检测算法,利用连续相邻两帧的车辆的运动信息对车辆进行检测和分割.光流的计算基于如下两个假设:a.任何物体点所观察到的亮度随时间是恒定不变的;b.图像平面内的邻近点以类似的方式进行移动.虽然这两个假设在实际应用中并不是总成立,但对于估算光流场,这些理想化的假设是必要的.由假设a可知运动实体的亮度沿运动路径保持恒定.于是设时变图像可表示为I(x,y,t),则根据该约束条件可以认为I(x,y,t)=I(x+Δx,y+Δy,t+Δt).(6)把式(6)等号的右边用泰勒级数展开可得 I(x+Δx,y+Δy,t+Δt)=I(x,y,t)+Ixdx+Iydy+Itdt+ε,式中:Ix,Iy,It是亮度函数I(x,y,t)关于x,y,t的偏导数;ε为高阶项.简单的忽略高阶项则有:Ixdx+Iydy+Itdt=0;Ixdx/dt+Iydy/dt+It=0,式中(dx/dt,dy/dt)=(u,v)为(x,y)点的光流.所以就得到了光流方程Ixu+Iyv+It=0.光流方程说明了这样一个现象:一个运动物体经过一固定摄像机,则图像上某一点灰度的时间变化率等于它的空间变化率乘上物体经过摄像机的运动速度.方程中包含两个未知量,若要求解光流(u,v),必须引入其他的约束条件.常用的约束光流方法是假设光流在这个图像上的变化平滑,即约束平滑条件:令光流速度的梯度幅度的平方值u2x+u2y和v2x+v2y最小[5].于是,可将光流场u,v的计算归结为如下的变分问题: (Ixu+Iyv+It)2+α2(u2x+u2y+v2x+v2y)=δmin. 利用迭代算法,可以得到(u,v)的递归解: u(k+1)=u(k)-Ix(Ixu(k)+Iyv(k)+It)/(α2+I2x+I2y); v(k+1)=v(k)-Iy(Ixu(k)+Iyv(k)+It)/(α2+I2x+I2y),式中u(k),v(k)是第k次迭代的速度估计(u(k),v(k))的邻域平均值.光流场检测结果如图4所示.
图4 车辆检测分析通过阈值化的分割方法,可以将感兴趣的车辆分割提取出来,并根据分割得到的区域的面积S、区域的长度H和宽度W进行识别,此区域满足车辆信息的条件是:Smin