CN107194386A - A kind of intersection electric bicycle travel speed acquisition methods based on video - Google Patents

A kind of intersection electric bicycle travel speed acquisition methods based on video Download PDF

Info

Publication number
CN107194386A
CN107194386A CN201710593862.3A CN201710593862A CN107194386A CN 107194386 A CN107194386 A CN 107194386A CN 201710593862 A CN201710593862 A CN 201710593862A CN 107194386 A CN107194386 A CN 107194386A
Authority
CN
China
Prior art keywords
mrow
msub
electric bicycle
video
frame 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.)
Pending
Application number
CN201710593862.3A
Other languages
Chinese (zh)
Inventor
何杰
张燕
吴冠鹤
成诚
刘子洋
卢文慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN201710593862.3A priority Critical patent/CN107194386A/en
Publication of CN107194386A publication Critical patent/CN107194386A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/40Scenes; Scene-specific elements in video content
    • G06V20/41Higher-level, semantic clustering, classification or understanding of video scenes, e.g. detection, labelling or Markovian modelling of sport events or news items
    • G06V20/42Higher-level, semantic clustering, classification or understanding of video scenes, e.g. detection, labelling or Markovian modelling of sport events or news items of sport video content
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/50Context or environment of the image
    • G06V20/52Surveillance or monitoring of activities, e.g. for recognising suspicious objects
    • G06V20/53Recognition of crowd images, e.g. recognition of crowd congestion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/50Context or environment of the image
    • G06V20/56Context or environment of the image exterior to a vehicle by using sensors mounted on the vehicle
    • G06V20/58Recognition of moving objects or obstacles, e.g. vehicles or pedestrians; Recognition of traffic objects, e.g. traffic signs, traffic lights or roads
    • G06V20/584Recognition of moving objects or obstacles, e.g. vehicles or pedestrians; Recognition of traffic objects, e.g. traffic signs, traffic lights or roads of vehicle lights or traffic lights

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Multimedia (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Software Systems (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开了一种基于视频的道路交叉口电动自行车行驶速度获取方法,包括如下步骤:1、采集道路交叉口的车辆行驶视频,设置视频中的检测区域为多条人行道围合的区域;2、对采集的视频进行处理,识别跟踪视频中道路交叉口的电动自行车,获取行驶轨迹坐标点;3、将电动自行车在视频图像上的质心坐标转换为实际二维平面坐标,计算其行驶速度。该方法能够准确获取道路交叉口内部电动自行车的行驶速度。

The invention discloses a video-based method for acquiring the driving speed of an electric bicycle at a road intersection, comprising the following steps: 1. Collecting a vehicle driving video at a road intersection, and setting the detection area in the video as an area surrounded by multiple sidewalks; 2. 1. Process the collected video, identify and track the electric bicycle at the intersection of the road in the video, and obtain the coordinate points of the driving trajectory; 3. Convert the centroid coordinates of the electric bicycle on the video image into actual two-dimensional plane coordinates, and calculate its driving speed. This method can accurately obtain the driving speed of the electric bicycle inside the road intersection.

Description

一种基于视频的道路交叉口电动自行车行驶速度获取方法A video-based method for obtaining the speed of electric bicycles at road intersections

技术领域technical field

本发明属于交通监控领域,具体涉及一种利用计算机视觉技术获取电动自行车行驶速度的方法。The invention belongs to the field of traffic monitoring, and in particular relates to a method for obtaining the running speed of an electric bicycle by using computer vision technology.

背景技术Background technique

近年来,电动自行车在交通出行结构中的占比日渐增长。电动自行车相较于步行出行更加节省体力,且速度更快;而与汽车相比,它具有电力驱动无排放、环境友好污染少、不怕拥堵且成本低的优点。电动自行车这种快捷、清洁且成本低廉的出行方式十分受到人们的青睐,也逐渐成为城市中短途出行的主要交通方式。In recent years, the proportion of electric bicycles in the transportation structure has been increasing. Compared with walking, electric bicycles are more energy-saving and faster; compared with cars, they have the advantages of electric drive, no emissions, environmental friendliness, less pollution, no fear of congestion and low cost. Electric bicycles, a fast, clean and low-cost way of travel, are very popular among people, and have gradually become the main mode of transportation for short-distance travel in cities.

然而,电动自行车这种出行方式在我国迅速发展与普及的同时,也给城市道路交叉口增添了拥堵、安全等问题。电动自行车骑行者不顾道路安全法律法规,如在交叉口内经常发生电动自行车驶入机动车辆的行驶空间、违反交通信号、超速抢行和电动车主私自改装电动车进行货物运输等行为,是导致电动自行车行驶安全性低的重要因素。如何合理分析、规范电动自行车在交叉口的驾驶行为,还需获取电动自行车的行驶参数,包括行驶速度、安全行驶区域和安全行驶速度。However, while the travel mode of electric bicycles is rapidly developing and popularizing in my country, it also brings problems such as congestion and safety to urban road intersections. Electric bicycle riders disregard road safety laws and regulations, such as electric bicycles often entering the driving space of motor vehicles at intersections, violating traffic signals, speeding and rushing, and electric vehicle owners modifying electric vehicles for cargo transportation without permission, which is the cause of electric bicycles. An important factor of low driving safety. How to reasonably analyze and standardize the driving behavior of electric bicycles at intersections requires obtaining the driving parameters of electric bicycles, including driving speed, safe driving area and safe driving speed.

目前国内外对电动自行车的微观交通数据的提取主要为模糊的数学推算,这种调查方式与人的主观能动性关系很大,调查人员对数据的观测具有不稳定性,推算出的微观交通数据往往精度较低。At present, the extraction of micro-traffic data of electric bicycles at home and abroad is mainly based on fuzzy mathematical calculations. This investigation method has a lot to do with people's subjective initiative. The investigators' observation of the data is unstable, and the calculated micro-traffic data is often The precision is lower.

发明内容Contents of the invention

发明目的:针对现有技术中存在的问题,本发明提供了一种基于视频的道路交叉口电动自行车行驶速度获取方法,该方法能够准确获取道路交叉口内部电动自行车的行驶速度。Purpose of the invention: Aiming at the problems existing in the prior art, the present invention provides a video-based method for obtaining the speed of electric bicycles at road intersections, which can accurately obtain the speed of electric bicycles inside road intersections.

技术方案:本发明采用如下技术方案:一种基于视频的道路交叉口电动自行车行驶速度获取方法,包括如下步骤:Technical solution: The present invention adopts the following technical solution: a video-based method for obtaining the driving speed of an electric bicycle at a road intersection, comprising the following steps:

(1)采集道路交叉口的车辆行驶视频,设置视频中的检测区域为多条人行道围合的区域;(1) Collect the vehicle driving video at the road intersection, and set the detection area in the video as the area surrounded by multiple sidewalks;

(2)对采集的视频进行处理,识别跟踪视频中道路交叉口的电动自行车,获取行驶轨迹坐标点;(2) Process the collected video, identify and track the electric bicycle at the road intersection in the video, and obtain the coordinate points of the driving track;

(3)将电动自行车在视频图像上的质心坐标转换为实际二维平面坐标,计算其行驶速度。(3) Convert the coordinates of the center of mass of the electric bicycle on the video image to the actual two-dimensional plane coordinates, and calculate its driving speed.

步骤(2)包括如下步骤:Step (2) comprises the following steps:

(2-1)建立所采集的道路交叉口视频的背景图像模型f(x,y);可以采用混合高斯背景建模法建立背景图像模型f(x,y);(2-1) set up the background image model f (x, y) of the collected road intersection video; Can adopt the Gaussian mixture background modeling method to set up the background image model f (x, y);

(2-2)将视频中第t帧原始图像与步骤(2-1)建立的背景图像模型作差,得到第t帧前景运动目标A(xt,yt);前景运动目标计算公式为:(2-2) Make a difference between the original image of the tth frame in the video and the background image model established in step (2-1), and obtain the foreground moving target A(x t , y t ) of the tth frame; the calculation formula of the foreground moving target is: :

A(xt,yt)=f(xt,yt)-f(x,y)A(x t ,y t )=f(x t ,y t )-f(x,y)

其中,f(xt,yt)为第t帧原始图像,f(x,y)为背景图像;Among them, f(x t , y t ) is the original image of frame t, and f(x, y) is the background image;

(2-3)选取视频中多帧图像,对每一帧图像中的前景运动目标为电动自行车的前景团块,计算其团块面积、团块外接矩形的长度、团块外接矩形的宽度;(2-3) select multiple frames of images in the video, and the foreground moving target in each frame of images is the foreground blob of an electric bicycle, calculate its blob area, the length of the blob circumscribed rectangle, the width of the blob circumscribed rectangle;

(2-4)建立电动自行车参数模型,所述模型的参数为电动自行车团块面积S,团块外接矩形的长度L、团块外接矩形的宽度W;对(2-3)中获取的参数进行统计,选取分布概率为η的参数值域作为电动自行车参数模型的取值范围;(2-4) set up electric bicycle parameter model, the parameter of described model is electric bicycle agglomerate area S, the length L of agglomerate circumscribed rectangle, the width W of agglomerate circumscribed rectangle; To the parameter that obtains in (2-3) Carry out statistics, select the parameter range that distribution probability is η as the value range of electric bicycle parameter model;

(2-5)对于视频中第i帧图像和第j帧图像,获取前景运动目标Ai和Aj,计算对应的前景团块面积、团块外接矩形的长度、团块外接矩形的宽度,根据电动自行车参数模型中参数取值范围,判断Ai和Aj是否为电动自行车;(2-5) For the i-th frame image and the j-th frame image in the video, obtain the foreground moving objects A i and A j , calculate the corresponding foreground blob area, the length of the blob circumscribed rectangle, and the blob circumscribed rectangle width, According to the value range of parameters in the electric bicycle parameter model, judge whether A i and A j are electric bicycles;

(2-6)如果Ai和Aj均为电动自行车,计算Ai和Aj的质心pi和pj,质心距离dij,如果满足条件:则Ai和Aj为同一辆电动自行车,pi和pj即为所述电动自行车在第i帧图像和第j帧图像中的行驶轨迹坐标点。(2-6) If both A i and A j are electric bicycles, calculate the centroids p i and p j of A i and A j , and the centroid distance d ij , if the conditions are met: Then A i and A j are the same electric bicycle, and p i and p j are the driving trajectory coordinate points of the electric bicycle in the i-th frame image and the j-th frame image.

步骤(3)包括如下步骤:Step (3) comprises the following steps:

(3-1)电动自行车在第a帧图像和第b帧图像中的实际行驶距离dab为:(3-1) The actual driving distance d ab of the electric bicycle in the first frame image and the second frame image is:

其中(ua,va)为第a帧图像中电动自行车的二维地面坐标;(ub,vb)为第b帧图像中电动自行车的二维地面坐标;Where (u a , v a ) is the two-dimensional ground coordinates of the electric bicycle in the a-th frame image; (u b , v b ) is the two-dimensional ground coordinates of the electric bicycle in the b-th frame image;

(3-2)电动自行车在第a帧图像和第b帧图像间的平均速度vab为:(3-2) The average speed v ab of the electric bicycle between the first frame image and the second frame image is:

其中tab为第a帧图像和第b帧图像之间的时间差。Where t ab is the time difference between the a-th frame image and the b-th frame image.

电动自行车在第f帧图像的瞬时速度Speedf为:The instantaneous speed Speed f of the image in frame f of the electric bicycle is:

其中r为视频的帧率,(uf,vf)为第f帧图像中电动自行车的二维地面坐标,(uf+1,vf+1)为第f+1帧图像中电动自行车的二维地面坐标。Where r is the frame rate of the video, (u f , v f ) is the two-dimensional ground coordinates of the electric bicycle in the fth frame image, (u f+1 , v f+1 ) is the electric bicycle in the f+1th frame image The two-dimensional ground coordinates of .

有益效果:与现有技术相比,本发明公开的基于视频的道路交叉口电动自行车行驶速度获取方法具有以下优点:1、能够准确获取到交叉口内部的电动自行车的行驶轨迹和行驶速度;2、对不同交叉口所拍摄的视频源有着检测方法可重复、准确率高和参数易修改的优良特性,对电动自行车在交叉口的行驶行为分析有着重要的意义。Beneficial effects: Compared with the prior art, the video-based method for obtaining the driving speed of electric bicycles at road intersections disclosed by the present invention has the following advantages: 1. It can accurately obtain the driving trajectory and driving speed of electric bicycles inside the intersection; 2. 1. The video sources taken at different intersections have the excellent characteristics of repeatable detection methods, high accuracy and easy modification of parameters, which is of great significance to the analysis of driving behavior of electric bicycles at intersections.

附图说明Description of drawings

图1为本发明提供的方法的流程图;Fig. 1 is the flowchart of the method provided by the present invention;

图2为电动自行车跟踪步骤的详细流程图;Figure 2 is a detailed flow chart of the electric bicycle tracking steps;

图3为原始图像与背景图像检测区域图;Fig. 3 is an original image and a background image detection area diagram;

图4为前景图像降噪前后前景图像对比图;Fig. 4 is a comparison diagram of the foreground image before and after noise reduction of the foreground image;

图5为电动自行车筛选效果图;Fig. 5 is the effect drawing of electric bicycle screening;

图6为道路交叉口电动车安全行驶范围图;Fig. 6 is a diagram of the safe driving range of electric vehicles at road intersections;

图7为道路交叉口电动自行车最高速度小时分布图。Fig. 7 is the hourly distribution diagram of the maximum speed of electric bicycles at road intersections.

具体实施方式detailed description

本实施例以南京市中山东路-太平北路交叉口的电动自行车行驶速度的获取为例,进一步阐明本发明。In this embodiment, the acquisition of the driving speed of an electric bicycle at the intersection of Zhongshan East Road-Taiping North Road in Nanjing is taken as an example to further illustrate the present invention.

如图1所示,本发明提供一种基于视频的道路交叉口电动自行车行驶速度获取方法,包括如下步骤:As shown in Figure 1, the present invention provides a kind of video-based method for obtaining the traveling speed of an electric bicycle at a road intersection, comprising the following steps:

(1)采集道路交叉口的车辆行驶视频,设置视频中的检测区域为多条人行道围合的区域;(1) Collect the vehicle driving video at the road intersection, and set the detection area in the video as the area surrounded by multiple sidewalks;

以俯视视角拍摄中山东路-太平北路十字道路交叉口;设置检测区域1为该交叉口多条人行道围合的区域,如图3(b)中的多边形区域。Shoot the Zhongshan East Road-Taiping North Road intersection from a top-down perspective; set the detection area 1 as the area surrounded by multiple sidewalks at the intersection, such as the polygonal area in Figure 3(b).

(2)对采集的视频进行处理,识别跟踪视频中道路交叉口的电动自行车,获取行驶轨迹坐标点;如图2所示,具体包括如下步骤:(2) process the collected video, identify and track the electric bicycle at the road intersection in the video, and obtain the driving trajectory coordinate points; as shown in Figure 2, specifically include the following steps:

(2-1)建立所采集的道路交叉口视频的背景图像模型f(x,y);本实施例采用混合高斯背景建模法建立背景图像模型f(x,y);(2-1) set up the background image model f (x, y) of the collected road intersection video; The present embodiment adopts the mixed Gaussian background modeling method to set up the background image model f (x, y);

(2-2)将视频中第t帧原始图像与步骤(2-1)建立的背景图像模型作差,得到第t帧前景运动目标A(xt,yt);前景运动目标计算公式为:(2-2) Make a difference between the original image of the tth frame in the video and the background image model established in step (2-1), and obtain the foreground moving target A(x t , y t ) of the tth frame; the calculation formula of the foreground moving target is: :

A(xt,yt)=f(xt,yt)-f(x,y)A(x t ,y t )=f(x t ,y t )-f(x,y)

其中,f(xt,yt)为第t帧原始图像,f(x,y)为背景图像;Among them, f(x t , y t ) is the original image of frame t, and f(x, y) is the background image;

得到前景运动目标A(xt,yt)后,对A(xt,yt)进行滤波操作,以去除电动自行车阴影部分和噪声,本实施例中对A(xt,yt)进行高斯滤波,滤波效果如图4所示,其中(a)为滤波前的前景图像;(b)为滤波后的前景图像。After obtaining the foreground moving target A(x t , y t ), filter operation is performed on A(x t , y t ) to remove the shadow part and noise of the electric bicycle. In this embodiment, A(x t , y t ) is filtered Gaussian filtering, the filtering effect is shown in Figure 4, where (a) is the foreground image before filtering; (b) is the foreground image after filtering.

(2-3)选取视频中多帧图像,对每一帧图像中的前景运动目标为电动自行车的前景团块,计算其团块面积、团块外接矩形的长度、团块外接矩形的宽度;(2-3) select multiple frames of images in the video, and the foreground moving target in each frame of images is the foreground blob of an electric bicycle, calculate its blob area, the length of the blob circumscribed rectangle, the width of the blob circumscribed rectangle;

(2-4)建立电动自行车参数模型,所述模型的参数为电动自行车团块面积S,团块外接矩形的长度L、团块外接矩形的宽度W;对步骤(2-3)中获取的参数进行统计,选取分布概率为80%的参数值域作为电动自行车参数模型的取值范围;(2-4) set up electric bicycle parameter model, the parameter of described model is electric bicycle agglomerate area S, the length L of agglomerate circumscribed rectangle, the width W of agglomerate circumscribed rectangle; To the step (2-3) that obtains The parameters are counted, and the parameter value range with a distribution probability of 80% is selected as the value range of the electric bicycle parameter model;

(2-5)对于视频中第i帧图像和第j帧图像,获取前景运动目标Ai和Aj,计算对应的前景团块面积、团块外接矩形的长度、团块外接矩形的宽度,根据电动自行车参数模型中参数取值范围,判断Ai和Aj是否为电动自行车;如图5所示,图中方框为判断为电动自行车的前景图像。(2-5) For the i-th frame image and the j-th frame image in the video, obtain the foreground moving objects A i and A j , calculate the corresponding foreground blob area, the length of the blob circumscribed rectangle, and the blob circumscribed rectangle width, According to the value range of the parameters in the electric bicycle parameter model, judge whether A i and A j are electric bicycles; as shown in Figure 5, the box in the figure is the foreground image judged to be electric bicycles.

(2-6)如果Ai和Aj均为电动自行车,计算Ai和Aj的质心pi和pj,质心距离dij,如果两帧图像中的电动自行车为同一辆,那么其在极短时间内所行驶的距离不会超过一个极小的数值,因此设置帧差阈值fth来控制时间,同时设置距离阈值dth来控制距离。本实施例中设置dth为20像素,fth为20帧。如果满足条件:则Ai和Aj为同一辆电动自行车,pi和pj即为所述电动自行车在第i帧图像和第j帧图像中的行驶轨迹坐标点。(2-6) If both A i and A j are electric bicycles, calculate the centroids p i and p j of A i and A j , and the centroid distance d ij , if the electric bicycles in the two frames of images are the same, then their The distance traveled in a very short time will not exceed a very small value, so set the frame difference threshold f th to control the time, and set the distance threshold d th to control the distance. In this embodiment, d th is set to 20 pixels, and f th is set to 20 frames. If the conditions are met: Then A i and A j are the same electric bicycle, and p i and p j are the driving trajectory coordinate points of the electric bicycle in the i-th frame image and the j-th frame image.

质心计算公式为:The formula for calculating the centroid is:

其中,mpq为图像的p+q阶矩,(xc,yc)为图像的质心,A(x,y)为检测到的电动自行车团块图像,N和M分别为图像长度和宽度值。Among them, m pq is the p+q order moment of the image, (x c , y c ) is the centroid of the image, A(x, y) is the detected electric bicycle clump image, N and M are the length and width of the image respectively value.

第i帧图像和第j帧图像中团块质心间的距离计算公式为:The formula for calculating the distance between the mass centroids of the i-th frame image and the j-th frame image is:

其中(xci,yci)为第i帧图像中电动自行车团块的质心坐标。Where (x ci , y ci ) is the coordinates of the center of mass of the electric bicycle clump in the i-th frame image.

(3)将电动自行车在视频图像上的质心坐标转换为实际二维平面坐标,计算其行驶速度。(3) Convert the coordinates of the center of mass of the electric bicycle on the video image to the actual two-dimensional plane coordinates, and calculate its driving speed.

本实施例中运用齐次坐标转换方法将像素坐标转换为二维地面坐标,转换距离与实际距离误差分析如表1所示:In this embodiment, the homogeneous coordinate conversion method is used to convert the pixel coordinates into two-dimensional ground coordinates. The error analysis between the conversion distance and the actual distance is shown in Table 1:

表1Table 1

在这四个样本中,检验的所有误差率均在5%以下,平均误差率为3%。该交叉口长为20米,电动自行车驶过该交叉口的总误差为60厘米,对于电动自行车的冲突分析可以忽略。Across the four samples, all error rates for the tests were below 5%, with an average error rate of 3%. The intersection is 20 meters long, and the total error of the electric bicycle driving through the intersection is 60 cm, which can be ignored for the conflict analysis of electric bicycles.

计算两帧图像间电动自行车平均速度的步骤为:The steps to calculate the average speed of the electric bicycle between two frames of images are:

(3-1)电动自行车在第a帧图像和第b帧图像中的实际行驶距离dab为:(3-1) The actual driving distance d ab of the electric bicycle in the first frame image and the second frame image is:

其中(ua,va)为第a帧图像中电动自行车的二维地面坐标;(ub,vb)为第b帧图像中电动自行车的二维地面坐标;Where (u a , v a ) is the two-dimensional ground coordinates of the electric bicycle in the a-th frame image; (u b , v b ) is the two-dimensional ground coordinates of the electric bicycle in the b-th frame image;

(3-2)电动自行车在第a帧图像和第b帧图像间的平均速度vab为:(3-2) The average speed v ab of the electric bicycle between the first frame image and the second frame image is:

其中tab为第a帧图像和第b帧图像之间的时间差。Where t ab is the time difference between the a-th frame image and the b-th frame image.

电动自行车在第f帧图像的瞬时速度Speedf为:The instantaneous speed Speed f of the image in frame f of the electric bicycle is:

其中r为视频的帧率,(uf,vf)为第f帧图像中电动自行车的二维地面坐标,(uf+1,vf+1)为第f+1帧图像中电动自行车的二维地面坐标。Where r is the frame rate of the video, (u f , v f ) is the two-dimensional ground coordinates of the electric bicycle in the fth frame image, (u f+1 , v f+1 ) is the electric bicycle in the f+1th frame image The two-dimensional ground coordinates of .

对通过道路交叉口的电动自行车行驶速度进行统计,可以得到在该道路交叉口内的安全行驶区域与安全行驶速度,步骤如下:The safe driving area and safe driving speed in the road intersection can be obtained by making statistics on the driving speed of electric bicycles passing through the road intersection. The steps are as follows:

统计电动自行车经过道路交叉口的最高速度,国家规定的最大行驶车速20km/s,即5.5m/s,对未超速行驶的电动自行车获取其行驶轨迹累计图,绘制交叉口内部安全行驶区域图,如图6所示。Statize the maximum speed of electric bicycles passing through road intersections. The maximum driving speed stipulated by the state is 20km/s, that is, 5.5m/s. For electric bicycles that are not speeding, obtain the cumulative trajectory of their driving trajectories, and draw a map of the safe driving area inside the intersection. As shown in Figure 6.

统计高峰小时内通过道路交叉口电动自行车的最高速度,最高速度分布情况如表2所示,最高速度分布图如图7所示。The maximum speed of electric bicycles passing through road intersections in peak hours is counted. The distribution of the maximum speed is shown in Table 2, and the distribution of the maximum speed is shown in Figure 7.

表2Table 2

由表2和图7可知:电动自行车最高速度分布在5.5-6.5m/s、6.5-7.5m/s、7.5-8.5m/s、8.5m/s以上这四个区间内均为超速行驶的电动车,其分布概率分别为5%、4%、2%和0%,超速行驶的电动车最大速度总体分布概率为11%。From Table 2 and Figure 7, it can be seen that the maximum speed of electric bicycles is distributed in the four intervals above 5.5-6.5m/s, 6.5-7.5m/s, 7.5-8.5m/s, and 8.5m/s. For electric vehicles, the distribution probabilities are 5%, 4%, 2% and 0% respectively, and the overall distribution probability of the maximum speed of the electric vehicles driving over the speed limit is 11%.

未超速行驶的电动自行车最高速度分布在2.5-3.5m/s之间最为集中,概率为46%。将未超速的速度数据由小到大排列,去除前15%的速度数据与后15%的速度数据,选取剩余速度参数分布区间作为安全行驶速度区间为[2.5,5.0],单位为m/s。The maximum speed distribution of electric bicycles without speeding is most concentrated between 2.5-3.5m/s, with a probability of 46%. Arrange the non-speeding speed data from small to large, remove the first 15% of the speed data and the last 15% of the speed data, and select the remaining speed parameter distribution interval as the safe driving speed interval [2.5, 5.0], the unit is m/s .

Claims (5)

1.一种基于视频的道路交叉口电动自行车行驶速度获取方法,其特征在于,包括如下步骤:1. A video-based method for obtaining the driving speed of an electric bicycle at a road intersection, is characterized in that, comprising the steps: (1)采集道路交叉口的车辆行驶视频,设置视频中的检测区域为多条人行道围合的区域;(1) Collect the vehicle driving video at the road intersection, and set the detection area in the video as the area surrounded by multiple sidewalks; (2)对采集的视频进行处理,识别跟踪视频中道路交叉口的电动自行车,获取行驶轨迹坐标点;(2) Process the collected video, identify and track the electric bicycle at the road intersection in the video, and obtain the coordinate points of the driving track; (3)将电动自行车在视频图像上的质心坐标转换为实际二维平面坐标,计算其行驶速度。(3) Convert the coordinates of the center of mass of the electric bicycle on the video image to the actual two-dimensional plane coordinates, and calculate its driving speed. 2.根据权利要求1所述的基于视频的道路交叉口电动自行车行驶速度获取方法,其特征在于,步骤(2)包括如下步骤:2. the video-based road crossing electric bicycle travel speed acquisition method according to claim 1, is characterized in that, step (2) comprises the steps: (2-1)建立所采集的道路交叉口视频的背景图像模型f(x,y);(2-1) set up the background image model f (x, y) of the road crossing video that is collected; (2-2)将视频中第t帧原始图像与步骤(2-1)建立的背景图像模型作差,得到第t帧前景运动目标A(xt,yt);前景运动目标计算公式为:(2-2) Make a difference between the original image of the tth frame in the video and the background image model established in step (2-1), and obtain the foreground moving target A(x t , y t ) of the tth frame; the calculation formula of the foreground moving target is: : A(xt,yt)=f(xt,yt)-f(x,y)A(x t ,y t )=f(x t ,y t )-f(x,y) 其中,f(xt,yt)为第t帧原始图像,f(x,y)为背景图像;Among them, f(x t , y t ) is the original image of frame t, and f(x, y) is the background image; (2-3)选取视频中多帧图像,对每一帧图像中的前景运动目标为电动自行车的前景团块,计算其团块面积、团块外接矩形的长度、团块外接矩形的宽度;(2-3) select multiple frames of images in the video, and the foreground moving target in each frame of images is the foreground blob of an electric bicycle, calculate its blob area, the length of the blob circumscribed rectangle, the width of the blob circumscribed rectangle; (2-4)建立电动自行车参数模型,所述模型的参数为电动自行车团块面积S,团块外接矩形的长度L、团块外接矩形的宽度W;对(2-3)中获取的参数进行统计,选取分布概率为η的参数值域作为电动自行车参数模型的取值范围;(2-4) set up electric bicycle parameter model, the parameter of described model is electric bicycle agglomerate area S, the length L of agglomerate circumscribed rectangle, the width W of agglomerate circumscribed rectangle; To the parameter that obtains in (2-3) Carry out statistics, select the parameter range that distribution probability is η as the value range of electric bicycle parameter model; (2-5)对于视频中第i帧图像和第j帧图像,获取前景运动目标Ai和Aj,计算对应的前景团块面积、团块外接矩形的长度、团块外接矩形的宽度,根据电动自行车参数模型中参数取值范围,判断Ai和Aj是否为电动自行车;(2-5) For the i-th frame image and the j-th frame image in the video, obtain the foreground moving objects A i and A j , calculate the corresponding foreground blob area, the length of the blob circumscribed rectangle, and the blob circumscribed rectangle width, According to the value range of parameters in the electric bicycle parameter model, judge whether A i and A j are electric bicycles; (2-6)如果Ai和Aj均为电动自行车,计算Ai和Aj的质心pi和pj,质心距离dij,如果满足条件:则Ai和Aj为同一辆电动自行车,pi和pj即为所述电动自行车在第i帧图像和第j帧图像中的行驶轨迹坐标点。(2-6) If both A i and A j are electric bicycles, calculate the centroids p i and p j of A i and A j , and the centroid distance d ij , if the conditions are met: Then A i and A j are the same electric bicycle, and p i and p j are the driving trajectory coordinate points of the electric bicycle in the i-th frame image and the j-th frame image. 3.根据权利要求1所述的基于视频的道路交叉口电动自行车行驶速度获取方法,其特征在于,步骤(3)包括如下步骤:3. the video-based road crossing electric bicycle traveling speed acquisition method according to claim 1, is characterized in that, step (3) comprises the steps: (3-1)电动自行车在第a帧图像和第b帧图像中的实际行驶距离dab为:(3-1) The actual driving distance d ab of the electric bicycle in the first frame image and the second frame image is: <mrow> <msub> <mi>d</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> <mrow> <msub> <mi>d</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>a</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>b</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> 其中(ua,va)为第a帧图像中电动自行车的二维地面坐标;(ub,vb)为第b帧图像中电动自行车的二维地面坐标;Where (u a , v a ) is the two-dimensional ground coordinates of the electric bicycle in the a-th frame image; (u b , v b ) is the two-dimensional ground coordinates of the electric bicycle in the b-th frame image; (3-2)电动自行车在第a帧图像和第b帧图像间的平均速度vab为:(3-2) The average speed v ab of the electric bicycle between the first frame image and the second frame image is: <mrow> <msub> <mi>v</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>d</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <msub> <mi>t</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mfrac> </mrow> <mrow> <msub> <mi>v</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <mo>=</mo> <mfrac> <msub> <mi>d</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> <msub> <mi>t</mi> <mrow> <mi>a</mi> <mi>b</mi> </mrow> </msub> </mfrac> </mrow> 其中tab为第a帧图像和第b帧图像之间的时间差。Where t ab is the time difference between the a-th frame image and the b-th frame image. 4.根据权利要求1所述的基于视频的道路交叉口电动自行车行驶速度获取方法,其特征在于,电动自行车在第f帧图像的瞬时速度Speedf为:4. the method for obtaining the speed of electric bicycles at road intersections based on video according to claim 1, wherein the instantaneous speed Speed f of electric bicycles at the f frame image is: <mrow> <msub> <mi>Speed</mi> <mi>f</mi> </msub> <mo>=</mo> <mi>r</mi> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>f</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mi>f</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> <mrow> <msub> <mi>Speed</mi> <mi>f</mi> </msub> <mo>=</mo> <mi>r</mi> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>u</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mrow> <mi>f</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>v</mi> <mi>f</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mrow> <mi>f</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> 其中r为视频的帧率,(uf,vf)为第f帧图像中电动自行车的二维地面坐标,(uf+1,vf+1)为第f+1帧图像中电动自行车的二维地面坐标。Where r is the frame rate of the video, (u f , v f ) is the two-dimensional ground coordinates of the electric bicycle in the fth frame image, (u f+1 , v f+1 ) is the electric bicycle in the f+1th frame image The two-dimensional ground coordinates of . 5.根据权利要求2所述的基于视频的道路交叉口电动自行车行驶速度获取方法,其特征在于,步骤(2-1)采用混合高斯背景建模法建立背景图像模型f(x,y)。5. The method for obtaining the speed of electric bicycles at road intersections based on video according to claim 2, wherein step (2-1) adopts a mixed Gaussian background modeling method to establish a background image model f (x, y).
CN201710593862.3A 2017-07-20 2017-07-20 A kind of intersection electric bicycle travel speed acquisition methods based on video Pending CN107194386A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710593862.3A CN107194386A (en) 2017-07-20 2017-07-20 A kind of intersection electric bicycle travel speed acquisition methods based on video

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710593862.3A CN107194386A (en) 2017-07-20 2017-07-20 A kind of intersection electric bicycle travel speed acquisition methods based on video

Publications (1)

Publication Number Publication Date
CN107194386A true CN107194386A (en) 2017-09-22

Family

ID=59884101

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710593862.3A Pending CN107194386A (en) 2017-07-20 2017-07-20 A kind of intersection electric bicycle travel speed acquisition methods based on video

Country Status (1)

Country Link
CN (1) CN107194386A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615862A (en) * 2018-12-29 2019-04-12 南京市城市与交通规划设计研究院股份有限公司 Road vehicle movement of traffic state parameter dynamic acquisition method and device
CN113380035A (en) * 2021-06-16 2021-09-10 山东省交通规划设计院集团有限公司 Road intersection traffic volume analysis method and system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102592454A (en) * 2012-02-29 2012-07-18 北京航空航天大学 Intersection vehicle movement parameter measuring method based on detection of vehicle side face and road intersection line
CN103971521A (en) * 2014-05-19 2014-08-06 清华大学 Method and device for detecting road traffic abnormal events in real time
CN105786895A (en) * 2014-12-25 2016-07-20 日本电气株式会社 Calculating method and device of discharge amount of road intersection

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102592454A (en) * 2012-02-29 2012-07-18 北京航空航天大学 Intersection vehicle movement parameter measuring method based on detection of vehicle side face and road intersection line
CN103971521A (en) * 2014-05-19 2014-08-06 清华大学 Method and device for detecting road traffic abnormal events in real time
CN105786895A (en) * 2014-12-25 2016-07-20 日本电气株式会社 Calculating method and device of discharge amount of road intersection

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
盛能: "混合交通流中的自行车识别及参数提取", 《计算机应用研究》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109615862A (en) * 2018-12-29 2019-04-12 南京市城市与交通规划设计研究院股份有限公司 Road vehicle movement of traffic state parameter dynamic acquisition method and device
CN113380035A (en) * 2021-06-16 2021-09-10 山东省交通规划设计院集团有限公司 Road intersection traffic volume analysis method and system
CN113380035B (en) * 2021-06-16 2022-11-11 山东省交通规划设计院集团有限公司 Road intersection traffic volume analysis method and system

Similar Documents

Publication Publication Date Title
Ashalatha et al. Critical gap through clearing behavior of drivers at unsignalised intersections
CN102324183B (en) Method for detecting and shooting vehicle based on composite virtual coil
Munigety et al. Semiautomated tool for extraction of microlevel traffic data from videographic survey
CN105513371B (en) A kind of highway parking offense detection method based on Density Estimator
CN107316010A (en) A kind of method for recognizing preceding vehicle tail lights and judging its state
CN106935035A (en) Parking offense vehicle real-time detection method based on SSD neutral nets
CN102073852B (en) Multiple vehicle segmentation method based on optimum threshold values and random labeling method for multiple vehicles
CN104835319A (en) Method for estimating vehicle import behavior on high-grade road bottleneck zone on-ramp
CN103500322A (en) Automatic lane line identification method based on low-altitude aerial images
CN104318781B (en) Based on the travel speed acquisition methods of RFID technique
CN105632186A (en) Method and device for detecting vehicle queue jumping behavior
CN105389996A (en) Traffic operation condition characteristic parameter extraction method based on big data
CN104933859A (en) Macroscopic fundamental diagram-based method for determining bearing capacity of network
CN105069859A (en) Vehicle driving state monitoring method and apparatus thereof
CN105118305A (en) Vehicle management platform of outdoor parking lot exit
CN102156989B (en) Vehicle blocking detection and segmentation method in video frame
CN107507433A (en) A kind of control method of big data analysis dynamic traffic signal system
CN110705484A (en) Method for recognizing illegal behavior of continuously changing lane by using driving track
CN111489555A (en) A traffic operation state prediction method, device and system
CN106548628A (en) The road condition analyzing method that a kind of view-based access control model space transition net is formatted
CN202422420U (en) Illegal parking detection system based on video monitoring
CN105844915A (en) Method for determining traffic flow fundamental diagram in variable speed limit control state
CN105869402A (en) Highway section speed correction method based on multiple types of floating car data
CN107578048A (en) A vehicle detection method in far-sighted scenes based on rough classification of vehicle types
CN107194386A (en) A kind of intersection electric bicycle travel speed acquisition methods based on video

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170922