CN111445699B - 基于实时车辆轨迹的交叉口交通冲突判别方法 - Google Patents
基于实时车辆轨迹的交叉口交通冲突判别方法 Download PDFInfo
- Publication number
- CN111445699B CN111445699B CN202010287744.1A CN202010287744A CN111445699B CN 111445699 B CN111445699 B CN 111445699B CN 202010287744 A CN202010287744 A CN 202010287744A CN 111445699 B CN111445699 B CN 111445699B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- conflict
- traffic
- time
- image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30236—Traffic on road, railway or crossing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30241—Trajectory
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30248—Vehicle exterior or interior
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Traffic Control Systems (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于实时车辆轨迹的交叉口交通冲突判别方法,包括如下步骤:步骤一、采集实时视频图像数据,对其进行图像灰度化、中值滤波和二值化的预处理;步骤二、根据混合高斯背景模型建立视频背景模型,并通过形态学滤波提取车辆目标;步骤三、根据扩展卡尔曼滤波方法的最小外接矩形框的质心坐标,获取车辆的运行轨迹,并跟踪运动目标;步骤四、采集交叉口交通流轨迹速度参数数据;步骤五、基于实时的交通流数据,创建车辆的交通冲突判别模型和冲突计数方法。本发明提供了基于视频的车辆运行轨迹的交通冲突判别和计数方法,可以实时检测交叉口的交通冲突并计数,更加的高效准确,通过与人工计数法获得的交通冲突数相比较,准确度较高。
Description
技术领域
本发明属于交通工程领域,具体涉及一种基于实时车辆轨迹的交叉口交通冲突判别方法。
背景技术
平面交叉口作为道路系统中的重要组成部分,由于交通量大、冲突点多及视线盲区大,并且是各向交通流的汇集点,因而存在着诸多交通安全风险。目前,交通冲突的采集方法主要是通过人工观察的手段,这种方法实时性差且耗时费力,采集的信息不够准确。随着智能交通技术的发展,利用视频检测技术对交通冲突进行判别计数成为了可能。但现有利用视频技术检测交通冲突系统的研究还存在着诸多问题,如目标检测不准确,跟踪效果不好等。
发明内容
本发明为解决上述问题,提出一种基于实时车辆轨迹的交叉口交通冲突判别方法。
本发明涉及一种基于实时车辆轨迹的交叉口交通冲突判别方法,包括如下步骤:
步骤一、采集实时视频图像数据,对其进行图像灰度化、中值滤波和二值化的预处理;
步骤二、根据混合高斯背景模型建立视频背景模型,并通过形态学滤波提取车辆目标;
步骤三、根据扩展卡尔曼滤波方法的最小外接矩形框的质心坐标,获取车辆的运行轨迹,并跟踪运动目标;
步骤四、采集交叉口交通流轨迹速度参数数据;
步骤五、基于实时的交通流数据,创建车辆的交通冲突判别模型和冲突计数方法。
进一步地,步骤一中,图像灰度化、中值滤波和二值化的具体过程为:
(1)采用移动平均法对视频图像灰度化,将彩色像素转换为灰度像素;
(2)对灰度化的图像进行中值进行滤波去噪处理;二维中值滤波的公式如下:
g(x,y)=med{f(x-k,y-l),(k,l∈W)}
其中,f(x,y),g(x,y)分别为原始图像和处理后图像;
(3)选择一个合适的阈值,利用图像二值化对图像中感兴趣的部分分离,把整个图像分割为两部分,灰度值大于这个基准值的部分灰度赋值为0,灰度值小于这个基准值的部分灰度赋值为255,最佳阈值由Otsu法来确定,Otsu法公式如下:
类间方差σ2=p1(μ1-μ)2+p2(μ2-μ)2
式中μ1、μ2——两部分各自的灰度均值;
p1、p2——两部分各自的概率;
μ——图像整体灰度均值;
图像阈值化处理的一般变换函数表达式为:
式中E(i,j)——像素的的灰度
t——选定的阈值。
进一步地,步骤二中,根据混合高斯背景模型建立视频背景模型的步骤如下:
(1)初始化K组高斯模型
当前帧中的像素X,分别与K个模型做差分,若第i组差分结果小于阂值,称像素Xt,与第i组高斯模型匹配,否则称为不匹配;
(2)更新权重,更新公式为:
ωi,t=(1-a)ωi,t-1+Mi,ta
式中a为更新速率,当像素Xt与第i组模型匹配时,Mi=1,否则Mi=0;
(3)更新与像素Xt匹配的模型,更新公式为:
μt=(1-ρ)μt-1+ρXt
式中,ρ为更新速率,α为更新速率。均值与标准差的更新速率不是固定值,像素Xt与背景模型的匹配度越高,更新速率越快;
(4)在二值图像中标记运动目标;
(5)返回步骤二处理下一帧。
进一步地,步骤二中,通过形态学滤波提取车辆目标,是对二值化及差分处理以后的图像存在的许多孤立的点、孤立的小区域、小间隙和孔洞进行数学形态学滤波处理,先腐蚀消除目标图像中的无用点,再膨胀将与目标接触的所有点合并到该目标,对车辆图像进行检测。
进一步地,步骤三中,首先提取外接矩形跟踪框,找到最小外接矩形的质心坐标,进而计算相关的交通流参数。
进一步地,步骤三中,提取外接矩形跟踪框时,运动区域的最小外接矩形坐标信息被定义为两个坐标点(xmin,ymin),(xmax,ymax),各个值具体求解公式为:
用边界框标识的运动区域外接矩形框的四个坐标为A(xmin,ymin),B(xmax,ymin),C(xmin,ymax),D(xmax,ymax)。
进一步地,所述外接矩形框的跟踪,即只有0和255的二值图像中,0代表背景像素,255代表运动区域像素;当第一次出现大于0.5的像素点时,记录这时的行号记为t1,当最后一次出现大于0.5的像素点时,记录这时的行号为t2;以相同的方式从左到右扫描经过处理的二值图像,对每一列的值为255的像素点进行计数,当第一次出现大于0.5的像素点时,记录这时的列号记为l1,当最后一次出现大于0.5的像素点时,记录这时的列号为l2。在记录了t1、t2、l1、l2这四个值后,行间距TΔ为TΔ=t1-t2,列间距LΔ为:LΔ=l1-l2,根据坐标即可画出最小外接矩形框。
进一步地,跟踪质心获取车辆轨迹时,根据外接矩形框的四个坐标A(xmin,ymin),B(xmax,ymin),C(xmin,ymax),D(xmax,ymax)得到最小外接矩阵的高度为:H=Xmax-Xmin;宽度为:L=Ymax-Ymin;质心的坐标为(Xmin+H/2,Ymin+L/2),对提取到的车辆记录其每一帧或其中若干帧图像中的矩形框质心用于数据统计,连续帧的矩形框质心即为车辆的行驶轨迹。
进一步地,步骤四中,通过获取的车辆相关数据进行交通流参数的计算,可得时间间隔t=∑Δt,Δt为视频的每帧图像时间间隔;车辆流率ht为车头时距,即连续车辆经过同一检测点的时间间隔;车辆速度d为车辆驶过的距离,t为时间间隔。
进一步地,步骤五中,具体操作步骤如下:
(1)对当前车辆的运动轨迹和前进方向为参考,考虑其它各个方向的交通流与当前车辆冲突的可能性,可以将冲突再细分为七类:直行与对向左转冲突、同向左转冲突、同向右转冲突、同向直行冲突、右转车流与对向左转车流的合流冲突、同向直行车流与右转车流的分流冲突和同向直行车流与左转车流的分流冲突;
(2)对交通冲突进行判别,判断两车的行驶轨迹是否为碰撞轨迹,如果两个运动物体在同一时刻的轨迹有交叉,则认为发生交通冲突,交通冲突判别的具体的步骤为:
步骤a:获取车辆速度参数和轨迹位置参数;
步骤b:计算两辆车的距离,判断是否车辆间的距离相近;
步骤c:计算最小距离,判断是否有潜在碰撞可能;
步骤d:若有潜在碰撞可能,计算临界冲突距离;
步骤e:根据当前的两辆车的距离,判断两辆车处于何种状态;
步骤f:若车辆小于临界冲突距离,则发生交通冲突车辆处于交通冲突状态;
(3)在交通冲突识别后,结合视频图像数据,对车辆在某段时间内发生的交通冲突进行计数,具体过程为:
步骤a:基于获取的车辆轨迹和速度参数,判断车辆是否足够接近;
步骤b:车辆足够接近,判断车辆是否处于平稳运行状态;
步骤c:车辆处于平稳运行状态,根据基于车辆轨迹进行交通冲突的识别方法判断车辆状态;
步骤d:两辆车处于冲突状态,对车辆进行标记;
步骤e:标记处于冲突状态的车辆后,对冲突状态的车辆进行跟踪;
步骤f:判断其是否有运动状态改变的行为;
步骤g:判断其是否重新回归到平稳运行状态,若车辆重新回归到平稳运行状态,则记为一次交通冲突,并对车辆解除冲突状态标记。
本发明提供了基于视频的车辆运行轨迹的交通冲突判别和计数方法,可以实时检测交叉口的交通冲突并计数,更加的高效准确,通过与人工计数法获得的交通冲突数相比较,准确度较高。
附图说明
图1为本发明的技术流程概图;
图2为本发明详细的技术流程图;
图3a为本发明的第99帧原图像;图3b本发明的灰度化处理的图像;图3c为本发明中值滤波后的图像;图3d为本发明二值化处理后的图像;
图4为本发明的运动目标提取例图;
图5为本发明的运动目标跟踪例图;
图6为本发明的摄像扫描线示意图;
图7为本发明跟踪车辆质心坐标获取的实际坐标、速度、方向角及加速度信息;
图7a为本发明中直行与对向左转冲突示意图;图7b为本发明中同向左转冲突示意图;图7c为本发明中同向右转冲突示意图;图7d为本发明中同向直行冲突示意图;图7e为本发明中右转车流与对向左转冲突示意图;图7f为本发明中同向直行与右转冲突示意图;图7g为本发明中同向直行与左转冲突冲突示意图;
图8为本发明中两车辆之间的冲突示意图;
图9为本发明与人工计数法在不同时段下的交通冲突数示意图。
具体实施方式
实施例1
本实施例涉及一种基于实时车辆轨迹的交叉口交通冲突判别方法的具体过程,如图1所示:
步骤一、采集实时视频图像数据,对其进行图像灰度化、中值滤波和二值化的预处理;
步骤二、根据混合高斯背景模型建立视频背景模型,并通过形态学滤波提取车辆目标;
步骤三、根据扩展卡尔曼滤波方法的最小外接矩形框的质心坐标,获取车辆的运行轨迹,并跟踪运动目标;
步骤四、采集交叉口交通流轨迹速度参数数据;
步骤五、基于实时的交通流数据,创建车辆的交通冲突判别模型和冲突计数方法。
实施例2
结合图2说明本实施方式,以哈尔滨市嵩山路与淮河路交叉口为例,该交叉口右转车不受信号限制,因此会出现直行与右转车发生交通冲突。
本实施例中,上述步骤一对图像灰度化、中值滤波和二值化进行预处理的具体过程为:
(1)采用移动平均法对视频图像灰度化,将彩色像素转换为灰度像素公式为
GRAY=0.299R+0.587G+0.114B
式中R、G、B——原图像的红、绿、蓝分量;
处理前后的的图像见图3a和3b。
(2)对灰度化的图像进行中值滤波去噪处理见图3c,二维中值滤波的公式如下:
g(x,y)=med{f(x-k,y-l),(k,l∈W)}
其中,f(x,y),g(x,y)分别为原始图像和处理后图像。
(3)利用图像二值化对图像中感兴趣的部分分离,选择一个合适的阈值,把整个图像分割为两部分,灰度值大于这个基准值的部分灰度赋值为0(黑色),灰度值小于这个基准值的部分灰度赋值为255(白色),最佳阈值由Otsu法(最大类间方差法)来确定,见图3d,Otsu法公式如下:
类间方差σ2=p1(μ1-μ)2+p2(μ2-μ)2
式中μ1、μ2——两部分各自的灰度均值;
p1、p2——两部分各自的概率;
μ——图像整体灰度均值。
图像阈值化处理的一般变换函数表达式为:
式中E(i,j)——像素的的灰度
t——选定的阈值
实施例3
结合图2和图4来说明此实施例,本实施例中,上述步骤二根据混合高斯背景模型建立视频背景模型,并通过形态学滤波提取车辆目标的具体过程为:
(1)利用混合高斯分布模型建立视频的背景模型,进行背景差分提取车辆图像,步骤如下:
步骤A、初始化K组高斯模型;
当前帧中的像素X,分别与K个模型做差分,若第i组差分结果小于阂值(一般取2.5σ),称像素Xt,与第i组高斯模型匹配,否则称为不匹配;
步骤B、更新权重,更新公式为:
ωi,t=(1-a)ωi,t-1+Mi,ta
式中a为更新速率,当像素Xt与第i组模型匹配时,Mi=1,否则Mi=0
步骤C、更新与像素Xt匹配的模型,更新公式为:
μt=(1-ρ)μt-1+ρXt
图中X是被处理的对象,B是结构元素。
实施例4
结合图2和图5来说明此实施例,在上述步骤三中,根据扩展卡尔曼滤波方法的最小外接矩形框的质心坐标来获取车辆的运行轨迹,进行运动目标的跟踪,具体过程为:
(1)将卡尔曼滤波扩展到了非线性的系统,得到EKF算法:关于系统的方程:
状态方程:xk+1=f(xk,uk)+wk
输出方程:yk=h(xk)+vk
在每一个步长里计算下面的导数矩阵和评价状态估计。
执行以下的卡尔曼滤波方程:
Kk=PkCT(CPkCT+R)-1
Pk+1=A(I-KkC)PkAT+Q
(2)外接矩形跟踪框的提取,运动区域的最小外接矩形坐标信息被定义为两个坐标点(xmin,ymin),(xmax,ymax),各个值具体求解公式为:
用边界框标识的运动区域外接矩形框的四个坐标为A(xmin,ymin),B(xmax,ymin),C(xmin,ymax),D(xmax,ymax)。
(3)外接矩形框的跟踪,即只有0和255的二值图像中,0代表背景像素,255代表运动区域像素。当第一次出现大于0.5的像素点时(0.5为本文中所选取的判定值,可选取别的数值)记录这时的行号记为t1,当最后一次出现大于0.5的像素点时,记录这时的行号为t2。以相同的方式从左到右扫描经过处理的二值图像,对每一列的值为255的像素点进行计数,当第一次出现大于0.5的像素点时,记录这时的列号记为l1,当最后一次出现大于0.5的像素点时,记录这时的列号为l2。在记录了t1、t2、l1、l2这四个值后,行间距TΔ为TΔ=t1-t2,列间距LΔ为:LΔ=l1-l2,根据坐标即可画出最小外接矩形框。
(4)跟踪质心获取车辆轨迹,根据外接矩形框的四个坐标为A(xmin,ymin),B(xmax,ymin),C(xmin,ymax),D(xmax,ymax)得最小外接矩阵的高度为:H=Xmax-Xmin;宽度为:L=Ymax-Ymin;所以质心的坐标为(Xmin+H/2,Ymin+L/2),对提取到的车辆记录其每一帧或其中若干帧图像中的矩形框质心用于数据统计,连续帧的矩形框质心即为车辆的行驶轨迹。
实施例5
结合图2、图6和表1来说明此实施例,本实施中,上述步骤四根据以上步骤,获取交叉口交通流轨迹速度参数,具体过程为:
(1)通过获取的车辆相关数据进行交通流参数的计算,可得时间间隔t=∑Δt,Δt为视频的每帧图像时间间隔。
设车辆识别特征点所处的像素位置坐标为(x1,y1),下一帧中在列上所处的位置坐标为(x2,y2),则像素距离如图6中,m处为摄像头,h为摄像头距离地面的垂直距离。假设特征点移动s2像素距离所对应的实际距离为d2,设置相互之间夹角相同的若干扫描线,夹角为B、C、D…角度为A的扫描线为最近的一条扫描线,所以可以求出每条扫描线对应的实际距离,每一条扫描线即代表了每一列像素的像素距离。我们可使用数学方法求出扫描线间的实际距离d1、d2、d3,实际距离与像素距离的关系如下:
d1=h×tanA
d1+d2=h×tan(A+B)
d1+d2+d3=h×tan(A+B+C)
……
最后得到跟踪的Car70质心坐标获取的实际坐标、速度、方向角及加速度信息如图7所示。
表1
实施例6
结合图2和图8来说明此实施例,上述步骤五基于实时的交通流数据提出了车辆的交通冲突判别模型和冲突计数方法,具体过程为:
(1)以当前车辆的运动轨迹和前进方向为参考,考虑其它各个方向的交通流与当前车辆冲突的可能性,可以将冲突再细分为七类:直行与对向左转冲突(如图7a所示)、同向左转冲突(如图7b所示)、同向右转冲突(如图7c所示)、同向直行冲突(如图7d所示)、右转车流与对向左转车流的合流冲突(如图7e所示)、同向直行车流与右转车流的分流冲突(如图7f所示)和同向直行车流与左转车流的分流冲突(如图7g所示)。
(2)对交通冲突进行判别,判断两车的行驶轨迹是否为碰撞轨迹,如果两个运动物体在同一时刻的轨迹有交叉,则认为发生交通冲突,具体的步骤为:
1)潜在冲突状态预测
当两辆车以当前状态行驶时,可以通过运动学公式判别其在未来是否会发生碰撞。其中车辆用质点表示,则两辆车不同时间的距离可以表示为:
式中,D(t)表示t时刻两点距离,x(t),y(t)表示t时刻目标在坐标系中的位置。若车辆做匀速运动,则有
xi(t)=xi(0)+vi(0)tcosθi
yi(t)=yi(0)+vi(0)tsinθi
其中,xi(0)、yi(0)表示目标0时刻的位置,xi(t)、yi(t)表示目标t时刻的位置,θi表示目标运动方向与x轴夹角,vi(0)表示目标的速度。
则两目标间的距离表示为:
则当t=tm时,D(t)取得最小值dmin
2)临界冲突距离计算
临界冲突距离就是驾驶员从发现危险开始并采取措施直到车辆停止,车辆间刚好不发生碰撞所得到的距离。计算公式为:
其中,S为临界冲突距离,sstop为制动距离,Sfeel为反应距离,v0为初始速度,a为制动减速度,tfeel为反应时间。
3)交通冲突判别
步骤a:获取车辆速度参数和轨迹位置参数;
步骤b:计算两辆车的距离,判断是否车辆间的距离相近;
步骤c:计算最小距离,判断是否有潜在碰撞可能;
步骤d:若有潜在碰撞可能,计算临界冲突距离;
步骤e:根据当前的两辆车的距离,判断两辆车处于何种状态;
步骤f:若车辆小于临界冲突距离,则发生交通冲突车辆处于交通冲突状态。
(3)在交通冲突识别后还要结合视频图像数据对车辆在一段时间内发生的交通冲突进行计数,具体过程为:
步骤a:基于获取的车辆轨迹和速度参数,判断车辆是否足够接近;
步骤b:车辆足够接近,判断车辆是否处于平稳运行状态;
步骤c:车辆处于平稳运行状态,根据基于车辆轨迹进行交通冲突的识别方法判断车辆状态;
步骤d:两辆车处于冲突状态,对车辆进行标记;
步骤e:标记处于冲突状态的车辆后,对冲突状态的车辆进行跟踪;
步骤f:判断其是否有运动状态改变的行为,如减速、加速;
步骤g:判断其是否重新回归到平稳运行状态,若车辆重新回归到平稳运行状态,则记为一次交通冲突,并对车辆解除冲突状态标记。
实施例7
结合图9和表2来说明此实施方式,本实施例中,以哈尔滨市嵩山路与淮河路交叉口为例,验证此方法的合理性,按照上述方法对该交叉口的交通冲突情况进行实时检测,识别交通冲突的个数并且判断交通冲突的类型,根据该交叉口的实际的相位方案,共得到五种交通类型,分别为同向左转冲突、同向右转冲突、同向直行冲突、同向直行与右转冲突和右转与对向左转冲突,并进行交通冲突计数。
以7:30-7:45时段下同向直行与右转的冲突为例,进行交通冲突判别。对选取的发生冲突的两辆车编号A和B,A车从西向南右转,B车从北向南直行。通过获取的两辆车的坐标,方向角以及速度等参数,对同一时刻两辆车相应的参数A车的坐标为(-107,26),方向角为169度,车速5.95m/s,B车的坐标为(-54,89),方向角为91度,车速为7.16m/s。由两辆车的运动状态,存在潜在的碰撞可能。由于右转车的通行权比直行车低,因此右转车会减速避让直行车。只需要计算A车的临界冲突距离即可,司机的反应时间为1s,A车的减速度为1.5m/s2,则临界冲突距离为计算两辆车潜在的最小冲突距离dmin和到达此状态所用的时间tm,若按照此状态继续行驶下去,车辆在
9.89s后存在最小距离为4.98m,两辆车位置与预测碰撞点的最小距离小于临界冲突距离,存在交通冲突。通过与实际相对比,右转车与直行车有明显的冲突行为,右转车减速避让直行车通过交叉口。通过对该车的运动轨迹跟踪,冲突结束后,该车恢复至原先的平稳行驶状态,至此交通冲突结束,视为一次交通冲突并计数。为验证本发明的准确性,将人工统计得到的交通冲突数和本发明的统计的交通冲突数相对比,最终验证了该发明的可行性。
表2
本发明的基于实时车辆轨迹的交叉口交通冲突判别方法,提出了一种基于混合高斯背景模型的交通流视频检测方法来对实时的车辆视频数据进行车辆识别,跟踪并提取以扩展卡尔曼滤波方法的最小外接矩形框的质心坐标来获取车辆的运行轨迹,来求取速度等交通流数据,提出了车辆的交通冲突判别模型和冲突计数方法,对交叉口的交通冲突类型进行判别。
首先对交叉口的视频图像进行灰度化处理,将将彩色像素转换为灰度像素,便于提取信息;对灰度化的图像进行中值滤波去噪处理,利用阈值分割法进行图像二值化,将车辆运动部分分离出来。
再利用混合高斯分布模型建立视频的背景模型,进行背景差分提取车辆图像;对二值化及差分处理以后的图像存在的许多孤立的点、孤立的小区域、小间隙和孔洞进行数学形态学滤波处理,先腐蚀消除目标图像中的无用点,再膨胀将与目标接触的所有点合并到该目标,对车辆图像进行检测。然后以扩展卡尔曼滤波方法的最小外接矩形框的质心坐标来获取车辆的运行轨迹的方法,提取外接矩形跟踪框,找到最小外接矩形的质心坐标,进而计算相关的交通流参数。
最后根据计算获得的交通流参数判断两车的行驶轨迹是否为碰撞轨迹,计算临界冲突距离并根据潜在冲突状态预测判断是否会发生冲突,并结合车辆的运行状态对判断交通冲突类型并计数。
Claims (1)
1.一种基于实时车辆轨迹的交叉口交通冲突判别方法,其特征在于,包括如下步骤:
步骤一、采集实时视频图像数据,对其进行图像灰度化、中值滤波和二值化的预处理;
步骤二、根据混合高斯背景模型建立视频背景模型,并通过形态学滤波提取车辆目标;
步骤三、根据扩展卡尔曼滤波方法的最小外接矩形框的质心坐标,获取车辆的运行轨迹,并跟踪运动目标;
步骤四、采集交叉口交通流轨迹速度参数数据;
步骤五、基于实时的交通流数据,创建车辆的交通冲突判别模型和冲突计数方法;
步骤一中,图像灰度化、中值滤波和二值化的具体过程为:
(1)采用移动平均法对视频图像灰度化,将彩色像素转换为灰度像素;
(2)对灰度化的图像进行中值进行滤波去噪处理;二维中值滤波的公式如下:
g(x,y)=med{f(x-k,y-l),(k,l∈W)}
其中,f(x,y),g(x,y)分别为原始图像和处理后图像;
(3)选择一个特定的阈值,利用图像二值化对图像中感兴趣的部分分离,把整个图像分割为两部分,灰度值大于这个基准值的部分灰度赋值为0,灰度值小于这个基准值的部分灰度赋值为255,最佳阈值由Otsu法来确定,Otsu法公式如下:
类间方差σ2=p1(μ1-μ)2+p2(μ2-μ)2
式中μ1、μ2——两部分各自的灰度均值;
p1、p2——两部分各自的概率;
μ——图像整体灰度均值;
图像阈值化处理的一般变换函数表达式为:
式中E(i,j)——像素的灰度
t——选定的阈值;
步骤二中,根据混合高斯背景模型建立视频背景模型的步骤如下:
(1)初始化K组高斯模型
当前帧中的像素X,分别与K个模型做差分,若第i组差分结果小于阂值,称像素Xt,与第i组高斯模型匹配,否则称为不匹配;
(2)更新权重,更新公式为:
ωi,t=(1-a)ωi,t-1+Mi,ta
式中a为更新速率,当像素Xt与第i组模型匹配时,Mi=1,否则Mi=0;
(3)更新与像素Xt匹配的模型,更新公式为:
μt=(1-ρ)μt-1+ρXt
式中,ρ为更新速率,α为更新速率;均值与标准差的更新速率不是固定值,像素Xt与背景模型的匹配度越高,更新速率越快;
(4)在二值图像中标记运动目标;
(5)返回步骤二处理下一帧;
步骤二中,通过形态学滤波提取车辆目标,是对二值化及差分处理以后的图像存在的许多孤立的点、孤立的小区域、小间隙和孔洞进行数学形态学滤波处理,先腐蚀消除目标图像中的无用点,再膨胀将与目标接触的所有点合并到该目标,对车辆图像进行检测;
步骤三中,首先提取外接矩形跟踪框,找到最小外接矩形的质心坐标,进而计算相关的交通流参数;
步骤三中,提取外接矩形跟踪框时,运动区域的最小外接矩形坐标信息被定义为两个坐标点(xmin,ymin),(xmax,ymax),各个值具体求解公式为:
用边界框标识的运动区域外接矩形框的四个坐标为A(xmin,ymin),B(xmax,ymin),C(xmin,ymax),D(xmax,ymax);
所述外接矩形框的跟踪,即只有0和255的二值图像中,0代表背景像素,255代表运动区域像素;当第一次出现大于0.5的像素点时,记录这时的行号记为t1,当最后一次出现大于0.5的像素点时,记录这时的行号为t2;以相同的方式从左到右扫描经过处理的二值图像,对每一列的值为255的像素点进行计数,当第一次出现大于0.5的像素点时,记录这时的列号记为l1,当最后一次出现大于0.5的像素点时,记录这时的列号为l2;在记录了t1、t2、l1、l2这四个值后,行间距TΔ为TΔ=t1-t2,列间距LΔ为:LΔ=l1-l2,根据坐标即可画出最小外接矩形框;
跟踪质心获取车辆轨迹时,根据外接矩形框的四个坐标A(xmin,ymin),B(xmax,ymin),C(xmin,ymax),D(xmax,ymax)得到最小外接矩阵的高度为:H=Xmax-Xmin;宽度为:L=Ymax-Ymin;质心的坐标为(Xmin+H/2,Ymin+L/2),对提取到的车辆记录其每一帧或其中若干帧图像中的矩形框质心用于数据统计,连续帧的矩形框质心即为车辆的行驶轨迹;
步骤四中,通过获取的车辆相关数据进行交通流参数的计算,可得时间间隔t=∑Δt,Δt为视频的每帧图像时间间隔;车辆流率ht为车头时距,即连续车辆经过同一检测点的时间间隔;车辆速度d为车辆驶过的距离,t为时间间隔;
步骤五中,具体操作步骤如下:
(1)对当前车辆的运动轨迹和前进方向为参考,考虑其它各个方向的交通流与当前车辆冲突的可能性,将冲突再细分为七类:直行与对向左转冲突、同向左转冲突、同向右转冲突、同向直行冲突、右转车流与对向左转车流的合流冲突、同向直行车流与右转车流的分流冲突和同向直行车流与左转车流的分流冲突;
(2)对交通冲突进行判别,判断两车的行驶轨迹是否为碰撞轨迹,如果两个运动物体在同一时刻的轨迹有交叉,则认为发生交通冲突;
具体的步骤为:
1)潜在冲突状态预测
当两辆车以当前状态行驶时,通过运动学公式判别其在未来是否会发生碰撞;其中车辆用质点表示,则两辆车不同时间的距离表示为:
式中,D(t)表示t时刻两点距离,x(t),y(t)表示t时刻目标在坐标系中的位置;若车辆做匀速运动,则有
xi(t)=xi(0)+vi(0)tcosθi
yi(t)=yi(0)+vi(0)tsinθi
其中,xi(0)、yi(0)表示目标0时刻的位置,xi(t)、yi(t)表示目标t时刻的位置,θi表示目标运动方向与x轴夹角,vi(0)表示目标的速度;
则两目标间的距离表示为:
则当t=tm时,D(t)取得最小值dmin
2)临界冲突距离计算
临界冲突距离就是驾驶员从发现危险开始并采取措施直到车辆停止,车辆间刚好不发生碰撞所得到的距离;计算公式为:
其中,S为临界冲突距离,Sstop为制动距离,Sfeel为反应距离,v0为初始速度,a为制动减速度,tfeel为反应时间;
3)交通冲突判别的具体的步骤为:
步骤a:获取车辆速度参数和轨迹位置参数;
步骤b:计算两辆车的距离,判断是否车辆间的距离相近;
步骤c:计算最小距离,判断是否有潜在碰撞可能;
步骤d:若有潜在碰撞可能,计算临界冲突距离;
步骤e:根据当前的两辆车的距离,判断两辆车处于何种状态;
步骤f:若车辆小于临界冲突距离,则发生交通冲突车辆处于交通冲突状态;
(3)在交通冲突识别后,结合视频图像数据,对车辆在某段时间内发生的交通冲突进行计数,具体过程为:
步骤a:基于获取的车辆轨迹和速度参数,判断车辆是否足够接近;
步骤b:车辆足够接近,判断车辆是否处于平稳运行状态;
步骤c:车辆处于平稳运行状态,根据基于车辆轨迹进行交通冲突的识别方法判断车辆状态;
步骤d:两辆车处于冲突状态,对车辆进行标记;
步骤e:标记处于冲突状态的车辆后,对冲突状态的车辆进行跟踪;
步骤f:判断其是否有运动状态改变的行为;
步骤g:判断其是否重新回归到平稳运行状态,若车辆重新回归到平稳运行状态,则记为一次交通冲突,并对车辆解除冲突状态标记。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010287744.1A CN111445699B (zh) | 2020-04-13 | 2020-04-13 | 基于实时车辆轨迹的交叉口交通冲突判别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010287744.1A CN111445699B (zh) | 2020-04-13 | 2020-04-13 | 基于实时车辆轨迹的交叉口交通冲突判别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111445699A CN111445699A (zh) | 2020-07-24 |
CN111445699B true CN111445699B (zh) | 2021-10-26 |
Family
ID=71654096
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010287744.1A Active CN111445699B (zh) | 2020-04-13 | 2020-04-13 | 基于实时车辆轨迹的交叉口交通冲突判别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111445699B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112200131A (zh) * | 2020-10-28 | 2021-01-08 | 鹏城实验室 | 一种基于视觉的车辆碰撞检测方法、智能终端及存储介质 |
CN112561940B (zh) * | 2020-12-08 | 2022-07-12 | 中国人民解放军陆军工程大学 | 密集多目标参数提取的方法、装置及终端设备 |
CN113192331B (zh) * | 2021-04-26 | 2022-04-15 | 吉林大学 | 一种网联环境下面向骑行安全的智能预警系统及预警方法 |
CN113553986B (zh) * | 2021-08-02 | 2022-02-08 | 浙江索思科技有限公司 | 一种船舶上运动目标检测方法及系统 |
CN113744529A (zh) * | 2021-09-06 | 2021-12-03 | 上海辉略企业发展集团有限公司 | 一种基于交通冲突点的交通安全事故黑点识别方法及系统 |
CN114267189B (zh) * | 2021-12-20 | 2022-11-18 | 黑龙江工程学院 | 一种快速路出口匝道与衔接交叉口联合控制方法 |
CN114155715B (zh) * | 2022-02-07 | 2022-05-06 | 北京图盟科技有限公司 | 一种冲突点检测方法、装置、设备及可读存储介质 |
CN114241398A (zh) * | 2022-02-23 | 2022-03-25 | 深圳壹账通科技服务有限公司 | 基于人工智能的车辆定损方法、装置、设备及存储介质 |
CN114815857B (zh) * | 2022-06-28 | 2022-10-11 | 广东邦盛北斗科技股份公司 | 基于北斗导航的智能化农机管理方法、系统及云平台 |
CN115273467B (zh) * | 2022-07-15 | 2023-07-04 | 南京莱斯信息技术股份有限公司 | 一种基于多区域检测数据的车流相互影响事件识别方法 |
CN115168810B (zh) * | 2022-09-08 | 2022-11-29 | 南京慧尔视智能科技有限公司 | 一种交通数据的生成方法、装置、电子设备和存储介质 |
CN115578862B (zh) * | 2022-10-26 | 2023-09-19 | 中国建设基础设施有限公司 | 交通流量转换方法、装置、计算设备和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002056494A (ja) * | 2000-08-10 | 2002-02-22 | Gen Tec:Kk | 道路上を走行する車両のトラッキング方法及び該方法を実施するプログラムを記録した記録媒体 |
CN102368351B (zh) * | 2011-10-19 | 2014-10-29 | 北京航空航天大学 | 一种消解无信号交叉口两车交通冲突的方法 |
CN104282020A (zh) * | 2014-09-22 | 2015-01-14 | 中海网络科技股份有限公司 | 一种基于目标运动轨迹的车辆速度检测方法 |
CN107103292A (zh) * | 2017-04-12 | 2017-08-29 | 湖南源信光电科技股份有限公司 | 一种基于运动车辆跟踪的车流量统计方法 |
CN109643367A (zh) * | 2016-07-21 | 2019-04-16 | 御眼视觉技术有限公司 | 用于自主车辆导航的众包和分发稀疏地图以及车道测量 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4480299B2 (ja) * | 2001-06-21 | 2010-06-16 | 富士通マイクロエレクトロニクス株式会社 | 移動体を含む画像の処理方法及び装置 |
US20070239809A1 (en) * | 2006-04-06 | 2007-10-11 | Michael Moseler | Method for calculating a local extremum, preferably a local minimum, of a multidimensional function E(x1, x2, ..., xn) |
JP4982472B2 (ja) * | 2008-12-15 | 2012-07-25 | 日本電信電話株式会社 | 映像監視システム、映像監視方法および映像監視プログラム |
CN108922176B (zh) * | 2018-06-27 | 2021-05-14 | 华南理工大学 | 一种混合交通冲突态势的快速判断方法 |
CN109598946B (zh) * | 2018-11-19 | 2022-01-07 | 南京理工大学 | 一种基于雷达体制的多车道测速方法 |
CN110298122B (zh) * | 2019-07-03 | 2021-04-27 | 北京理工大学 | 基于冲突消解的无人驾驶车辆城市交叉口左转决策方法 |
-
2020
- 2020-04-13 CN CN202010287744.1A patent/CN111445699B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002056494A (ja) * | 2000-08-10 | 2002-02-22 | Gen Tec:Kk | 道路上を走行する車両のトラッキング方法及び該方法を実施するプログラムを記録した記録媒体 |
CN102368351B (zh) * | 2011-10-19 | 2014-10-29 | 北京航空航天大学 | 一种消解无信号交叉口两车交通冲突的方法 |
CN104282020A (zh) * | 2014-09-22 | 2015-01-14 | 中海网络科技股份有限公司 | 一种基于目标运动轨迹的车辆速度检测方法 |
CN109643367A (zh) * | 2016-07-21 | 2019-04-16 | 御眼视觉技术有限公司 | 用于自主车辆导航的众包和分发稀疏地图以及车道测量 |
CN107103292A (zh) * | 2017-04-12 | 2017-08-29 | 湖南源信光电科技股份有限公司 | 一种基于运动车辆跟踪的车流量统计方法 |
Non-Patent Citations (3)
Title |
---|
Investigate moped-car conflicts in China using a naturalistic driving;Yi G. Glaser et al.;《Journal of Safety Research》;20171021;第63卷;171-175 * |
基于交通冲突理论的信号交叉口安全评价技术;郭延永;《中国博士学位论文全文数据库工程科技Ⅱ辑》;20170215(第2期);C034-100 * |
基于视频的平面交叉口机动车交通冲突检测技术研究;张方方;《中国优秀硕士学位论文全文数据库工程科技II辑》;20080815(第8期);C034-266 * |
Also Published As
Publication number | Publication date |
---|---|
CN111445699A (zh) | 2020-07-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111445699B (zh) | 基于实时车辆轨迹的交叉口交通冲突判别方法 | |
Datondji et al. | A survey of vision-based traffic monitoring of road intersections | |
Kanhere et al. | Real-time incremental segmentation and tracking of vehicles at low camera angles using stable features | |
JP6519262B2 (ja) | 立体物検出装置、立体物検出方法、立体物検出プログラム、及び移動体機器制御システム | |
Hadi et al. | Vehicle detection and tracking techniques: a concise review | |
Malik | Fast vehicle detection with probabilistic feature grouping and its application to vehicle tracking | |
EP1780675B1 (en) | Object detector | |
Nieto et al. | Road environment modeling using robust perspective analysis and recursive Bayesian segmentation | |
CN110379168B (zh) | 一种基于Mask R-CNN的交通车辆信息获取方法 | |
KR100969995B1 (ko) | 영상처리기법을 이용한 신호교차로의 교통상충 판단 시스템및 방법 | |
JPH06337938A (ja) | 画像処理装置及び距離測定装置 | |
CN113516853B (zh) | 一种针对复杂监控场景的多车道车流量检测方法 | |
CN115113206B (zh) | 井下有轨车辅助驾驶的行人与障碍物检测方法 | |
Li et al. | Robust vehicle tracking for urban traffic videos at intersections | |
Fries et al. | Autonomous convoy driving by night: The vehicle tracking system | |
Javadi et al. | A robust vision-based lane boundaries detection approach for intelligent vehicles | |
Jaiswal et al. | Comparative analysis of CCTV video image processing techniques and application: a survey | |
CN115240471A (zh) | 一种基于图像采集的智慧厂区避撞预警方法和系统 | |
Liu et al. | Lane-level and full-cycle multivehicle tracking using low-channel roadside LiDAR | |
Kini | Real time moving vehicle congestion detection and tracking using OpenCV | |
Khan | Estimating Speeds and Directions of Pedestrians in Real-Time Videos: A solution to Road-Safety Problem. | |
Kanhere et al. | Real-time detection and tracking of vehicle base fronts for measuring traffic counts and speeds on highways | |
CN113177509A (zh) | 一种倒车行为识别方法及装置 | |
Bhaskar et al. | Enhanced and effective parallel optical flow method for vehicle detection and tracking | |
CN114581863A (zh) | 车辆危险状态识别方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |