发明名称 基于序列图像和基准图匹配的飞行器全参数导航方法
摘要 本发明涉及基于序列图像和基准图匹配的飞行器全参数导航方法。该方法在飞行器上安装摄像机、惯导系统、数据处理器件和存储器件,且数据存储器上存储有预定区域的地面基准图数据,构成基于景象匹配的飞行器全参数导航系统。通过实时图像与基准图的匹配,并与惯导进行组合导航的卡尔曼滤波,即可实现飞行器导航的高精度全参数测量,具有广泛、重要的应用前景。
申请公布号 CN102829785B 申请公布日期 2014.12.31
申请号 CN201210313150.9 申请日期 2012.08.30
申请人 中国人民解放军国防科学技术大学 发明人 于起峰;雷志辉;尚洋;刘晓春;李鑫;李强;李想;苏昂
分类号 G01C21/20(2006.01)I;G01C11/00(2006.01)I 主分类号 G01C21/20(2006.01)I
代理机构 湖南省国防科技工业局专利中心 43102 代理人 冯青
主权项 基于序列图像和基准图匹配的飞行器全参数导航方法,其特征在于,该方法实施过程是:1) 在飞行器上安装摄像机、惯导系统、数据处理器件和存储器件,且数据存储器上存储有预定区域的地面基准图数据,构成基于景象匹配的飞行器全参数导航系统;2) 飞行器进入预定区域后,开启基于景象匹配的飞行器全参数导航系统,摄像机对地面进行连续拍摄,采用飞行器位姿参数估计方法或飞行器全参数测量的组合导航方法,完成飞行器全参数的测量,具体按以下步骤进行:2.1) 在每帧实时图像上选取多个特征点,与带有大地水平坐标和高程的基准图进行匹配,得到特征点在大地绝对坐标系下的三维坐标;2.2) 利用匹配得到特征点图像坐标和三维大地绝对坐标,以及已知的摄像机内参数和安装角,求解摄像机成像模型,得到飞行器三维位置和姿态参数;2.3) 将解算出的飞行器位姿参数拟合出飞行器的速度、加速度、角速度、角加速度导航参数,或者将解算出的飞行器位姿参数和惯导进行组合导航,得到高精度的飞行器全参数导航测量结果;飞行器位姿参数估计方法具体为:在实时图上选取多个特征点与基准图匹配,得到特征点图像坐标以及对应的三维大地绝对坐标,根据摄像测量学知识,设地面上有<img file="2012103131509100001dest_path_image001.GIF" wi="14" he="16" />个点,其大地绝对坐标为<img file="16430dest_path_image002.GIF" wi="73" he="28" />,图像坐标为<img file="2012103131509100001dest_path_image003.GIF" wi="50" he="28" />,<img file="83000dest_path_image004.GIF" wi="20" he="25" />、<img file="dest_path_image005.GIF" wi="21" he="26" />为摄像机等效焦距,<img file="237005dest_path_image006.GIF" wi="21" he="25" />、<img file="dest_path_image007.GIF" wi="22" he="26" />为摄像机像主点坐标,<img file="566355dest_path_image008.GIF" wi="76" he="25" />为摄像机光心在大地绝对坐标系中的位置坐标,由中心投影成像模型,建立以下方程:<img file="dest_path_image009.GIF" wi="336" he="105" />(1)其中,大地绝对坐标系到摄像机坐标系的旋转矩阵为<img file="17453dest_path_image010.GIF" wi="18" he="18" />,则<img file="60233dest_path_image010.GIF" wi="18" he="18" />可由飞行器姿态角以及摄像机安装角表示,如下:<img file="dest_path_image011.GIF" wi="118" he="76" />(2)当有至少三个对应特征点对时,根据以上成像方程可以利用非线性迭代最优化的方法求解得到飞行器三个位置量和三个姿态角;飞行器全参数测量的组合导航方法包括方法,具体为:考虑姿态角误差、速度误差、位置误差、陀螺仪测量误差、加速度计测量误差,由惯导误差方程得到组合导航的误差状态方程如下:<img file="362251dest_path_image012.GIF" wi="158" he="25" />(5)其中:<img file="dest_path_image013.GIF" wi="553" he="31" /><img file="178898dest_path_image014.GIF" wi="258" he="33" />式中<img file="dest_path_image015.GIF" wi="14" he="22" />表示姿态角;<img file="230424dest_path_image016.GIF" wi="13" he="16" />表示飞行器速度;<img file="dest_path_image017.GIF" wi="14" he="16" />表示飞行器位置;<img file="3076dest_path_image018.GIF" wi="14" he="16" />表示陀螺仪测量误差;<img file="dest_path_image019.GIF" wi="17" he="20" />表示加速度计测量误差;下标<i>E</i><i>、</i><i>N</i><i>、</i><i>U</i>分别表示直角坐标系的东向、北向和天向,下标g表示陀螺仪,下标a表示加速度计;以解算出的飞行器位姿参数与惯导输出的位姿参数之差作为量测值,构建组合导航的量测方程;<img file="616942dest_path_image020.GIF" wi="104" he="25" />(6)其中:<img file="dest_path_image021.GIF" wi="452" he="30" />式中下标<i>I</i>表示惯性导航输出;下标<i>V</i>表示基于景象匹配的视觉导航输出,基于式(5)、(6)的组合导航为卡尔曼滤波方法,根据组合导航的滤波输出,得到高精度的飞行器导航的全参数测量结果如下:<img file="107835dest_path_image022.GIF" wi="92" he="100" />(7)其中:<img file="dest_path_image023.GIF" wi="122" he="28" />式中<img file="635156dest_path_image024.GIF" wi="17" he="16" />表示飞行器姿态角速度。
地址 410073 湖南省长沙市砚瓦池正街47号