CN108537101A - 一种基于状态识别的行人定位方法 - Google Patents
一种基于状态识别的行人定位方法 Download PDFInfo
- Publication number
- CN108537101A CN108537101A CN201810010536.XA CN201810010536A CN108537101A CN 108537101 A CN108537101 A CN 108537101A CN 201810010536 A CN201810010536 A CN 201810010536A CN 108537101 A CN108537101 A CN 108537101A
- Authority
- CN
- China
- Prior art keywords
- data
- data segment
- pedestrian
- new
- state
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/20—Movements or behaviour, e.g. gesture recognition
-
- 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
- G01C21/18—Stabilised platforms, e.g. by gyroscope
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Engineering & Computer Science (AREA)
- Psychiatry (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- Social Psychology (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
Abstract
一种基于状态识别的行人定位方法,包括以下步骤:步骤1,行人执行多状态的运动,通过移动设备内置传感器采集多状态的原始数据,训练多状态分类器C;步骤2,行人执行待定位的运动,通过移动设备内置传感器采集待定位的原始数据;步骤3,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段;步骤4,利用多状态分类器C,识别出该数据段的状态A;步骤5,在该数据段内,基于识别的状态A,采用合适的航位推算模型,更新行人的三维定位。本发明提供了一种稳定性较好、定位精度较高的基于状态识别的行人定位方法。
Description
技术领域
本发明涉及导航定位技术领域,具体来说涉及一种基于状态识别的行人定位方法。
背景技术
行人定位技术主要用于在室内外环境中,利用或不利用基础设施帮助行人进行定位,能够在诸如老人小孩的安全监测、地震火灾中救援人员的准确定位、大型建筑如医院机场内的导航、基于位置的社交服务等领域提供技术支持,因此具有重要的应用价值和研究意义。
在室外环境中,全球定位系统(Global Positioning System,GPS)是一种主要的定位设备,但在室内环境和部分室外城市环境中,受建筑物遮挡的影响,GPS信号会受到干扰甚至阻断,无法提供稳定、准确的定位。
行人航位推算是一种基于已知的初始位置,通过步伐检测、步长估计、方位估计,逐步更新定位的技术。实际情况中,地面条件较为复杂,除平地外,楼梯、斜坡也很常见。相关研究表明,在不同的地面条件下,行人的运动姿势具有较大的差异性。因此当不区别状态采用同一种行人航位推算模型时,会导致定位精度较低。
发明内容
为了克服已有行人定位方法在多场景下定位精度较低的不足,本发明提供了一种适应多场景的定位精度较高的基于状态识别的行人定位方法,能够识别行人运动的状态,包括平地行走、上坡、下坡、上楼梯、下楼梯、站立等状态,根据状态采用合适的航位推算模型,获得定位,以得到室内外均适用的高精度的行人定位技术。
本发明解决其技术问题所采用的技术方案是:
一种基于状态识别的行人定位方法,包括以下步骤:
步骤1,行人执行多状态的运动,通过移动设备内置传感器采集多状态的原始数据,训练多状态分类器C;
步骤2,行人执行待定位的运动,通过移动设备内置传感器采集待定位的原始数据;
步骤3,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段;
步骤4,利用多状态分类器C,识别出该数据段的状态A;
步骤5,在该数据段内,基于识别的状态A,采用合适的航位推算模型,更新行人的三维定位。
进一步,在所述步骤1中,多状态分类器的训练步骤如下:
步骤1-1,将移动设备放置于裤子前侧口袋,利用移动设备内置的加速度计、重力计、陀螺仪和磁力计,分别采集平地行走、上坡、下坡、上楼梯、下楼梯、站立状态下的传感器数据,得到训练集原始数据的集合采样点数为Nraw,并标记真实的状态,其中,ax,ay,az分别为加速度在设备坐标系x,y,z轴上的分量,gx,gy,gz分别为重力加速度在设备坐标系x,y,z轴上的分量,分别为移动设备的俯仰角、偏航角、滚转角,为移动设备磁力计读取的方位角;
步骤1-2,将训练集原始数据的集合中的数据,利用基于步频的滑动窗口做数据分割,窗口长度为取决于行人步频的介于一个周期和两个周期之间的采样点数Nwin,相邻窗口的覆盖率为77,则每类新数据被分割成M个数据段:
得到训练集原始数据的数据段集合
步骤1-3,对训练集原始数据的数据段集合中的数据段进行预处理,得到预处理数据的数据段,预处理数据包括设备坐标系中加速度的模、世界坐标系中的垂直加速度和水平加速度、腿部张角4类数据,预处理方式如下:
设备坐标系中加速度的模aN:
为了获得世界坐标系中的垂直加速度和水平加速度,首先,定义四元数:
标准化的四元数:
由标准化的四元数计算从设备坐标系到世界坐标系的旋转矩阵R:
从而得到加速度在世界坐标系x,y,z轴上的分量
则在世界坐标系中的垂直加速度aV:
在世界坐标系中的水平加速度aH:
腿部张角θ:
预处理数据的数据段与原始数据的数据段共同构成训练集新数据段的集合
步骤1-4,对训练集新数据段的集合内的每个数据段进行特征提取,获得所有数据段对应的M×N的特征矩阵Ft=[f1,f2,...,fN],其中,N为特征的种类数,fi,(i=1,2,...,N)表示第i种特征的M×1的列向量;
步骤1-5,通过主成分分析法,得到特征矩阵Ft从高维到低维映射的转移矩阵T,进而低维特征矩阵
得到主成分表示的M×N′的低维特征矩阵其中,N′为低维特征的种类数,f′i,(i=1,2,...,N′)表示第i种低维特征的M×1的列向量;
步骤1-6,把低维特征矩阵中列向量f′i,(i=1,2,...,N′)的数值属性作为标准化参数P,如最大值和最小值,或均值和方差等,将f′i,(i=1,2,...,N′)标准化到区间[a,b],得到标准化的特征矩阵其中,f′i,(i=1,2,...,N′)表示第i种标准化的低维特征的M×1的列向量;
步骤1-7,将标准化的特征矩阵及真实的状态作为输入,利用多状态支持向量机算法,将多状态分类器拆分成一系列二元分类器,利用序列最小最优化算法计算二元分类器的最优解,训练得到多状态分类器C={C1,C2,...CI},I表示二元分类器的个数。
再进一步,在所述步骤2中,行人执行待定位的运动,将移动设备放置于裤子前侧口袋,采集移动设备内置的加速度计、重力计、陀螺仪和磁力计传感器的数据,得到待定位的原始数据的集合
更进一步,在所述步骤3中,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段的步骤如下:
步骤3-1,当待定位的原始数据的集合Sraw中的每类数据达到Nwin个采样点,利用长度为Nwin的滑动窗口截断,得到待定位的原始数据的数据段Swin;否则,回到步骤2;
步骤3-2,对待定位的原始数据的数据段Swin执行步骤1-3,得到待定位的新数据段Snew。
在所述步骤4中,利用多状态分类器C识别出该数据段的状态A的步骤如下:
步骤4-1,对步骤3获得的待定位的新数据段Snew,执行步骤1-4,得到特征矩阵F;
步骤4-2,利用步骤1-5中的转移矩阵T对特征矩阵F进行降维,得到低维的特征矩阵Fpca:
Fpca=F×T;
步骤4-3,利用步骤1-6中的标准化参数P,对Fpca进行标准化,得到标准化的特征矩阵Fnorm;
步骤4-4,利用步骤1-7得到的多状态分类器C,输入标准化的特征矩阵Fnorm,输出该数据段的状态A。
在所述步骤5中,在该数据段内,基于识别的状态A,采用航位推算模型,更新行人的三维定位的步骤如下:
步骤5-1,该步骤实现步伐检测,在待定位的新数据段Snew内,从未被定位过的采样点开始遍历,对腿部张角θ进行极值检测,获得极值及其采样点索引j表示当前极值的序号,利用极大值阈值TMAX、极小值阈值TMIN和最小间隔ΔN排除不符合要求的极值,获得有效极值及其采样点索引s表示当前有效极值的序号:当大于TMAX且距离最近一个有效极值的采样点间隔大于ΔN,则视为有效极大值,或当小于TMIN且距离最近一个有效极值的采样点间隔大于ΔN,则视为有效极小值,此时更新有效极值及其采样点索引当检测到一个腿部张角θ的有效极值时,即检测到一个步伐,新数据段Snew内的步数Nstep=Nstep+1;直到新数据段Snew内的采样点遍历结束,若检测到有效极值及其采样点索引(s=1,2,...,Nstep)时,执行步骤5-2,否则,执行步骤5-6;
步骤5-2,该步骤实现步长估计,在新数据段Snew内,选择适用于状态A的步长模型估计步长:当A={平地行走、上坡、下坡}时,利用行人腿长l、第s个有效极值和上一个有效极值及该步对应的地面坡度αs估计步长;当A={上楼梯、下楼梯}时,利用台阶的宽度wstair估计步长;当A={站立}时,步长为零;则第s步步长Ls:
,(s=1,2,...,Nstep)
其中,αs=0表示平地,αs>0表不上坡,αs<0表示下坡;
步骤5-3,该步骤实现步高估计,在新数据段Snew内,选择适用于识别状态A的步高模型估计高度:当A={平地行走、上坡、下坡}时,利用第s步对应的地面坡度αs、步长Ls估计步高;当A={上楼梯、下楼梯}时,利用台阶的高度hstair估计步高,上楼梯时为正,下楼梯时为负;当A={站立}时,步高为零;则第s步步高Hs:
步骤5-4,该步骤实现方位角估计,在新数据段Snew内,结合陀螺仪和磁力计确定方位,利用初始时刻磁力计获取的方位角和磁偏角校正陀螺仪的偏航角得到在世界坐标系中相对于真北方向的绝对方位
第s步的方位角取有效极值前后个采样点的绝对方位的均值:
步骤5-5,该步骤实现行人定位的更新,选择适用于识别状态A的三维位置更新模型,在世界坐标系中,基于位置更新前的行人定位(x0,y0,z0),新数据段Snew最终更新的位置即为当前行人的定位:
其中,系数ks:
步骤5-6,为了更新行人在下一个时刻的定位,保留原始数据的数据段Sraw中在采样点[Nwin×(1-η),Nwin]之间的原始数据,作为与下一个时刻的原始数据的数据段重合的部分,回到步骤2。
所述步骤1-6中,所述数值属性为最大值和最小值或均值和方差。
本发明的有益效果主要表现在:能够识别行人运动的状态,包括平地行走、上坡、下坡、上楼梯、下楼梯、站立等状态,根据状态采用航位推算模型,获得定位,以得到室内外均适用的高精度的行人定位技术。
附图说明
图1为本发明的训练多状态分类器的流程图;
图2为本发明的基于状态识别的行人定位方法的流程图;
图3为本发明实施例中室内的实验结果。
图4为本发明实施例中室外的实验结果。
具体实施方式
下面结合附图对本发明作进一步描述。
实施例1
参照图1~图4,一种基于状态识别的行人定位方法,包括以下步骤:
步骤1,行人执行多状态的运动,通过移动设备内置传感器采集多状态的原始数据,训练多状态分类器C;
步骤2,行人执行待定位的运动,通过移动设备内置传感器采集待定位的原始数据;
步骤3,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段;
步骤4,利用多状态分类器C,识别出该数据段的状态A;
步骤5,在该数据段内,基于识别的状态A,采用合适的航位推算模型,更新行人的三维定位。
其中,所述步骤1的流程图如图1所示,多状态分类器的训练步骤如下:
步骤1-1,将移动设备放置于裤子前侧口袋,利用移动设备内置的加速度计、重力计、陀螺仪、磁力计等,分别采集平地行走、上坡、下坡、上楼梯、下楼梯、站立等状态下的传感器数据,得到训练集原始数据的集合 采样点数为Nraw,并标记真实的状态,其中,ax,ay,az分别为加速度在设备坐标系x,y,z轴上的分量,gx,gy,gz分别为重力加速度在设备坐标系x,y,z轴上的分量,分别为移动设备的俯仰角、偏航角、滚转角,为移动设备磁力计读取的方位角;
步骤1-2,将训练集原始数据的集合中的数据,利用基于步频的滑动窗口做数据分割,窗口长度为取决于行人步频的介于一个周期和两个周期之间的采样点数Nwin,本实施例中,对于快速、正常、慢速三种步频的运动,Nw分别等于50、70、90,相邻窗口的覆盖率为η=50%,则每类新数据被分割成M个数据段:
得到训练集原始数据的数据段集合
步骤1-3,对训练集原始数据的数据段集合中的数据段进行预处理,得到预处理数据的数据段,预处理数据包括设备坐标系中加速度的模、世界坐标系中的垂直加速度和水平加速度、腿部张角等4类数据,具体预处理方式如下:
设备坐标系中加速度的模aN:
为了获得世界坐标系中的垂直加速度和水平加速度,首先,定义四元数:
标准化的四元数:
由标准化的四元数计算从设备坐标系到世界坐标系的旋转矩阵R:
从而得到加速度在世界坐标系x,y,z轴上的分量
则在世界坐标系中的垂直加速度aV:
世界坐标系中的水平加速度aH:
腿部张角θ:
预处理数据的数据段与原始数据的数据段共同构成训练集新数据段的集合
步骤1-4,对训练集新数据段的集合内的每个数据段进行特征提取,本实施例中,对预处理数据的数据段提取时域上极大值、极小值、均值、方差、能量值,及频域上快速傅立叶变换后的主频率幅度,对原始数据的数据段提取加速度运动单元aAU:
其中,分别表示原始数据的数据段内第i个采样点的加速度在设备坐标系x,y,z轴上的分量,分别表示原始数据的数据段内加速度在设备坐标系x,y,z轴上的分量的均值;
从而获得所有数据段对应的M×N的特征矩阵Ft=[f1,f2,...,fN],其中,N为特征的种类数,本实施例中N=25,fi,(i=1,2,...,N)表示第i种特征的M×1的列向量;
步骤1-5,通过主成分分析法,得到特征矩阵Ft从高维到低维映射的转移矩阵T,进而低维特征矩阵
得到主成分表示的M×N′的低维特征矩阵其中,N′为低维特征的种类数,f′i,(i=1,2,...,N′)表示第i种低维特征的M×1的列向量,本实施例中取N′=10;
步骤1-6,将低维特征矩阵中列向量f′i,(i=1,2,...,N′)的数值属性作为标准化参数P,本实施例采用f′i,(i=1,2,...,N′)的最大值fmaxi和最小值fmini,使f′i,(i=1,2,...,N′)标准化到区间[a,b],
步骤1-7,将标准化的特征矩阵及真实的状态作为输入,利用多状态支持向量机算法,将多状态分类器拆分成一系列二元分类器,本实施例采用一对一的方法,即任意两两状态之间构建一个二元分类器,六个状态共构建15个二元分类器,利用序列最小最优化算法计算二元分类器的最优解,训练得到多状态分类器C={C1,C2,...CI},I=15表示二元分类器的个数。
步骤2至步骤5的流程图如图2所示:
在所述步骤2中,行人执行待定位的运动,将移动设备放置于裤子前侧口袋,采集移动设备内置的加速度计、重力计、陀螺仪、磁力计等传感器的数据,得到待定位的原始数据的集合
更进一步,在所述步骤3中,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段的步骤如下:
步骤3-1,当待定位的原始数据的集合Sraw中的每类数据达到Nwin个采样点,利用长度为Nwin的滑动窗口截断,得到待定位的原始数据的数据段Swin;否则,回到步骤2;
步骤3-2,对待定位的原始数据的数据段Swin执行步骤1-3,得到待定位的新数据段Snew;
在步骤4中,利用多状态分类器C识别出该数据段的状态A的步骤如下:
步骤4-1,对步骤3获得的待定位的新数据段Snew,执行步骤1-4,得到特征矩阵F;
步骤4-2,利用步骤1-5中的转移矩阵T对特征矩阵F进行降维,得到低维的特征矩阵Fpca:
Fpca=F×T;
步骤4-3,利用步骤1-6中的标准化参数P,本实施例中为fmaxi,fmini,(i=1,2,...,N′),对Fpca进行标准化,得到标准化的特征矩阵Fnorm;
步骤4-4,利用步骤1-7得到的多状态分类器C,输入标准化的特征矩阵Fnorm,输出该数据段的状态A。
在所述步骤5中,在该数据段内,基于识别的状态A,采用航位推算模型,更新行人的三维定位的步骤如下:
步骤5-1,该步骤实现步伐检测,在待定位的新数据段Snew内,从未被定位过的采样点开始遍历,对腿部张角θ进行极值检测,获得极值及其采样点索引j表示当前极值的序号,利用极大值阈值TMAX、极小值阈值TMIN和最小间隔ΔN排除不符合要求的极值,获得有效极值及其采样点索引s表示当前有效极值的序号:当大于TMAX且距离最近一个有效极值的采样点间隔大于ΔN,则视为有效极大值,或当小于TMIN且距离最近一个有效极值的采样点间隔大于ΔN,则视为有效极小值,此时更新有效极值及其采样点索引当检测到一个腿部张角θ的有效极值时,即检测到一个步伐,新数据段Snew内的步数Nstep=Nstep+1;直到新数据段Snew内的采样点遍历结束,若检测到有效极值及其采样点索引(s=1,2,...,Nstep)时,执行步骤5-2,否则,执行步骤5-6。
步骤5-2,该步骤实现步长估计,在新数据段Snew内,选择适用于状态A的步长模型估计步长:当A={平地行走、上坡、下坡}时,利用行人腿长l、第s个有效极值和上一个有效极值及该步对应的地面坡度αs估计步长;当A={上楼梯、下楼梯}时,利用台阶的宽度wstair估计步长;当A={站立}时,步长为零;则第s步步长Ls:
,(s=1,2,...,Nstep)
其中,αs=0表示平地,αs>0表示上坡,αs<0表示下坡;
步骤5-3,该步骤实现步高估计,在新数据段Snew内,选择适用于识别状态A的步高模型估计高度:当A={平地行走、上坡、下坡}时,利用第s步对应的地面坡度αs、步长Ls估计步高;当A={上楼梯、下楼梯}时,利用台阶的高度hstair估计步高,上楼梯时为正,下楼梯时为负;当A={站立}时,步高为零;则第s步步高Hs:
步骤5-4,该步骤实现方位角估计,在新数据段Snew内,结合陀螺仪和磁力计确定方位,利用初始时刻磁力计获取的方位角和磁偏角校正陀螺仪的偏航角得到在世界坐标系中相对于真北方向的绝对方位
第s步的方位角取有效极值前后个采样点的绝对方位的均值:
步骤5-5,该步骤实现行人定位的更新,选择适用于识别状态A的三维位置更新模型,在世界坐标系中,基于位置更新前的行人定位(x0,y0,z0),新数据段Snew最终更新的位置即为当前行人的定位:
其中,系数ks:
步骤5-6,为了更新行人在下一个时刻的定位,保留原始数据的数据段Sraw中在采样点[Nwin×(1-η),Nwin]之间的原始数据,作为与下一个时刻的原始数据的数据段重合的部分,回到步骤2。
实施例2
为了评估“基于状态识别的行人定位方法”在实际环境中的定位性能,基于实施例1中的技术流程,进行实验,使用设备为苹果公司的iPhone 6手机。利用基于iPhone 6开发的应用,采集手机内置的加速度计、重力计、陀螺仪、磁力计(GPS)的数据,采样频率为50Hz。
首先,验证多状态分类器。实验者按照正常、慢速、快速三种步频,分别进行平地行走、上坡、下坡、上楼梯、下楼梯、站立六种运动,标记真实的状态,组成训练集,用于训练多状态分类器;为了评估多状态分类器的性能,重复上述运动,标记真实的状态,组成测试集。利用训练出的多状态分类器,对测试集进行分类,将分类状态与真实状态进行比较,得到分类的准确度,如表1所示。
表1
其次,验证多状态的行人定位方法,实验环境包括室内和室外。室内实验在浙江大学某建筑内进行,测试长度为101.3米,一段完整的运动包含状态:上楼梯-平地-下坡-平地-上坡-平地-下楼梯-平地,涉及7次状态转换、3次180度转弯、1次90度转弯,起点与终点重合,视为原点。台阶宽度、高度、坡度均为测量多次取均值,分别为0.5米,0.06米,6.8度。基于本方法的实验路径和真实路径如图3所示。
室外实验在浙江大学某建筑外进行,测试长度为107.7米,一段完整的运动包含状态:上坡-平地-下楼梯-平地-下楼梯-平地,涉及5次状态转换,7次90度转弯,起点与终点重合,视为原点。台阶宽度、高度、坡度均为测量多次取均值,分别为0.35米,0.13米,5.1度。基于本方法的实验路径和真实路径如图4所示。
本实验利用步数误差、距离误差、终点误差来描述系统的误差,其中,终点误差表示估计终点到真实终点的欧氏距离与总路程的比值,其他的误差表示估计值与真实值的差的绝对值与真实值的比值。多次实验的平均结果,如表2所示,利用GPS的定位作为对比。
表2
从表2中结果可知,本技术方案在室内和室外环境,均能提供高精度的定位信息,具体来说,在室内环境中,本技术方案能够代替GPS进行定位,在室外且建筑物环绕的环境中,本技术方案优于GPS的定位结果。
以上对本发明所提供的一种基于状态识别的行人定位系统进行了详尽的介绍,实施例只是用于帮助阐述本发明的方法,对本发明的保护范围不构成任何限制,凡采用等同变化或者等效替换形成的技术方案,均落在本发明的权利保护范围之内。
Claims (7)
1.一种基于状态识别的行人定位方法,其特征在于:所述行人定位方法包括以下步骤:
步骤1,行人执行多状态的运动,通过移动设备内置传感器采集多状态的原始数据,训练多状态分类器C;
步骤2,行人执行待定位的运动,通过移动设备内置传感器采集待定位的原始数据;
步骤3,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段;
步骤4,利用多状态分类器C,识别出该数据段的状态A;
步骤5,在该数据段内,基于识别的状态A,采用合适的航位推算模型,更新行人的三维定位。
2.如权利要求1所述的一种基于状态识别的行人定位方法,其特征在于:在所述步骤1中,多状态分类器的训练步骤如下:
步骤1-1,将移动设备放置于裤子前侧口袋,利用移动设备内置的加速度计、重力计、陀螺仪和磁力计,分别采集平地行走、上坡、下坡、上楼梯、下楼梯、站立状态下的传感器数据,得到训练集原始数据的集合采样点数为Nraw,并标记真实的状态,其中,ax,ay,az分别为加速度在设备坐标系x,y,z轴上的分量,gx,gy,gz分别为重力加速度在设备坐标系x,y,z轴上的分量,分别为移动设备的俯仰角、偏航角、滚转角,为移动设备磁力计读取的方位角;
步骤1-2,将训练集原始数据的集合中的数据,利用基于步频的滑动窗口做数据分割,窗口长度为取决于行人步频的介于一个周期和两个周期之间的采样点数Nwin,相邻窗口的覆盖率为η,则每类新数据被分割成M个数据段:
得到训练集原始数据的数据段集合
步骤1-3,对训练集原始数据的数据段集合中的数据段进行预处理,得到预处理数据的数据段,预处理数据包括设备坐标系中加速度的模、世界坐标系中的垂直加速度和水平加速度、腿部张角4类数据,预处理方式如下:
设备坐标系中加速度的模aN:
为了获得世界坐标系中的垂直加速度和水平加速度,首先,定义四元数:
标准化的四元数:
由标准化的四元数计算从设备坐标系到世界坐标系的旋转矩阵R:
从而得到加速度在世界坐标系x,y,z轴上的分量
则在世界坐标系中的垂直加速度aV:
在世界坐标系中的水平加速度aH:
腿部张角θ:
预处理数据的数据段与原始数据的数据段共同构成训练集新数据段的集合
步骤1-4,对训练集新数据段的集合内的每个数据段进行特征提取,获得所有数据段对应的M×N的特征矩阵Ft=[f1,f2,…,fN],其中,N为特征的种类数,fi,(i=1,2,…,N)表示第i种特征的M×1的列向量;
步骤1-5,通过主成分分析法,得到特征矩阵Ft从高维到低维映射的转移矩阵T,进而低维特征矩阵
得到主成分表示的M×N′的低维特征矩阵其中,N′为低维特征的种类数,f′i,(i=1,2,…,N′)表示第i种低维特征的M×1的列向量;
步骤1-6,把低维特征矩阵中列向量f′i,(i=1,2,…,N′)的数值属性作为标准化参数P,将f′i,(i=1,2,…,N′)标准化到区间[a,b],得到标准化的特征矩阵其中,f″i,(i=1,2,…,N′)表示第i种标准化的低维特征的M×1的列向量;
步骤1-7,将标准化的特征矩阵及真实的状态作为输入,利用多状态支持向量机算法,将多状态分类器拆分成一系列二元分类器,利用序列最小最优化算法计算二元分类器的最优解,训练得到多状态分类器C={C1,C2,…CI},I表示二元分类器的个数。
3.如权利要求1或2所述的一种基于状态识别的行人定位方法,其特征在于:在所述步骤2中,行人执行待定位的运动,将移动设备放置于裤子前侧口袋,采集移动设备内置的加速度计、重力计、陀螺仪和磁力计传感器的数据,得到待定位的原始数据的集合
4.如权利要求1或2一种基于状态识别的行人定位方法,其特征在于:在所述步骤3中,对移动设备内置传感器采集的待定位的原始数据进行数据分割、数据预处理,获得数据段的步骤如下:
步骤3-1,当待定位的原始数据的集合Sraw中的每类数据达到Nwin个采样点,利用长度为Nwin的滑动窗口截断,得到待定位的原始数据的数据段Swin;否则,回到步骤2;
步骤3-2,对待定位的原始数据的数据段Swin执行步骤1-3,得到待定位的新数据段Snew。
5.如权利要求1或2所述的一种基于状态识别的行人定位方法,其特征在于:在所述步骤4中,利用多状态分类器C识别出该数据段的状态A的步骤如下:
步骤4-1,对步骤3获得的待定位的新数据段Snew,执行步骤1-4,得到特征矩阵F;
步骤4-2,利用步骤1-5中的转移矩阵T对特征矩阵F进行降维,得到低维的特征矩阵Fpca:
Fpca=F×T;
步骤4-3,利用步骤1-6中的标准化参数P,对Fpca进行标准化,得到标准化的特征矩阵Fnorm;
步骤4-4,利用步骤1-7得到的多状态分类器C,输入标准化的特征矩阵Fnorm,输出该数据段的状态A。
6.如权利要求1或2所述的一种基于状态识别的行人定位方法,其特征在于:在所述步骤5中,在该数据段内,基于识别的状态A,采用航位推算模型,更新行人的三维定位的步骤如下:
步骤5-1,该步骤实现步伐检测,在待定位的新数据段Snew内,从未被定位过的采样点开始遍历,对腿部张角θ进行极值检测,获得极值及其采样点索引j表示当前极值的序号,利用极大值阈值TMAX、极小值阈值TMIN和最小间隔ΔN排除不符合要求的极值,获得有效极值及其采样点索引s表示当前有效极值的序号:当大于TMAX且距离最近一个有效极值的采样点间隔大于ΔN,则视为有效极大值,或当小于TMIN且距离最近一个有效极值的采样点间隔大于ΔN,则视为有效极小值,此时更新有效极值及其采样点索引当检测到一个腿部张角θ的有效极值时,即检测到一个步伐,新数据段Snew内的步数Nstep=Nstep+1;直到新数据段Snew内的采样点遍历结束,若检测到有效极值及其采样点索引时,执行步骤5-2,否则,执行步骤5-6;
步骤5-2,该步骤实现步长估计,在新数据段Snew内,选择适用于状态A的步长模型估计步长:当A={平地行走、上坡、下坡}时,利用行人腿长l、第s个有效极值和上一个有效极值及该步对应的地面坡度αs估计步长;当A={上楼梯、下楼梯}时,利用台阶的宽度wstair估计步长;当A={站立}时,步长为零;则第s步步长Ls:
其中,αs=0表示平地,αs>0表示上坡,αs<0表示下坡;
步骤5-3,该步骤实现步高估计,在新数据段Snew内,选择适用于识别状态A的步高模型估计高度:当A={平地行走、上坡、下坡}时,利用第s步对应的地面坡度αs、步长Ls估计步高;当A={上楼梯、下楼梯}时,利用台阶的高度hstair估计步高,上楼梯时为正,下楼梯时为负;当A={站立}时,步高为零;则第s步步高Hs:
步骤5-4,该步骤实现方位角估计,在新数据段Snew内,结合陀螺仪和磁力计确定方位,利用初始时刻磁力计获取的方位角和磁偏角校正陀螺仪的偏航角得到在世界坐标系中相对于真北方向的绝对方位
第s步的方位角取有效极值前后个采样点的绝对方位的均值:
步骤5-5,该步骤实现行人定位的更新,选择适用于识别状态A的三维位置更新模型,在世界坐标系中,基于位置更新前的行人定位(x0,y0,z0),新数据段Snew最终更新的位置即为当前行人的定位:
其中,系数ks:
步骤5-6,为了更新行人在下一个时刻的定位,保留原始数据的数据段Sraw中在采样点[Nwin×(1-η),Nwin]之间的原始数据,作为与下一个时刻的原始数据的数据段重合的部分,回到步骤2。
7.如权利要求2所述的一种基于状态识别的行人定位方法,其特征在于:所述步骤1-6中,所述数值属性为最大值和最小值或均值和方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810010536.XA CN108537101B (zh) | 2018-01-05 | 2018-01-05 | 一种基于状态识别的行人定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810010536.XA CN108537101B (zh) | 2018-01-05 | 2018-01-05 | 一种基于状态识别的行人定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108537101A true CN108537101A (zh) | 2018-09-14 |
CN108537101B CN108537101B (zh) | 2022-03-25 |
Family
ID=63489843
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810010536.XA Active CN108537101B (zh) | 2018-01-05 | 2018-01-05 | 一种基于状态识别的行人定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108537101B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109597031A (zh) * | 2018-11-19 | 2019-04-09 | 桂林电子科技大学 | 一种基于svm和步行者航位推算融合的室内定位方法 |
CN109827568A (zh) * | 2019-01-29 | 2019-05-31 | 东北大学秦皇岛分校 | 基于mems传感器的多层建筑中行人高度位置估计方法 |
CN110132276A (zh) * | 2019-04-30 | 2019-08-16 | 杭州电子科技大学 | 一种基于行人运动状态的自适应步长估计方法 |
CN110501008A (zh) * | 2019-07-29 | 2019-11-26 | 西安电子科技大学 | 人体运动模型的自主进化方法 |
CN112833889A (zh) * | 2019-11-22 | 2021-05-25 | 北京初速度科技有限公司 | 一种车辆的定位方法及装置 |
CN113709860A (zh) * | 2021-08-12 | 2021-11-26 | 电子科技大学 | 一种层间过渡区域的室内定位装置和方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103471589A (zh) * | 2013-09-25 | 2013-12-25 | 武汉大学 | 一种室内行人行走模式识别和轨迹追踪的方法 |
CN104215238A (zh) * | 2014-08-21 | 2014-12-17 | 北京空间飞行器总体设计部 | 一种智能手机室内定位方法 |
CN104757976A (zh) * | 2015-04-16 | 2015-07-08 | 大连理工大学 | 一种基于多传感器融合的人体步态分析方法和系统 |
CN104931049A (zh) * | 2015-06-05 | 2015-09-23 | 北京信息科技大学 | 一种基于运动分类的行人自主定位方法 |
CN105547298A (zh) * | 2015-12-25 | 2016-05-04 | 北京京元智慧应急技术有限公司 | 智能手机内置传感器与Wi-Fi融合的室内动态连续定位方法 |
US9565521B1 (en) * | 2015-08-14 | 2017-02-07 | Samsung Electronics Co., Ltd. | Automatic semantic labeling based on activity recognition |
CN106574838A (zh) * | 2014-07-03 | 2017-04-19 | 德州仪器公司 | 行人导航装置及方法 |
CN106705968A (zh) * | 2016-12-09 | 2017-05-24 | 北京工业大学 | 基于姿态识别和步长模型的室内惯性导航算法 |
CN106767810A (zh) * | 2016-11-23 | 2017-05-31 | 武汉理工大学 | 一种基于移动终端的wifi和视觉信息的室内定位方法及系统 |
US20170241787A1 (en) * | 2014-09-15 | 2017-08-24 | Oxford University Innovation Limited | Determining the position of a mobile device in a geographical area |
-
2018
- 2018-01-05 CN CN201810010536.XA patent/CN108537101B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103471589A (zh) * | 2013-09-25 | 2013-12-25 | 武汉大学 | 一种室内行人行走模式识别和轨迹追踪的方法 |
CN106574838A (zh) * | 2014-07-03 | 2017-04-19 | 德州仪器公司 | 行人导航装置及方法 |
CN104215238A (zh) * | 2014-08-21 | 2014-12-17 | 北京空间飞行器总体设计部 | 一种智能手机室内定位方法 |
US20170241787A1 (en) * | 2014-09-15 | 2017-08-24 | Oxford University Innovation Limited | Determining the position of a mobile device in a geographical area |
CN104757976A (zh) * | 2015-04-16 | 2015-07-08 | 大连理工大学 | 一种基于多传感器融合的人体步态分析方法和系统 |
CN104931049A (zh) * | 2015-06-05 | 2015-09-23 | 北京信息科技大学 | 一种基于运动分类的行人自主定位方法 |
US9565521B1 (en) * | 2015-08-14 | 2017-02-07 | Samsung Electronics Co., Ltd. | Automatic semantic labeling based on activity recognition |
CN105547298A (zh) * | 2015-12-25 | 2016-05-04 | 北京京元智慧应急技术有限公司 | 智能手机内置传感器与Wi-Fi融合的室内动态连续定位方法 |
CN106767810A (zh) * | 2016-11-23 | 2017-05-31 | 武汉理工大学 | 一种基于移动终端的wifi和视觉信息的室内定位方法及系统 |
CN106705968A (zh) * | 2016-12-09 | 2017-05-24 | 北京工业大学 | 基于姿态识别和步长模型的室内惯性导航算法 |
Non-Patent Citations (4)
Title |
---|
QINGLIN TIAN等: "A Multi-Mode Dead Reckoning System for Pedestrian Tracking Using Smartphones", 《IEEE SENSORS JOURNAL》 * |
SHENG GUO等: "Activity Recognition and Semantic Description for Indoor Mobile Localization", 《SENSORS (BASEL)》 * |
张玲: "基于移动设备的多状态行人定位系统研究", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
韦丹: "基于移动终端的室内关键定位技术研究", 《中国优秀博硕士学位论文全文数据库(硕士) 信息科技辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109597031A (zh) * | 2018-11-19 | 2019-04-09 | 桂林电子科技大学 | 一种基于svm和步行者航位推算融合的室内定位方法 |
CN109827568A (zh) * | 2019-01-29 | 2019-05-31 | 东北大学秦皇岛分校 | 基于mems传感器的多层建筑中行人高度位置估计方法 |
CN110132276A (zh) * | 2019-04-30 | 2019-08-16 | 杭州电子科技大学 | 一种基于行人运动状态的自适应步长估计方法 |
CN110132276B (zh) * | 2019-04-30 | 2020-12-08 | 杭州电子科技大学 | 一种基于行人运动状态的自适应步长估计方法 |
CN110501008A (zh) * | 2019-07-29 | 2019-11-26 | 西安电子科技大学 | 人体运动模型的自主进化方法 |
CN110501008B (zh) * | 2019-07-29 | 2021-03-26 | 西安电子科技大学 | 人体运动模型的自主进化方法 |
CN112833889A (zh) * | 2019-11-22 | 2021-05-25 | 北京初速度科技有限公司 | 一种车辆的定位方法及装置 |
CN112833889B (zh) * | 2019-11-22 | 2022-12-23 | 北京魔门塔科技有限公司 | 一种车辆的定位方法及装置 |
CN113709860A (zh) * | 2021-08-12 | 2021-11-26 | 电子科技大学 | 一种层间过渡区域的室内定位装置和方法 |
CN113709860B (zh) * | 2021-08-12 | 2022-06-03 | 电子科技大学 | 一种层间过渡区域的室内定位装置和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108537101B (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108537101A (zh) | 一种基于状态识别的行人定位方法 | |
CN103968827B (zh) | 一种可穿戴式人体步态检测的自主定位方法 | |
CN104897154B (zh) | 基于手机的室内地磁导航方法及装置 | |
CN104897161B (zh) | 基于激光测距的室内平面地图制图方法 | |
CN104713554B (zh) | 一种基于mems惯性器件与安卓智能手机融合的室内定位方法 | |
CN105004340B (zh) | 结合惯性导航技术和指纹定位技术的定位误差修正方法 | |
CN102128625B (zh) | 重力辅助惯性导航系统中重力图匹配的初始匹配方法 | |
CN107830862A (zh) | 一种基于智能手机的室内定位行人追踪的方法 | |
CN108225304A (zh) | 基于多源传感器室内快速定位方法与系统 | |
CN108844533A (zh) | 一种基于多传感器融合和姿态解算的自由姿态pdr定位方法 | |
CN107218938A (zh) | 基于人体运动模型辅助的穿戴式行人导航定位方法和设备 | |
CN103791902B (zh) | 适用于高机动载体的星敏感器自主导航方法 | |
CN105716605B (zh) | 一种重力辅助惯性导航系统匹配方法 | |
CN107421535A (zh) | 一种基于磁场特征和加速度信息计步的室内行人定位系统 | |
CN103217154B (zh) | 一种煤矿井下人员定位方法及装置 | |
CN107179079A (zh) | 基于pdr与地磁融合的室内定位方法 | |
CN105043380A (zh) | 基于微机电传感器、WiFi定位、磁场匹配的室内导航方法 | |
CN106767828A (zh) | 一种手机室内定位解决方法 | |
CN104380043B (zh) | 在用户携带的便携式电子设备中的信息确定 | |
CN106705959A (zh) | 检测移动终端航向的方法和装置 | |
CN109029450A (zh) | 一种室内定位方法 | |
CN109211229A (zh) | 一种基于手机传感器及WiFi特征的人员室内定位方法 | |
CN106352879B (zh) | 一种基于图优化的uwb定位与编码器融合的位姿估计方法 | |
Wu et al. | Indoor positioning system based on inertial MEMS sensors: Design and realization | |
CN109637126A (zh) | 一种基于v2x终端的交通对象识别系统及其方法 |
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 |