CN107808524A - 一种基于无人机的道路交叉口车辆检测方法 - Google Patents
一种基于无人机的道路交叉口车辆检测方法 Download PDFInfo
- Publication number
- CN107808524A CN107808524A CN201711239604.1A CN201711239604A CN107808524A CN 107808524 A CN107808524 A CN 107808524A CN 201711239604 A CN201711239604 A CN 201711239604A CN 107808524 A CN107808524 A CN 107808524A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- video
- intersection
- frame
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000012545 processing Methods 0.000 claims abstract description 10
- 241000283070 Equus zebra Species 0.000 claims description 21
- 238000001514 detection method Methods 0.000 claims description 14
- 230000006870 function Effects 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 7
- 238000002187 spin decoupling employing ultra-broadband-inversion sequences generated via simulated annealing Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000003287 optical effect Effects 0.000 claims description 4
- 238000009416 shuttering Methods 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims description 2
- 238000005260 corrosion Methods 0.000 claims description 2
- 230000007797 corrosion Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 238000001914 filtration Methods 0.000 claims description 2
- 238000009499 grossing Methods 0.000 claims description 2
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims description 2
- 230000001629 suppression Effects 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims description 2
- 238000002620 method output Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000011835 investigation Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
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/017—Detecting movement of traffic to be counted or controlled identifying vehicles
- G08G1/0175—Detecting movement of traffic to be counted or controlled identifying vehicles by photographing vehicles, e.g. when violating traffic rules
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/25—Determination of region of interest [ROI] or a volume of interest [VOI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/08—Detecting or categorising vehicles
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及一种基于无人机的道路交叉口车辆检测方法，包括有以下步骤：(1)将无人机悬停在道路交叉口中心的上方，对交叉口行车区域进行垂直拍摄，获得交叉口行车区域的视频；(2)对获得的视频进行稳像处理，去除视频各帧视频帧中的非正常偏移；(3)对于视频中的各帧视频帧，执行以下步骤对其中的交叉口行车区域中的车辆进行检测：S1.提取视频帧的背景图像；S2.识别道路交叉口的行车区域；S3.通过背景差分法获取交叉口行车区域的前景图像；S4.排除交叉口行车区域前景图像中的非车辆目标；S5.得到道路交叉口行车区域中车辆的影像。
Description
技术领域
本发明涉及交叉口车辆检测领域，更具体地，涉及一种基于无人机的道路交叉口车辆检测方法。
背景技术
随着汽车保有量的日益增长，道路上的车辆数量迅速增加，交通压力倍增。平面交叉口是道路交通的枢纽，城市的交通问题往往突出地表现在交叉口上，因此了解平面交叉口的车流量是十分重要的。道路交叉口车流量调查的目的是通过长期连续性观测或短期间隙性和临时性观测，搜集交通量资料，获得有关交通量的实况、通行能力、流向分布、交通量变化及高峰小时交通量和交通组成等方面的资料，了解交通量在时间、空间上的变化和分布规律，为交通规划、道路建设、交通控制与管理、工程经济分析等提供必要的数据，以便对交叉口的运行效能做出准确的评价，提出交通管理、控制措施或改建、扩建方案。
目前常见的交叉口车流量调查方法有：人工计数法、机械计数法、录像法和航摄法。这些常见的方法存在适用性差、费用高等问题，近年来，随着视频处理技术与无人机的发展，越来越多的道路交通调查采用多旋翼无人机进行，特别是道路交叉口的车流量调查。要想通过视频完成交通调查。其中车辆检测是关键，后续的跟踪和流量统计都是建立在高精度的车辆检测基础上的。
常用的车辆检测算法分为两类，一类是将运动前景和背景分离的方法，获取车辆前景区域，如背景差法和帧差法，这种方法依赖于高质量的背景图像和车辆的速度等运动特征。由于无人机存在轻微的抖动、旋转等现象，导致各视频帧的拍摄范围和角度存在一定差别，很难直接从视频得到高质量的背景图像；同时交叉口附近的车辆存在大量的等红灯、转向等运动，或者车身颜色与路面很相近，导致车辆检测的精度不高。第二种是采用机器学习方法，根据车辆的颜色、角点和尺度特征，对车辆样本和非车辆样本进行计算和学习，如基于HOG特征方法、基于HAAR特征方法等，然后直接对视频帧上的车辆进行识别。这种方法需要大量的训练才能达到高精度，对车辆周围环境也很敏感，而各个交叉口的环境恰恰各不相同，使这类算法很难用于车辆检测。
发明内容
本发明为克服目前道路交叉口车辆检测方法存在的适用性差、精度低等问题，针对道路交叉口附近车辆的运动特征以及车辆检测的普适性、高精度等实际要求，提出了一种基于无人机的道路交叉口车辆检测方法。
为实现以上发明目的，采用的技术方案是：
一种基于无人机的道路交叉口车辆检测方法，包括有以下步骤：
(1)将无人机悬停在道路交叉口中心的上方，对交叉口行车区域进行垂直拍摄，获得交叉口行车区域的视频；
(2)对获得的视频进行稳像处理，去除视频各帧视频帧中的非正常偏移；
(3)对于视频中的各帧视频帧，执行以下步骤对其中的交叉口行车区域中的车辆进行检测：
S1.提取视频帧的背景图像；
S2.识别道路交叉口的行车区域；
S3.通过背景差分法获取交叉口行车区域的前景图像；
S4.排除交叉口行车区域前景图像中的非车辆目标；
S5.得到道路交叉口行车区域中车辆的影像。
与现有技术相比，本发明的有益效果是：
本发明通过无人机在道路交叉口上方采集得到的视频来对道路交叉口的车辆进行检测，不需要大量的人力和高成本，只需要设计合适的算法能在无人机拍摄的到的视频中提取出车辆的前景影像，实现车辆检测。本发明有效降低检测成本和操作复杂度，是一种便捷、适用性强的检测方法。
附图说明
图1为方法的流程示意图。
图2为方法的具体实施示意图。
图3为经历过稳像处理的视频图像的示意图。
图4为使用背景差法进行前景提取和车辆识别的示意图。
图5为对前景图像进行去噪处理，以及根据车辆图斑大小阈值进行筛选，得到的车辆识别结果图。
图6为实际被检测出车辆的示意图。
具体实施方式
附图仅用于示例性说明，不能理解为对本专利的限制；
以下结合附图和实施例对本发明做进一步的阐述。
实施例1
如图1所示，本发明涉及一种基于无人机的道路交叉口车辆检测方法，包括有以下步骤：
(1)将无人机悬停在道路交叉口中心的上方，对交叉口行车区域进行垂直拍摄，获得交叉口行车区域的视频；
(2)对获得的视频进行稳像处理，去除视频各帧视频帧中的非正常偏移；
(3)对于视频中的各帧视频帧，执行以下步骤对其中的交叉口行车区域中的车辆进行检测：
S1.提取视频帧的背景图像；
S2.识别道路交叉口的行车区域；
S3.通过背景差分法获取交叉口行车区域的前景图像；
S4.排除交叉口行车区域前景图像中的非车辆目标；
S5.得到道路交叉口行车区域中车辆的影像。
实施例2
本实施例提供了一种基于无人机的道路交叉口车辆检测方法，如图2所示，其包括有以下步骤：
(1)利用无人机悬停在道路交叉口的中心上方，对交叉口范围进行垂直拍摄，获取交叉口车辆等运动目标的视频。
如上所述的无人机的拍摄高度应在能够拍到清晰路面信息的范围内，以便于得到清晰的交叉口信息，从而提取出道路交叉口车流量信息。
(2)利用SUSAN快速角点检测算法对视频各帧视频帧中的特征点进行检测，然后利用光流法基于检测到的特征点进行运动矢量的计算，利用得到的相邻邻视频帧的偏移量x,y对相邻的视频帧进行运动补偿，获得稳定的输出序列。采用光流法消除视频的旋转、平移等现象，实现视频稳像。
其中，所述利用SUSAN快速角点检测算法对视频各帧视频帧中的特征点进行检测的具体过程如下：
S21.以核点像素(x,y)为中心做半径为r的圆，用圆形模板遍历图像；
S22.使用相似比较函数计算模板中各像素点与核点像素的相似度：
其中C(r,r0)为相似比较函数；I(r0)为核点像素的灰度值；I(r)为圆形模板区域内其他像素点的灰度值；t为SUSAN的阈值；
S23.由相似比较函数的计算可以得到每个像素的USAN区域大小，其表达式为
S24.以不同的核点像素为中心做半径为r的圆，对视频帧图像进行遍历，然后按照步骤S22～23的方式计算每个核点像素的USAN区域大小；
S25.得到每个核点像素的USAN区域大小后，再由角点响应函数产生候选角点，最后使用非极大值抑制得到特征点：
式中，g为几何阈值，
而光流法计算运动矢量并进行运动补偿的具体过程如下：
S12.对于视频中相邻的视频帧，以前一帧作为参考帧，后一帧为当前帧，建立光流约束方程：
Ixu+Ixv+It＝0
式中，u,v为x,y方向的2个速度分量，Ix,Iy,It为灰度对于x方向、y方向和时间的偏导，
S13.根据L－K局部平滑假设的条件可知在小范围内光流是一致的，从而得到光流的能量函数为:
Ω为光流一致的邻域范围；W(x)为该邻域范围内不同特征点的权值，令V＝(u,v)T，通过最小二乘法求得方程解如下:
V＝[ATW2A]-1ATW2b
公式中，W＝diag[W(X1),…W(Xn)]，b＝-[It(X1),…It(Xn)]T；
对求解到的u,v去模，即可得到特征点在相邻视频帧的偏移量x,y；
S14.利用得到的相邻邻视频帧的偏移量x,y对相邻的视频帧进行运动补偿，获得稳定的输出序列。
视频拍摄期间无人机存在旋转抖动等现象。如图3所示，由于进行了稳像处理，那些没有像素值的区域全部补0，所以图像周围存在一定的黑边现象，但不影响后续车辆识别。
(3)进行背景建模，利用建立的高斯混合模型提取视频帧的背景图像。
所述高斯混合模型中，视频帧中的每一个像素点按不同的权值的多个高斯分布的叠加来建模，每种高斯分布对应像素点一个可能产生的颜色状态，对于随机变量X的观测数据集{x1,x2,…,xn}，xt为t时刻像素点的灰度值，xt出现的概率为：
ωi,t表示在t时刻像素点的混合高斯模型中第i个高斯分布的权重，且是像素在t时刻的GB彩色向量；n取3；K表示高斯分布的数量，和分别表示t时刻第i个高斯分布的均值向量和协方差矩阵，σ表示标准差，In表示单位矩阵。
视频内容变化会导致像素点发生变化，高斯混合模型需要根据新的像素值不断地学习更新。更新的方法是将混合高斯模型中的K个高斯分布按照由大到小排序，当有新的图像帧到来时，将其中每个像素的彩色向量值Xt与该像素已存在的K个高斯分布逐一匹配检验，如果满足则定义Xt与高斯混合模型中第i个高斯分布相匹配，反之则不匹配；对匹配的第i个高斯参数的参数按下式更新：
μi，t＝(1-ρi，t)μi，t-1+ρi，tXt
其中α是用户自己定义的学习速率，ρi,t为参数学习速率，且如果构建该像素的所有高斯分布中不存在与Xt相匹配的高斯分布，那么将Xt作为均值赋给权重最小的高斯分布，并对该高斯分布重新赋以较小的权重和较大的方差；匹配检验后，按下式更新同一个像素的各个高斯分布的权重：ωi,t＝(1-α)ωi,t-1+αMi,t
上式中，如果第i个高斯分布与Xt匹配，则Mi,t取1，否则取0；更新完成后，各高斯分布的权值被归一化。
全部像素点的高斯混合模型建成后，模型中每个高斯分布按照的值进行排序，然后将前B个高斯分布被选作背景模型，即：
其中T是背景高斯分布在整个概率分布中所占的最小比例，需要不断调整T的值才能获得更好的分布。如果选择一个数值较小的T，则背景模型通常是单峰。如果是这种情况，只使用最可能的分布可以减少处理过程。如果T数值比较大，通过重复的背景运动引起的多模态分布可能会导致一个以上的颜色被包括在背景模型。这将导致透明度效果，允许背景接受两个或多个单独的颜色。
(4)利用等高线法对交叉口的斑马线进行识别，因为道路交叉口一般有斑马线标志的人行横道，各个斑马线首尾连接，围成的区域就是交叉口行车区域。在交叉口行车区域，车辆一般是运动的，只考虑行车区域内的情形，这对背景差法识别车辆非常有利。
等高线法识别斑马线的过程如下：等高线法是一种边缘检测算法，其输入是二值图像，即高斯混合模型背景建模得到的背景图像的二值化。输出是各种地物的轮廓多边形，其中包括斑马线各个条带的多边形。然后根据斑马线的面积特征、形状特征和等间距特征，将其他多边形过滤掉，剩下的基本是斑马线条带的多边形，这样斑马线就识别出来了。由于车辆或行道树遮挡等因素，不能正确识别所有的斑马线条带，需要对斑马线识别结果进一步处理。一般地，斑马线朝向交叉口部分遮挡因素少，多边形形状清晰整齐，所以可以提取各个多边形最靠近交叉口的坐标，再分别对各条斑马线的坐标采用霍夫变换变换进行直线拟合，再对斜率相近的交叉直线进行融合，使每条斑马线只有一条直线。最后对各个直线顺序求交，则以交点为顶点围成的凸多边形就是交叉口行车区域。
(5)通过背景差分法获取交叉口行车区域的前景图像，示意图如图4所示。
所述背景差分法获取前景图像的具体过程如下：
其中fb(x,y)表示视频帧的背景图像，fk(x,y)表示视频帧，Dk(x,y)为前景图像，T为阈值。
(6)使用腐蚀膨胀操作去除交叉口行车区域前景图像里面的细小密集的噪声点，同时保证车辆形状较为完好，使用中值滤波进一步处理去除剩余的噪声，然后使用二值化操作加强目标效果。具体如图5所示。
(7)设定需要识别的目标大小的范围，将非车辆目标和其他因素排除，非机动车、人力三轮车、行人等干扰，得到的目标就是车辆的影像，将它们的范围用矩形标示出来，得到车辆识别结果，如图6所示。
显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.一种基于无人机的道路交叉口车辆检测方法，其特征在于：包括有以下步骤：
(1)将无人机悬停在道路交叉口中心的上方，对交叉口行车区域进行垂直拍摄，获得交叉口行车区域的视频；
(2)对获得的视频进行稳像处理，去除视频各帧视频帧中的非正常偏移；
(3)对于视频中的各帧视频帧，执行以下步骤对其中的交叉口行车区域中的车辆进行检测：
S1.提取视频帧的背景图像；
S2.识别道路交叉口的行车区域；
S3.通过背景差分法获取交叉口行车区域的前景图像；
S4.排除交叉口行车区域前景图像中的非车辆目标；
S5.得到道路交叉口行车区域中车辆的影像。
2.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述步骤(2)对获得的视频进行稳像处理的具体过程如下：
S11.利用SUSAN快速角点检测算法对视频各帧视频帧中的特征点进行检测；
S12.对于视频中相邻的视频帧，以前一帧作为参考帧，后一帧为当前帧，建立光流约束方程：
Ixu+Ixv+It＝0
式中，u,v为x,y方向的2个速度分量，Ix,Iy,It为灰度对于x方向、y方向和时间的偏导，
S13.根据L－K局部平滑假设的条件可知在小范围内光流是一致的，从而得到光流的能量函数为:
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
<mo>&Element;</mo>
<mi>&Omega;</mi>
</mrow>
</munder>
<msup>
<mi>W</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mi>x</mi>
</msub>
<mi>u</mi>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>y</mi>
</msub>
<mi>v</mi>
<mo>+</mo>
<msub>
<mi>I</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
Ω为光流一致的邻域范围；W(x)为该邻域范围内不同特征点的权值，令V＝(u,v)T，通过最小二乘法求得方程解如下:
V＝[ATW2A]-1ATW2b
公式中，W＝diag[W(X1),…W(Xn)]，b＝-[It(X1),…It(Xn)]T；
对求解到的u,v去模，即可得到特征点在相邻视频帧的偏移量x,y；
S14.利用得到的相邻邻视频帧的偏移量x,y对相邻的视频帧进行运动补偿，获得稳定的输出序列。
3.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述利用SUSAN快速角点检测算法对视频各帧视频帧中的特征点进行检测的具体过程如下：
S21.以核点像素(x,y)为中心做半径为r的圆，用圆形模板遍历图像；
S22.使用相似比较函数计算模板中各像素点与核点像素的相似度：
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>exp</mi>
<mo>{</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&lsqb;</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>I</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mi>t</mi>
</mfrac>
<mo>}</mo>
</mrow>
其中C(r,r0)为相似比较函数；I(r0)为核点像素的灰度值；I(r)为圆形模板区域内其他像素点的灰度值；t为SUSAN的阈值；
S23.由相似比较函数的计算可以得到每个像素的USAN区域大小，其表达式为
<mrow>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>r</mi>
<mo>&Element;</mo>
<mi>C</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</munder>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
S24.以不同的核点像素为中心做半径为r的圆，对视频帧图像进行遍历，然后按照步骤S22～23的方式计算每个核点像素的USAN区域大小；
S25.得到每个核点像素的USAN区域大小后，再由角点响应函数产生候选角点，最后使用非极大值抑制得到特征点：
式中，g为几何阈值，
4.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述步骤S1通过高斯混合模型提取出视频帧的背景图像。
5.根据权利要求4所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述高斯混合模型中，视频帧中的每一个像素点按不同的权值的多个高斯分布的叠加来建模，每种高斯分布对应像素点一个可能产生的颜色状态，对于随机变量X的观测数据集{x1,x2,…,xn}，xt为t时刻像素点的灰度值，xt出现的概率为：
<mrow>
<mi>P</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<mfrac>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mi>n</mi>
<mn>2</mn>
</mfrac>
</msup>
<mo>|</mo>
<msub>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msup>
<mo>|</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>X</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>u</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msup>
</mrow>
ωi,t表示在t时刻像素点的混合高斯模型中第i个高斯分布的权重，且 是像素在t时刻的GB彩色向量；n取3；K表示高斯分布的数量，和分别表示t时刻第i个高斯分布的均值向量和协方差矩阵，σ表示标准差，In表示单位矩阵。
6.根据权利要求5所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述高斯混合模型在提取每帧视频帧的背景图像之间，会进行更新，具体如下：
将混合高斯模型中的K个高斯分布按照由大到小排序，当有新的图像帧到来时，将其中每个像素的彩色向量值Xt与该像素已存在的K个高斯分布逐一匹配检验，如果满足则定义Xt与高斯混合模型中第i个高斯分布相匹配，反之则不匹配；对匹配的第i个高斯参数的参数按下式更新：
μi,t＝(1-ρi,t)μi,t-1+ρi,tXt
<mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>&sigma;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中α是用户自己定义的学习速率，ρi,t为参数学习速率，且如果构建该像素的所有高斯分布中不存在与Xt相匹配的高斯分布，那么将Xt作为均值赋给权重最小的高斯分布，并对该高斯分布重新赋以较小的权重和较大的方差；匹配检验后，按下式更新同一个像素的各个高斯分布的权重：ωi,t＝(1-α)ωi,t-1+αMi,t
上式中，如果第i个高斯分布与Xt匹配，则Mi,t取1，否则取0；更新完成后，各高斯分布的权值被归一化。
7.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述步骤S2使用等高线法识别交叉口的斑马线，等高线法输出各种地物的轮廓多边形，包括斑马线各个条带的多边形；根据斑马线的面积特征、形状特征和等间距特征，将其他多边形过滤掉，剩下是斑马线条带的多边形，完成对斑马线的识别；提取各条斑马线多边形最靠近交叉口的坐标，然后分别对各条斑马线的坐标采用霍夫变换变换进行直线拟合，再对斜率相近的交叉直线进行融合；多条斑马线所覆盖的区域即为交叉口行车区域。
8.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述背景差分法获取前景图像的具体过程如下：
其中fb(x,y)表示视频帧的背景图像，fk(x,y)表示视频帧，Dk(x,y)为前景图像，T为阈值。
9.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述步骤S3得到前景图像后，使用腐蚀膨胀操作去除前景图像中的噪声点，使用中值滤波进一步处理去除前景图像中剩余的噪声，然后使用二值化操作加强目标效果。
10.根据权利要求1所述的基于无人机的道路交叉口车辆检测方法，其特征在于：所述步骤S4排除非车辆目标的具体原理如下：设定需要识别的车辆目标的大小范围，将大小范围不符合的目标进行删除。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711239604.1A CN107808524B (zh) | 2017-11-30 | 2017-11-30 | 一种基于无人机的道路交叉口车辆检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711239604.1A CN107808524B (zh) | 2017-11-30 | 2017-11-30 | 一种基于无人机的道路交叉口车辆检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107808524A true CN107808524A (zh) | 2018-03-16 |
| CN107808524B CN107808524B (zh) | 2021-09-24 |
Family
ID=61582090
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711239604.1A Expired - Fee Related CN107808524B (zh) | 2017-11-30 | 2017-11-30 | 一种基于无人机的道路交叉口车辆检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107808524B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109003451A (zh) * | 2018-07-04 | 2018-12-14 | 北京航空航天大学 | 一种基于低空无人机的交叉口od矩阵估算方法 |
| CN109919957A (zh) * | 2019-01-08 | 2019-06-21 | 同济大学 | 一种基于动态视觉传感器的角点检测方法 |
| CN109993134A (zh) * | 2019-04-04 | 2019-07-09 | 中山大学 | 一种基于hog和svm分类器的道路交叉口车辆检测方法 |
| CN111523527A (zh) * | 2020-07-02 | 2020-08-11 | 平安国际智慧城市科技股份有限公司 | 特种运输车辆监控方法、装置、介质及电子设备 |
| CN112542050A (zh) * | 2020-12-23 | 2021-03-23 | 重庆市市政设计研究院有限公司 | 基于无人倾斜摄影的复杂立交交通分析装置及方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103413444A (zh) * | 2013-08-26 | 2013-11-27 | 深圳市川大智胜科技发展有限公司 | 一种基于无人机高清视频的交通流调查处理方法 |
| CN105321342A (zh) * | 2015-12-07 | 2016-02-10 | 北京航空航天大学 | 一种基于航拍视频的交叉口车辆排队长度检测方法 |
| CN105739523A (zh) * | 2015-12-07 | 2016-07-06 | 北京航空航天大学 | 一种警用车载无人机监测系统及控制方法 |
| CN106548633A (zh) * | 2016-10-20 | 2017-03-29 | 中国科学院深圳先进技术研究院 | 一种路网潮汐交通流可变导向车道控制方法 |
-
2017
- 2017-11-30 CN CN201711239604.1A patent/CN107808524B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103413444A (zh) * | 2013-08-26 | 2013-11-27 | 深圳市川大智胜科技发展有限公司 | 一种基于无人机高清视频的交通流调查处理方法 |
| CN105321342A (zh) * | 2015-12-07 | 2016-02-10 | 北京航空航天大学 | 一种基于航拍视频的交叉口车辆排队长度检测方法 |
| CN105739523A (zh) * | 2015-12-07 | 2016-07-06 | 北京航空航天大学 | 一种警用车载无人机监测系统及控制方法 |
| CN106548633A (zh) * | 2016-10-20 | 2017-03-29 | 中国科学院深圳先进技术研究院 | 一种路网潮汐交通流可变导向车道控制方法 |
Non-Patent Citations (2)
| Title |
|---|
| 冯仰敏: "运动目标图像检测与跟踪问题研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 程德强 等: "一种基于改进光流法的电子稳像算法", 《煤炭学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109003451A (zh) * | 2018-07-04 | 2018-12-14 | 北京航空航天大学 | 一种基于低空无人机的交叉口od矩阵估算方法 |
| CN109919957A (zh) * | 2019-01-08 | 2019-06-21 | 同济大学 | 一种基于动态视觉传感器的角点检测方法 |
| CN109919957B (zh) * | 2019-01-08 | 2020-11-27 | 同济大学 | 一种基于动态视觉传感器的角点检测方法 |
| CN109993134A (zh) * | 2019-04-04 | 2019-07-09 | 中山大学 | 一种基于hog和svm分类器的道路交叉口车辆检测方法 |
| CN111523527A (zh) * | 2020-07-02 | 2020-08-11 | 平安国际智慧城市科技股份有限公司 | 特种运输车辆监控方法、装置、介质及电子设备 |
| CN111523527B (zh) * | 2020-07-02 | 2020-10-27 | 平安国际智慧城市科技股份有限公司 | 特种运输车辆监控方法、装置、介质及电子设备 |
| CN112542050A (zh) * | 2020-12-23 | 2021-03-23 | 重庆市市政设计研究院有限公司 | 基于无人倾斜摄影的复杂立交交通分析装置及方法 |
| CN112542050B (zh) * | 2020-12-23 | 2024-03-22 | 重庆市市政设计研究院有限公司 | 基于无人倾斜摄影的复杂立交交通分析装置及方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107808524B (zh) | 2021-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110175576B (zh) | 一种结合激光点云数据的行驶车辆视觉检测方法 | |
| CN107808524A (zh) | 一种基于无人机的道路交叉口车辆检测方法 | |
| Hadi et al. | Vehicle detection and tracking techniques: a concise review | |
| CN105373135B (zh) | 一种基于机器视觉的飞机入坞引导和机型识别的方法及系统 | |
| WO2017041396A1 (zh) | 一种车道线数据的处理方法、装置、存储介质及设备 | |
| CN108364466B (zh) | 一种基于无人机交通视频的车流量统计方法 | |
| CN110379168B (zh) | 一种基于Mask R-CNN的交通车辆信息获取方法 | |
| Azevedo et al. | Automatic vehicle trajectory extraction by aerial remote sensing | |
| CN102073846B (zh) | 基于航拍图像的交通信息获取方法 | |
| CN103632158B (zh) | 森林防火监控方法及森林防火监控系统 | |
| CN104318263A (zh) | 一种实时高精度人流计数方法 | |
| CN103093198B (zh) | 一种人群密度监测方法及装置 | |
| CN111105495A (zh) | 一种融合视觉语义信息的激光雷达建图方法及系统 | |
| CN103426179B (zh) | 一种基于均值偏移多特征融合的目标跟踪方法及装置 | |
| CN105404857A (zh) | 一种基于红外的夜间智能车前方行人检测方法 | |
| CN109584294A (zh) | 一种基于激光点云的路面点云提取方法和装置 | |
| CN107139666A (zh) | 越障识别系统及方法 | |
| CN110348332A (zh) | 一种交通视频场景下机非人多目标实时轨迹提取方法 | |
| CN112798811A (zh) | 速度测量方法、装置和设备 | |
| CN106228570A (zh) | 一种真值数据确定方法和装置 | |
| CN109740609A (zh) | 一种轨距检测方法及装置 | |
| CN116188417A (zh) | 基于slam和图像处理的裂缝检测及其三维定位方法 | |
| Mijić et al. | Traffic sign detection using YOLOv3 | |
| CN116643291A (zh) | 一种视觉与激光雷达联合剔除动态目标的slam方法 | |
| CN103456029A (zh) | 一种抗相似颜色和光照变化干扰的Mean Shift跟踪方法 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210924 |