基于惯性传感器的行人室内位置跟踪方法,本发明专利技术涉及一种行人室内位置跟踪方法,具体涉及基于惯性传感器的PDR方法和PF融合地图信息算法。本发明专利技术是要解决单纯使用PDR对行人位置跟踪时,因为惯性传感器自身所具有的长时漂移性而造成位置估计误差大甚至错误估计等情况。一、根据加速度传感器进行检步及步长估计;二、根据陀螺仪测量数据中的三轴角速度变化进行航向角估计,对航向角进行校正,然后根据校正后的航向角与第一步步长进行航迹推算,最后根据步长、航向角通过PDR方法估计位置:三、通过粒子滤波将地图信息和PDR估计结果融合,即完成了基于惯性传感器的行人室内位置跟踪方法。本发明专利技术应用于室内定位技术领域。
【技术实现步骤摘要】
基于惯性传感器的行人室内位置跟踪方法
本专利技术涉及一种行人室内位置跟踪方法,具体涉及基于惯性传感器的PDR方法和PF融合地图信息算法。
技术介绍
近几年,微电子机械系统(Micro-ElectroMechanicalSystems,MEMS)的发展,使得惯性传感器在智能移动终端的应用得到普及。基于惯性传感器的行人室内导航系统,以其无需铺设外部设施的低成本优势,逐渐成为研究热门。该系统的基本原理是采用行人航迹推算法(PDR,PedestrianDead-Reckoning),根据惯性传感器(如加速度计、陀螺仪)的测量数据,主要涉及到检步、估计步长、估计航向角等技术,因此这种系统有时也称为SHSs(Step-and-HeadingSystems)。早期采用的是可穿戴设备,将传感器安装于鞋子、头盔、置于口袋、腰部等,行走中脚部能更多的反映运动特征,所以基于脚部传感器能更好的检步,然而需要额外购置专门设备且不便携带,不适用于普通的行人室内导航系统。在行走过程中,加速度传感器会输出一定的行走特征,通过对测量数据进行峰值检测、过零检测、自相关匹配、频谱分析等,其中一种或几种相结合来识别每一步,后两种因为耗时长而导致实时性不强。对于步长估计,最粗略的是直接设定步长为一个常数,因为行人在常速行走下有一个平均步长,但是现实中,步长受一个人的身高、体型、换步频率等因素影响,不可以一概为均值。研究者们普遍接受一个最大最小加速度差值四次根方的公式计算,但也有人提出其它如基于步频和步长有线性关系,穿戴脚部的传感器采用零速更新(ZUPT,Zero-velocityUPdaTe),基于行人身高等步长估计方法。而航向估计中,最简单的估计方法是将传感器固定于行人,使传感器坐标系与地球坐标系(指南针坐标系,N-E坐标系)三轴重合,这样传感器测得的方位角即是行人的航向角,则经常通过借助重力在加速度传感器三轴上的分量以及磁强计(电子罗盘)来计算求得,而一般情况下,传感器的放置是无固定姿态的(即姿态任意性),这是需要借助其它方法来估计航向,比如主成份分析法(PCA,PrimaryComponentsAnalyze)估计航向角。粒子滤波(PF,ParticleFilter)是在Kalman滤波的状态模型的高斯、线性假设条件不满足的情况下兴起的滤波跟踪算法,它通过使用蒙特卡洛法抽样最大后验估计,从而逼近贝叶斯最优估计。目前各种室内定位如WLAN室内定位、超宽带室内定位、基于视觉室内定位、惯性传感器室内定位等都采用粒子滤波来改善定位精度以及解决多目标跟踪问题,并且使用粒子滤波来融合不同定位系统的数据,从而相互辅助提高定位结果。经常采用粒子滤波来将地图信息加入到位置估计中,根据地图信息来调整粒子滤波中的参数,从而不仅提高精度,也防止了跟踪目标出现穿墙等特殊情况的发生。
技术实现思路
本专利技术是要解决单纯使用PDR对行人位置跟踪时,因为惯性传感器自身所具有的长时漂移性而造成位置估计误差大甚至错误估计等情况,而提供了基于惯性传感器的行人室内位置跟踪方法。基于惯性传感器的行人室内位置跟踪方法按以下步骤实现:一、根据加速度传感器进行检步及步长估计,具体过程为:采用峰值-零值-谷值-时间间隔的方法进行检步,即每行走一步包含1个最大加速度、2个零值、1个最小加速度,并且满足时间间隔合理,按照成人常速行走的速度为2~4步每秒,设置时间间隔下限为250毫秒,S0表示起步,Si,i=1…9表示第i步,步长估计通过下式来获得其中sLen是估计的步长,acci、accAve分别代表每一步中的加速度值、平均加速度值,N表示每步中采集到的数据个数;二、根据陀螺仪测量数据中的三轴角速度变化进行航向角估计,对航向角进行校正,然后根据校正后的航向角与第一步步长进行航迹推算,最后根据步长、航向角通过PDR方法估计位置:其中,表示第k步的PDR估计位置hk代表第k步估计的航向角,sLenk代表第k步的步长;三、通过粒子滤波将地图信息和PDR估计结果融合,即完成了基于惯性传感器的行人室内位置跟踪方法。专利技术效果:基于惯性传感器的行人位置跟踪方法,借助智能手机中的惯性传感器模块,行人航迹推算法初步估计位置,然后采用粒子滤波算法融合地图信息,为用户提供实时位置信息。在行人航迹推算法实现中,提出一种基于加速度值和时间的PZVT检步算法,通过步长估计,针对陀螺仪的角度积分有长时漂移特点,利用三轴角度变化互相矫正,在实验环境中进行了测试,位置估计误差相对减小。使用粒子滤波结合地图信息对估计位置进行滤波,在一定程度上提高了位置估算的准确度。附图说明图1是基于惯性传感器的PDR原理框图;图2是手机传感器坐标系、电子罗盘坐标系示图;图3是基于惯性传感器行人位置跟踪实验环境示意图;图4是Z轴加速度值以及检步结果;图5是未校正及校正后航向角对比图;图6是使用PF和只使用DR的位置跟踪对比图。具体实施方式具体实施方式一:本实施方式的基于惯性传感器的行人室内位置跟踪方法按以下步骤实现:一、根据加速度传感器进行检步及步长估计;二、根据陀螺仪测量数据中的三轴角速度变化进行航向角估计,对航向角进行校正,然后根据校正后的航向角与第一步步长进行航迹推算,最后根据步长、航向角通过PDR方法估计位置:其中,表示第k步的PDR估计位置hk代表第k步估计的航向角,sLenk代表第k步的步长;三、通过粒子滤波将地图信息和PDR估计结果融合,即完成了基于惯性传感器的行人室内位置跟踪方法。具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一具体过程如下:采用峰值-零值-谷值-时间间隔的方法进行检步,即每行走一步包含1个最大加速度、2个零值、1个最小加速度,并且满足时间间隔合理,按照成人常速行走的速度为2~4步每秒,设置时间间隔下限为250毫秒,S0表示起步,Si,i=1…9表示第i步,步长估计通过下式来获得其中sLen是估计的步长,acci、accAve分别代表每一步中的加速度值、平均加速度值,N表示每步中采集到的数据个数。其它步骤及参数与具体实施方式一相同。具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中航向角估计方法为:首先对围绕x、y、z的三轴角速度进行时间积分,分别得到俯仰角、滚动角、方位角,记作Pitch,Roll,Azimuth;采用(3)对航向角进行第一步校正:heading=c1·Pitch+c2·Roll+c3·Azimuth(3)其中,俯仰角Pitch表示围绕着x轴的旋转量,滚动角Roll表示围绕着y轴的旋转量,方位角Azimuth表示围绕着z三轴的旋转量,其中c1,c2,c3为对应的加权系数;对第一步校正后的航向角进行第二步校正,即进行逐次平滑,具体通过下式计算:mean(·)表示取平均值。其它步骤及参数与具体实施方式一或二相同。具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三的通过粒子滤波将地图信息和PDR估计结果融合方法如下:粒子滤波是一种采用蒙特卡洛方法来逼近后验概率的贝叶斯滤波方法,不要求运动系统是线性或者高斯的;一、采用的粒子滤波的状态方程为:其中,xk,yk表示第k步时的位置坐标,xk-1,yk-1表示第k-1步的位置坐标,hk是估计出的航向角,本文档来自技高网...
【技术保护点】
基于惯性传感器的行人室内位置跟踪方法,其特征在于它按以下步骤实现:一、根据加速度传感器进行检步及步长估计;二、根据陀螺仪测量数据中的三轴角速度变化进行航向角估计,对航向角进行校正,然后根据校正后的航向角与第一步步长进行航迹推算,最后根据步长、航向角通过PDR方法估计位置:poskPDR=xkyk=xk-1yk-1+sLenk·cos(hk)sin(hk)---(1)]]>其中,表示第k步的PDR估计位置hk代表第k步估计的航向角,sLenk代表第k步的步长;三、通过粒子滤波将地图信息和PDR估计结果融合,即完成了基于惯性传感器的行人室内位置跟踪方法。
【技术特征摘要】
1.基于惯性传感器的行人室内位置跟踪方法,其特征在于它按以下步骤实现:一、根据加速度传感器进行检步及步长估计,具体过程为:采用峰值-零值-谷值-时间间隔的方法进行检步,即每行走一步包含1个最大加速度、2个零值、1个最小加速度,并且满足时间间隔合理,按照成人常速行走的速度为2~4步每秒,设置时间间隔下限为250毫秒,S0表示起步,Si,i=1…9表示第i步,步长估计通过下式来获得其中sLen是估计的步长,acci、accAve分别代表每一步中的加速度值、平均加速度值,N表示每步中采集到的数据个数;二、根据陀螺仪测量数据中的三轴角速度变化进行航向角估计,对航向角进行校正,然后根据校正后的航向角与第一步步长进行航迹推算,最后根据步长、航向角通过PDR方法估计位置:其中,表示第k步的PDR估计位置hk代表第k步估计的航向角,sLenk代表第k步的步长;三、通过粒子滤波将地图信息和PDR估计结果融合,即完成了基于惯性传感器的行人室内位置跟踪方法。2.根据权利要求1中所述的基于惯性传感器的行人室内位置跟踪方法,其特征在于所述步骤二中航向角估计方法为:首先对围绕x、y、z的三轴角速度进行时间积分,分别得到俯仰角、滚动角、方位角,记作Pitch,Roll,Azimuth;采用(3)对航向角进行第一步校正:heading=c1·Pitch+c2·Roll+c3·Azimuth(3)其中,俯仰角Pitch表示围绕着x轴的旋转量,滚动角Roll表示围绕着y轴的旋转量,方位角Azimuth表示围绕着z三轴的旋转量,其中c1,c2,c3为对应的加权系数;对第一步校正后的航向角进行第二步校正,即进行逐次平滑,具体通过下式计算:
【专利技术属性】
技术研发人员:马琳,邓仲哲,秦丹阳,何晨光,徐玉滨,崔扬,
申请(专利权)人:哈尔滨工业大学,
类型:发明
国别省市:黑龙江;23
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。