一种利用隐马尔可夫地图匹配算法增强手机基站定位精度的方法技术

本发明专利技术属于手机基站定位领域，涉及一种利用隐马尔可夫地图匹配算法增强手机基站定位精度的方法。利用隐马尔可夫地图匹配算法来提高手机基站定位的精度，算法将定位点匹配到候选匹配道路上的概率称作初始匹配概率，将候选匹配道路发生转换的概率称作道路转移概率，算法使用动态编程来快速找到初始匹配概率和道路转移概率乘积最大的路径。考虑到车辆在实际行驶过程中不可能频繁地转换道路，我们在算法计算定位点道路转移概率时引入了约束限制进行改进，从而提高了算法性能。

【技术实现步骤摘要】
一种增强手机基站定位精度的方法
本专利技术属于地图匹配修正
，涉及一种利用隐马尔可夫地图匹配算法增强手机基站定位精度的方法。地图匹配是目前广泛应用于导航系统定位的一种修正方法，是一种以计算机软件技术为基础，通过特定的模型和算法，对导航、定位中的误差进行修正的技术。就是将车辆获取的带有误差的车辆轨迹信息，运用一定方法和算法匹配到数字交通地图上的正确位置。通过使用地图匹配算法，可以较好的弥补定位点存在的精度误差，因此地图匹配技术也是提升导航、定位系统性能的一项重要技术。这种匹配技术前提是假定车辆必须行驶在路上，使用模式识别理论综合分析定位点形成的轨迹和电子地图中的道路网信息，将测得的车辆位置信息和数字交通地图中的相关数据进行对比，并使用算法对定位点进行匹配，最终得出车辆所在路段，并相应找出车辆在该路段上的具体位置，然后将车辆定位点由该道路段之外的位置投影到该道路段上，从而校正了定位点的误差。
技术介绍
当前常用的地图匹配算法有点到点、点到曲线(最短距离法)、曲线到曲线、基于模糊逻辑、模式识别、误差代价函数等。在实际运用中，无论采用哪一种大都是通过比较定位轨迹与侯选道路数据形状的相似性、定位轨迹与候选道路方向的相似性、定位轨迹到侯选道路距离的远近，以及利用路网的拓扑性质来确定车辆行驶道路。点到点和点到曲线的算法实现比较简单，不需利用路网的连通性和历史轨迹，具有实时性的优点，但此方法对定位精度要求较高，并且路段间距离较近时无法完成相关匹配。曲线到曲线的算法则利用了道路网的连通性和历史轨迹，匹配率较高。但这种方法需要知道车辆行驶起始节点，并且要求定位轨迹与道路具有相当的相似性才能实现算法。误差代价函数算法综合利用了各种定位信息和路网信息，匹配准确率高，适合于复杂的道路网络，但比较复杂。通过对手机定位误差的实验分析可知，手机定位点随机分布在路段周围，手机定位误差较GPS大，且定位采样频率较低，因此定位轨迹的形状同真实行驶路线的相似性不明显，而现有算法大都是通过比较定位轨迹与侯选道路数据形状、方向的相似性，定位轨迹到侯选道路距离的远近，以及利用路网的拓扑性质来确定车辆行驶道路。若仅用形状相似性规则判断行驶路段，极易导致误匹配错误。此外，现有算法大都以定位信息精度和定位采样频率较高的GPS定位数据为基础，由于手机定位的定位信息精度和采样频率较低，这些算法均不能直接用于以手机基站定位数据为定位信息的匹配算法。
技术实现思路
隐马尔可夫模型能够有效地平滑整合定位误差数据和路径约束。利用隐马尔可夫地图匹配算法来提高手机基站定位的精度。定位点匹配到候选匹配道路上的概率称作初始匹配概率，候选匹配道路发生转换的概率称作道路转移概率，使用动态编程来快速找到初始匹配概率和道路转移概率乘积最大的路径。考虑到车辆在实际行驶过程中不可能频繁地转换道路，计算定位点道路转移概率时引入了约束限制进行改进，从而提高了性能。本专利技术用模拟数据检验时表明，即使采样周期长达30秒，误差也仅为0.11％，在较长的采样周期中，该方法对测量值噪声误差的鲁棒性高达50米的标准偏差。在隐马尔可夫地图匹配算法中，路段之间的转换是通过路网的连通来管理。如图1所示，HMM离散状态是Nr个道路段，用ri表示，其中i＝1,…,Nr，不同道路段之间的转换发生在道路相接的交叉路口。zt表示某一时刻测量的坐标值，t＝1时，临近z1的三条道路用三个黑点在第一列中表示。t＝2时，分别从这3条路上最近的点指向靠近z2的两条道路。t＝3时做同样的工作。每个t时刻选取的道路段组成一个道路网格，目标是在道路网格中找到最可能的路径。本专利技术的方法利用隐马尔可夫地图匹配算法来提高手机基站定位的精度，将定位点匹配到候选匹配道路上的概率称作初始匹配概率，将候选匹配道路发生转换的概率称作道路转移概率，使用动态编程来快速找到初始匹配概率和道路转移概率乘积最大的路径。考虑到车辆在实际行驶过程中不可能频繁地转换道路，在计算定位点道路转移概率时引入了约束限制进行改进，从而提高性能。附图说明图1是zt候选匹配道路段ri及道路段间转换示意图。图2是测量值zt在候选匹配道路段ri上的候选匹配点示意图。图3是本专利技术的示例图。图4是Δs修正示例图。具体实施方式以下结合技术方案和附图详细叙述本专利技术的具体实施例。1、初始匹配概率点zt与每一条临近的道路ri存在一个匹配概率p(zt|ri)。xt,i表示点zt匹配到道路ri上的点。由图2可以发现，点z1匹配到候选匹配道路r1，r2，r3上的候选匹配点为x1,1，x1,2，x1,3。坐标点和候选匹配点的大圆距离为||zt-xt,i||greatcircle。相对于t+1时的点zt+1来说，zt+1匹配到道路rj上的点为xt+1,j。点xt,i和点xt+1,j之间车辆行驶的距离称为“路径距离”，记做||Xt,i-xt+1,j||route。对于正确匹配来说，坐标点和匹配点之间的这种误差是由定位误差形成的。算法模拟的定位误差为零均值的高斯误差，这意味着：σz在式中代表定位测量值的标准偏差。初始匹配概率πi，i＝1,…,Nr，其表示在地图匹配开始运行前，从所有路段中找出车辆所在路段的概率。当部分HMM的构想值均匀分布于πi上时，假设没有测量值被采用，则算法在第一个测量值上开始，这时有πi＝p(zt|ri)，即使用第一测量值z1。2、道路转移概率，将候选匹配道路发生转换同每个测量值zt一样，下一个测量zt+1也有一个可能匹配道路的列表。转移概率是指在这两个时间(t到t+1)内，车辆在候选匹配道路之间的转移概率。对于测量值zt和其候选匹配道路段ri来说，在ri上最可能匹配点是xt,i。同样，测量值zt+1在其候选匹配道路段rj上的候选匹配点是xt+1,j。我们比较测量点之间的“路径距离”||xt,i-xt+1,j||route和大圆距离||zt-zt+1||greatcircle，正确匹配时这两个值相当接近，通过比较这两个值就可以判断道路是否发生转换。举例说明，在图3中有三个道路段r1，r2和r3，两个测量点zt和zt+1。测量点zt有xt,1和xt,3两个候选匹配点，测量点zt+1的候选匹配点为xt+1,2。||zt-zt+1||greatcircle是两个测量点的大圆距离，路径距离有两个，分别为||xt,1-xt+1,2||route和||xt,3-xt+1,2||route。与不正确匹配相比，正确匹配时的路径距离更接近大圆距离。该算法通过对正确匹配的大圆距离和路径距离差的绝对值进行分析，得出指数概率分布：式中：dt＝|||zt-zt+1||greatcircle-||xt,i*-xt+1,j*||route|其中，i*和j*指的是车辆的真实行驶道路段。3、优化考虑到正常行驶时，驾驶人员在到达目的地之前会尽量沿同一条道路上行驶，不会轻意改变行驶道路，更不会频繁地转换道路。因此频繁转换的候选匹配道路可能是由于定位误差产生的错误匹配。为此，我们引入一个参数Δs来对这一行为进行约束。当定位点的候选匹配道路发生改变，候选匹配点离开原行1驶的道路，其路径距离就远离大圆距离50米进行惩罚，就是当路径距离大于大圆距离时，路径距离加上50米，路径距离小于大圆距离时，路径距离减去50米。不正本文档来自技高网...

【技术保护点】
一种利用隐马尔可夫地图匹配算法增强手机基站定位精度的方法，其特征在于如下步骤，(1)初始匹配概率点zt与每一条临近的道路ri存在一个匹配概率p(zt|ri)；xt,i表示点zt匹配到道路ri上的点；点z1匹配到候选匹配道路r1，r2，r3上的候选匹配点为x1,1，x1,2，x1,3；坐标点和候选匹配点的大圆距离为‖zt‑xt,i‖great circle；相对于t+1时的点zt+1来说，zt+1匹配到道路rj上的点为xt+1,j；点xt,i和点xt+1,j之间车辆行驶的距离称为“路径距离”，记做‖xt,i‑xt+1,j‖route；模拟的定位误差为零均值的高斯误差，即：p(zt|ri)=12πσze-0.5(||zt-xt,i||great circle/σz)2---(1)]]>σz在式中代表定位测量值的标准偏差；初始匹配概率πi，i＝1…Nr，表示在地图匹配开始运行前，从所有路段中找出车辆所在路段的概率；当部分HMM的构想值均匀分布于πi上时，假设没有测量值被采用，则在第一个测量值上开始，有πi＝p(zt|ri)，使用第一测量值z1；(2)道路转移概率转移概率是指在这两个时间(t到t+1)内，车辆在候选匹配道路之间的转移概率；对于测量值zt和其候选匹配道路段ri来说，在ri上最可能匹配点是xt,i；测量值zt+1在其候选匹配道路段rj上的候选匹配点是xt+1,j；测量点之间的“路径距离”||xt,i‑xt+1,j||route和大圆距离||zt‑zt+1||great circle，正确匹配时这两个值相当接近，通过比较这两个值判断道路是否发生转换；通过对正确匹配的大圆距离和路径距离差的绝对值进行分析，得出指数概率分布：p(dt)=1βe-dt/β---(2)]]>式中：dt=|||zt-zt+1||great circle-||xt,i*-xt+1,j*||route|]]>β=1ln(2)mediant(|||zt-zt+1||great circle-||xt,i*-xt+1,j*||route|)]]>其中，i*和j*指的是车辆的真实行驶道路段；(3)优化引入一个参数Δs来对这一行为进行约束；当定位点的候选匹配道路发生改变，候选匹配点离开原行1驶的道路，其路径距离就远离大圆距离50米进行惩罚，就是当路径距离大于大圆距离时，路径距离加上50米，路径距离小于大圆距离时，路径距离减去50米；这样不正确匹配点的路径距离就会远离大圆距离，更容易找出正确匹配点；(4)输出最佳匹配路径通过用公式(1)中的匹配概率和公式(2)中的转换概率，使用维特比算法在HMM点阵中来计算出最佳匹配路径；维特比算法使用动态编程来快速找到匹配概率和转移概率乘积最大的路径，从而推断出坐标点的正确匹配道路段；如果对于从t时刻的点zt匹配到道路ri上的点xt,i到达t+1时刻的点zt+1匹配到rj道路上的点xt+1,j发生道路转换时，要考虑点zt是通过哪条道路段rx到达点zt+1，不能只单独考虑道路匹配概率p(zt\ri)和p(zt+1\rj)，要考虑两点之间发生的道路转换概率，即匹配概率乘以相应的道路转换概率，从中选出最大值即为此时点的最佳匹配道路段，并记下该点在此道路段上的最佳匹配值是由哪条道路段转换而来，以此类推；算法对定位点的处理区分为对第一个定位点的处理和对普通定位点的处理；具体如下：对于一个点zt，它匹配到每条候选道路段上的状态概率为p(zt|ri)；如果该点是车辆定位点的第一个点，则不需要对此概率进行修正，如果该点不是该车辆定位点的第一个点，仅为分组后该组点的第一个点，则要对状态概率p(zt|ri)进行修正；需要将状态概率乘以上时刻点的匹配道路段到本时刻点的匹配路段的转移概率，即p(zt|ri)＝p(zt|ri)×p(dt‑1)；对于其余点的处理，只需考虑每个点zt匹配到每条可行弧ri的状态概率；设p(zt|ri)表示第t个点zt匹配到弧ri的最大概率，设从弧ri转移到弧rj的转移概率为p(ri|rj),第t+1个点定位到弧rj的状态概率为p(zt+1|rj)，则其转移方程为p(zt+1\rj)＝max(p(zt+1\rj),p(zt\ri)×p(zt+1\rj)×p(ri\rj))；如果p(zt+1|rj)取值为p(zt\ri)×p(zt+1\rj)×p(ri\rj)，则需要记录转换路径，在输出结果时以此反推最佳匹配道路段；当计算完全部点的状态概率后，确定最后一个点的最佳匹配道路段，并以此道路段反向推出从第一个点到最后一个点的车辆实际行驶轨迹的道路段；整个行驶轨迹计算完毕后，最终输出车辆的最佳匹配路径。...

【技术特征摘要】
1.一种增强手机基站定位精度的方法，其特征在于如下步骤，(1)初始匹配概率点zt与每一条临近的道路ri存在一个匹配概率p(zt|ri)；xt,i表示点zt匹配到道路ri上的点；点z1匹配到候选匹配道路r1，r2，r3上的候选匹配点为x1,1，x1,2，x1,3；坐标点和候选匹配点的大圆距离为‖zt-xt,i‖greatcircle；相对于t+1时的点zt+1来说，zt+1匹配到道路rj上的点为xt+1,j；点xt,i和点xt+1,j之间车辆行驶的距离称为“路径距离”，记做‖xt,i-xt+1,j‖route；模拟的定位误差为零均值的高斯误差，即：σz在式中代表定位测量值的标准偏差；初始匹配概率πi，i＝1,…,Nr，表示在地图匹配开始运行前，从所有路段中找出车辆所在路段的概率；当部分HMM的构想值均匀分布于πi上时，假设没有测量值被采用，则在第一个测量值上开始，有πi＝p(zt|ri)，使用第一测量值z1；(2)道路转移概率转移概率是指在这两个时间t到t+1内，车辆在候选匹配道路之间的转移概率；对于测量值zt和其候选匹配道路段ri来说，在ri上最可能匹配点是xt,i；测量值zt+1在其候选匹配道路段rj上的候选匹配点是xt+1,j；测量点之间的“路径距离”||xt,i-xt+1,j||route和大圆距离||zt-zt+1||greatcircle，正确匹配时这两个值相当接近，通过比较这两个值判断道路是否发生转换；通过对正确匹配的大圆距离和路径距离差的绝对值进行分析，得出指数概率分布：式中：其中，i*和j*指的是车辆的真实行驶道路段；(3)优化引入一个参数Δs来对这一行为进行约束；当定位点的候选匹配道路发生改变，候选匹配点离...

【专利技术属性】
技术研发人员：申彦明，张鹏飞，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21