一种基于低采样率浮动车数据的全局投票地图匹配方法
【专利摘要】一种基于低采样率浮动车数据的全局投票地图匹配方法，该方法在浮动车GPS轨迹数据的基础上，考虑道路网络的拓扑结构和不同距离的空间位置相邻的GPS轨迹数据对地图匹配过程的影响。定义了新的地图匹配函数指标，综合考虑道路网络的几何学特性和拓扑结构对匹配的影响，得出作为初始地图匹配结果的静态匹配矩阵(SMM)；在此基础上，定义反映节点之间距离影响的距离加权函数对SMM进行修正，得到动态匹配矩阵(DMM)，最后，对DMM采用全局投票的方法，找出最优轨迹作为地图匹配的结果。本发明能在低采样率的情况下，以较低的计算时间复杂度实现较好的地图匹配效果。
【专利说明】一种基于低采样率浮动车数据的全局投票地图匹配方法
【技术领域】
[0001]本发明涉及交通领域，特别是指一种基于低采样率浮动车数据的全局投票地图匹配方法。
【背景技术】
[0002]浮动车一般是指安装了车载GPS定位装置并行驶在城市主干道上的公交汽车或出租车。浮动车技术是近年来兴起的一种动态交通信息检测的技术，是近年来智能交通系统中所采用的获取道路交通信息的先进技术手段之一。浮动车可以定期将车辆编号、时刻、方向、经纬度坐标等数据传送到调度中心，经过信息处理，可以方便地获得整个城市路网的实时交通信息。地图匹配方法是浮动车数据处理的关键技术之一，它可以最大程度地校正GPS卫星定位所产生的误差，从而矫正车辆轨迹偏离道路的现象，并进一步将GPS数据转化为道路的交通信息。
[0003]传统的地图匹配方法，GPS采样率往往都很高，匹配的正确率才得以保证。但实际在浮动车技术中，为了节约能源以及浮动车本身特性，返回给调度中心的数据往往都是低采样率的。在低采样率的情况下(比如2分钟)，在车辆的速度仅为30km/h，两个采样点之间的距离就达到了 1000米，两个GPS点之间的关联信息会大部分丢失。因此，传统的地图匹配方法不能很好运用在基于低采样率的浮动车数据的匹配过程。以此为背景，为了提高匹配方法在低采样率GPS数据输入情况下的正确匹配率，本文提出的全局投票方法考虑道路网络的拓扑结构以及GPS点之间的相互影响。
[0004]考虑几何学特性的地图匹配方法，根据匹配过程，这类方法又可以分为点对点的匹配，点对曲线的匹配，曲线对曲线的匹配。由于缺乏对整体道路网络拓扑结构的考虑，在复杂的道路环境和低采样率的浮动车数据情况下，这些方法很容易出现匹配误差。同时这些匹配方法不考虑全局GPS点之间的相互影响，在某个GPS点的匹配过程中，只利用该点本身的信息，或者仅仅考虑其相邻节点对其匹配过程的影响，但针对节点间关联信息大量丢失的低采样率浮动车数据时，很难得到正确的匹配结果。
【发明内容】
[0005]为了克服现有的全区投票地图匹配方法的正确率较低的不足，本发明提供一种一种基于低采样率浮动车数据的全局投票地图匹配方法，在道路网络的模型基础上，首先计算当前GPS轨迹点可能匹配的候选集合，考虑集合中候选点的近距概率和连通概率，得出静态分析矩阵。再考虑节点之间距离的不同导致影响程度的不同，得出加权匹配矩阵，在此基础上，各个GPS轨迹点之间相互投票，完成匹配方法。
[0006]为了解决上述技术问题采用的技术方案如下:
[0007]—种基于低采样率浮动车数据的全局投票地图匹配方法，所述匹配方法包括以下步骤:
[0008]步骤一:构建有向道路网络G (V，E)，其中V为道路的交叉点、起点或终点，E为被两个相邻交叉口分隔出的路段，定义一条路径为:在道路网络中，选取两个节点Vi, Vj,找到一条互相连接的路段集合S1 — S2 — s3...— sn,并且S1.start = Vi, sn.end = Vj,即路段S1的起点为Vi,路段Sn的终点为Vj ;在道路网络的基础上，处理浮动车GPS数据，得到浮动车GPS轨迹数据T ；
[0009]步骤二:针对浮动车GPS轨迹数据T % — p2 — p3...— Pn中的轨迹点Pi,
选取在半径r范围内的所有路段为候选路段<，k表示点Pi的第k个候选路段；
通过投影得到对应的候选节点:如果轨迹点Pi在路段的范围内存在垂点，则选取该垂点作为轨迹点Pi的候选节点，记为cf，并且选取垂线的长度为该轨迹点与该路段的距离；否则，选取该路段离轨迹点较近的起点或者终点为候选节点，选取轨迹点与该路段起点或者终点的连线长度为该轨迹点与该路段的距离；根据路段投影过程，得出每个GPS轨迹点的候选点集合；
[0010]步骤三:静态分析过程，定义匹配函数指标+4Cf)，其中
c\力点Pg的第j个候选节点，Cki为点Pi的第k个候选节点，i≥2 ;对每个候选点进行匹配函数计算，得出静态分析矩阵；
[0011]步骤四:加权得分过程，定义与两个轨迹点Pi，Pj之间的距离有关的函数的距离加权函数WJ = f (dist (Pi, Pj)), (Iist(PpPj)为轨迹点Pi, Pj之间的欧式距离，I≤i≤η,I≤j≤n,i古j，距离加权函数代表节点之间距离越远，节点之间影响越弱，对每一个GPS点，经过全局计算，得出加权分析矩阵；
[0012]步骤五:全局投票，对每个轨迹点Pi，在其对应的加权分析矩阵Gi中，寻找一条最优的通过每个轨迹点的候选节点的路径，定义为单点投票过程；某个候选节点被选取，对这个候选节点进行记录加一操作，表示一次投票选中过程，针对全部候选节点集合，依次重复上述操作，会得到全部候选节点的得票情况；最后，对每个轨迹点Pi,在其候选节点集合中选择得票数最多的候选点作为全局候选节点结果，即得到了全局匹配结果。
[0013]进一步，所述步骤一中，一个GPS点P的记录包括经度、纬度、时间，所有的GPS点的集合定义为GPS记录L，GPS轨迹T定义为一组采样间隔为Λ t的GPS点的序列；一个路段s为一条道路被两个相邻交叉口分隔出的一段路，每个路段包含起始点、结束点、路段长度。
[0014]再进一步,所述步骤三中，根据步骤二得到的候选节点集合，定义候选节点的近距概率MP为候选卢c是正确的匹配结果的概率，I≤k≤n，即:
【权利要求】
1.一种基于低采样率浮动车数据的全局投票地图匹配方法，其特征在于:所述匹配方法包括以下步骤: 步骤一:构建有向道路网络G(V，E)，其中V为道路的交叉点、起点或终点，E为被两个相邻交叉口分隔出的路段，定义一条路径为:在道路网络中，选取两个节点Vi, Vj，找到一条互相连接的路段集合S1 — S2 — s3...— sn,并且S1.start = Vi, sn.end = Vj,即路段S1的起点为Vi，路段Sn的终点为Vj ;在道路网络的基础上，处理浮动车GPS数据，得到浮动车GPS轨迹数据T ； 步骤二:针对浮动车GPS轨迹数据T -.P1 — p2 — p3...— pn中的轨迹点Pi, I≤i≤η,选取在半径r范围内的所有路段为候选路段<，k表示点Pi的第k个候选路段；通过投影得到对应的候选节点:如果轨迹点Pi在路段4的范围内存在垂点，则选取该垂点作为轨迹点Pi的候选节点，记为Cf，并且选取垂线的长度为该轨迹点与该路段的距离；否则，选取该路段离轨迹点较近的起点或者终点为候选节点，选取轨迹点与该路段起点或者终点的连线长度为该轨迹点与该路段的距离；根据路段投影过程，得出每个GPS轨迹点的候选点集合；步骤三:静态分析过程，定义匹配函数指标巧匕一 =MPic1;:)+ CP(cU ,其中&为点Pm的第j个候选节点，cf为点Pi的第k个候选节点，i 2 ;对每个候选点进行匹配函数计算，得出静态分析矩阵；步骤四:加权得分过程，定义与两个轨迹点Pi，Pj之间的距离有关的函数的距离加权函数沙V' =.Rdisi(J)卜Pi)), dist (Pi, Pj)为轨迹点Pi, P」之间的欧式距离，I≤i≤η, I≤j≤n，i幸j，距离加权函数代表节点之间距离越远，节点之间影响越弱，对每一个GPS点，经过全局计算，得出加权分析矩阵； 步骤五:全局投票，对每个轨迹Api，在其对应的加权分析矩阵Gi中，寻找一条最优的通过每个轨迹点的候选节点的路径，定义为单点投票过程；某个候选节点被选取，对这个候选节点进行记录加一操作，表示一次投票选中过程，针对全部候选节点集合，依次重复上述操作，会得到全部候选节点的得票情况；最后，对每个轨迹点Pi，在其候选节点集合中选择得票数最多的候选点作为全局候选节点结果，即得到了全局匹配结果。
2.如权利要求1所述的一种基于低采样率浮动车数据的全局投票地图匹配方法，其特征在于:所述步骤一中，一个GPS点P的记录包括经度、纬度、时间，所有的GPS点的集合定义为GPS记录L，GPS轨迹T定义为一组采样间隔为Λ t的GPS点的序列；一个路段s为一条道路被两个相邻交叉口分隔出的一段路，每个路段包含起始点、结束点、路段长度。
3.如权利要求1或2所述的一种基于低采样率浮动车数据的全局投票地图匹配方法，其特征在于:所述步骤三中，根据步骤二得到的候选节点集合，定义候选节点的近距概率MP为候选点cf是正确的匹配结果的概率 ，即:
其中为该候选卢I与轨迹点Pi的欧式距离；考虑道路网络的拓扑结构，定义两个候选节点CU、cf之间的连通概率CP为:
其中，l1-1-1为两个相邻轨迹点P1-1和轨迹点Pi之间的欧氏距离，s(i_m，k)为通过候选节点cU和cf的最短路径长度； 选取相邻轨迹点Pm和Pi相应的候选节点和Cf，经过匹配函数计算，得出静态分析矩阵M，针对每条路径cη≥i≥2，计算出匹配函数的值，得出静态分析矩阵M:
其中Mi为路径cL经过匹配函数计算，得到的矩阵:
其中，1≤j≤m, 1≤k≤n, m, η分别为。卜” Ci的候选节点总个数。
4.如权利要求3所述的一种基于低采样率浮动车数据的全局投票地图匹配方法，其特征在于:所述步骤四中，针对每个轨迹点Pi,定义一个η-l维距离影响矩阵W1:
其中，距离加权函数Wij采用
,其中δ为阈值，接下 来，加权分析矩阵Gi, 1 ≤ i≤ η:
【文档编号】G01C21/30GK104197945SQ201410427108
【公开日】2014年12月10日 申请日期:2014年8月27日 优先权日:2014年8月27日
【发明者】杨旭华, 赵久强, 汪向飞 申请人:浙江工业大学