CN108426582A - 行人室内三维地图匹配方法 - Google Patents
行人室内三维地图匹配方法 Download PDFInfo
- Publication number
- CN108426582A CN108426582A CN201810176554.5A CN201810176554A CN108426582A CN 108426582 A CN108426582 A CN 108426582A CN 201810176554 A CN201810176554 A CN 201810176554A CN 108426582 A CN108426582 A CN 108426582A
- Authority
- CN
- China
- Prior art keywords
- point
- pedestrian
- height
- state point
- indicate
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/28—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network with correlation of data from several navigational instruments
- G01C21/30—Map- or contour-matching
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/26—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for navigation in a road network
- G01C21/34—Route searching; Route guidance
- G01C21/3446—Details of route searching algorithms, e.g. Dijkstra, A*, arc-flags, using precalculated routes
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了行人室内三维地图匹配方法,属于室内行人定位技术领域。利用捆绑于行人脚面上的MEMS‑INS传感器采集行人室内的运动信息并解算行人的速度、位置与航向,分析室内结构并创建状态点,根据导航输出位置信息与状态点位置信息建立条件随机场模型。以定长行走距离提取室内行人的水平二维位置信息,零速时刻提取行人高度信息,获取CRF模型的观测点,并分别记录室内行人的二维位置与高度信息的采样时刻。建立室内电子地图,并根据室内结构信息创建状态点,存储状态点坐标。本发明能够实现行人三维位置定位,算法精度高;采用二维位置与高度信息分开匹配的方法,简化算法复杂性;采用临近时刻对二维位置与高度信息进行融合,提高了地图匹配的正确率。
Description
技术领域
本发明属于室内行人定位技术领域,涉及一种室内环境已知条件下基于微机电***(Micro-Electro-Mechanical System,MEMS)技术的行人惯性导航***(InertialNavigation System,INS)同步定位、匹配算法结构设计与位置信息融合的三维地图匹配方法。
背景技术
近些年,随着物流、智慧病房、新概念超市的出现与发展,室内导航受到了学术界及工程界的广泛关注。MEMS-INS的集成化与小型化促使其成为导航领域的主导技术。
然而,惯性器件误差会随时间积累,如果不对惯性器件误差进行有效的修正,最终会导致行人的位置轨迹发散。为了解决惯性导航位置误差的问题,不同的地图匹配算法被提出。例如粒子滤波算法、基于主航向的地图匹配算法、基于隐形马尔科夫的地图匹配算法等。但这些算法更多的关注行人的二维信息,行人室内行走不仅包括二维平面坐标的移动,还可以通过楼梯进行高度上的改变。因此,行人室内三维地图匹配方法具有重要的理论意义及应用价值。
Y.M′eneroux等人提出通过平均Hausdorff距离和面积差两种常见的度量方法来解决匹配精度与网络质量指数之间的关系问题,同时还提供了参考网络的影响上限;而Xiao Z等人提出采用条件随机场(conditional random fields,CRF)算法对导航输出位置与室内状态点进行数学建模来解决惯性导航位置误差发散的问题。然而,前者作为一种室外地图匹配所考虑的更多是为车辆进行路径匹配,行人室内行走比较随意,路径概念相对较弱,因此,该方法不适用于室内环境;而后者虽然考虑了室内的特殊环境,以状态点辅助地图匹配,但忽视了室内的高度信息,只是对单一的楼层进行实验研究,从而导致了只能呈现行人二维的位置信息。
为了有效实现行人在室内环境条件下三维地图匹配,将状态点与传感器采集的导航输出位置信息进行条件随机场模型进行联立,实现行人室内三维定位。通过三维位置定位,提高算法精度与地图匹配的正确率,并简化算法的复杂性是非常重要的。
发明内容
针对佩戴MEMS-INS传感器的行人在室内环境条件下导航位置导致的发散问题,本发明提供了一种面向室内行人导航的三维地图匹配方法。利用捆绑于行人脚面上的MEMS-INS传感器采集行人室内的运动信息并解算行人的速度、位置与航向,分析室内结构并创建状态点,根据导航输出位置信息与状态点位置信息建立条件随机场模型,实现行人室内三维定位。
为实现上述技术目的,本发明采用的技术方案为行人室内三维地图匹配方法,包括以下步骤:
步骤1:数据采集,初步解算室内行人的三维位置和航向。
步骤1.1,行人通过佩戴MEMS-INS传感器,采集行人运动数据,行人运动数据包括:三轴加速度数据和三轴陀螺数据。
步骤1.2,应用捷联解算算法对采集的行人运动数据进行三维位置和航向信息进行求解。
步骤2:提取条件随机场模型的观测点。
步骤2.1,以定长行走距离提取室内行人的水平二维位置信息,零速时刻提取行人高度信息,获取CRF模型的观测点,并分别记录室内行人的二维位置与高度信息的采样时刻。
步骤3:建立室内电子地图,并根据室内结构信息创建状态点,存储状态点坐标。
步骤3.1,利用已知的室内地图创建室内电子地图,并将室内电子地图存储于导航计算机中。
步骤3.2,求取电子地图中数值最小与最大的坐标点作为状态点的取值范围,然后以等间距的状态点覆盖整个室内范围,存储状态点信息。
步骤3.3,以台阶高度为标准,添加高度信息上的状态点。
步骤4,基于条件随机场算法的二维位置地图匹配算法。
步骤4.1,根据二维位置观测点坐标与各个状态点坐标之间关系建立特征方程;
步骤4.2,根据观测点的方位信息与与之对应时刻状态点之间的方位角建立特征方程。
步骤4.3,建立基于条件随机场的二维地图匹配数学模型,并求取以二维位置为观测序列条件下状态序列的最大概率,最大概率序列即为位置的最佳状态匹配。
步骤5,基于条件随机场算法的高度信息地图匹配算法。
步骤5.1,依据台阶的高度和行人行走的台阶的极限值,将行人每一步行走的高度划分为不同状态。
步骤5.2,建立以高度为观测点和与该观测点对应的状态点之间的特征方程。
步骤5.3,求取之前所有相邻观测点高度信息与匹配的状态点高度之间的均方差。
步骤5.4,以均方差作为另一个特征,根据相邻状态点高度与各个状态点的高度之差的关系建立特征方程。
步骤5.5,建立基于条件随机场的高度地图匹配数学模型,并求取以高度为观测序列条件下状态序列的最大概率,最大概率序列即为高度的最佳状态匹配。
步骤6:二维位置与高度信息融合
步骤6.1,利用二维位置最佳匹配的状态序列查询相应的观测序列的采样时间进行存储,同时利用高度最佳匹配状态序列查询相应观测序列的采样时间进行存储。
步骤6.2,利用临近时间的方法对二维最佳匹配点与高度最佳匹配点进行组合。
步骤6.3,将匹配的三维位置对惯性导航***输出的三维位置信息进行修正。
与现有技术相比,本发明具有以下有益效果:
第一,实现行人三维位置定位,算法精度高;第二,采用二维位置与高度信息分开匹配的方法,简化算法的复杂性;第三,采用临近时刻对二维位置与高度信息进行融合,提高了地图匹配的正确率。
附图说明
图1是本发明所涉及方法的框架结构图;
图2是行人佩戴传感器的方式;
图3是惯性导航解算的***流程图;
图4是经过处理前后室内结构对比图;
图5是经过处理后的电子地图和状态点;
图6是室内三维定位***的流程图;
图7是惯性导航输出的三维位置信息;
图8是三维室内地图匹配输出行人的匹配轨迹。
具体实施方式
以下结合附图对本发明作进一步详细说明。
本发明所涉及方法的框架结构如图1所示,包括以下步骤:
步骤1:数据采集,初步解算室内行人的三维位置和航向。
步骤1.1,行人通过佩戴MEMS-INS传感器,采集行人运动数据,行人佩戴传感器的方式如图2所示。行人运动数据数据包括:三轴加速度数据、三轴陀螺数据。
步骤1.2,采用捷联解算算法对采集的行人运动数据进行三维位置航向信息进行求解,惯性导航解算的***流程图如图3所示。
基于物理学,在采样间隔很短的条件下,速度位置关系满足:
其中,pn表示行人的三维位置坐标;t表示采样时刻,dn表示行人行走的距离;T为采样间隔;vn表示行人的速度信息,an表示行人加速度;表示由载体坐标系b到导航坐标系n的坐标转换矩阵;gn为重力加速度。
步骤2:提取条件随机场模型的观测点。
步骤2.1,以定长行走距离提取水平二维位置信息,零速时刻提取行人高度信息,获取CRF模型的观测点,并分别记录行人水平二维位置与高度信息的采样时刻。
(1)水平二维位置观测点提取模型如下,当满足如下条件:
记录此刻位置点坐标作为二维位置观测点,同时记录下相应的时间:
time1(tob1)=t
其中Distance表示当前时刻与上一时刻采样点坐标的欧氏距离;Pob1(tob1)表示tob1观测点的二维坐标;表示t时刻导航输出的位置坐标;Threshold为设置的行走距离阈值;time1表示所有观测点所对应的时刻。
(2)高度观测点提取模型如下,当满足条件:
if v(t)==0
记录此刻高度信息作为高度观测点,同时记录下相应的时间:
time2(tob2)=t
其中Hob2(tob2)表示第tob2个高度观测点;time2表示高度观测点对应的时刻;表示t时刻的导航输出的高度信息。
步骤3:建立室内电子地图,并根据室内结构信息创建状态点,存储状态点坐标。
步骤3.1,利用已知的室内地图创建室内电子地图,电子地图如图4所示,并将电子地图存储于导航计算机中。
步骤3.2,求取电子地图中数值最小与最大的坐标点作为状态点的取值范围,然后以等间距的状态点覆盖整个取值范围,存储所有状态点坐标信息(X,Y)。
(1)求解电子地图中数值最小与最大的坐标点:
pmin=(xmin,ymin)=min(X,Y)
pmax=(xmax,ymax)=max(X,Y)
其中<pmin,pmax>表示状态点所覆盖范围的最小位置点、最大位置点。
(2)根据最大位置点、最小位置点得到状态点的分布范围为:
r1=pmin=(xmin,ymin)
r2=(xmin,ymax)
r3=(xmax,ymin)
r4=pmax=(xmax,ymax)
其中,<r1,r2,r3,r4>表示状态点范围矩阵的四个顶点坐标。
(3)状态点坐标求解,选取最小的坐标点作为第一个状态点,选取Threshold长度作为状态点之间的间隔,状态点模型如下:
state1(0,0)=r1
state1=(is,js)=r1+(is×Threshold,js×Threshold)<r4
其中state1为所有的状态点总称;(is,js)表示状态点存储的状态。
步骤3.3,由于行人行走楼道时每一步高度变化量为台阶高度的整数倍,因此,添加楼梯区域内的状态点高度是以台阶高度为标准,设台阶高度stair_high结合二维位置的状态点分布的规律,台阶高度的运算模型如下:
State2(N)=0+N×stair_high
其中,State2为所有台阶高度状态点集合,N表示个数,State2(N)由台阶个数决定。结合state1与State2得状态点三维坐标,存储所有状态点坐标信息。则三维状态点与数字地图如图5所示。
步骤4,基于条件随机场算法的二维位置地图匹配算法,匹配及融合方法流程图如图6所示。
步骤4.1,根据二维位置观测点坐标与各个状态点坐标之间关系建立特征方程;二维位置对应的状态点被表示为Sp。
其中,fc表示观测点坐标与状态点坐标之间的关系;Sp(tob1)表示tob1时刻状态点的(x,y)坐标,并且Pob1(tob1)表示tob1时刻观测点的(x,y)坐标,并且σc表示状态点与观测点距离误差的协方差。
步骤4.2,计算观测点的方位信息,并根据观测点的方位信息与该观测点的方位信息所对应时刻状态点之间的方位角建立特征方程。
其中σθ表示观测方位角误差的协方差;b(Sp(tob1-1),Sp(tob1))表示tob1-1与tob1时刻状态点之间的函数;θ(Sp(tob1-1),Sp(tob1)表示tob1-1与tob1时刻状态点之间方位角函数,该函数以地图坐标系X轴正方向为参考。
步骤4.3,建立基于条件随机场的二维地图匹配数学模型,并求取以二维位置为观测序列条件下状态序列的最大概率,最大概率序列即为位置的最佳状态匹配。
采用维特比算法计算最大概率状态点序列SP *。
其中,λp,μp分别表示二维地图匹配数学模型中每个特征对应的权值,在此将权值全部设置为1;i,l表示特征函数的个数;Zob1是规范化因子。
步骤5,基于条件随机场算法的高度信息地图匹配算法。
步骤5.1,依据台阶的高度和行人行走的台阶的极限值,估计行人每一步行走所对应的状态。
利用行人每次行走所跨越的台阶的极限,假设行人每次最多跨过NT个台阶,则每次行人行走的高度变化量应该在(-NT,NT)这个范围内。因此,行人每次行走的所对应的高度的状态为:
Sh=(((-NT)×stair_high),((1-NT)×stair_high),…,(NT×stair_high))
步骤5.2,建立以高度为观测点和与之对应的状态点之间的特征方程。
其中,Sh(tob2)表示tob2时刻状态点高度;Hob2(tob2)表示tob2时刻观测点高度;g表示观测点与状态点高度之间的函数关系;σh表示状态点与观测点之间高度之间的协方差;stair_high表示每个台阶的高度;b(Sh(tob2-1),Sh(tob2))表示tob2-1与tob2时刻状态点之间的函数;h(Sh(tob2-1),Sh(tob2))表示tob2-1与tob2时刻状态点之间的相对高度函数。
步骤5.3,求取之前所有观测点高度信息与匹配的状态点高度之间的均方差。
δHob2(tob2)=Hob2(tob2)-Sh *(tob2)
δHob2观测点高度和匹配状态点高度之间的误差矢量,表示平均误差矢量。
步骤5.4,以均方差作为另一个特征,根据相邻状态点高度与各个状态点的高度之差的关系建立特征方程。
gc(Sh(tob2),Hob2(tob2),S_Hob2(tob2))=(Hob2(tob2)-Sh(tob2))-S_Hob2(tob2)
其中,σs表示高度误差协方差;S_Hob2(tob2)表示tob2时刻之前所有观测误差的协方差。
步骤5.5,建立基于条件随机场的高度地图匹配数学模型,并求取以高度为观测序列条件下状态序列的最大概率,最大概率序列即为高度的最佳状态匹配。
采用维特比算法计算最大概率状态点序列Sh*。
其中,λh,μh表示特征对应的权值,在此将权值全部设置为1。i,l表示特征函数的个数;Zob2是规范化因子。
步骤6:二维位置与高度信息融合
步骤6.1,利用二维位置最佳匹配的状态序列查询相应的观测序列的采样时间进行存储,同时利用高度最佳匹配状态序列查询相应观测序列的采样时间进行存储。
步骤6.2,利用临近时间的方法对二维最佳匹配点与高度最佳匹配点进行组合。
time=|time1(tob1)-time2(kob2)|kob2=1…tob2
其中要求:0<tob1<tob2。找出time中最小的值,并将对应的kob2标记为
S*(tob1)=<SP *(tob1),S(tob1)>
其中,S(tob1)表示tob1时刻采用临近点融合得到的行人高度信息;S*为最终行人的轨迹信息。
步骤6.3,将匹配的三维位置对惯性导航***输出的三维位置信息进行修正反馈。修正反馈的数学模型如下:
pn(t)=S*(tob1)
为了验证算法的有效性,进行了实验验证。以某室内办公环境为例,实验场所包括走廊,楼道两种室内环境。惯性导航***解算出的三维位置信息如图7所示。可见,导航输出无论二维轨迹还是高度信息都存在误差,定位精度不准确。基于条件随机场的三维室内地图匹配算法如图8所示。实验结果显示该方法匹配的结果具有很高的准确性和有效性。
Claims (2)
1.行人室内三维地图匹配方法,其特征在于:该方法包括以下步骤,
步骤1:数据采集,初步解算室内行人的三维位置和航向;
步骤1.1,行人通过佩戴MEMS-INS传感器,采集行人运动数据,行人运动数据包括:三轴加速度数据和三轴陀螺数据;
步骤1.2,应用捷联解算算法对采集的行人运动数据进行三维位置和航向信息进行求解;
步骤2:提取条件随机场模型的观测点;
步骤2.1,以定长行走距离提取室内行人的水平二维位置信息,零速时刻提取行人高度信息,获取CRF模型的观测点,并分别记录室内行人的二维位置与高度信息的采样时刻;
步骤3:建立室内电子地图,并根据室内结构信息创建状态点,存储状态点坐标;
步骤3.1,利用已知的室内地图创建室内电子地图,并将室内电子地图存储于导航计算机中;
步骤3.2,求取电子地图中数值最小与最大的坐标点作为状态点的取值范围,然后以等间距的状态点覆盖整个室内范围,存储状态点信息;
步骤3.3,以台阶高度为标准,添加高度信息上的状态点;
步骤4,基于条件随机场算法的二维位置地图匹配算法;
步骤4.1,根据二维位置观测点坐标与各个状态点坐标之间关系建立特征方程;
步骤4.2,根据观测点的方位信息与与之对应时刻状态点之间的方位角建立特征方程;
步骤4.3,建立基于条件随机场的二维地图匹配数学模型,并求取以二维位置为观测序列条件下状态序列的最大概率,最大概率序列即为位置的最佳状态匹配;
步骤5,基于条件随机场算法的高度信息地图匹配算法;
步骤5.1,依据台阶的高度和行人行走的台阶的极限值,将行人每一步行走的高度划分为不同状态;
步骤5.2,建立以高度为观测点和与该观测点对应的状态点之间的特征方程;
步骤5.3,求取之前所有相邻观测点高度信息与匹配的状态点高度之间的均方差;
步骤5.4,以均方差作为另一个特征,根据相邻状态点高度与各个状态点的高度之差的关系建立特征方程;
步骤5.5,建立基于条件随机场的高度地图匹配数学模型,并求取以高度为观测序列条件下状态序列的最大概率,最大概率序列即为高度的最佳状态匹配;
步骤6:二维位置与高度信息融合
步骤6.1,利用二维位置最佳匹配的状态序列查询相应的观测序列的采样时间进行存储,同时利用高度最佳匹配状态序列查询相应观测序列的采样时间进行存储;
步骤6.2,利用临近时间的方法对二维最佳匹配点与高度最佳匹配点进行组合;
步骤6.3,将匹配的三维位置对惯性导航***输出的三维位置信息进行修正。
2.根据权利要求1所述的行人室内三维地图匹配方法,其特征在于:包括以下步骤:
步骤1:数据采集,初步解算室内行人的三维位置和航向;
步骤1.1,行人通过佩戴MEMS-INS传感器,采集行人运动数据;行人运动数据数据包括:三轴加速度数据、三轴陀螺数据;
步骤1.2,采用捷联解算算法对采集的行人运动数据进行三维位置航向信息进行求解;
基于物理学,在采样间隔很短的条件下,速度位置关系满足:
其中,pn表示行人的三维位置坐标;t表示采样时刻,dn表示行人行走的距离;T为采样间隔;vn表示行人的速度信息,an表示行人加速度;表示由载体坐标系b到导航坐标系n的坐标转换矩阵;gn为重力加速度;
步骤2:提取条件随机场模型的观测点;
步骤2.1,以定长行走距离提取水平二维位置信息,零速时刻提取行人高度信息,获取CRF模型的观测点,并分别记录行人水平二维位置与高度信息的采样时刻;
(1)水平二维位置观测点提取模型如下,当满足如下条件:
记录此刻位置点坐标作为二维位置观测点,同时记录下相应的时间:
time1(tob1)=t
其中Distance表示当前时刻与上一时刻采样点坐标的欧氏距离;Pob1(tob1)表示tob1观测点的二维坐标;表示t时刻导航输出的位置坐标;Threshold为设置的行走距离阈值;time1表示所有观测点所对应的时刻;
(2)高度观测点提取模型如下,当满足条件:
if v(t)==0
记录此刻高度信息作为高度观测点,同时记录下相应的时间:
time2(tob2)=t
其中Hob2(tob2)表示第tob2个高度观测点;time2表示高度观测点对应的时刻;表示t时刻的导航输出的高度信息;
步骤3:建立室内电子地图,并根据室内结构信息创建状态点,存储状态点坐标;
步骤3.1,利用已知的室内地图创建室内电子地图,并将电子地图存储于导航计算机中;
步骤3.2,求取电子地图中数值最小与最大的坐标点作为状态点的取值范围,然后以等间距的状态点覆盖整个取值范围,存储所有状态点坐标信息(X,Y);
(1)求解电子地图中数值最小与最大的坐标点:
pmin=(xmin,ymin)=min(X,Y)
pmax=(xmax,ymax)=max(X,Y)
其中<pmin,pmax>表示状态点所覆盖范围的最小位置点、最大位置点;
(2)根据最大位置点、最小位置点得到状态点的分布范围为:
r1=pmin=(xmin,ymin)
r2=(xmin,ymax)
r3=(xmax,ymin)
r4=pmax=(xmax,ymax)
其中,<r1,r2,r3,r4>表示状态点范围矩阵的四个顶点坐标;
(3)状态点坐标求解,选取最小的坐标点作为第一个状态点,选取Threshold长度作为状态点之间的间隔,状态点模型如下:
state1(0,0)=r1
state1=(is,js)=r1+(is×Threshold,js×Threshold)<r4
其中state1为所有的状态点总称;(is,js)表示状态点存储的状态;
步骤3.3,由于行人行走楼道时每一步高度变化量为台阶高度的整数倍,因此,添加楼梯区域内的状态点高度是以台阶高度为标准,设台阶高度stair_high结合二维位置的状态点分布的规律,台阶高度的运算模型如下:
State2(N)=0+N×stair_high
其中,State2为所有台阶高度状态点集合,N表示个数,State2(N)由台阶个数决定;结合state1与State2得状态点三维坐标,存储所有状态点坐标信息;
步骤4,基于条件随机场算法的二维位置地图匹配算法;
步骤4.1,根据二维位置观测点坐标与各个状态点坐标之间关系建立特征方程;二维位置对应的状态点被表示为Sp;
其中,fc表示观测点坐标与状态点坐标之间的关系;Sp(tob1)表示tob1时刻状态点的(x,y)坐标,并且Pob1(tob1)表示tob1时刻观测点的(x,y)坐标,并且σc表示状态点与观测点距离误差的协方差;
步骤4.2,计算观测点的方位信息,并根据观测点的方位信息与该观测点的方位信息所对应时刻状态点之间的方位角建立特征方程;
其中σθ表示观测方位角误差的协方差;b(Sp(tob1-1),Sp(tob1))表示tob1-1与tob1时刻状态点之间的函数;θ(Sp(tob1-1),Sp(tob1)表示tob1-1与tob1时刻状态点之间方位角函数,该函数以地图坐标系X轴正方向为参考;
步骤4.3,建立基于条件随机场的二维地图匹配数学模型,并求取以二维位置为观测序列条件下状态序列的最大概率,最大概率序列即为位置的最佳状态匹配;
采用维特比算法计算最大概率状态点序列SP *;
其中,λp,μp分别表示二维地图匹配数学模型中每个特征对应的权值,在此将权值全部设置为1;i,l表示特征函数的个数;Zob1是规范化因子;
步骤5,基于条件随机场算法的高度信息地图匹配算法;
步骤5.1,依据台阶的高度和行人行走的台阶的极限值,估计行人每一步行走所对应的状态;
利用行人每次行走所跨越的台阶的极限,假设行人每次最多跨过NT个台阶,则每次行人行走的高度变化量应该在(-NT,NT)这个范围内;因此,行人每次行走的所对应的高度的状态为:
Sh=(((-NT)×stair_high),((1-NT)×stair_high),…,(NT×stair_high))
步骤5.2,建立以高度为观测点和与之对应的状态点之间的特征方程;
其中,Sh(tob2)表示tob2时刻状态点高度;Hob2(tob2)表示tob2时刻观测点高度;g表示观测点与状态点高度之间的函数关系;σh表示状态点与观测点之间高度之间的协方差;stair_high表示每个台阶的高度;b(Sh(tob2-1),Sh(tob2))表示tob2-1与tob2时刻状态点之间的函数;h(Sh(tob2-1),Sh(tob2))表示tob2-1与tob2时刻状态点之间的相对高度函数;
步骤5.3,求取之前所有观测点高度信息与匹配的状态点高度之间的均方差;
δHob2(tob2)=Hob2(tob2)-Sh *(tob2)
δHob2观测点高度和匹配状态点高度之间的误差矢量,表示平均误差矢量;
步骤5.4,以均方差作为另一个特征,根据相邻状态点高度与各个状态点的高度之差的关系建立特征方程;
gc(Sh(tob2),Hob2(tob2),S_Hob2(tob2))=(Hob2(tob2)-Sh(tob2))-S_Hob2(tob2)
其中,σs表示高度误差协方差;S_Hob2(tob2)表示tob2时刻之前所有观测误差的协方差;
步骤5.5,建立基于条件随机场的高度地图匹配数学模型,并求取以高度为观测序列条件下状态序列的最大概率,最大概率序列即为高度的最佳状态匹配;
采用维特比算法计算最大概率状态点序列Sh *;
其中,λh,μh表示特征对应的权值,在此将权值全部设置为1;i,l表示特征函数的个数;Zob2是规范化因子;
步骤6:二维位置与高度信息融合
步骤6.1,利用二维位置最佳匹配的状态序列查询相应的观测序列的采样时间进行存储,同时利用高度最佳匹配状态序列查询相应观测序列的采样时间进行存储;
步骤6.2,利用临近时间的方法对二维最佳匹配点与高度最佳匹配点进行组合;
time=|time1(tob1)-time2(kob2)|kob2=1…tob2
其中要求:0<tob1<tob2;找出time中最小的值,并将对应的kob2标记为
其中,S(tob1)表示tob1时刻采用临近点融合得到的行人高度信息;S*为最终行人的轨迹信息;
步骤6.3,将匹配的三维位置对惯性导航***输出的三维位置信息进行修正反馈;修正反馈的数学模型如下:
pn(t)=S*(tob1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810176554.5A CN108426582B (zh) | 2018-03-03 | 2018-03-03 | 行人室内三维地图匹配方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810176554.5A CN108426582B (zh) | 2018-03-03 | 2018-03-03 | 行人室内三维地图匹配方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108426582A true CN108426582A (zh) | 2018-08-21 |
CN108426582B CN108426582B (zh) | 2021-07-30 |
Family
ID=63157697
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810176554.5A Active CN108426582B (zh) | 2018-03-03 | 2018-03-03 | 行人室内三维地图匹配方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108426582B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110032709A (zh) * | 2019-01-24 | 2019-07-19 | 太原理工大学 | 一种用于地理坐标转换中异常点的定位与估值方法 |
CN110337065A (zh) * | 2019-05-09 | 2019-10-15 | 南京工程学院 | 一种基于三维地图的智能枢纽人员定位监测预警***及方法 |
CN110543917A (zh) * | 2019-09-06 | 2019-12-06 | 电子科技大学 | 一种利用行人惯导轨迹与视频信息的室内地图匹配方法 |
CN111982132A (zh) * | 2019-05-22 | 2020-11-24 | 合肥四维图新科技有限公司 | 数据处理方法、装置和存储介质 |
CN113720332A (zh) * | 2021-06-30 | 2021-11-30 | 北京航空航天大学 | 一种基于层高模型的楼层自主识别方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8605998B2 (en) * | 2011-05-06 | 2013-12-10 | Toyota Motor Engineering & Manufacturing North America, Inc. | Real-time 3D point cloud obstacle discriminator apparatus and associated methodology for training a classifier via bootstrapping |
CN104023228A (zh) * | 2014-06-12 | 2014-09-03 | 北京工业大学 | 一种基于全局运动估计的自适应室内视觉定位方法 |
CN106871894A (zh) * | 2017-03-23 | 2017-06-20 | 北京工业大学 | 一种基于条件随机场的地图匹配方法 |
CN107179085A (zh) * | 2016-03-10 | 2017-09-19 | 中国科学院地理科学与资源研究所 | 一种面向稀疏浮动车数据的条件随机场地图匹配方法 |
CN107635204A (zh) * | 2017-09-27 | 2018-01-26 | 深圳大学 | 一种运动行为辅助的室内融合定位方法及装置、存储介质 |
CN108322889A (zh) * | 2018-02-01 | 2018-07-24 | 深圳市交投科技有限公司 | 一种多源数据融合室内定位的方法、存储介质及智能终端 |
-
2018
- 2018-03-03 CN CN201810176554.5A patent/CN108426582B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8605998B2 (en) * | 2011-05-06 | 2013-12-10 | Toyota Motor Engineering & Manufacturing North America, Inc. | Real-time 3D point cloud obstacle discriminator apparatus and associated methodology for training a classifier via bootstrapping |
CN104023228A (zh) * | 2014-06-12 | 2014-09-03 | 北京工业大学 | 一种基于全局运动估计的自适应室内视觉定位方法 |
CN107179085A (zh) * | 2016-03-10 | 2017-09-19 | 中国科学院地理科学与资源研究所 | 一种面向稀疏浮动车数据的条件随机场地图匹配方法 |
CN106871894A (zh) * | 2017-03-23 | 2017-06-20 | 北京工业大学 | 一种基于条件随机场的地图匹配方法 |
CN107635204A (zh) * | 2017-09-27 | 2018-01-26 | 深圳大学 | 一种运动行为辅助的室内融合定位方法及装置、存储介质 |
CN108322889A (zh) * | 2018-02-01 | 2018-07-24 | 深圳市交投科技有限公司 | 一种多源数据融合室内定位的方法、存储介质及智能终端 |
Non-Patent Citations (2)
Title |
---|
MINGRONG REN 等: "Indoor Pedestrian Navigation Based on Conditional Random Field Algorithm", 《MICROMACHINES》 * |
张晓军 等: "一种基于语音识别的室内定位方法", 《小型微型计算机***》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110032709A (zh) * | 2019-01-24 | 2019-07-19 | 太原理工大学 | 一种用于地理坐标转换中异常点的定位与估值方法 |
CN110032709B (zh) * | 2019-01-24 | 2023-04-14 | 太原理工大学 | 一种用于地理坐标转换中异常点的定位与估值方法 |
CN110337065A (zh) * | 2019-05-09 | 2019-10-15 | 南京工程学院 | 一种基于三维地图的智能枢纽人员定位监测预警***及方法 |
CN111982132A (zh) * | 2019-05-22 | 2020-11-24 | 合肥四维图新科技有限公司 | 数据处理方法、装置和存储介质 |
CN110543917A (zh) * | 2019-09-06 | 2019-12-06 | 电子科技大学 | 一种利用行人惯导轨迹与视频信息的室内地图匹配方法 |
CN110543917B (zh) * | 2019-09-06 | 2021-09-28 | 电子科技大学 | 一种利用行人惯导轨迹与视频信息的室内地图匹配方法 |
CN113720332A (zh) * | 2021-06-30 | 2021-11-30 | 北京航空航天大学 | 一种基于层高模型的楼层自主识别方法 |
CN113720332B (zh) * | 2021-06-30 | 2022-06-07 | 北京航空航天大学 | 一种基于层高模型的楼层自主识别方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108426582B (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108426582A (zh) | 行人室内三维地图匹配方法 | |
Zhou et al. | Activity sequence-based indoor pedestrian localization using smartphones | |
Renaudin et al. | Evaluating indoor positioning systems in a shopping mall: The lessons learned from the IPIN 2018 competition | |
Dong et al. | ViNav: A vision-based indoor navigation system for smartphones | |
CN107635204B (zh) | 一种运动行为辅助的室内融合定位方法及装置、存储介质 | |
CN104897161B (zh) | 基于激光测距的室内平面地图制图方法 | |
Badino et al. | Visual topometric localization | |
CN109945856A (zh) | 基于惯性/雷达的无人机自主定位与建图方法 | |
KR102226846B1 (ko) | Imu 센서와 카메라를 이용한 하이브리드 실내 측위 시스템 | |
CN108051002A (zh) | 基于惯性测量辅助视觉的运输车空间定位方法及*** | |
CN105516929B (zh) | 建立室内地图数据的方法、室内定位方法及对应装置 | |
KR20100072590A (ko) | 동적 환경에서 모바일 플랫폼의 지도 작성방법 | |
CN110207704B (zh) | 一种基于建筑楼梯场景智能识别的行人导航方法 | |
CA3125501A1 (en) | Method and system for tracking a mobile device | |
Lu et al. | Hybrid navigation method of INS/PDR based on action recognition | |
Wu et al. | A pedestrian dead-reckoning system for walking and marking time mixed movement using an SHSs scheme and a foot-mounted IMU | |
CN112729301A (zh) | 一种基于多源数据融合的室内定位方法 | |
CN110412596A (zh) | 一种基于图像信息和激光点云的机器人定位方法 | |
Ishihara et al. | Beacon-guided structure from motion for smartphone-based navigation | |
CN115574816A (zh) | 仿生视觉多源信息智能感知无人平台 | |
CN113971438A (zh) | 荒漠环境下多传感器融合的定位与建图方法 | |
Kessler et al. | Multi-Sensor indoor pedestrian navigation system with vision aiding | |
CN112487234A (zh) | 基于***型层次聚类的轨迹分析方法、设备及存储介质 | |
CN116429112A (zh) | 多机器人协同定位方法和装置、设备及存储介质 | |
Saito et al. | Pre-driving needless system for autonomous mobile robots navigation in real world robot challenge 2013 |
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 |