一种有立交桥遮挡的车辆轨迹匹配方法
【技术领域】
[0001] 本发明涉及一种有立交桥遮挡的车辆轨迹匹配方法，属于计算机视觉与数字图像 处理技术领域。
【背景技术】
[0002] 在城市交通问题的研究中，道路车流量数据、车辆时空轨迹数据对于探索车辆在 城市道路中的运行规律，解决城市交通问题具有重要意义。随着无人机在民用市场的开放， 利用无人机监控城市道路交通运行状况具有很大优势。无人机高空拍摄覆盖范围广，直观 呈现交通运行状态，且采集并保存的道路交通运行视频对于提取车辆时空轨迹数据具有重 大价值。但在城市路网中，立交桥的存在阻断了车辆在航拍视频中跟踪的连续性。车辆在驶 入立交桥下后消失在画面中，此时目标丢失，当驶出立交桥后，车辆又重新出现在画面中。 如何解决有立交桥遮挡的车辆轨迹匹配算法，对于精确统计车流量、保持车辆跟踪的连续 性、提取连续的车辆时空轨迹数据至关重要。
【发明内容】
[0003] 本发明所要解决的技术问题是提供一种能够解决有立交桥遮挡情况下的道路航 拍视频车辆轨迹匹配方法。
[0004] 一种针对航拍视频的有立交桥遮挡的车辆轨迹匹配方法，包括以下步骤：
[0005] (1)检测立交桥边界，确定立交桥区域，利用Sobel边缘检测结合霍夫变换确定立 交桥边缘；
[0006] (2)对跟踪车辆提取梯度方向直方图（H0G)特征，作为该车辆的模板特征；
[0007] (3)根据当前跟踪车辆的行驶速度，估算出其行驶出立交桥的时间区间；
[0008] (4)在立交桥的另一侧检测对应的跟踪车辆。在立交桥的另一侧设定搜索区域，当 视频处于步骤（3)估算出的车辆驶出时间区间时，在搜索区域内利用滑动窗口法结合模板 匹配技术检测对应的跟踪车辆，将检测到的车辆坐标与之前驶入立交桥的坐标串联起来， 实现车辆轨迹的匹配。
[0009] 步骤（1)中所述的检测立交桥边界的具体步骤如下：
[0010] la)根据车辆行驶方向确定Sobel边缘检测卷积模板。由于立交桥边缘方向与车 辆行驶方向垂直，故若在航拍视频中车辆行驶方向为水平方向时（水平方向是指与航拍视 频图像水平边缘平行的车辆行驶方向），Sobel算子的卷积模板Gx为：
[0012] 若航拍视频中车辆行驶方向为竖直方向时，Sobel算子的卷积模板Gy为：
[0014] lb)利用Sobel算子检测立交桥边缘，确定立交桥区域。对视频的一帧f (x, y)利 用步骤la)确定的卷积模板进行卷积运算，设定边缘阈值T得到边缘灰度图G :
[0015] G = Gx*f or G = Gy*f
[0016] lc)将边缘灰度图利用Otsu阈值分割法转化为二值图像，采用霍夫变换检测直 线，计算检测到直线的长度^，若^大于图像宽度的1/2,则该直线处于立交桥的边缘。将 所有满足条件的直线延长，最终确定立交桥区域的边缘。
[0017] 步骤（2)中所述的提取梯度方向直方图（H0G)特征的具体步骤如下：
[0018] 2a)对待跟踪车辆图像的颜色空间进行归一化，并采用中心对称的一维点模版 [-1，0, 1]进行卷积运算，得到每个像素点的X、Y方向的梯度幅值Gx、Gy，则该像素点的梯度 大小
，梯度方向为D = arctan(Gy/Gx);
[0019] 2b)将待跟踪车辆图像整体划分成相同大小的单元格（Cell)，单元格数量依单元 格的尺寸而定；在每一个单元格内按照梯度方向统计一维梯度方向直方图向量，方法为：
[0020] 每个单元格内的每个像素均有梯度大小值和梯度方向值与之对应，将单元格内每 个像素的梯度方向0°~180°按照20°的间隔划分为9个区间，即：0°~20°，20°~ 40°，依次类推；如：一个单元格中的某个像素的梯度方向是10°，那么它属于0°~20° 区间，在该区间上加上像素的梯度值，统计单元格内的所有像素，从而构成梯度方向直方图 向量，是9维的，9个区间。随后将像素的梯度方向D映射到9个区间上，并以梯度幅值作 为权值，构成单元格的梯度方向直方图向量H1= [h u h2,. . . hs，h9]，其中匕为每个区间的权 值；
[0021] 2c)将待跟踪车辆图像整体划分成包含多个相邻单元格的相同大小的矩形区块 (Block)，区块间可以重叠。将区块内所有单元格的一维梯度方向直方图向量私串联起来 构成串联向量[氏，氏，...，札]，η为区块内单元格的数量，采用L2-Norm方式对串联向量 [4? · · ·，Hn]进行归一化处理构成区块特征向量Fi= [f ;，f2, · · ·，fj，其中m = ηΧ9·
[0022] 2d)将待跟踪车辆图像中所有的区块特征向量Fi串联起来构成梯度方向直方图特 征向量，即所跟踪车辆的模板特征向量Vf;
[0023] 步骤（3)中所述的估算车辆行驶出立交桥的时间区间的具体方法如下：
[0024] 由车辆的跟踪数据可以实时估算出车辆在驶入桥下时的行驶速度v(piXel S/ frame)与加速度a(pixels/frame2)。假设车辆在立交桥下做勾加速运动，即以初始速度v、 加速度a驶过桥下，根据步骤（1)得到的立交桥边缘直线，可以计算出桥宽为s，那么依据s =vt+0. 5at2可计算出驶出时刻t。将区间[t-fps, t+fps]作为车辆行驶出立交桥的时间 区间，其中fps(frames/s)为视频的帧率。
[0025] 步骤⑷中所述的检测对应跟踪车辆的具体步骤如下：
[0026] 4a)在立交桥车辆驶出的一侧设定搜索区域；
[0027] 4b)利用滑动窗口法遍历搜索区域，利用步骤（2)所述方法提取每个窗口的梯度 方向直方图（H0G)特征向量
[0028] 4c)计算模板特征向量Vf与滑动窗口特征向量#的相似度d :
[0030] 其中
：表示对向量Vf与V；'的差采用L2-Norm的方式进行归一化，而 I I Vf I I表示对向量v f本身进行L2-Norm归一化。
[0031] 相似度最大的窗口即为所跟踪的车辆，将检测到的车辆坐标与驶入立交桥时的坐 标串联起来实现跟踪车辆轨迹的匹配。
[0032] 本发明的优点在于：
[0033] (1)本发明针对道路航拍视频中出现的立交桥遮挡影响车辆跟踪连续性的问题， 通过在立交桥的另一侧检测对应的车辆，最后将坐标串联起来，实现了车辆轨迹的精确匹 配，解决了立交桥的遮挡，保证了车辆跟踪的连续性；
[0034] (2)本发明采用基于梯度方向直方图特征（H0G)的模板匹配技术，通过提取滑动 窗口内图像的H0G特征向量，与对应的车辆模板H0G特征向量进行匹配比对实现了对特定 车辆的精确定位，提高了轨迹匹配的准确性与稳定性；
[0035] (3)本发明通过立交桥区域的检测，明确了车辆的驶出区域，进而缩小了驶出立交 桥车辆的检测范围，降低了检测时间，同时得到的立交桥宽度数据可以用来估算车辆驶出 立交桥的时间；
[0036] (4)本发明通过估算车辆驶出立交桥的时间，避免了相同车型的干扰，保证了车辆 轨迹的正确匹配，提高了算法整体的鲁棒性。
【附图说明】
[0037] 图1是本发明实施例提供的有立交桥遮挡的车辆轨迹匹配方法流程图；
[0038] 图2是本发明实施例中检测立交桥边界的示意图；
[0039] 图3是本发明实施例中提取车辆梯度方向直方图特征示意图；
[0040] 图4是本发明实施例中在搜索区域内进行特征模板匹配的示意图；
【具体实施方式】
[0041] 下面将结合附图和实施例对本发明作进一步的详细说明。
[0042] 图1展示了本发明提供的一种针对道路航拍视频的有立交桥遮挡的车辆轨迹匹 配算法的执行流程图。在本发明中，车辆轨迹匹配指的是把同一车辆在驶入立交桥前的坐 标和驶出立交桥的坐标连接起来，保证车辆轨迹跟踪的连续性，同时也确保车流量统计的 正确。其中的关键技术在于如何能够在立交桥驶出的一侧准确地匹配到之前的车辆。本 发明通过检测立交桥边界，确定立交桥区域，针对性地把图像区域分为两部分。在车辆驶 出区域利用提取的H0G特征模板进行车辆模板匹配检测，从而将跟踪车辆对应起来。具体 地，在本实施中，选用一段存在立交桥遮挡的道路航拍视频，拍摄高度为150m，视频分辨率 为1920*1080像素。本实施例提供的一种针对有立交桥遮挡的车辆轨迹匹配算法的具体实 施步骤如下：
[0043] 步骤1 :检测立交桥边界，确定立交桥区域
[0044] 通过检测立交桥的边缘确定立交桥区域、立交桥的宽度，能够大体确定车辆驶出 立交桥的位置和时间。检测立交桥边界的具体步骤如下：
[0045] la)根据车辆行驶方向确定Sobel边缘检测卷积模板。由于立交桥边缘方向与车 辆行驶方向垂直，故若车辆行驶方向为水平方向时，Sobel算子的卷积模板6：(为：
[0047] 若车辆行驶方向为竖直方向时，Sobel算子的卷积模板Gy为：
[0049] lb)利用Sobel算子检测立交桥边缘，确定立交桥区域。对原图像f (X，y)利用步 骤la)确定的卷积模板进行卷积运算，设定边缘阈值T得到边缘灰度图G :
[0050] G = Gx*f or G = Gx*f
[0051] 本实施例选用的道路航拍视频中，车辆沿水平方向行驶，立交桥的边缘为竖直方 向，故选取Gx卷积模板过滤掉水平边缘，得到边缘灰度图像。
[0052] lc)将边缘灰度图利用Otsu阈值分割法转化为二值图像（图2 (a)所示