CN113514055B - 一种地基星敏感器大气折射与对地姿态的联合估计方法 - Google Patents

一种地基星敏感器大气折射与对地姿态的联合估计方法 Download PDF

Info

Publication number
CN113514055B
CN113514055B CN202110778625.0A CN202110778625A CN113514055B CN 113514055 B CN113514055 B CN 113514055B CN 202110778625 A CN202110778625 A CN 202110778625A CN 113514055 B CN113514055 B CN 113514055B
Authority
CN
China
Prior art keywords
star
difference
gas
observation
coordinate system
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
Application number
CN202110778625.0A
Other languages
English (en)
Other versions
CN113514055A (zh
Inventor
江洁
杨季三
田凌峰
张广军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN202110778625.0A priority Critical patent/CN113514055B/zh
Priority to PCT/CN2021/118406 priority patent/WO2023279529A1/zh
Publication of CN113514055A publication Critical patent/CN113514055A/zh
Application granted granted Critical
Publication of CN113514055B publication Critical patent/CN113514055B/zh
Priority to US18/190,164 priority patent/US20230243655A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
    • 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
    • G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
    • 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
    • G01C21/20—Instruments for performing navigational calculations

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • General Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Navigation (AREA)

Abstract

本发明提供了一种地基星敏感器大气折射与对地姿态联合估计方法。步骤包括：星敏感器拍摄星图，经过星图识别得到观测星与参考星匹配关系；根据观测时间和地点将参考星转换到地理坐标系下，得到参考星真天顶距；采用简化的大气折射模型和真天顶距以及初始蒙气差系数计算估计蒙气差，对参考星进行蒙气差补偿；计算星敏感器对地姿态，将观测星重投影到地理坐标系下，计算观测蒙气差和蒙气差误差；如果观测蒙气差与估计蒙气差不重合，则根据蒙气差误差调整蒙气差系数，重新补偿参考星并计算姿态，直到观测蒙气差与估计蒙气差重合，获得星敏感器大气折射与对地姿态联合估计结果。本发明实现了星敏感器大气折射的实时自主估计和消除，不依赖外部传感器，不受经验公式约束。

Description

一种地基星敏感器大气折射与对地姿态的联合估计方法
技术领域
本发明涉及星敏感器技术领域，具体涉及一种地基星敏感器大气折射与对地姿态的联合估计方法。
背景技术
星敏感器是一种以恒星作为参考，经星图摄影、星点质心定位、星图识别/跟踪、姿态解算等处理步骤，获得空间三轴姿态的嵌入式传感设备。作为目前姿态测量精度最高的仪器，星敏感器已在各类航行载体上得到广泛应用。对于近地面空间载体使用的星敏感器，如弹载、机载、船载等星敏感器，统称为地基星敏感器(terrestrial star tracker)。这类星敏感器为了实现高精度的三轴姿态测量，必须要克服星光大气折射(蒙气差)导致的测量误差问题。
目前常用的蒙气差校正方法主要是经验公式法，该方法将恒星天顶距(视天顶距/真天顶距)和测站大气参数(气温、气压、湿度等)代入经验蒙气差多项展开式，计算得到经验蒙气差大小，使用经验蒙气差校正恒星天顶距，可以在一定程度上减小三轴姿态测量误差。然而该方法依赖大气参数传感器，如温度、湿度、气压传感器等，并受传感器误差及经验公式误差的影响。如果使用的经验蒙气差公式是关于视天顶距的，则还需要依赖星敏感器的先验姿态，会进一步降低蒙气差估计的准确性。
此外还有一些利用星光大气折射中的不变量进行蒙气差估计的方法，如利用星光大气折射前后方位角不变的特点，可以先估计天顶方向，再估计出蒙气差；还有利用星间角距不变原则构建卡尔曼滤波方程估计蒙气差系数的方法。但是这些方法都不能同时满足实时性与最优性的要求，所以难以投入实用。
发明内容
本发明的目的在于克服以往蒙气差修正方法的不足，提出了一种同时满足实时性与最优性要求的基于大气折射系数反馈优化的地基星敏感器姿态估计方法。该方法依据星光大气折射简化模型、星图识别结果、天文坐标系转换链、QUEST姿态解算特点以及蒙气差与最优姿态的关系，构造了大气折射系数反馈机制，可以使大气折射系数反馈寻优并实时获得地基星敏感器的最优对地姿态。
本发明的技术方案如下：
本发明公开了一种地基星敏感器大气折射与对地姿态的联合估计方法，包括以下步骤：
步骤S100，星敏感器拍摄星图，经过星图识别得到观测星与参考星匹配关系；
步骤S200，根据观测时间和地点将所述参考星转换到地理坐标系下，得到所述参考星的真天顶距；
步骤S300，得到所述地理坐标系下蒙气差补偿后的参考星矢量；
步骤S400，得到星敏感器坐标系到所述地理坐标系的姿态余弦矩阵；
步骤S500，基于所述姿态余弦矩阵，得到蒙气差误差；
步骤S600，对所述蒙气差误差进行线性最小二乘拟合，得到拟合斜率和截距，若拟合斜率绝对值大于截止条件，则使大气折射系数估计值叠加反馈增益，反馈增益与拟合斜率正相关，并返回步骤S300；若拟合斜率绝对值小于截止条件，则终止迭代，得到最优大气折射系数和补偿过大气折射误差的星敏感器对地姿态。
进一步的，步骤S100中，所述星敏感器拍摄观测星图，获得成像面上星点坐标，得到星敏感器坐标系下所有观测星矢量，再进行星图识别，得到观测星与导航星表上参考星的匹配关系。
进一步的，步骤S200中，将步骤S100中得到的参考星从导航星表采用的天球坐标系转换到星敏感器观测位置的地理坐标系，得到所述地理坐标系下无蒙气差的参考星矢量，进而得到参考星的真天顶距。
进一步的，步骤S300中，采用简化的大气折射模型和真天顶距以及初始蒙气差系数计算估计蒙气差，对参考星进行蒙气差补偿，得到所述地理坐标系下蒙气差补偿后的参考星矢量。
进一步的，步骤S400中，利用QUEST算法代入星敏感器坐标系下观测星矢量和地理坐标系下蒙气差补偿后的参考星矢量，计算星敏感器坐标系到地理坐标系的姿态四元数，进而可得到星敏感器坐标系到地理坐标系的姿态余弦矩阵。
进一步的，所述对参考星进行蒙气差补偿是指将参考星天顶距减去对应的蒙气差估计值后得到新的参考星矢量。
进一步的，步骤S500中，将所述星敏感器坐标系下观测星矢量经所述姿态余弦矩阵转换到所述地理坐标系下，并和地理坐标系下参考星比较，计算观测蒙气差大小，再将所述观测蒙气差与估计蒙气差求差，得到蒙气差误差。
进一步的，步骤S600中，对蒙气差误差进行线性最小二乘拟合，得到拟合斜率和截距。
与现有技术相比，本发明有如下有益的技术效果：
(1)、本发明实现了大气折射系数的实时自主估计而不依赖其他外部传感器，如温度、气压、湿度、姿态传感器等，不受传感器及经验模型误差影响，降低了地基星敏感器消除蒙气差的难度及系统成本。
(2)、本发明实现了大气折射系数反馈寻优机制，相比以往蒙气差估计方法，能同时满足实时性与最优性要求，具有较高实用性。
(3)、本发明的方法消除了由大气折射引起的星敏感器对地姿态测量误差。该方法采用了简化的星光大气折射模型，在星图识别算法和天文坐标系转换链的基础上，根据QUEST姿态解算特点以及蒙气差与最优姿态的关系，构造了大气折射系数反馈机制。利用大气折射系数反馈机制，能够实时、最优地估计出大气折射系数并同时得到消除大气折射影响的地基星敏感器对地姿态。
附图说明
图1为本发明的一种地基星敏感器大气折射与对地姿态的联合估计方法流程图；
图2为本发明的星光大气折射简化模型示意图；
图3为本发明的大气折射系数反馈机制示意图；
图4为本发明的转台倾角20°时星敏感器连续1000帧观星测试数据图；
图5为本发明的转台倾角20°时星敏感器单帧数据蒙气差估计效果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式并参照附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
本发明提供了一种地基星敏感器大气折射与对地姿态的联合估计方法，具体的，包括以下步骤：
步骤S100，星敏感器拍摄星图，经过星图识别得到观测星与参考星匹配关系。
具体的，星敏感器拍摄观测星图，获得成像面上星点坐标，得到星敏感器坐标系下所有观测星矢量，再进行星图识别，得到观测星与导航星表上参考星的匹配关系。
具体的，根据所述星敏感器内参、畸变模型及小孔成像模型，得到星敏感器坐标系下所有观测星矢量。
具体的，所述星敏感器拍摄观测星图，得到一系列星敏感器坐标系下观测矢量Ws，Ws＝[ws,1,ws,2,ws,3,…,ws,n]。
所述星图识别方法可采用子图同构类算法(如三角形算法、金字塔算法等)，也可采用模式识别类算法(如栅格算法、采用径向和环向特征的识别算法等)，还可以采用其他类型的识别算法(如基于神经网络的星图识别算法等)。所述星图识别算法不唯一，需结合实际条件选择合适的算法，但需保证识别结果正确。所述参考星，其信息包含赤经、赤纬、自行、视差、径向速度。
具体的，星敏感器拍摄观测星图，获得成像面上星点坐标[X,Y]，X＝[X1,X2,X3,…,Xn]T，Y＝[Y1,Y2,Y3,…,Yn]T。根据镜头畸变模型、小孔成像模型以及星敏感器内参(主点X0,Y0、焦距f、畸变参数p1,p2,q1,q2、像元尺寸dx,dy)，获得星敏感器坐标系下的观测星矢量Ws＝[ws,1,ws,2,ws,3,…,ws,n]，
其中：
Figure BDA0003155299980000051
其中：
Figure BDA0003155299980000052
再使用改进三角形算法进行星图识别得到导航星表上与观测星匹配的参考星Vr＝[vr,1,vr,2,vr,3,…,vr,n]，其中vr,i＝[RAi,DEi,pmRAi,pmDEi,Plxi,rvi]T，其中，RAi表示赤经，DEi表示赤纬，pmRAi是自行的经度分量,pmDEi是自行的纬度分量，Plxi表示视差，rvi表示径向速度。
步骤S200，根据观测时间和地点将所述参考星转换到地理坐标系下，得到所述参考星的真天顶距。
具体的，将步骤S100中得到的参考星从导航星表采用的天球坐标系转换到星敏感器观测位置的地理坐标系，得到所述地理坐标系下无蒙气差的参考星矢量，进而得到参考星的真天顶距。
具体的，根据星敏感器拍摄时间和地点以及地球自转信息，基于高精度的天文坐标系转换链，将所述步骤S100中得到的参考星从导航星表采用的天球坐标系转换到星敏感器观测位置的地理坐标系。
具体的，参考星坐标系转换时，星敏感器拍摄时间为曝光中心时刻，表示方法可选星敏感器曝光中心时刻UTC时间，精确到毫秒；地点为曝光中心时刻星敏感器的经度、纬度、高度，表示方法可选WGS84椭球参考系，精确到秒和米。
具体的，地球自转信息包括地球当日极移和DUT1时间，可由国际地球自转服务(International Earth Rotation Service,IERS)提供。在国际天文联合会(International Astronomical Union,IAU)基础天文学标准(Standards Of FundamentalAstronomy,SOFA)提供的高精度坐标系转换链基础上，进行参考星坐标系转换，还需获得导航星表上参考星的赤经、赤纬、自行、视差、径向速度等参数。
具体的，观测点地理坐标系即为东北天坐标系，得到地理坐标系下参考矢量Vg＝[vg,1,vg,2,vg,3,…,vg,n]，天顶矢量U＝[0,0,1]T，真天顶距为Z＝[Z1,Z2,Z3,…,Zn]T，其中：
Zi＝arccos(U·vg,i) (3)
具体的，获得星敏感器曝光中心时刻UTC时间T、星敏感器所在经纬高坐标P(WGS84椭球参考系)、地球当日极移Pm＝[px,py]以及dut1时间，其中Pm和dut1可由国际地球自转服务(IERS,https://www.iers.org/)提供。再使用国际天文联合会(IAU)基础天文学标准(SOFA,http://www.iausofa.org/)提供的高精度坐标系转换链，将参考星从导航星表采用的天球坐标系转换到星敏感器观测位置的地理坐标系，得到地理坐标系下参考矢量：
Vg＝[vg,1,vg,2,vg,3,…,vg,n]，
Vg＝iauCST(Vr,T,P,Pm,dut1) (4)
其中，iau(*)表示采用iauSOFA函数完成参考星坐标从天球系转换到地理坐标系，Vr是参考星天球系下坐标，T是观测时刻，P是观测位置，Pm是地球当日极移，dut1是UT1时间与UTC时间差值常数。
已知天顶矢量U＝[0,0,1]T，计算所有参考星真天顶距Z＝[Z1,Z2,Z3,…,Zn]T：
Zi＝arccos(U·vg,i) (5)
步骤S300，得到所述地理坐标系下蒙气差补偿后的参考星矢量。
具体的，根据星光大气折射简化模型，如图2所示，得到蒙气差计算公式，代入大气折射系数估计值和所有参考星的所述真天顶距，计算所有参考星的估计蒙气差大小并对参考星矢量沿仰角增加的方向做蒙气差修正，得到地理坐标系下蒙气差补偿后的参考星矢量。
具体的，星光大气折射简化模型是指，将地球大气层看作由无数同心球层组成，大气密度和折射率随高度下降而增加，并且对于星敏感器视场内的大气，可忽略大气球层曲率而认为从观测点到星际空间是1～n的平面层，大气外星光的入射角为真天顶距Z，星光传播路径在层与层之间发生折射现象，光线是一条凸向观测点铅垂线的曲线。设最外层大气折射率为nn，根据折射定律有：
sin(Z)＝nnsin(Zn) (6)
Zn为星光矢量在最外层大气的真天顶距，对于相邻的下一层有：
nn sin(Zn)＝nn-1sin(Zn-1) (7)
直到最下一层有：
n2 sin(Z2)＝n1sin(Z1) (8)
综合上式可得：
sin(Z)＝n1sin(Z1) (9)
其中Z1即为星光矢量传播到最底层时的真天顶距，对于观测点来说即为视天顶距，用Z′表示，Z′<Z。蒙气差ρ定义为真天顶距Z与视天顶距Z′之差：
ρ＝Z-Z′ (10)
所以式(9)可写为：
sin(Z)＝n1 sin(Z-ρ)＝n1 sin(Z)cos(ρ)-n1cos(Z)sin(ρ) (11)
由于ρ极小，取cos(ρ)＝1，sin(ρ)＝ρ，可得蒙气差ρ为：
Figure BDA0003155299980000071
其中，N为大气折射系数，以角秒为单位。当n1取标准大气折射率1.000292时，N0＝60.21"。大气折射系数
Figure BDA0003155299980000072
初值可设为标准大气折射系数
Figure BDA0003155299980000073
根据蒙气差计算公式(12)，可代入大气折射系数估计值
Figure BDA0003155299980000074
和所有参考星的真天顶距Z，计算所有参考星矢量的估计蒙气差大小
Figure BDA0003155299980000075
其中：
Figure BDA0003155299980000076
再对参考星矢量沿仰角增加的方向做蒙气差修正，得到地理坐标系下蒙气差补偿后的参考星矢量Vgc＝[vgc,1,vgc,2,vgc,3,…,vgc,n]，其中：
Figure BDA0003155299980000077
其中，I表示单位矩阵，ni表示参考星矢量vg,i增加仰角
Figure BDA0003155299980000078
时的转轴矢量：
ni＝vg,i×U (15)
ni ∧表示由矢量元素构成的反对称矩阵，ni(1)、ni(2)、ni(3)分别表示矢量ni的第1、2、3个元素：
Figure BDA0003155299980000081
步骤S400，得到星敏感器坐标系到地理坐标系的姿态余弦矩阵。
具体的，利用QUEST算法代入星敏感器坐标系下观测星矢量和地理坐标系下蒙气差补偿后的参考星矢量，计算星敏感器坐标系到地理坐标系的姿态四元数，进而可得到星敏感器坐标系到地理坐标系的姿态余弦矩阵。
具体的，利用QUEST算法代入星敏感器坐标系下观测星矢量和地理坐标系下蒙气差补偿后的参考星矢量，计算得到星敏感器坐标系到地理坐标系的姿态四元数：
Qout＝QUEST(Ws,Vgc) (17)
其中，QUEST(*)指QUEST算法，Ws是星敏系下观测星矢量，Vgc是地理系下蒙气差补偿后的参考星矢量；
式中，Qout＝[v,q4]，其中v为四元数虚部v＝[q1,q2,q3]，q1、q2、q3为四元数虚部各元素，q4为四元数实部，进而得到星敏感器对地姿态余弦矩阵为：
Figure BDA0003155299980000082
步骤S500，基于所述姿态余弦矩阵，得到蒙气差误差；
具体的，将星敏感器坐标系下观测星矢量经姿态余弦矩阵转换到地理坐标系下，并和地理坐标系下参考星比较，计算观测蒙气差大小，再将观测蒙气差与估计蒙气差求差，得到蒙气差误差。
其中，将星敏感器坐标系下观测星矢量经姿态余弦矩阵转换到地理坐标系下得到Wg＝[wg,1,wg,2,wg,3,…,wg,n]：
Figure BDA0003155299980000083
再和地理坐标系下参考星比较，计算观测蒙气差大小
Figure BDA0003155299980000084
其中：
Figure BDA0003155299980000085
再将观测蒙气差与估计蒙气差求差，得到蒙气差误差：
Figure BDA0003155299980000091
步骤S600，对蒙气差误差进行线性最小二乘拟合，得到拟合斜率和截距，若拟合斜率绝对值大于截止条件，则使大气折射系数估计值叠加反馈增益，反馈增益与拟合斜率正相关，并返回步骤S300；若拟合斜率绝对值小于截止条件，则终止迭代，得到最优大气折射系数和补偿过大气折射误差的星敏感器对地姿态。
具体的，根据大气折射系数反馈机制，对蒙气差误差进行线性最小二乘拟合，得到拟合斜率和截距。
其中，大气折射系数反馈机制是指，根据星光大气折射简化模型和QUEST算法最小二乘的特点，当
Figure BDA0003155299980000092
(N*为真值)时，蒙气差补偿后的参考星仰角小于真实位置，仰角分布范围大于观测星，不能与观测星重合，经过QUEST姿态解算后，视轴指向低于真实指向，姿态估计有误差，但是蒙气差误差与天顶距成正比，斜率大于零,且当
Figure BDA0003155299980000093
越远离N*则蒙气差误差斜率越大；当
Figure BDA0003155299980000094
时，蒙气差补偿后的参考星仰角等于真实位置，仰角分布范围与观测星重合，经过QUEST姿态解算后，视轴指向与真实指向重合，姿态估计无误差，蒙气差误差等于零，斜率等于零；当
Figure BDA0003155299980000095
时，蒙气差补偿后的参考星仰角大于真实位置，仰角分布范围小于观测星，不能与观测星重合，经过QUEST姿态解算后，视轴指向高于真实指向，姿态估计有误差，但是蒙气差误差与天顶距成反比，斜率小于零，且当
Figure BDA0003155299980000096
越远离N*则蒙气差误差斜率越小。据此构造大气折射系数反馈机制，用最小二乘法拟合蒙气差误差直线参数：
[k,b]T＝(MTM)-1MTΔP (22)
其中M＝[Z,L]，L＝[1,1,1,…,1]1×n T。如果|k|>ε，(ε为终止条件，接近0)，则说明
Figure BDA0003155299980000097
没有达到最优，需要更新
Figure BDA0003155299980000098
给
Figure BDA0003155299980000099
叠加增益G(k)，并返回步骤S300。其中，G(k)与k正相关，为使
Figure BDA00031552999800000910
快速收敛，可采用PID增益方式：
Figure BDA00031552999800000911
其中，gP、gI、gD分别为比例、积分、微分系数，ki是当前蒙气差误差斜率，ki-1是前一帧蒙气差误差斜率。如果|k|<ε，则可以认为
Figure BDA00031552999800000912
达到最优，此时得到最优大气折射系数，同时得到补偿过大气折射误差的星敏感器对地姿态。
具体实施效果
为验证本发明提出的星敏感器姿态估计方法的有效性，进行了外场观星实验。实验前将星敏感器安装在转台上，转台初始倾角为0°，安装面法线大约朝向天顶方向，转动转台6次，每次转动大约5°，当转台稳定后开始采集数据，共计7组星图数据各1000帧，采用本发明提出的方法进行姿态估计，再采用地理坐标系下观测矢量投影与参考矢量的平均角距误差(重投影误差)作为评价指标：
Figure BDA0003155299980000101
Figure BDA0003155299980000102
b表示一帧星图中所有观测星个数；B表示所有星图的个数。
得到校正前后姿态测量对比结果如表1所示，表1为转台倾角0°,5°,10°,……,30°时各组姿态估计效果统计，其中以转台倾角20°时为例，1000帧测试对比结果如图4所示，取其中一帧数据观察蒙气差估计情况如图5所示，表2为转台倾角20°时星敏感器单帧数据大气折射系数迭代过程。
表1
Figure BDA0003155299980000103
Figure BDA0003155299980000111
表2
Figure BDA0003155299980000112
综上所述，本发明提供了一种地基星敏感器大气折射与对地姿态联合估计方法，包括：星敏感器拍摄星图，经过星图识别得到观测星与参考星匹配关系；根据观测时间和地点将参考星转换到地理坐标系下，得到参考星真天顶距；采用简化的大气折射模型和真天顶距以及初始蒙气差系数计算估计蒙气差，对参考星进行蒙气差补偿；计算星敏感器对地姿态，将观测星重投影到地理坐标系下，计算观测蒙气差和蒙气差误差；如果观测蒙气差与估计蒙气差不重合，则根据蒙气差误差调整蒙气差系数，重新补偿参考星并计算姿态，直到观测蒙气差与估计蒙气差重合，获得星敏感器大气折射与对地姿态联合估计结果。本发明实现了星敏感器大气折射的实时自主估计和消除，不受经验公式约束，也不依赖外部传感器，如温度、气压、湿度、姿态传感器等，降低了地基星敏感器消除蒙气差的难度及系统成本；本方法采用了简化的星光大气折射模型，在星图识别算法和天文坐标系转换链的基础上，根据QUEST姿态解算特点以及蒙气差与最优姿态的关系，构造了大气折射系数反馈机制。利用大气折射系数反馈机制，能够实时、最优地估计出大气折射系数并同时得到消除大气折射影响的地基星敏感器对地姿态。
应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (2)

1.一种地基星敏感器大气折射与对地姿态的联合估计方法，其特征在于，所述方法包括以下步骤：
步骤S100，星敏感器拍摄星图，获得成像面上星点坐标，得到星敏感器坐标系下所有观测星矢量，再进行星图识别，得到观测星与参考星匹配关系；
步骤S200，根据观测时间和地点将步骤S100中得到的参考星从导航星表采用的天球坐标系转换到星敏感器观测位置的地理坐标系，得到所述地理坐标系下无蒙气差的参考星矢量，进而得到参考星的真天顶距；
步骤S300，采用简化的大气折射模型和真天顶距以及初始蒙气差系数计算估计蒙气差，对参考星进行蒙气差补偿，得到所述地理坐标系下蒙气差补偿后的参考星矢量；对参考星进行蒙气差补偿是将参考星天顶距减去对应的蒙气差估计值后得到新的参考星矢量：简化的大气折射模型是指，将地球大气层看作由无数同心球层组成，大气密度和折射率随高度下降而增加，并且对于星敏感器视场内的大气，忽略大气球层曲率而认为从观测点到星际空间是
Figure 846949DEST_PATH_IMAGE001
的平面层，大气外星光的入射角为真天顶距
Figure 608975DEST_PATH_IMAGE002
，星光传播路径在层与层之间发生折射现象，光线是一条凸向观测点铅垂线的曲线，设最外层大气折射率为
Figure 644933DEST_PATH_IMAGE003
，根据折射定律有：
Figure 520748DEST_PATH_IMAGE004
（6）
Figure 709284DEST_PATH_IMAGE005
为星光矢量在最外层大气的真天顶距，对于相邻的下一层有：
Figure 346938DEST_PATH_IMAGE006
（7）
直到最下一层有：
Figure 368728DEST_PATH_IMAGE007
（8）
综合上式可得：
Figure 625396DEST_PATH_IMAGE008
（9）
其中
Figure 58652DEST_PATH_IMAGE009
即为星光矢量传播到最底层时的真天顶距，对于观测点来说即为视天顶距，用
Figure DEST_PATH_IMAGE010
表示，
Figure 165410DEST_PATH_IMAGE011
，蒙气差
Figure DEST_PATH_IMAGE012
定义为真天顶距
Figure 582485DEST_PATH_IMAGE013
与视天顶距
Figure 125069DEST_PATH_IMAGE010
之差：
Figure DEST_PATH_IMAGE014
（10）
所以式（9）可写为：
Figure 740727DEST_PATH_IMAGE015
（11）
由于
Figure 329971DEST_PATH_IMAGE012
极小，取
Figure 63703DEST_PATH_IMAGE016
，
Figure 786808DEST_PATH_IMAGE017
，可得蒙气差
Figure 538864DEST_PATH_IMAGE012
为：
Figure 324110DEST_PATH_IMAGE018
（12）
其中,
Figure 263247DEST_PATH_IMAGE019
为大气折射系数，以角秒为单位，当
Figure 383519DEST_PATH_IMAGE020
取标准大气折射率1.000292时，
Figure 990081DEST_PATH_IMAGE021
，大气折射系数估计值
Figure 937439DEST_PATH_IMAGE022
初值设为标准大气折射系数
Figure 629452DEST_PATH_IMAGE023
，根据蒙气差计算公式（12），代入大气折射系数估计值
Figure 818993DEST_PATH_IMAGE022
和所有参考星的真天顶距
Figure 14482DEST_PATH_IMAGE013
，计算所有参考星矢量的估计蒙气差大小
Figure DEST_PATH_IMAGE024
，其中：
Figure 801916DEST_PATH_IMAGE025
（13）
再对参考星矢量沿仰角增加的方向做蒙气差修正，得到地理坐标系下蒙气差补偿后的参考星矢量
Figure DEST_PATH_IMAGE026
，其中：
Figure 433755DEST_PATH_IMAGE027
（14）
其中，
Figure DEST_PATH_IMAGE028
表示单位矩阵，
Figure 131715DEST_PATH_IMAGE029
表示参考星矢量
Figure 696557DEST_PATH_IMAGE030
增加仰角
Figure 123734DEST_PATH_IMAGE031
时的转轴矢量：
Figure DEST_PATH_IMAGE032
（15）
Figure 649393DEST_PATH_IMAGE033
表示由矢量元素构成的反对称矩阵，
Figure 682203DEST_PATH_IMAGE034
、
Figure 976918DEST_PATH_IMAGE035
、
Figure 420668DEST_PATH_IMAGE036
分别表示矢量
Figure 823837DEST_PATH_IMAGE037
的第1、2、3个元素：
Figure 644025DEST_PATH_IMAGE038
（16）
步骤S400，利用QUEST算法代入星敏感器坐标系下观测星矢量和地理坐标系下蒙气差补偿后的参考星矢量，计算星敏感器坐标系到地理坐标系的姿态四元数，进而可得到星敏感器坐标系到地理坐标系的姿态余弦矩阵：利用QUEST算法代入星敏感器坐标系下观测星矢量和地理坐标系下蒙气差补偿后的参考星矢量，计算得到星敏感器坐标系到地理坐标系的姿态四元数：
Figure 822941DEST_PATH_IMAGE039
其中，QUEST(*)指QUEST算法，Ws是星敏系下观测星矢量；Vgc是地理系下蒙气差补偿后的参考星矢量；
式中，
Figure 686860DEST_PATH_IMAGE040
，其中
Figure 62478DEST_PATH_IMAGE041
为四元数虚部
Figure 702669DEST_PATH_IMAGE042
，q1、 q2 、q3为四元数虚部各元素，
Figure 847342DEST_PATH_IMAGE043
为四元数实部，进而得到星敏感器对地姿态余弦矩阵为：
Figure 147743DEST_PATH_IMAGE044
步骤S500，将所述星敏感器坐标系下观测星矢量经所述姿态余弦矩阵转换到所述地理坐标系下，并和地理坐标系下参考星比较，计算观测蒙气差大小，再将所述观测蒙气差与估计蒙气差求差，得到蒙气差误差：其中，将星敏感器坐标系下观测星矢量经姿态余弦矩阵转换到地理坐标系下得到
Figure 10656DEST_PATH_IMAGE045
：
Figure 197748DEST_PATH_IMAGE046
再和地理坐标系下参考星比较，计算观测蒙气差大小
Figure 196928DEST_PATH_IMAGE047
，其中：
Figure 933809DEST_PATH_IMAGE048
再将观测蒙气差与估计蒙气差求差，得到蒙气差误差：
Figure 752860DEST_PATH_IMAGE049
步骤S600，对所述蒙气差误差进行线性最小二乘拟合，得到拟合斜率和截距，若拟合斜率绝对值大于截止条件，则使大气折射系数估计值叠加反馈增益，反馈增益与拟合斜率正相关，并返回步骤S300；若拟合斜率绝对值小于截止条件，则终止迭代，得到最优大气折射系数和补偿过大气折射误差的星敏感器对地姿态。
2.根据权利要求1所述的一种地基星敏感器大气折射与对地姿态的联合估计方法，其特征在于，步骤S600中，根据大气折射系数反馈机制，对蒙气差误差进行线性最小二乘拟合，得到拟合斜率和截距；
其中，大气折射系数反馈机制是指，根据星光大气折射简化模型和QUEST算法最小二乘的特点，当
Figure DEST_PATH_IMAGE050
时，
Figure 406958DEST_PATH_IMAGE051
为真值，蒙气差补偿后的参考星仰角小于真实位置，仰角分布范围大于观测星，不能与观测星重合，经过QUEST姿态解算后，视轴指向低于真实指向，姿态估计有误差，但是蒙气差误差与天顶距成正比，斜率大于零,且当
Figure DEST_PATH_IMAGE052
越远离
Figure 713174DEST_PATH_IMAGE051
则蒙气差误差斜率越大；当
Figure 853912DEST_PATH_IMAGE053
时，蒙气差补偿后的参考星仰角等于真实位置，仰角分布范围与观测星重合，经过QUEST姿态解算后，视轴指向与真实指向重合，姿态估计无误差，蒙气差误差等于零，斜率等于零；当
Figure DEST_PATH_IMAGE054
时，蒙气差补偿后的参考星仰角大于真实位置，仰角分布范围小于观测星，不能与观测星重合，经过QUEST姿态解算后，视轴指向高于真实指向，姿态估计有误差，但是蒙气差误差与天顶距成反比，斜率小于零，且当
Figure 409527DEST_PATH_IMAGE052
越远离
Figure 70578DEST_PATH_IMAGE051
则蒙气差误差斜率越小，据此构造大气折射系数反馈机制，用最小二乘法拟合蒙气差误差直线参数：
Figure 44350DEST_PATH_IMAGE055
其中
Figure DEST_PATH_IMAGE056
，
Figure 60717DEST_PATH_IMAGE057
，如果
Figure DEST_PATH_IMAGE058
，
Figure 992376DEST_PATH_IMAGE059
为终止条件，则说明
Figure 299861DEST_PATH_IMAGE052
没有达到最优，需要更新
Figure 144451DEST_PATH_IMAGE052
，给
Figure 879189DEST_PATH_IMAGE052
叠加增益
Figure DEST_PATH_IMAGE060
，并返回步骤S300，其中，
Figure 347079DEST_PATH_IMAGE060
与
Figure 612582DEST_PATH_IMAGE061
正相关，为使
Figure 75793DEST_PATH_IMAGE052
快速收敛，可采用PID增益方式：
Figure 981433DEST_PATH_IMAGE062
其中，
Figure 765980DEST_PATH_IMAGE063
、
Figure 415267DEST_PATH_IMAGE064
、
Figure 998564DEST_PATH_IMAGE065
分别为比例、积分、微分系数，
Figure 75105DEST_PATH_IMAGE066
是当前蒙气差误差斜率，
Figure 90158DEST_PATH_IMAGE067
是前一帧蒙气差误差斜率，如果
Figure 543136DEST_PATH_IMAGE068
，则认为
Figure 121885DEST_PATH_IMAGE052
达到最优，此时得到最优大气折射系数，同时得到补偿过大气折射误差的星敏感器对地姿态。
CN202110778625.0A 2021-07-09 2021-07-09 一种地基星敏感器大气折射与对地姿态的联合估计方法 Active CN113514055B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202110778625.0A CN113514055B (zh) 2021-07-09 2021-07-09 一种地基星敏感器大气折射与对地姿态的联合估计方法
PCT/CN2021/118406 WO2023279529A1 (zh) 2021-07-09 2021-09-15 一种地基星敏感器大气折射与对地姿态的联合估计方法
US18/190,164 US20230243655A1 (en) 2021-07-09 2023-03-27 Method for joint estimation of stellar atmospheric refraction and star tracker atitude

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110778625.0A CN113514055B (zh) 2021-07-09 2021-07-09 一种地基星敏感器大气折射与对地姿态的联合估计方法

Publications (2)

Publication Number Publication Date
CN113514055A CN113514055A (zh) 2021-10-19
CN113514055B true CN113514055B (zh) 2022-10-28

Family

ID=78066583

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110778625.0A Active CN113514055B (zh) 2021-07-09 2021-07-09 一种地基星敏感器大气折射与对地姿态的联合估计方法

Country Status (3)

Country Link
US (1) US20230243655A1 (zh)
CN (1) CN113514055B (zh)
WO (1) WO2023279529A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116952206B (zh) * 2023-09-19 2024-02-02 中国科学院光电技术研究所 一种基于星光测姿的空间目标检测方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107449444A (zh) * 2017-07-17 2017-12-08 中国人民解放军国防科学技术大学 一种多星图姿态关联的星敏感器内参数标定方法及其装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8355868B2 (en) * 2009-09-30 2013-01-15 Ge Aviation Systems Llc Method and system for spectral image celestial navigation
CN102261921B (zh) * 2011-06-14 2013-01-30 哈尔滨工业大学 一种修正大气折射对星敏感器精度影响的方法
CN104977021B (zh) * 2014-04-02 2018-06-12 北京自动化控制设备研究所 一种星光导航环境误差补偿方法
RU2566379C1 (ru) * 2014-04-14 2015-10-27 Открытое акционерное общество "Ракетно-космическая корпорация "Энергия" имени С.П. Королева" Способ определения величины атмосферной рефракции в условиях космического полета
CN104236546B (zh) * 2014-09-10 2017-01-11 中国空间技术研究院 一种卫星星光折射导航误差确定与补偿方法
CN104280049B (zh) * 2014-10-20 2017-04-05 北京控制工程研究所 一种高精度星敏感器外场精度测试方法
GB2575971A (en) * 2018-07-23 2020-02-05 Atlantic Inertial Systems Ltd A navigation system
CN111537003B (zh) * 2020-06-19 2021-09-07 北京航空航天大学 一种基于折射面共线的星光大气折射测量校正方法
CN112050832B (zh) * 2020-08-21 2023-04-18 中国科学院国家空间科学中心 一种航空星敏感器大气折射误差修正方法
CN113049006B (zh) * 2021-03-25 2024-01-16 北京航空航天大学 一种基于模型参数估计的星光矢量蒙气差修正方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107449444A (zh) * 2017-07-17 2017-12-08 中国人民解放军国防科学技术大学 一种多星图姿态关联的星敏感器内参数标定方法及其装置

Also Published As

Publication number Publication date
CN113514055A (zh) 2021-10-19
US20230243655A1 (en) 2023-08-03
WO2023279529A1 (zh) 2023-01-12

Similar Documents

Publication Publication Date Title
CN111537003B (zh) 一种基于折射面共线的星光大气折射测量校正方法
CN106871927B (zh) 一种无人机光电吊舱安装误差标校方法
US7751645B2 (en) Precision optical systems with performance characterization and uses thereof
CN104913780B (zh) 集成gnss和ccd天顶筒的高精度垂线偏差快速测量方法
US20120257792A1 (en) Method for Geo-Referencing An Imaged Area
CN102261921B (zh) 一种修正大气折射对星敏感器精度影响的方法
JP2008506167A (ja) 画像に関連するロケーションを確定する方法および装置
CN110631567B (zh) 一种差分天空偏振罗盘大气折射误差的反演及修正方法
CN107490364A (zh) 一种大角度倾斜成像航空相机对地目标定位方法
CN106873004B (zh) 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法
CN110030978A (zh) 一种全链路光学卫星几何成像模型构建方法及系统
CN113514055B (zh) 一种地基星敏感器大气折射与对地姿态的联合估计方法
CN108362268A (zh) 一种基于视频测量的自动天文测量方法和测量系统
Yan et al. Correction of atmospheric refraction geolocation error for high resolution optical satellite pushbroom images
KR20230040921A (ko) 티베트고원(靑藏高原) 지역에서 광학적 원격탐지 위성영상의 대기굴절 위치고정 오차 수정방법
CN114838740A (zh) 一种考虑不同经纬度区域的卫星图像几何定标方法
CN115597586A (zh) 一种基于大气偏振模式对称性的航向角提取方法
Speyerer et al. In-flight geometric calibration of the lunar reconnaissance orbiter camera
CN106996779A (zh) 基于gnss的紫外敏感器系统误差在轨标定方法
JP2009509125A (ja) 画像に関連する位置を決定するための方法および装置
CN110046430A (zh) 基于椭球分层大气折射的光学卫星影像精确对地定位方法
CN113049006B (zh) 一种基于模型参数估计的星光矢量蒙气差修正方法
CN113720352B (zh) 一种具有蒙气差效应的星图模拟方法
Guienko Geometric accuracy of Ikonos: zoom in
Vasilyuk et al. Calibration of the intrinsic parameters of the digital camera of a star tracker based on ground-based observations of stars, taking atmospheric refraction and aberration of light into account

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