一种采用面阵光学载荷的单星自主定位方法技术

本发明专利技术公开了一种采用面阵光学载荷的单星自主定位方法，包含以下步骤：S1、对地观测后获得一幅遥感图像，与地面标志点信息库匹配并提取出遥感图像中包含的所有地面标志点，从所有地面标志点中选取三个用于定位的地面标志点；S2、对三个用于定位的地面标志点分别进行空间几何位置建模，结合地面标志点坐标系、相机坐标系、像平面坐标系和地心惯性坐标系之间的几何关系，得到航天器定位模型；S3、采用预设处理算法对航天器定位模型进行求解,得到单星自主定位输出量,以完成对单星自主定位。本发明专利技术仅使用航天器对地观测载荷便可完成自主定位，从而使航天器容易地拟合出航天器轨道参数，使航天器测定轨摆脱外界条件，抗干扰性和可靠性得到极大增强。

【技术实现步骤摘要】

本专利技术航天器定轨研究
，具体涉及航天器测定轨技术，尤其涉及一种采用面阵光学载荷的单星自主定位方法。
技术介绍
随着各国航天技术的发展，在轨卫星数量急剧增加，使地面测控站负担过重，测控系统易发生信息阻塞，从而导致在轨航天器安全运行受到威胁。为保证航天器安全可靠并减轻地面负担，航天器应具有自主运行能力，这就要求航天器具有在轨高精度自主测定轨能力。考察现有的航天器测定轨方法，大多是以地基测定轨、天基测定轨、天文测定轨为技术实现途径。地基测定轨的精度最高，但自主性缺失，且地面测控站受限的布设直接影响了定位精度和轨道覆盖率；天基测定轨精度较高，有一定的自主性，但天基系统本身的轨道仍然需要地面测控站的观测和支持，且天基系统极易受到人为干扰；天文测定轨是通过观测自然天体获取航天器轨道信息的，其自主性极高，但观测精度受到地平仪等测量设备的影响，误差在几百米～几公里之间。
技术实现思路
本专利技术的目的在于提供一种采用面阵光学载荷的单星自主定位方法，仅使用航天器对地观测载荷便可完成自主定位，从而使航天器容易地拟合出航天器轨道参数，使航天器测定轨摆脱外界条件，抗干扰性和可靠性得到极大增强。为了达到上述目的，本专利技术通过以下技术方案实现：一种采用面阵光学载荷的单星自主定位方法，其特点是，包含以下步骤：S1、航天器通过面阵光学载荷对地观测后获得一幅遥感图像，经过图像预处理后，根据地面标志点信息库匹配并提取出遥感图像中包含的所有地面标志点，采用一预设提取算法，从所有地面标志点中选取三个用于定位的地面标志点；S2、建立地面标志点坐标系、相机坐标系、像平面坐标系及地心惯性坐标系，并对三个用于定位的地面标志点分别进行空间几何位置建模，结合地面标志点坐标系、相机坐标系、像平面坐标系和地心惯性坐标系之间的几何关系，推导可得基于面阵图像的航天器定位模型；S3、采用预设处理算法对航天器定位模型进行求解,得到单星自主定位输出量,以完成对单星自主定位。所述的步骤S1中预设提取算法包含：对所有地面标志点进行排列组合，形成若干个不同的观测组合，每一观测组合包含三个用于定位的地面标志点；对每种观测组合进行几何表述，获得几何表示矩阵；通过分析几何表述矩阵，对每种观测组合的几何分布进行分析和计算，并根据分析和计算结果选择最佳观测组合用于定位解算。所述的步骤S1中预设提取算法包含：记所有地面标志点的个数为N，则从中选取三个用于定位的地面标志点，则有种组合方式，共有M个观测组合；记任一观测组合中三个地面标志点为pi,pj,pk，则该观测组合对应的几何表述矩阵如式（1）所示；G=-cosθisinαi-cosθicosαi-sinθi-cosθjsinαj-cosθjcosαj-sinθj-cosθksinαk-cosθkcosαk-sinθk---(1)]]>式中，G表示观测组合对应的几何表述矩阵；θi,θj,θk为在航天器轨道坐标系下地面标志点与航天器的连线同轨道坐标系+Z轴的夹角，即高低角；αi,αj,αk为在航天器轨道坐标系下地面标志点与航天器连线的投影与轨道坐标系+X轴的夹角，即方位角；假设定位过程中，三个地面标志点pi,pj,pk的测量误差均呈零均值正态分布，则该观测组合对应的误差权系数矩阵如式（2）所示；H=(GTG)-1=h11h22h33---(2)]]>式中，H表示误差权系数矩阵；hii表示误差权系数矩阵H的对角元素，i＝1,2,3；则三个标志点pi,pj,pk组成的观测组合，其几何分布系数可按照（3）式计算：Jijk=tr(H)=h11+h22+h33---(3)]]>式中，Jijk表示几何分布系数；hii表示误差权系数矩阵H的对角元素，i＝1,2,3；对M个观测组合，分别计算其几何分布系数Jijk，取M个Jijk中的最小值min(Jijk)对应的观测组合作为本幅面阵光学遥感图像中用于航天器自主定位的最佳观测组合，对应的地面标志点为选择出来的用于定位的地面标志点。所述的步骤S2中包含：三个用于定位的地面标志点分别表示为Tp、Tm和Tn；建立相机坐标系Oc-XcYcZc、像平面坐标系O-UV、地面标志点Tp坐标系Op-XpYpZp及地心惯性坐标系OI-XIYIZI，其中，地面标志点Tp坐标系Op-XpYpZp定义为，原点在地面标志点Tp，将地面标志点Tp与星下点连线，转过β1角后定为Xp轴方向，Xp轴与地面标志点星下点连线组成平面，垂直于该平面为Zp轴，与Xp/Zp轴垂直为Yp轴；焦距为f的面阵光学载荷在成像时，地面标志点Tp在像平面上的坐标可表示为其在相机坐标系中的坐标可表示为则该标志点对应的高低角和方位角可按式（4）和（5）得出；αTp=argtg(vTp/uTp)---(4)]]>βTp=argtg(-f(uTp)2+(vTp)2)---(5)]]>记地面标志点Tp在地心惯性坐标系下的坐标为记成像时刻航天器的星下点在地心惯性坐标系下的坐标为记航天器在地心惯性坐标系下的位置为星下点与航天器间的向量在地心惯性系下记为星下点与地面标志点Tp的向量在地心惯性系下记为地面标志点Tp与航天器连线间的距离记为则根据星地几何关系，可得方程组（6）：r→=R→Tp+AI-GAG-T1(xISp,yISp,zISp)|ρTp|cosβTpcosαTpcosβTpsinαTpsinβTp---(6)]]>式（6）中包含共7个未知量，其中表示星下点坐标由标志点坐标系转到地心地固系，AI-G表示地心地固系转到地心惯性系；联立地面标志点Tm和地面标志点Tn的几何关系方程组，可得；r→=R→Tp+AI-GAG-T1(xISp,yISp,zISp)|ρTp|cosβTpcosαTpcosβTpsinαTpsinβTpr→=R→Tm+AI-GAG-T1(xISp,yISp,zISp)|ρTm|cosβTmcosαTmcosβTmsinαTmsinβTmr→=R→Tn+AI-GAG-T1(xISp,yISp,zISp)|ρTn|cosβTncosαTncosβTnsinαTnsinβTn---(7)]]>式（7）中包含9个方程，可求解其中9个未知量，该模型即为基于面阵图像的航天器定位模型。所述的步骤S3中预设处理算法为无迹卡尔曼滤波算法。所述的步骤S3中包含：采用地心惯性坐标系下的航天器定位模型建立状态方程和观测方程；利用无迹卡尔曼滤波实现状态量的实时非线性滤波估计，将地心惯性坐标系下的状态量转换为轨道坐标系下的状态量输出。本专利技术一种采用面本文档来自技高网...

【技术保护点】
一种采用面阵光学载荷的单星自主定位方法，其特征在于，包含以下步骤：S1、航天器通过面阵光学载荷对地观测后获得一幅遥感图像，经过图像预处理后，根据地面标志点信息库匹配并提取出遥感图像中包含的所有地面标志点，采用一预设提取算法，从所有地面标志点中选取三个用于定位的地面标志点；S2、建立地面标志点坐标系、相机坐标系、像平面坐标系及地心惯性坐标系，并对三个用于定位的地面标志点分别进行空间几何位置建模，结合地面标志点坐标系、相机坐标系、像平面坐标系和地心惯性坐标系之间的几何关系，推导可得基于面阵图像的航天器定位模型；S3、采用预设处理算法对航天器定位模型进行求解,得到单星自主定位输出量,以完成对单星自主定位。

【技术特征摘要】
1.一种采用面阵光学载荷的单星自主定位方法，其特征在于，包含以下步骤：S1、航天器通过面阵光学载荷对地观测后获得一幅遥感图像，经过图像预处理后，根据地面标志点信息库匹配并提取出遥感图像中包含的所有地面标志点，采用一预设提取算法，从所有地面标志点中选取三个用于定位的地面标志点；S2、建立地面标志点坐标系、相机坐标系、像平面坐标系及地心惯性坐标系，并对三个用于定位的地面标志点分别进行空间几何位置建模，结合地面标志点坐标系、相机坐标系、像平面坐标系和地心惯性坐标系之间的几何关系，推导可得基于面阵图像的航天器定位模型；S3、采用预设处理算法对航天器定位模型进行求解,得到单星自主定位输出量,以完成对单星自主定位。2.如权利要求1所述的单星自主定位方法，其特征在于，所述的步骤S1中预设提取算法包含：对所有地面标志点进行排列组合，形成若干个不同的观测组合，每一观测组合包含三个用于定位的地面标志点；对每种观测组合进行几何表述，获得几何表示矩阵；通过分析几何表述矩阵，对每种观测组合的几何分布进行分析和计算，并根据分析和计算结果选择最佳观测组合用于定位解算。3.如权利要求1或2所述的单星自主定位方法，其特征在于，所述的步骤S1中预设提取算法包含：记所有地面标志点的个数为N，则从中选取三个用于定位的地面标志点，则有种组合方式，共有M个观测组合；记任一观测组合中三个地面标志点为pi,pj,pk，则该观测组合对应的几何表述矩阵如式(1)所示；G=-cosθisinαi-cosθicosαi-sinθi-cosθjsinαj-cosθjcosαj-sinθj-cosθksinαk-cosθkcosαk-sinθk---(1)]]>式中，G表示观测组合对应的几何表述矩阵；θi,θj,θk为在航天器轨道坐标系下地面标志点与航天器的连线同轨道坐标系+Z轴的夹角，即高低角；αi,αj,αk为在航天器轨道坐标系下地面标志点与航天器连线的投影与轨道坐标系+X轴的夹角，即方位角；假设定位过程中，三个地面标志点pi,pj,pk的测量误差均呈零均值正态分布，则该观测组合对应的误差权系数矩阵如式(2)所示；H=(GTG)-1=h11h22h33---(2)]]>式中，H表示误差权系数矩阵；hii表示误差权系数矩阵H的对角元素，i＝1,2,3；则三个标志点pi,pj,pk组成的观测组合，其几何分布系数可按照(3)式计算：Jijk=tr(H)=h11+h22+h33---(3)]]>式中，Jijk表示几何分布系数；hii表示误差权系数矩阵H的对角元素，i＝1,2,3；对M个观测组合，分别计算其几何分布系数Jijk，取M个Jijk中的最小值min(Jijk)对应的观测组合作为本幅面阵光学遥感图像中用于航天器自主定位的最佳观测组合，对应的地面标志点为选择出来的用于定位的地面标志点。4.如权利要求1所述的单星自主定位方法，其特征在于，所述的步骤S...

【专利技术属性】
技术研发人员：王兆龙，孙俊，曹姝清，田少雄，
申请(专利权)人：上海航天控制技术研究所，
类型：发明
国别省市：上海;31