发明名称 一种基于建筑物规则结构信息的室内定位方法
摘要 本发明公开了一种基于建筑物规则结构信息的室内定位方法,首先通过一般惯性传感器定位算法解算得到初级定位结果,然后根据矩形结构建筑物内的八个先验方向在粒子滤波算法框架内调整粒子权值,校正航向漂移,最后经过重采样可以得到轨迹修正后的定位结果,其中八个先验方向包含矩形四个边的朝向与四个经插值得到的辅助方向。本发明可利用建筑物的矩形规则结构对定位结果进行航向校正,显著提高了定位精度。本发明所提出的方法稳健性好,在行走轨迹复杂的条件下,仍能准确估计航向,抑制了航向漂移,提高了长时间室内定位的精确度。
申请公布号 CN104197936A 申请公布日期 2014.12.10
申请号 CN201410457534.7 申请日期 2014.09.10
申请人 长沙格致电子科技有限公司 发明人 王雅敏
分类号 G01C21/16(2006.01)I;G01C25/00(2006.01)I;G06F19/00(2011.01)I 主分类号 G01C21/16(2006.01)I
代理机构 长沙正奇专利事务所有限责任公司 43113 代理人 郭立中
主权项 一种基于建筑物规则结构信息的室内定位方法,包括以下步骤:1)确定导航坐标系,确定室内运动目标的初始位置与初始朝向;2)根据建筑结构的朝向,确定矩形四个边在导航坐标系下的初始角度ψ<sub>1</sub>、ψ<sub>3</sub>、ψ<sub>5</sub>与ψ<sub>7</sub>,将所述初始角度两两进行线性插值,得到四个辅助方向的角度ψ<sub>2</sub>、ψ<sub>4</sub>、ψ<sub>6</sub>与ψ<sub>8</sub>;3)利用惯性传感器测量目标的运动数据,并根据运动数据解算得到目标位移和角度的变化量,所述变化量逐次累加得到目标的相对方位,相对方位与初始方位相加得到初级定位结果,在所述导航坐标系下表示为{x<sub>k</sub>,y<sub>k</sub>},下坐标k表示第k步;4)在初级定位结果中,通过当前位置{x<sub>k</sub>,y<sub>k</sub>}与前一位置{x<sub>k‑1</sub>,y<sub>k‑1</sub>}得到目标当前位置步长l<sub>k</sub>与转向角△θ<sub>k</sub>:△θ<sub>k</sub>=θ<sub>k</sub>‑θ<sub>k‑1</sub><maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>l</mi><mi>k</mi></msub><mo>=</mo><msqrt><msup><msub><mi>&Delta;x</mi><mi>k</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>&Delta;y</mi><mi>x</mi></msub><mn>2</mn></msup></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000568047030000011.GIF" wi="390" he="99" /></maths>其中,<img file="FDA0000568047030000012.GIF" wi="454" he="142" />△x<sub>k</sub>=x<sub>k</sub>‑x<sub>k‑1</sub>;△y<sub>k</sub>=y<sub>k</sub>‑y<sub>k‑1</sub>;5)将当前位置步长l<sub>k</sub>与转向角△θ<sub>k</sub>分别加上高斯噪声,得到当前位置步长与转向角的噪声粒子l<sub>noise</sub><sup>[1:N]</sup>与△θ<sub>noise</sub><sup>[1:N]</sup>;其中,N为粒子数目;6)利用当前位置步长与转向角的噪声粒子l<sub>noise</sub><sup>[1:N]</sup>与△θ<sub>noise</sub><sup>[1:N]</sup>和前一位置的粒子集{x<sub>k‑1</sub><sup>[1:N]</sup>,y<sub>k‑1</sub><sup>[1:N]</sup>,l<sub>k‑1</sub><sup>[1:N]</sup>,θ<sub>k‑1</sub><sup>[1:N]</sup>},得到新的粒子集{x<sub>k</sub><sup>[1:N]</sup>,y<sub>k</sub><sup>[1:N]</sup>,l<sub>k</sub><sup>[1:N]</sup>,θ<sub>k</sub><sup>[1:N]</sup>}:l<sub>k</sub><sup>[i]</sup>=l<sub>k</sub>+l<sub>noise</sub><sup>[i]</sup>θ<sub>k</sub><sup>[i]</sup>=θ<sub>k‑1</sub><sup>[i]</sup>+△θ<sub>k</sub>+△θ<sub>noise</sub><sup>[i]</sup>x<sub>k</sub><sup>[i]</sup>=x<sub>k‑1</sub><sup>[i]</sup>+l<sub>k</sub><sup>[i]</sup>cos(θ<sub>k</sub><sup>[i]</sup>);y<sub>k</sub><sup>[i]</sup>=y<sub>k‑1</sub><sup>[i]</sup>+l<sub>k</sub><sup>[i]</sup>sin(θ<sub>k</sub><sup>[i]</sup>)其中,1≤i≤N,表示第i个粒子;θ<sub>k</sub><sup>[i]</sup>则表示第k步时,第i个粒子的朝向;△θ<sub>k</sub>表示通过运动传感器计算得到的转角,即朝向变化量;△θ<sub>noise</sub>表示高斯噪声模型下带来的随机转角;x<sub>k‑1</sub><sup>[i]</sup>、y<sub>k‑1</sub><sup>[i]</sup>分别表示第k‑1个位置下第i个粒子的横坐标和纵坐标;7)根据上述步骤2)确定的矩形建筑内的八个方向,按照下式更新每个粒子的权值:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><msub><mi>W</mi><mi>k</mi></msub><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>8</mn></munderover><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><msup><mrow><mo>(</mo><msup><msub><mi>&theta;</mi><mi>k</mi></msub><mrow><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msup><mo>-</mo><msub><mi>&psi;</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>2</mn><msup><mi>&pi;&sigma;</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000568047030000021.GIF" wi="614" he="156" /></maths>其中,ψ<sub>j</sub>为建筑物的引导向、σ<sup>2</sup>为选用的高斯噪声模型的方差、W<sub>k</sub><sup>[i]</sup>为第i个粒子在第k步的权值;8)对每个粒子的权值分别进行归一化处理,得到归一化权值w<sub>k</sub><sup>[1:N]</sup>;9)根据归一化权值w<sub>k</sub><sup>[1:N]</sup>,按照剩余重采样方法对所有粒子进行重采样,得到位置的粒子{x<sub>k</sub><sup>'[1:N]</sup>,y<sub>k</sub><sup>'[1:N]</sup>};10)对上述位置的粒子{x<sub>k</sub><sup>'[1:N]</sup>,y<sub>k</sub><sup>'[1:N]</sup>}取均值,得到航向校正后的位置:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mover><msub><mi>x</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><msub><mi>x</mi><mi>k</mi></msub><mrow><mo>&prime;</mo><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msup><mo>.</mo></mrow>]]></math><img file="FDA0000568047030000022.GIF" wi="351" he="153" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><msub><mi>y</mi><mi>k</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><msub><mi>y</mi><mi>k</mi></msub><mrow><mo>&prime;</mo><mo>[</mo><mi>i</mi><mo>]</mo></mrow></msup></mrow>]]></math><img file="FDA0000568047030000023.GIF" wi="316" he="134" /></maths>
地址 410008 湖南省长沙市开福区三一大道万煦园C9栋三门405