CN111366148B - 适用于机载光电观瞄系统多次观察的目标定位方法 - Google Patents
适用于机载光电观瞄系统多次观察的目标定位方法 Download PDFInfo
- Publication number
- CN111366148B CN111366148B CN202010230494.8A CN202010230494A CN111366148B CN 111366148 B CN111366148 B CN 111366148B CN 202010230494 A CN202010230494 A CN 202010230494A CN 111366148 B CN111366148 B CN 111366148B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- axis
- target
- geocentric
- photoelectric
- 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.)
- Active
Links
Images
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明属于机载光电侦察技术领域,具体涉及一种适用于机载光电观瞄系统多次观察的目标定位方法,该方法基于机载光电观瞄系统的预先校靶误差参数,机载光电观瞄系统位置和姿态测量单元的位置误差和姿态误差统计信息,机载光电观瞄系统的姿态误差统计信息,机载光电观瞄系统激光测距单元的测距误差统计信息,通过建立合适的卡尔曼滤波模型,可用于实时修正目标定位过程中上述各项参数误差对目标定位结果的影响,该方法可实现对静止目标的高精度地理定位。本发明对目标定位结果实时滤波,在现有光电系统的基础上不需要增加任何硬件资源,只需要增加适用于机载光电观瞄系统多次观察的目标定位软件包便可提升机载光电观瞄系统对静止目标的定位精度。
Description
技术领域
本发明属于机载光电侦察技术领域,具体涉及一种适用于机载光电观瞄系统多次观察的目标定位方法,该方法基于机载光电观瞄系统预先校靶参数建立健全的卡尔曼滤波模型,通过对目标多次观察过程中各项输入参数的随机误差进行修正,可实时结算滤波后的目标地理位置信息。
背景技术
当前的信息化战争对目标侦察提出了更高要求,在对感兴趣目标进行自动跟踪的同时,还需提供目标高精度的位置信息,如何能够准确的求解目标信息成为实施精确打击的首要解决问题。
机载光电观瞄系统目标定位过程中存在如下误差源:1)现有机载光电观瞄系统大多配有位置和姿态测量单元,并采取位置和姿态测量单元与光电平台固连的方式进行联接,这种方式与传统光电观瞄系统(未配有位置和姿态测量单元)通过减震器与载机组合导航系统联接相比,可消除载机组合导航系统与光电观瞄系统数据通信间的非同步误差,并有效减小有减震器带来的随机安装误差角,但同时引入了位置和姿态测量单元与光电平台的固定安装误差角;2)机载光电观瞄系统多采用两轴多框架的系统结构,光机状态无法完全保证两个转轴的垂直度,其将导致系统存在固定的转轴非正交误差角;3)系统对目标多次观察过程中存在各项输入参数的随机测量误差。
以上各误差源是机载光电观瞄系统目标定位误差的主要来源。
中国专利CN201410590112《系统无人机搜索目标快速定位方法》中公开了一种无人机搜索目标快速定位方法。该方法采用一架搜索无人机搜查旷野中目标,通过自身无线系统传输可视图像进行目标锁定,通过特定航行轨迹配合无人机多点视角切换和所处的姿态信息、高度信息和GPS定位信息进行多点定位迭代计算来快速获得旷野目标的具体定位信息。该专利可快速获取目标点定位信息,但并未考虑修正定位过程中存在的各种误差源,其降导致目标定位结果存在固定偏差和随机偏差,即最终获取的定位信息精度较低。
期刊《基于EKF的机载光电吊舱目标定位研究》为提高机载光电吊舱目标定位的精度,研究了光电吊舱测距、高低角和方位角的测量误差对目标定位精度的影响,推导并建立了目标定位的扩展卡尔曼滤波模型,通过该滤波模型抑制光电吊舱测距误差、高低角误差和方位角误差对目标点定位结果的影响来实现对目标定位精度的提升。该文献考虑了光电吊舱测距、高低角和方位角的测量误差对目标定位结果的影响,切定位滤波结果收敛速度较快,但其并未考虑载机姿态和位置的测量误差,同时忽略了光电吊舱与载机之间的安装误差和光电吊舱高低角与方位角之间的非正交误差,这些因素都想导致该文献提出的滤波模型与实际使用环境差异较大,模型的不准确性将直接导致滤波后的定位结果极易发散。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何提出一种适用于机载光电观瞄系统多次观察的目标定位方法。
(二)技术方案
为解决上述技术问题,本发明提供一种适用于机载光电观瞄系统多次观察的目标定位方法,所述方法包括以下步骤:
步骤1:建立地心坐标系e、地理坐标系g、光电平台坐标系b、瞄线坐标系s、瞄线俯仰零位坐标系s’;
步骤2:实时采集当前测距点数据组(αP,βP,γP,λP,LP,hP, R);其中:αP、βP、γP分别为在光电平台坐标系b下,位置和姿态测量单元输出的当前测距点载机俯仰角、横滚角和航向角;λP、LP、hP分别为在地心坐标系下,位置和姿态测量单元输出的当前测距点载机经度、纬度和高度信息;分别为在瞄线坐标系s下,机载光电观瞄系统在当前测距点输出的目标瞄线俯仰角和横滚角;R为机载光电观瞄系统在当前测距点输出的目标至光电观瞄系统的激光测距值;
其中:
上式中,为载机P在地心坐标系下的坐标,为目标T在瞄线坐标系下的坐标;RN为机载光电观瞄系统所在点的地理垂线与WGS-84参考旋转椭球体交点处的主曲率半径,Re为WGS-84参考旋转椭球体的半长轴,e为WGS-84参考旋转椭球体的扁率,e=1/298.257;
上式中, 分别为利用步骤2采集的当前时刻数据组数据计算所得的各旋转矩阵和载机P在地心坐标系下的坐标,为当前时刻采用的目标T在地心坐标系下的坐标;为目标T在瞄线俯仰零位坐标系s’下的估算值,xPT、yPT、zPT分别为向量的坐标值;
步骤6:采用以下公式定义卡尔曼滤波测量状态方程,选取地面点目标在地球直角坐标系下的位置误差作为状态变量,则对于地面点静止目标,其状态方程可表示为:
Xk=Φk/k-1Xk-1+Γk-1Wk-1
其中:
Φk/k-1=I3×3,Γk-1=ΔtI3×3
上式中,Xk、Xk-1分别为tk和tk-1时刻的目标点在地心坐标系下的位置误差;Δt为tk和tk-1时刻的时间差,即滤波周期;Φk/k-1为tk-1时刻至tk时刻的一步转移阵;Γk-1为系统噪声驱动阵;Wk-1为系统激励噪声序列,其方差用Qk表示;
步骤7:采用以下公式定义卡尔曼滤波测量等式:
步骤8:采用以下公式定义如下变量M1、M2、M3、M4,并利用步骤2采集的当前时刻数据组数据计算结果:
其中:
步骤9:采用以下公式对步骤7进行泰勒级数展开,并去掉高此项:
其中:
Gk=[G1 G2 G3 G4]T
上式中,Hk为量测阵;Gk为量测噪声驱动阵;Vk为量测噪声序列,其方差阵用Rk表示;ΔαP、ΔβP、ΔγP分别为位置和姿态测量单元输出的俯仰误差角、横滚误差角和航向误差角;ΔλP、ΔLP、ΔhP分别为位置和姿态测量单元输出的经度误差、纬度误差和高度误差; 分别为目标瞄线俯仰误差角和横滚误差角;ΔR为激光测距误差值;
步骤10:采用以下公式对状态变量进行递推估计:
Xk/k-1=Φk/k-1Xk-1
Pk=(I-KkHk)Pk/k-1
Xk=Xk/k-1+Kk(Zk-HkXk/k-1)
上式中,Xk/k-1为状态一步预测值;Xk为tk时刻的状态估计值;Xk-1为tk-1时刻的状态估计值;Kk为卡尔曼滤波增益;Pk/k-1为一步预测均方误差;Pk为估计均方误差;
步骤11:采用以下公式对目标T在地心坐标系下的坐标进行更新:
其中,所述地心坐标系根据描述的不同分为地心直角坐标系和地心球面坐标系。
其中,地心直角坐标系用o-xeyeze表示,原点o位于地心,xe、ye轴在地球赤道平面内,xe指向零度本初子午线,ze轴沿地球自转轴,ye轴穿越东经90°子午线与赤道的交点。
其中,地心球面坐标系的原点位于旋转椭球体中心,位置用经度、纬度和高度表示,即地球球面坐标系内坐标的表示方法为(λ,L,h)。
其中,所述地理坐标系g:用o-xgygzg表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xg轴指向正东方向,yg轴指向正北方向,zg轴指向天向。
其中,所述光电平台坐标系b:用o-xbybzb表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xb轴沿位置和姿态测量单元横轴向右,yb轴沿位置和姿态测量单元纵轴向前,zb轴沿位置和姿态测量单元立轴向上。
其中,所述瞄线坐标系s:用o-xsyszs表示,原点o位于光电平台回转中心,xs轴沿光电平台俯仰机构转轴向前,zs轴沿瞄线轴向上,ys轴与xs轴和zs轴构成右手坐标系。
其中,所述瞄线俯仰零位坐标系s’:用o-xs’ys’zs’表示,原点o位于光电平台回转中心,ys’轴沿光电平台俯仰机构转轴向前,zs’轴垂直于光电平台俯仰机构转轴和横滚机构转轴所在平面向上,xs’轴与ys’轴和zs’轴构成右手坐标系。
(三)有益效果
与现有技术相比较,本发明具备如下有益效果:
(1)本发明基于机载光电观瞄系统预先校靶所得的位置和姿态测量单元与光电平台的固定安装误差角、转轴非正交误差角,并基于机载光电观瞄系统位置和姿态测量单元的位置随机误差和姿态随机误差统计信息,机载光电观瞄系统的姿态随机误差统计信息,机载光电观瞄系统激光测距单元的测距随机误差统计信息,通过建立完善的卡尔曼滤波模型,可用于修正校靶所得的各固定误差对目标定位结果的影响,同时减弱上述各项随机误差对目标定位结果的影响,向用户实时提供高精度的目标定位结果。
(2)本发明通过软件算法对目标定位结果实时滤波,在现有光电系统的基础上不需要增加任何硬件资源,只需要增加适用于机载光电观瞄系统多次观察的目标定位软件包便可提升机载光电观瞄系统对静止目标的定位精度。
附图说明
图1是本发明目标定位方法工作流程图。
图2a-图2c是本发明滤波前后目标定位误差对比图。
具体实施方式
为使本发明的目的、内容、和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
为解决现有技术存在的问题,本发明提出一种适用于机载光电观瞄系统多次观察的目标定位方法,由目标定位软件模块实现,当机载光电观瞄系统跟踪上地面静止目标,且收到上级系统发出的目标定位指令时,该模块执行具体操作,并输出目标定位结果。
具体而言,本发明提供一种适用于机载光电观瞄系统多次观察的目标定位方法,如图1所示,所述方法包括以下步骤:
步骤1:建立地心坐标系e、地理坐标系g、光电平台坐标系b、瞄线坐标系s、瞄线俯仰零位坐标系s’;
步骤2:实时采集当前测距点数据组(αP,βP,γP,λP,LP,hP, R);其中:αP、βP、γP分别为在光电平台坐标系b下,位置和姿态测量单元输出的当前测距点载机俯仰角、横滚角和航向角;λP、LP、hP分别为在地心坐标系下,位置和姿态测量单元输出的当前测距点载机经度、纬度和高度信息;分别为在瞄线坐标系s下,机载光电观瞄系统在当前测距点输出的目标瞄线俯仰角和横滚角;R为机载光电观瞄系统在当前测距点输出的目标至光电观瞄系统的激光测距值;
其中:
上式中,为载机P在地心坐标系下的坐标,为目标T在瞄线坐标系下的坐标;RN为机载光电观瞄系统所在点的地理垂线与WGS-84参考旋转椭球体交点处的主曲率半径,Re为WGS-84参考旋转椭球体的半长轴,e为WGS-84参考旋转椭球体的扁率,e=1/298.257;
上式中,分别为利用步骤2采集的当前时刻数据组数据计算所得的各旋转矩阵和载机P在地心坐标系下的坐标,为当前时刻采用的目标T在地心坐标系下的坐标;为目标T在瞄线俯仰零位坐标系s’下的估算值,xPT、yPT、zPT分别为向量的坐标值;
步骤6:采用以下公式定义卡尔曼滤波测量状态方程,选取地面点目标在地球直角坐标系下的位置误差作为状态变量,则对于地面点静止目标,其状态方程可表示为:
Xk=Φk/k-1Xk-1+Γk-1Wk-1
其中:
Φk/k-1=I3×3,Γk-1=ΔtI3×3
上式中,Xk、Xk-1分别为tk和tk-1时刻的目标点在地心坐标系下的位置误差;Δt为tk和tk-1时刻的时间差,即滤波周期;Φk/k-1为tk-1时刻至tk时刻的一步转移阵;Γk-1为系统噪声驱动阵;Wk-1为系统激励噪声序列,其方差用Qk表示;
步骤7:采用以下公式定义卡尔曼滤波测量等式:
步骤8:采用以下公式定义如下变量M1、M2、M3、M4,并利用步骤2采集的当前时刻数据组数据计算结果:
其中:
步骤9:采用以下公式对步骤7进行泰勒级数展开,并去掉高此项:
其中:
Gk=[G1 G2 G3 G4]T
上式中,Hk为量测阵;Gk为量测噪声驱动阵;Vk为量测噪声序列,其方差阵用Rk表示;ΔαP、ΔβP、ΔγP分别为位置和姿态测量单元输出的俯仰误差角、横滚误差角和航向误差角;ΔλP、ΔLP、ΔhP分别为位置和姿态测量单元输出的经度误差、纬度误差和高度误差; 分别为目标瞄线俯仰误差角和横滚误差角;ΔR为激光测距误差值;
步骤10:采用以下公式对状态变量进行递推估计:
Xk/k-1=Φk/k-1Xk-1
Pk=(I-KkHk)Pk/k-1
Xk=Xk/k-1+Kk(Zk-HkXk/k-1)
上式中,Xk/k-1为状态一步预测值;Xk为tk时刻的状态估计值;Xk-1为tk-1时刻的状态估计值;Kk为卡尔曼滤波增益;Pk/k-1为一步预测均方误差;Pk为估计均方误差;
步骤11:采用以下公式对目标T在地心坐标系下的坐标进行更新:
其中,所述地心坐标系根据描述的不同分为地心直角坐标系和地心球面坐标系。
其中,地心直角坐标系用o-xeyeze表示,原点o位于地心,xe、ye轴在地球赤道平面内,xe指向零度本初子午线,ze轴沿地球自转轴,ye轴穿越东经90°子午线与赤道的交点。
其中,地心球面坐标系的原点位于旋转椭球体中心,位置用经度、纬度和高度表示,即地球球面坐标系内坐标的表示方法为(λ,L,h)。
其中,所述地理坐标系g:用o-xgygzg表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xg轴指向正东方向,yg轴指向正北方向,zg轴指向天向。
其中,所述光电平台坐标系b:用o-xbybzb表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xb轴沿位置和姿态测量单元横轴向右,yb轴沿位置和姿态测量单元纵轴向前,zb轴沿位置和姿态测量单元立轴向上。
其中,所述瞄线坐标系s:用o-xsyszs表示,原点o位于光电平台回转中心,xs轴沿光电平台俯仰机构转轴向前,zs轴沿瞄线轴向上,ys轴与xs轴和zs轴构成右手坐标系。
其中,所述瞄线俯仰零位坐标系s’:用o-xs’ys’zs’表示,原点o位于光电平台回转中心,ys’轴沿光电平台俯仰机构转轴向前,zs’轴垂直于光电平台俯仰机构转轴和横滚机构转轴所在平面向上,xs’轴与ys’轴和zs’轴构成右手坐标系。
实施例1
本实施例是针对机载光电观瞄系统多次观察的目标定位方法,该方法通过机载光电观瞄系统中的计算机实现。计算机带有存储器并装有目标定位软件模块。在载机飞行时,当机载光电观瞄系统跟踪上地面静止目标且收到上级系统发出的目标定位指令时,目标定位软件模块根据图1所示的流程执行以下操作步骤。
步骤1:建立地球坐标系e:所述地心坐标系根据描述的不同分为地心直角坐标系和地心球面坐标系。其中地心直角坐标系用o-xeyeze表示,原点o位于地心,xe、ye轴在地球赤道平面内,xe指向零度本初子午线,ze轴沿地球自转轴,ye轴穿越东经90°子午线与赤道的交点;地心球面坐标系的原点位于旋转椭球体中心,位置用经度、纬度和高度表示,即地球球面坐标系内坐标的表示方法为(λ,L,h);地理坐标系g:用o-xgygzg表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xg轴指向正东方向,yg轴指向正北方向,zg轴指向天向。光电平台坐标系b:用o-xbybzb表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xb轴沿位置和姿态测量单元横轴向右,yb轴沿位置和姿态测量单元纵轴向前,zb轴沿位置和姿态测量单元立轴向上。瞄线坐标系s:用o-xsyszs表示,原点o位于光电平台回转中心,xs轴沿光电平台俯仰机构转轴向前,zs轴沿瞄线轴向上,ys轴与xs轴和zs轴构成右手坐标系。瞄线俯仰零位坐标系s’:用o-xs’ys’zs’表示,原点o位于光电平台回转中心,ys’轴沿光电平台俯仰机构转轴向前,zs’轴垂直于光电平台俯仰机构转轴和横滚机构转轴所在平面向上,xs’轴与ys’轴和zs’轴构成右手坐标系。
步骤2:实时采集当前测距点数据组(αP,βP,γP,λP,LP,hP, R);其中:αP、βP、γP分别为在光电平台坐标系b下,位置和姿态测量单元输出的当前测距点载机俯仰角、横滚角和航向角;λP、LP、hP分别为在地球球面坐标系下,位置和姿态测量单元输出的当前测距点载机经度、纬度和高度信息;分别为在瞄线坐标系s下,机载光电观瞄系统在当前测距点输出的目标瞄线俯仰角和横滚角;R为机载光电观瞄系统在当前测距点输出的目标至光电观瞄系统的激光测距值。
在本优选实施例中,设定当前测距点真值为:
设定上述各项参数误差序列均为白噪声,其均值为0,均方根误差分别如下:
滤波次数:100次。
在本优选实施例中,设定当前系统校靶参数为:
其中:
上式中,为载机P在地球直角坐标系下的坐标,为目标T在瞄线坐标系下的坐标;RN为机载光电观瞄系统所在点的地理垂线与WGS-84参考旋转椭球体交点处的主曲率半径,Re为WGS-84参考旋转椭球体的半长轴,e为WGS-84参考旋转椭球体的扁率,e=1/298.257;
在本优选实施例中,计算得到的目标T在地球直角坐标系下的坐标真值为:
XT_real=[-1589994.96m 5050247.22m 3544416.88m]T,其对应地球球面坐标系下的坐标为:(107.4758303°,33.9779274°,0.62m);
将步骤二中的各项参数真值分别叠加其误差序列后,选取序列的第一组参数所求得的目标T在地球直角坐标系下的坐标作为卡尔曼滤波的首次输入值。
上式中,分别为利用步骤2采集的当前时刻数据组数据计算所得的各旋转矩阵和载机P在地球直角坐标系下的坐标,对于首次滤波,为选取各项输入参数叠加误差序列的第一组参数所求得的目标T在地球直角坐标系下的坐标,后续滤波中,为上一次滤波更新后的目标T在地球直角坐标系下的坐标值。
步骤6:采用以下算法定义卡尔曼滤波测量状态方程,选取地面点目标在地球直角坐标系下的位置误差作为状态变量,则对于地面点静止目标,其状态方程可表示为:
Xk=Φk/k-1Xk-1+Γk-1Wk-1
其中:
Φk/k-1=I3×3,Γk-1=ΔtI3×3
上式中,Xk、Xk-1分别为tk和tk-1时刻的目标点在地心坐标系下的位置误差;Δt为tk和tk-1时刻的时间差,即滤波周期;Φk/k-1为tk-1时刻至tk时刻的一步转移阵;Γk-1为系统噪声驱动阵;Wk-1为系统激励噪声序列,其方差用Qk表示;
步骤7:采用以下算法定义卡尔曼滤波测量等式输入:
步骤8:采用以下算法定义如下变量M1、M2、M3、M4,并利用步骤2采集的当前时刻数据组数据计算结果:
其中:
步骤9:采用以下算法对步骤7进行泰勒级数展开,并去掉高此项:
其中:
Gk=[G1 G2 G3 G4]T
上式中,Hk为量测阵;Gk为量测噪声驱动阵;Vk为量测噪声序列,其方差阵用Rk表示;ΔαP、ΔβP、ΔγP分别为位置和姿态测量单元输出的俯仰误差角、横滚误差角和航向误差角;ΔλP、ΔLP、ΔhP分别为位置和姿态测量单元输出的经度误差、纬度误差和高度误差; 分别为目标瞄线俯仰误差角和横滚误差角;ΔR为激光测距误差值;
步骤10:采用以下算法对状态变量进行递推估计:
Xk/k-1=Φk/k-1Xk-1
Pk=(I-KkHk)Pk/k-1
Xk=Xk/k-1+Kk(Zk-HkXk/k-1)
上式中,Xk/k-1为状态一步预测值;Xk为tk时刻的状态估计值;Xk-1为tk-1时刻的状态估计值;Kk为卡尔曼滤波增益;Pk/k-1为一步预测均方误差;Pk为估计均方误差;
在本优选实施例中,设定滤波估计均方误差初值P0为:
步骤11:采用以下算法在每一次滤波完成后,利用状态变量的估计值对目标T在地球直角坐标系下的坐标进行更新:
在本优选实施例中,未经本方法滤波前的目标点水平圆概率误差CEPT0和三维定位均方根误差σT0分别为:
CEPT0=59.84m;
σT0=93.03m;
经100次滤波后,目标点水平圆概率误差CEPT1和三维定位均方根误差σT1分别收敛至:
CEPT1=44.78m;
σT1=66.52m。
在本优选实施例中,滤波前后目标点在经度方向、纬度方向和高度方向误差对比图如图2a-图2c所示。由滤波结果和滤波前后定位误差对比图可知,本方法可有效修正目标定位过程中各项输入参数的随机误差项导致的定位误差,模型准确可靠。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。
Claims (6)
1.一种适用于机载光电观瞄系统多次观察的目标定位方法,其特征在于,所述方法包括以下步骤:
步骤1:建立地心坐标系e、地理坐标系g、光电平台坐标系b、瞄线坐标系s、瞄线俯仰零位坐标系s’;
步骤2:实时采集当前测距点数据组 其中:αP、βP、γP分别为在光电平台坐标系b下,位置和姿态测量单元输出的当前测距点载机俯仰角、横滚角和航向角;λP、LP、hP分别为在地心坐标系下,位置和姿态测量单元输出的当前测距点载机经度、纬度和高度信息;分别为在瞄线坐标系s下,机载光电观瞄系统在当前测距点输出的目标瞄线俯仰角和横滚角;R为机载光电观瞄系统在当前测距点输出的目标至光电观瞄系统的激光测距值;
其中:
上式中,为载机P在地心坐标系下的坐标,为目标T在瞄线坐标系下的坐标;RN为机载光电观瞄系统所在点的地理垂线与WGS-84参考旋转椭球体交点处的主曲率半径,Re为WGS-84参考旋转椭球体的半长轴,e为WGS-84参考旋转椭球体的扁率,e=1/298.257;
上式中,分别为利用步骤2采集的当前时刻数据组数据计算所得的各旋转矩阵和载机P在地心坐标系下的坐标,为当前时刻采用的目标T在地心坐标系下的坐标;为目标T在瞄线俯仰零位坐标系s’下的估算值,xPT、yPT、zPT分别为向量的坐标值;
步骤6:采用以下公式定义卡尔曼滤波测量状态方程,选取地面点目标在地球直角坐标系下的位置误差作为状态变量,则对于地面点静止目标,其状态方程可表示为:
Xk=Φk/k-1Xk-1+Γk-1Wk-1
其中:
Φk/k-1=I3×3,Γk-1=ΔtI3×3
上式中,Xk、Xk-1分别为tk和tk-1时刻的目标点在地心坐标系下的位置误差;Δt为tk和tk-1时刻的时间差,即滤波周期;Φk/k-1为tk-1时刻至tk时刻的一步转移阵;Γk-1为系统噪声驱动阵;Wk-1为系统激励噪声序列,其方差用Qk表示;
步骤7:采用以下公式定义卡尔曼滤波测量等式:
步骤8:采用以下公式定义如下变量M1、M2、M3、M4,并利用步骤2采集的当前时刻数据组数据计算结果:
其中:
步骤9:采用以下公式对步骤7进行泰勒级数展开,并去掉高此项:
其中:
Gk=[G1 G2 G3 G4]T
上式中,Hk为量测阵;Gk为量测噪声驱动阵;Vk为量测噪声序列,其方差阵用Rk表示;ΔαP、ΔβP、ΔγP分别为位置和姿态测量单元输出的俯仰误差角、横滚误差角和航向误差角;ΔλP、ΔLP、ΔhP分别为位置和姿态测量单元输出的经度误差、纬度误差和高度误差; 分别为目标瞄线俯仰误差角和横滚误差角;ΔR为激光测距误差值;
步骤10:采用以下公式对状态变量进行递推估计:
Xk/k-1=Φk/k-1Xk-1
Pk=(I-KkHk)Pk/k-1
Xk=Xk/k-1+Kk(Zk-HkXk/k-1)
上式中,Xk/k-1为状态一步预测值;Xk为tk时刻的状态估计值;Xk-1为tk-1时刻的状态估计值;Kk为卡尔曼滤波增益;Pk/k-1为一步预测均方误差;Pk为估计均方误差;
步骤11:采用以下公式对目标T在地心坐标系下的坐标进行更新:
其中,所述瞄线坐标系s:用o-xsyszs表示,原点o位于光电平台回转中心,xs轴沿光电平台俯仰机构转轴向前,zs轴沿瞄线轴向上,ys轴与xs轴和zs轴构成右手坐标系;
所述瞄线俯仰零位坐标系s’:用o-xs’ys’zs’表示,原点o位于光电平台回转中心,ys’轴沿光电平台俯仰机构转轴向前,zs’轴垂直于光电平台俯仰机构转轴和横滚机构转轴所在平面向上,xs’轴与ys’轴和zs’轴构成右手坐标系。
2.如权利要求1所述的适用于机载光电观瞄系统多次观察的目标定位方法,其特征在于,所述地心坐标系根据描述的不同分为地心直角坐标系和地心球面坐标系。
3.如权利要求2所述的适用于机载光电观瞄系统多次观察的目标定位方法,其特征在于,地心直角坐标系用o-xeyeze表示,原点o位于地心,xe、ye轴在地球赤道平面内,xe指向零度本初子午线,ze轴沿地球自转轴,ye轴穿越东经90°子午线与赤道的交点。
4.如权利要求3所述的适用于机载光电观瞄系统多次观察的目标定位方法,其特征在于,地心球面坐标系的原点位于旋转椭球体中心,位置用经度、纬度和高度表示,即地球球面坐标系内坐标的表示方法为(λ,L,h)。
5.如权利要求4所述的适用于机载光电观瞄系统多次观察的目标定位方法,其特征在于,所述地理坐标系g:用o-xgygzg表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xg轴指向正东方向,yg轴指向正北方向,zg轴指向天向。
6.如权利要求5所述的适用于机载光电观瞄系统多次观察的目标定位方法,其特征在于,所述光电平台坐标系b:用o-xbybzb表示,原点o位于机载光电观瞄系统位置和姿态测量单元质心,xb轴沿位置和姿态测量单元横轴向右,yb轴沿位置和姿态测量单元纵轴向前,zb轴沿位置和姿态测量单元立轴向上。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010230494.8A CN111366148B (zh) | 2020-03-27 | 2020-03-27 | 适用于机载光电观瞄系统多次观察的目标定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010230494.8A CN111366148B (zh) | 2020-03-27 | 2020-03-27 | 适用于机载光电观瞄系统多次观察的目标定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111366148A CN111366148A (zh) | 2020-07-03 |
CN111366148B true CN111366148B (zh) | 2022-11-22 |
Family
ID=71204851
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010230494.8A Active CN111366148B (zh) | 2020-03-27 | 2020-03-27 | 适用于机载光电观瞄系统多次观察的目标定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111366148B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112268558B (zh) * | 2020-10-15 | 2022-10-14 | 西北工业大学 | 一种基于无人机自主式光电平台的有源目标定位测速方法 |
CN112326199B (zh) * | 2020-10-29 | 2022-09-23 | 中国航空工业集团公司洛阳电光设备研究所 | 一种机载光学设备安装的视轴误差验证方法及验证系统 |
CN112672281B (zh) * | 2020-12-24 | 2022-11-01 | 中航贵州飞机有限责任公司 | 一种利用无线电链路装订光电侦察设备安装误差的方法 |
CN113238072B (zh) * | 2021-01-20 | 2022-11-08 | 西安应用光学研究所 | 一种适用于车载光电平台的运动目标解算方法 |
CN113252073A (zh) * | 2021-05-12 | 2021-08-13 | 无锡商业职业技术学院 | 一种应用于目标定位系统的现场标定方法及装置 |
CN113409400A (zh) * | 2021-06-17 | 2021-09-17 | 西安应用光学研究所 | 一种基于自动跟踪的机载光电系统目标地理定位方法 |
CN113686299B (zh) * | 2021-08-19 | 2023-05-02 | 中国科学院西安光学精密机械研究所 | 一种海上动态目标定位与移速预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5852792A (en) * | 1996-10-03 | 1998-12-22 | Lockheed Martin Corporation | Spacecraft boresight calibration filter |
US6064942A (en) * | 1997-05-30 | 2000-05-16 | Rockwell Collins, Inc. | Enhanced precision forward observation system and method |
CN104316949A (zh) * | 2014-10-29 | 2015-01-28 | 上海大学 | 一种无人机搜索目标快速定位方法 |
CN108535715A (zh) * | 2018-04-12 | 2018-09-14 | 西安应用光学研究所 | 一种适用于机载光电观瞄系统的大气折射下目标定位方法 |
CN109564285A (zh) * | 2016-07-29 | 2019-04-02 | 大众汽车有限公司 | 用于检测在移动单元的交通环境下的地面标志的方法和系统 |
CN110657808A (zh) * | 2019-10-16 | 2020-01-07 | 中国人民解放军火箭军工程大学 | 一种机载光电吊舱有源目标定位方法及系统 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3057656B1 (fr) * | 2016-10-14 | 2019-04-05 | Sofresud | Dispositif visuel pour la designation d’objectifs et procede de designation d’objectif utilisant ledit dispositif |
US10690525B2 (en) * | 2018-01-03 | 2020-06-23 | General Electric Company | Systems and methods associated with unmanned aerial vehicle targeting accuracy |
-
2020
- 2020-03-27 CN CN202010230494.8A patent/CN111366148B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5852792A (en) * | 1996-10-03 | 1998-12-22 | Lockheed Martin Corporation | Spacecraft boresight calibration filter |
US6064942A (en) * | 1997-05-30 | 2000-05-16 | Rockwell Collins, Inc. | Enhanced precision forward observation system and method |
CN104316949A (zh) * | 2014-10-29 | 2015-01-28 | 上海大学 | 一种无人机搜索目标快速定位方法 |
CN109564285A (zh) * | 2016-07-29 | 2019-04-02 | 大众汽车有限公司 | 用于检测在移动单元的交通环境下的地面标志的方法和系统 |
CN108535715A (zh) * | 2018-04-12 | 2018-09-14 | 西安应用光学研究所 | 一种适用于机载光电观瞄系统的大气折射下目标定位方法 |
CN110657808A (zh) * | 2019-10-16 | 2020-01-07 | 中国人民解放军火箭军工程大学 | 一种机载光电吊舱有源目标定位方法及系统 |
Non-Patent Citations (1)
Title |
---|
机载光电系统的地面多目标定位算法;闫明等;《应用光学》;20120731(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111366148A (zh) | 2020-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111366148B (zh) | 适用于机载光电观瞄系统多次观察的目标定位方法 | |
CN106871927B (zh) | 一种无人机光电吊舱安装误差标校方法 | |
CN112197761B (zh) | 一种高精度多旋翼机协同定位方法及系统 | |
US8666661B2 (en) | Video navigation | |
CN105698762A (zh) | 一种单机航迹上基于不同时刻观测点的目标快速定位方法 | |
CN106595668A (zh) | 一种用于光电吊舱的无源定位算法 | |
CN111102981B (zh) | 一种基于ukf的高精度卫星相对导航方法 | |
CN113252038B (zh) | 基于粒子群算法的航迹规划地形辅助导航方法 | |
CN111998855B (zh) | 光学望远镜共视观测确定空间目标初轨的几何方法及系统 | |
CN109633720B (zh) | 基于视频卫星的地面运动目标测量方法及其装置 | |
CN104360362A (zh) | 利用飞行器对观测目标进行定位的方法及系统 | |
CN112346104A (zh) | 一种无人机信息融合定位方法 | |
CN117455960B (zh) | 时变观测噪声条件下机载光电系统对地无源定位滤波方法 | |
CN112731281A (zh) | 一种空间碎片测角数据仿真方法 | |
CN116105730A (zh) | 基于合作目标卫星甚短弧观测的仅测角光学组合导航方法 | |
Mostafa et al. | Optical flow based approach for vision aided inertial navigation using regression trees | |
CN112833878B (zh) | 一种近地面多源天文自主导航方法 | |
CN113074715A (zh) | 无人机高精度目标定位方法 | |
CN112229421A (zh) | 基于李群最优估计的捷联惯性导航晃动基座粗对准方法 | |
CN114485624B (zh) | 一种基于恒星与卫星联合的全天时光学导航方法和装置 | |
CN115388890A (zh) | 基于视觉的多无人机协同对地目标定位方法 | |
CN113252039B (zh) | 面向地形辅助导航的粒子群快速匹配方法 | |
CN116309798A (zh) | 一种无人机成像定位方法 | |
CN112729305B (zh) | 一种基于单飞行器导引头图像信息的多目标定位方法 | |
CN112268558B (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 |