发明名称 一种飞行器前视末制导导航定位方法
摘要 本发明公开了一种飞行器前视导航定位方法,包括两个阶段:远距离导航阶段,首先识别地标形心在实时图中的位置,再结合目标形心与地标形心在大地坐标系下的位置偏差以及当前飞行器成像姿态参数,计算目标形心在实时图中的估计位置,进而计算飞行器的相关位置偏差值,并以目标形心在实时图中的估计位置为中心选取目标搜索区域;距离导航阶段,在目标搜索区域内对目标直接定位,依据直接定位结果计算飞行器的相关位置偏差值。本发明通过目标周边区域的地标识别,达到对目标的间接导航定位,保证了飞行器导航定位的鲁棒性和自适应性,克服了因误差引起的飞行器偏离规划航迹,导致飞行过程中目标不易直接识别而无法导航定位的缺陷。
申请公布号 CN101629827B 申请公布日期 2010.12.01
申请号 CN200910063620.9 申请日期 2009.08.14
申请人 华中科技大学 发明人 张天序;万美君;杨效余;王登位;汪小平;敖国庆;李成;毛麟;孙思远;郑璐璐;王曼
分类号 G01C21/24(2006.01)I 主分类号 G01C21/24(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 李智
主权项 1.一种飞行器前视末制导导航定位方法,包括远距离导航阶段和近距离导航阶段,具体为:A远距离导航阶段(A1)识别地标形心在实时图中的位置(X<sub>2</sub>,Y<sub>2</sub>);(A2)依据目标形心与地标形心在大地坐标系下的位置偏差(dx,dy)、地标形心在实时图中的位置(X<sub>2</sub>,Y<sub>2</sub>)和当前飞行器成像姿态参数,计算目标形心在实时图中的估计位置(X<sub>1</sub>,Y<sub>1</sub>);(A3)依据目标形心在实时图中的估计位置(X<sub>1</sub>,Y<sub>1</sub>)和当前飞行器成像姿态参数,计算飞行器的俯仰偏差角度dθ、方位偏差角度<img file="FSB00000192190900011.GIF" wi="61" he="50" />和飞行器位置偏差(ΔX,ΔY,ΔZ);(A4)在实时图中确定目标的搜索区域,搜索区域的中心点位置为(X<sub>1</sub>,Y<sub>1</sub>),搜索区域的行数=目标成像行数Y_s+Δ像素,搜索区域的列数=目标成像列数X_s+Δ像素,8≤Δ≤12;B近距离导航阶段(B1)目标识别:对于第一帧图像:在远距离识别阶段得到的搜索区域内进行目标直接识别;对于其它帧图像:采用历史帧的目标识别结果指导当前帧图像的目标识别;(B2)依据步骤(B1)得到的目标识别位置和当前飞行器成像姿态参数,计算飞行器的俯仰偏差角度dθ、方位偏差角度<img file="FSB00000192190900012.GIF" wi="63" he="50" />和飞行器位置偏差(ΔX,ΔY,ΔZ);所述步骤(A2)具体为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>Y</mi><mn>1</mn></msub><mo>=</mo><mfrac><mi>ROW</mi><mn>2</mn></mfrac><mo>+</mo><mfrac><mi>ROW</mi><mi>&phi;</mi></mfrac><mo>&times;</mo><mrow><mo>(</mo><mi>arctan</mi><mfrac><mi>h</mi><mrow><mi>U</mi><mo>+</mo><mi>cos</mi><mi>&alpha;</mi><mo>&times;</mo><mi>dy</mi><mo>-</mo><mi>sin</mi><mi>&alpha;</mi><mo>&times;</mo><mi>dx</mi></mrow></mfrac><mo>-</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></math>]]></maths><img file="FSB00000192190900014.GIF" wi="1668" he="210" />其中<maths num="0002"><![CDATA[<math><mrow><mi>U</mi><mo>=</mo><mfrac><mi>h</mi><mrow><mi>tan</mi><mo>[</mo><mi>&theta;</mi><mo>+</mo><mfrac><mrow><mi>&phi;</mi><mo>&times;</mo><mrow><mo>(</mo><msub><mi>Y</mi><mn>2</mn></msub><mo>-</mo><mfrac><mi>ROW</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mrow><mi>ROW</mi></mfrac><mo>]</mo></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mi>V</mi><mo>=</mo><mfrac><mi>h</mi><mrow><mi>sin</mi><mo>[</mo><mi>&theta;</mi><mo>+</mo><mfrac><mrow><mi>&phi;</mi><mo>&times;</mo><mrow><mo>(</mo><msub><mi>Y</mi><mn>2</mn></msub><mo>-</mo><mfrac><mi>ROW</mi><mn>2</mn></mfrac><mo>)</mo></mrow></mrow><mi>ROW</mi></mfrac><mo>]</mo></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><img file="FSB00000192190900023.GIF" wi="603" he="170" />R=V×W,φ为成像仪纵向成像视场角,<img file="FSB00000192190900024.GIF" wi="32" he="39" />为成像仪横向成像角,ROW为实时成像行数,COL为实时成像列数,α为成像方位角,θ为成像俯仰角,h为成像高度;所述步骤(A3)具体为:dθ=φ×(Y<sub>1</sub>-ROW/2)/ROW,<img file="FSB00000192190900025.GIF" wi="684" he="90" />(ΔX,ΔY,ΔZ)=(x-x<sub>0</sub>,y-y<sub>0</sub>,z-z<sub>0</sub>),其中(x′,y′,z′)=(y<sub>1</sub>×sinα-x<sub>1</sub>×cosα,-y<sub>1</sub>×cosα-x<sub>1</sub>×sinα,h),y<sub>1</sub>=h/tan[θ+φ×(Y<sub>1</sub>-ROW/2)/ROW],<img file="FSB00000192190900026.GIF" wi="905" he="122" />(x<sub>0</sub>,y<sub>0</sub>,z<sub>0</sub>)为实际航迹曲线f(x,y,z)到当前飞行器位置点(x′,y′,z′)距离最短的点的位置;所述步骤(A4)中,对于平面目标,目标成像行数Y_s以及目标成像列数X_s的计算公式为:<maths num="0004"><![CDATA[<math><mrow><mi>Y</mi><mo>_</mo><mi>s</mi><mo>=</mo><mn>2</mn><mo>&times;</mo><mfrac><mi>ROW</mi><mi>&phi;</mi></mfrac><mo>&times;</mo><mrow><mo>(</mo><mi>arctan</mi><mfrac><mi>h</mi><mrow><mfrac><mi>h</mi><mrow><mi>tan</mi><mi>&theta;</mi></mrow></mfrac><mo>-</mo><mfrac><msup><mi>Width</mi><mo>&prime;</mo></msup><mn>2</mn></mfrac></mrow></mfrac><mo>-</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow></math>]]></maths><img file="FSB00000192190900032.GIF" wi="1088" he="272" />对于立体目标,目标成像行数Y_s以及目标成像列数X_s的计算公式为:<maths num="0005"><![CDATA[<math><mrow><mi>Y</mi><mo>_</mo><mi>s</mi><mo>=</mo><mn>2</mn><mo>&times;</mo><mfrac><mi>ROW</mi><mi>&phi;</mi></mfrac><mo>&times;</mo><mo>[</mo><mi>arctan</mi><mfrac><mi>h</mi><mrow><mfrac><mi>h</mi><mrow><mi>tan</mi><mi>&theta;</mi></mrow></mfrac><mo>-</mo><mfrac><msup><mi>Width</mi><mo>&prime;</mo></msup><mn>2</mn></mfrac></mrow></mfrac><mo>-</mo><mi>arctan</mi><mfrac><mrow><mrow><mo>(</mo><mi>h</mi><mo>-</mo><mi>Height</mi><mo>)</mo></mrow><mo>&times;</mo><mi>tan</mi><mi>&theta;</mi></mrow><mi>h</mi></mfrac><mo>]</mo></mrow></math>]]></maths><img file="FSB00000192190900034.GIF" wi="1032" he="298" />其中Length′=Width×sinα-Length×cosαWidth′=Length×sinα-Width×cosα;Length为正北方向下目标的长,Width为正北方向下目标的宽,Height为目标实际高度。
地址 430074 湖北省武汉市洪山区珞喻路1037号