CN106767807A - 一种基于身高和运动特征的行人步长综合测量方法 - Google Patents
一种基于身高和运动特征的行人步长综合测量方法 Download PDFInfo
- Publication number
- CN106767807A CN106767807A CN201510810132.5A CN201510810132A CN106767807A CN 106767807 A CN106767807 A CN 106767807A CN 201510810132 A CN201510810132 A CN 201510810132A CN 106767807 A CN106767807 A CN 106767807A
- Authority
- CN
- China
- Prior art keywords
- walking
- length
- pedestrian
- height
- acceleration
- 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/20—Instruments for performing navigational calculations
- G01C21/206—Instruments for performing navigational calculations specially adapted for indoor navigation
-
- 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/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
Abstract
一种基于身高和运动特征的行人步长综合测量方法,该方法综合考虑不同个体身高、行走频率和每一行走周期内加速度方差三个因素,主要利用人体垂直向加速度信息进行步长测量。首先,利用穿戴在人体腰部的惯性测量单元(Inertial Measurement Unit,IMU)采集人体行走过程中身体重心处的加速度,然后由加速度数据计算行走频率和每一行走周期内的加速度方差,再利用行走频率、加速度方差和人体身高与步长的关系实现步长测量。本发明综合考虑了不同个体的身体特征和运动特征,具有可靠性强、精度高的特点,适用于室内自主式行人导航中的步长测量和导航定位,也可用于步态分析、运动训练和健康监测等领域。
Description
技术领域
本发明涉及一种基于身高和运动特征的行人步长综合测量方法,属于惯性测试领域,可提高行人的步长计算精度和行人导航的定位精度。适用于室内自主式行人导航中的步长测量和导航定位,也可用于步态分析、运动训练和健康监测等领域。
背景技术
随着城市化的不断发展和建筑水平的不断提高,现在城市中出现了越来越多大型且结构复杂,会严重遮挡卫星定位信号的建筑,比如大型商场、图书馆、机场航站楼、大型超市以及地下停车场等。当这些复杂室内环境中发生火灾、地震等灾难时,武警和消防人员抢救伤员时的安全尤为重要,因此亟需一套能够自主定位的室内行人导航系统。步长计算在行人导航领域是一个很重要的方面,步长测量精度的高低将直接影响基于步长和姿态角进行解算的行人导航系统定位精度。而且步长测量在其他领域,比如步态分析、运动训练和健康监测等也有着很大的作用和研究价值。以上需求背景要求得到行走过程中的高精度步长信息。现有的步长测量方法主要有基于运动肢体的几何模型估计步长和基于运动特征计算步长的方法。在行走过程中,行人肢体的几何模型并不规则,现有的方法均是对其进行近似以得到较为规则的几何模型,再利用勾股定理对步长进行计算。由于该方法忽略的影响因素较多且几何参数测量时会涉及到积分运算,故会导致误差增大,步长的测量精度降低。基于运动特征的计算方法较多,但是现有的方法对与步长相关的因素均不同程度有所忽略,而且没有考虑不同的行人个体之间的差异,所以会导致步长测量适应性差,精度低。
发明内容
本发明的目的:克服现有技术的不足,提出一种基于身高和运动特征的行人步长综合测量方法。该方法综合考虑行人身高、行走频率和一步之内的加速度方差三个信息,对行走过程中的步长进行计算,具有可靠性强、精度高的特点
本发明的技术解决方案为:一种基于身高和运动特征的行人步长综合测量方法,采用置于腰部后侧中间位置的IMU测量人体行走过程中身体重心处的加速度,由加速度信息检测和计算每一步的频率和每一行走周期内身体重心处垂直向加速度的方差,输入当前使用者的身高,从而对步长进行高精度测量,包括以下步骤:
(1)建立IMU的误差模型,并在转台上进行六位置标定实验,计算出误差模型的参数;
(2)将IMU置于行人身体腰部后侧中间位置,在行走过程中,对身体重心处附近的加速度进行测量;
(3)输入当前携带该IMU的行人身高h;
(4)根据行走过程中行人运动的周期性特征和加速度计垂直向输出,采用双窗口峰值检测法对行走过程中的步伐进行检测,得到每一步加速度峰值时刻时间ti;
(5)计算每一迈步周期(ti-1,ti)内的行走频率fi和加速度方差vi;
(6)根据身高h、行走频率fi和加速度方差vi三个因素与步长之间的关系,计算并输出行走过程中的步长SLi;
(7)重复步骤(4)至(6),直到测量过程结束。
本发明的原理:本装置采用IMU中的加速度计测量行人身体重心处的加速度,利用行走过程中的身体周期性运动特征来测量步长,核心器件是加速度计。考虑到不同个体之间的差异,引入行人身高因素h。
加速度计的输出为:
UA(a)=UA(0)+kAa+fA(a)+εA (1)
式中UA(a)是载体加速度a对应的加速度计输出电压,UA(0)是a为零时的加速度计输出电压,kA是加速度计标度因数,fA(a)是与加速度有关的误差项,εA是加速度计的噪声。
根据加速度信息检测出的每一步峰值时刻ti,得到每一步的行走频率fi:
与此同时,可以由下式计算得到每一迈步周期内身体重心处垂直向加速度方差vi:
式中at和分别为t时刻行人的垂直向加速度和该迈步周期内的垂直向加速度均值,ti-1和ti为第i步前后两个迈步周期的加速度峰值时刻时间,N为该迈步周期内传感器输出的采样数。
当前携带该IMU的行人身高为h,当不同的个体使用该装置测量步长时,输入当前使用者的身高信息。只要使用者不变,该值不需要重新输入。
在实际行走过程中,根据使用者身高h、行走频率fi和一步之内加速度方差vi与行走步长之间的关系,可以按照下式计算步长SLi:
SLi=h*(a*fi+b*vi+c)+d (4)
式中a、b、c、d为该步长测量方法中的四个参数,其值可预先通过多元非线性最小二乘拟合得出。
本发明与现有技术相比的优点在于:
(1)IMU放置于腰部,携带方便,测量准确;
(2)采用双窗口峰值检测法对步伐进行检测,漏检率低;
(3)综合考虑了身高h、行走频率fi和加速度方差vi多个因素进行步长测量,适应性更强,精度高。
附图说明
图1为本发明涉及步长测量方式示意图;
图2为本发明涉及步长测量计算流程图。
具体实施方式
图1是本发明涉及步长测量方式示意图,IMU佩戴在人体腰部后侧中间位置,测量行走过程中身体重心处的加速度。
图2是本发明涉及步长测量计算流程图,i表示第i个迈步周期,h是当前使用者的身高,fi为当前行走频率,vi是当前迈步周期身体重心处垂直向加速度的方差,SLi是当前测量的步长。
步长测量具体分为以下几个步骤:
(1)建立IMU的误差模型,并在转台上进行六位置标定实验,计算出误差模型的参数。
采取在转台上标定的方法标定出公式(1)中的模型误差系数UA(0)j、kAj、fA(a)j(j=1,2,3,分别代表X、Y、Z轴)。根据fA(a)j中取的与a有关项系数的不同可采取不同的标定方案,本装置采用六位置标定方案。
六位置方案是将惯性测量单元的x轴、y轴、z轴分别指向转台坐标系zp轴的正向和负向,共六个位置,在每一个位置列写相应的方程。利用不同位置下加速度计的输出值与标称值之间的关系列写方程组,其中惯性器件的输出是含有待定误差系数的,通过求解方程组解出这些误差系数就实现了对误差系数的标定。标定实验是在转台上进行的,记转台坐标系为oxpypzp,其中xp、yp、zp分别指向东、北、天,惯性测量单元的坐标系为oxyz。
(2)将IMU置于行人身体腰部后侧中间位置,在行走过程中,对身体重心处的加速度进行测量。本方法主要根据行走过程中身体重心处加速度信息进行检测和计算步长,腰部后侧中间位置最靠近身体重心,因此将IMU佩戴在腰部后侧中间位置。
(3)输入当前携带该IMU的行人身高h。当不同的个体使用该装置测量步长时,输入当前使用者的身高信息,只要使用者不变,该值不需要重新输入。
(4)根据行走过程中行人运动的周期性特征和加速度计垂直向输出,采用双窗口峰值检测法对行走过程中的步伐进行检测,得到每一步加速度峰值时刻时间ti。
在行人正常行走的一个迈步周期之内,身体重心处的垂直向加速度会先增大到一个峰值再减小到一个谷值依次循环。在本测量方法中,取两个相邻相等长度的窗口对垂直向加速度的输出数据进行检测,当检测到两个窗口中的最大值恰好各自位于左右窗口的相邻边缘时,即该迈步周期的峰值位于双窗口的正中央时,确定该点为这一次迈步周期的峰值,检测峰值的同时也将同一个迈步周期的谷值按照相同的原理检测标记出。同样的道理,依次对于采集到的垂直向加速度的数据进行检测,每一次峰值和谷值交替出现时确定为一个迈步,由此便可以检测出行走过程中的每一个迈步过程。
(5)计算每一迈步周期(ti-1,ti)内的行走频率fi和加速度方差vi。
由于已经检测出了行走时的每个迈步过程,根据式(2)和式(3)便可以计算得到每一步的行走频率fi和该迈步周期内身体重心处垂直向加速度方差vi。
(6)根据身高h、行走频率fi和加速度方差vi三个因素与步长之间的关系:正常行走情况下,行走步长与个体的身高成一定的比例关系,不同个体相同条件下行走步长不同;行走步长还与行走频率和一步之内的加速度方差存在线性的关系,步频越大,加速度方差越大,则步长越长。因此由(3)和(5)得到身高h、行走频率fi和加速度方差vi后,便可通过式(4)计算并输出行走过程中当前的步长SLi。
(7)重复(4)至(6),直到测量过程结束,便可以得到实际行走过程中每一步的步长。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (4)
1.一种基于身高和运动特征的行人步长综合测量方法,其特征在于包括以下步骤:
(1)建立IMU的误差模型,并在转台上进行六位置标定实验,计算出误差模型的参数;
(2)将IMU置于行人身体腰部后侧中间位置,在行走过程中,对身体重心处附近的加速度进行测量;
(3)输入当前携带该IMU的行人身高h;
(4)根据行走过程中行人运动的周期性特征和垂直向加速度计输出,采用双窗口峰值检测法对行走过程中的步伐进行检测,得到每一步加速度峰值时刻时间ti;
(5)计算每一行走周期(ti-1,ti)内的行走频率fi和加速度方差vi;
(6)根据身高h、行走频率fi和加速度方差vi三个因素与步长之间的关系,计算并输出行走过程中的步长SLi;
(7)重复步骤(4)至(6),直到测量过程结束。
2.根据权利要求1所述的基于身高和运动特征的行人步长综合测量方法,其特征在于:所述步骤(2)中,在实际行走过程中,该IMU放置于行人腰部后侧中间位置进行测量。
3.根据权利要求1所述的基于身高和运动特征的行人步长综合测量方法,其特征在于:所述步骤(4)中,采用双窗口峰值检测法对步伐进行检测。
4.根据权利要求1所述的基于身高和运动特征的行人步长综合测量方法,其特征在于:所述步骤(6)中,综合考虑了身高h、行走频率fi和加速度方差vi三个因素测量行走步长SLi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510810132.5A CN106767807B (zh) | 2015-11-20 | 2015-11-20 | 一种基于身高和运动特征的行人步长综合测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510810132.5A CN106767807B (zh) | 2015-11-20 | 2015-11-20 | 一种基于身高和运动特征的行人步长综合测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106767807A true CN106767807A (zh) | 2017-05-31 |
CN106767807B CN106767807B (zh) | 2020-07-03 |
Family
ID=58886231
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510810132.5A Active CN106767807B (zh) | 2015-11-20 | 2015-11-20 | 一种基于身高和运动特征的行人步长综合测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106767807B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115216A (zh) * | 2017-06-23 | 2019-01-01 | 北京方位捷讯科技有限公司 | 行人步长检测方法、装置及系统 |
CN110058265A (zh) * | 2018-01-18 | 2019-07-26 | 孙宏民 | 全球定位系统的分析方法 |
CN113790722A (zh) * | 2021-08-20 | 2021-12-14 | 北京自动化控制设备研究所 | 一种基于惯性数据时频域特征提取的行人步长建模方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101250215B1 (ko) * | 2012-05-31 | 2013-04-03 | 삼성탈레스 주식회사 | 칼만 필터와 보행 상태 추정 알고리즘을 이용한 보행자 관성항법 시스템 및 보행자 관성항법 시스템의 높이 추정 방법 |
US20130085677A1 (en) * | 2011-09-30 | 2013-04-04 | Apple Inc. | Techniques for improved pedometer readings |
CN104457781A (zh) * | 2014-12-22 | 2015-03-25 | 北京航空航天大学 | 一种基于单轴加速度计的自适应步数检测方法 |
CN104729524A (zh) * | 2013-12-18 | 2015-06-24 | 中国移动通信集团公司 | 一种步长估计方法、计步器及计步系统 |
CN104757976A (zh) * | 2015-04-16 | 2015-07-08 | 大连理工大学 | 一种基于多传感器融合的人体步态分析方法和系统 |
CN105022055A (zh) * | 2015-07-05 | 2015-11-04 | 吉林大学 | 一种imu室内定位方法 |
-
2015
- 2015-11-20 CN CN201510810132.5A patent/CN106767807B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130085677A1 (en) * | 2011-09-30 | 2013-04-04 | Apple Inc. | Techniques for improved pedometer readings |
KR101250215B1 (ko) * | 2012-05-31 | 2013-04-03 | 삼성탈레스 주식회사 | 칼만 필터와 보행 상태 추정 알고리즘을 이용한 보행자 관성항법 시스템 및 보행자 관성항법 시스템의 높이 추정 방법 |
CN104729524A (zh) * | 2013-12-18 | 2015-06-24 | 中国移动通信集团公司 | 一种步长估计方法、计步器及计步系统 |
CN104457781A (zh) * | 2014-12-22 | 2015-03-25 | 北京航空航天大学 | 一种基于单轴加速度计的自适应步数检测方法 |
CN104757976A (zh) * | 2015-04-16 | 2015-07-08 | 大连理工大学 | 一种基于多传感器融合的人体步态分析方法和系统 |
CN105022055A (zh) * | 2015-07-05 | 2015-11-04 | 吉林大学 | 一种imu室内定位方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109115216A (zh) * | 2017-06-23 | 2019-01-01 | 北京方位捷讯科技有限公司 | 行人步长检测方法、装置及系统 |
CN109115216B (zh) * | 2017-06-23 | 2021-06-15 | 北京方位捷讯科技有限公司 | 行人步长检测方法、装置及系统 |
CN110058265A (zh) * | 2018-01-18 | 2019-07-26 | 孙宏民 | 全球定位系统的分析方法 |
CN113790722A (zh) * | 2021-08-20 | 2021-12-14 | 北京自动化控制设备研究所 | 一种基于惯性数据时频域特征提取的行人步长建模方法 |
CN113790722B (zh) * | 2021-08-20 | 2023-09-12 | 北京自动化控制设备研究所 | 一种基于惯性数据时频域特征提取的行人步长建模方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106767807B (zh) | 2020-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104061934B (zh) | 基于惯性传感器的行人室内位置跟踪方法 | |
CN104296750B (zh) | 一种零速检测方法和装置以及行人导航方法和系统 | |
CN109579853B (zh) | 基于bp神经网络的惯性导航室内定位方法 | |
CN104406586B (zh) | 基于惯性传感器的行人导航装置和方法 | |
CN107289941B (zh) | 一种基于惯导的室内定位方法与装置 | |
CN206479647U (zh) | 定位系统及汽车 | |
CN101151508B (zh) | 行进方向计测装置和行进方向计测方法 | |
EP2850392B1 (en) | Method for step detection and gait direction estimation | |
CN106705968A (zh) | 基于姿态识别和步长模型的室内惯性导航算法 | |
CN102128625B (zh) | 重力辅助惯性导航系统中重力图匹配的初始匹配方法 | |
CN106908060A (zh) | 一种基于mems惯性传感器的高精度室内定位方法 | |
CN104757976A (zh) | 一种基于多传感器融合的人体步态分析方法和系统 | |
CN104197935B (zh) | 一种基于移动智能终端的室内定位方法 | |
CN106644208B (zh) | 一种骑行能力分析系统及分析方法 | |
CN108537101B (zh) | 一种基于状态识别的行人定位方法 | |
CN103175502A (zh) | 一种基于数据手套低速运动的姿态角检测方法 | |
CN104864873B (zh) | 一种利用人体运动特征辅助地图定位的方法 | |
CN102159920A (zh) | 用于处理来自加速度计的测量结果的方法 | |
CN109646009B (zh) | 基于便携式步态分析系统的步态时空参数的计算方法 | |
CN105771224B (zh) | 一种基于多传感器的运动评估系统 | |
CN106225786A (zh) | 一种自适应的行人导航系统零速区间检测方法 | |
CN107990895A (zh) | 一种基于穿戴式imu的建筑物楼层间行人轨迹跟踪方法及系统 | |
CN106643711A (zh) | 一种基于手持设备的室内定位方法及系统 | |
CN106767807A (zh) | 一种基于身高和运动特征的行人步长综合测量方法 | |
CN109141407A (zh) | 一种长输地下管道定位系统的误差补偿系统及方法 |
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 |