发明名称 多GPS/INS与交通图像融合定位方法
摘要 本发明公开了一种多GPS/INS与交通图像融合定位方法,属于信息融合技术与智能信息处理技术领域,用于对车辆进行精确定位。本发明将多GPS/INS以及交通图像视频等路况信息相结合,经融合处理得到等效测量值,带入系统的状态方程进行滤波,最终得到车辆的精确定位信息。当图像采集到的基准点信息时,直接标定该点位置的测量值,这样系统误差即可归零,减小了系统误差,提高了定位精确度,同时通过估计补充一些特殊路段(例如隧道)的数据信息。
申请公布号 CN101739840A 申请公布日期 2010.06.16
申请号 CN200910219169.5 申请日期 2009.11.26
申请人 西北工业大学 发明人 史忠科;王慧丽
分类号 G08G1/123(2006.01)I;G08G1/017(2006.01)I 主分类号 G08G1/123(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 黄毅新
主权项 1.一种多GPS/INS与交通图像融合定位方法,其特征在于包括下述步骤:(a)利用GPS、INS测量装置以及图像采集系统,对路况信息进行采集并记录,选取经度、纬度和高度作为状态变量,得到状态方程<img file="F2009102191695C0000011.GIF" wi="800" he="216" />式中,λ,φ,h分别是经度、纬度和高度,V<sub>E</sub>,V<sub>N</sub>,V<sub>U</sub>分别是GPS所测量到的东、北、天的速度,<img file="F2009102191695C0000012.GIF" wi="800" he="93" />是子午圈的曲率半径,<img file="F2009102191695C0000013.GIF" wi="800" he="98" />是卯酉圈的曲率半径;R<sub>e</sub>是地球的长半轴,ε为偏率;以及观测方程<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>X</mi><mi>G</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mi>h</mi><mo>)</mo></mrow><mi>cos</mi><mi></mi><mi>&phi;</mi><mi>cos</mi><mi>&lambda;</mi></mtd></mtr><mtr><mtd><mrow><msub><mi>Y</mi><mi>G</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>N</mi><mo>+</mo><mi>h</mi><mo>)</mo></mrow><mi>cos</mi><mi></mi><mi>&phi;</mi><mi>sin</mi><mi>&lambda;</mi></mrow></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>G</mi></msub><mo>=</mo><mo>[</mo><mi>N</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>e</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mi>h</mi><mo>]</mo><mi>sin</mi><mi>&phi;</mi></mtd></mtr></mtable></mfenced></math>]]></maths>式中,X,Y,Z是GPS接收机天线所对应的WGS-84坐标,<img file="F2009102191695C0000015.GIF" wi="363" he="122" />是当地卯酉圈半径,<img file="F2009102191695C0000016.GIF" wi="257" he="134" />表示地球的第一偏心率,a=6378137.000m,b=6356752.314m分别为地球长、短半轴长;(b)对不同噪声进行信息融合处理,其融合估计值是<maths num="0002"><![CDATA[<math><mrow><mover><mi>z</mi><mo>^</mo></mover><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&Lambda;</mi><mi>i</mi></msub><msub><mi>z</mi><mi>i</mi></msub></mrow></math>]]></maths>z<sub>i</sub>(i=1,2,…,m)为m个不同测量值,v<sub>i</sub>是测量误差,服从零均值方差阵R<sub>i</sub>(i=1,2,…,m)的白噪声,对于高斯白噪声<maths num="0003"><![CDATA[<math><mrow><msub><mi>&Lambda;</mi><mi>k</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>R</mi><mi>i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><msubsup><mi>R</mi><mi>k</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>R</mi><mi>k</mi></msub><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msubsup><mi>R</mi><mi>i</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>,</mo><mrow><mo>(</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msub><mi>&Lambda;</mi><mi>k</mi></msub><mo>=</mo><mi>I</mi><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&Lambda;</mi><mi>i</mi></msub></mrow></math>]]></maths>对于相关噪声,v<sub>i</sub>与v<sub>j</sub>的协方差阵为R<sub>ij</sub>,(i,j=1,2,…m)Λ=G<sup>-1</sup>U<maths num="0005"><![CDATA[<math><mrow><mi>&Lambda;</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>&Lambda;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&Lambda;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>&Lambda;</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><img file="F2009102191695C0000022.GIF" wi="742" he="296" /><maths num="0006"><![CDATA[<math><mrow><mi>U</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>-</mo><msub><mi>R</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>-</mo><msub><mi>R</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>G<sub>ik</sub>=R<sub>ik</sub>-R<sub>im</sub>-R<sub>km</sub>+R<sub>m</sub>,(i,k=1,2,…m-1)对于有色噪声,v<sub>i</sub>与v<sub>j</sub>的协方差阵为R<sub>ij</sub>Λ=(G′)<sup>-1</sup>U′<maths num="0007"><![CDATA[<math><mrow><mi>&Lambda;</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>&Lambda;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&Lambda;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>&Lambda;</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths><img file="F2009102191695C0000025.GIF" wi="784" he="346" /><maths num="0008"><![CDATA[<math><mrow><msup><mi>U</mi><mo>&CenterDot;</mo></msup><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>-</mo><msub><mi>R</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>-</mo><msub><mi>R</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>R</mi><mi>m</mi></msub><mo>-</mo><msub><mi>R</mi><mrow><mi>m</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow></math>]]></maths>G<sub>ik</sub>′=P<sub>ik</sub>-P<sub>im</sub>-P<sub>km</sub>+P<sub>mm</sub>(i,k=1,2,…m-1)<maths num="0009"><![CDATA[<math><mrow><msub><mi>P</mi><mi>ij</mi></msub><mo>=</mo><mfrac><mrow><msup><msub><mi>&lambda;</mi><mi>ij</mi></msub><mn>2</mn></msup><mo>+</mo><msup><msub><mi>&mu;</mi><mi>ij</mi></msub><mn>2</mn></msup></mrow><mn>2</mn></mfrac><msub><mi>A</mi><mi>ii</mi></msub><mo>+</mo><mfrac><mrow><msubsup><mi>&lambda;</mi><mi>ij</mi><mn>2</mn></msubsup><mo>+</mo><mn>1</mn><mo>/</mo><msup><msub><mi>&mu;</mi><mi>ij</mi></msub><mn>2</mn></msup></mrow><mn>2</mn></mfrac><msub><mi>A</mi><mi>jj</mi></msub><mo>+</mo><msub><mi>R</mi><mi>ij</mi></msub><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><mfrac><msub><mi>R</mi><mi>i</mi></msub><msubsup><mi>&lambda;</mi><mi>ji</mi><mn>2</mn></msubsup></mfrac><mo>+</mo><mfrac><msub><mi>R</mi><mi>j</mi></msub><msup><msub><mi>&lambda;</mi><mi>ij</mi></msub><mn>2</mn></msup></mfrac><mo>)</mo></mrow></mrow></math>]]></maths>式中,λ<sub>ji</sub>,μ<sub>ij</sub>为任意非零实数;(c)GPS、INS的测量信息与交通图像信息融合处理后,对交通图像检测到的基准点估计值进行标定,将所检测到的基准点的精确位置作为所估计位置的准确值;(d)将步骤(b)所得等效测量值与步骤(c)基准点的标定信息代入状态方程进行滤波,得到车辆的精确定位信息,实现对车辆的精确定位。
地址 710072 陕西省西安市友谊西路127号