CN114234962A - 多星敏感器在轨热变形修正方法以及存储介质和电子设备 - Google Patents

多星敏感器在轨热变形修正方法以及存储介质和电子设备 Download PDF

Info

Publication number
CN114234962A
CN114234962A CN202111326115.6A CN202111326115A CN114234962A CN 114234962 A CN114234962 A CN 114234962A CN 202111326115 A CN202111326115 A CN 202111326115A CN 114234962 A CN114234962 A CN 114234962A
Authority
CN
China
Prior art keywords
star sensor
star
optical axis
main optical
thermal deformation
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.)
Granted
Application number
CN202111326115.6A
Other languages
English (en)
Other versions
CN114234962B (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.)
Shanghai Aerospace Control Technology Institute
Original Assignee
Shanghai Aerospace Control Technology Institute
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 Shanghai Aerospace Control Technology Institute filed Critical Shanghai Aerospace Control Technology Institute
Priority to CN202111326115.6A priority Critical patent/CN114234962B/zh
Publication of CN114234962A publication Critical patent/CN114234962A/zh
Application granted granted Critical
Publication of CN114234962B publication Critical patent/CN114234962B/zh
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/20—Instruments for performing navigational calculations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Manufacturing & Machinery (AREA)
  • Astronomy & Astrophysics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Navigation (AREA)

Abstract

本发明公开了一种多星敏感器在轨热变形修正方法以及存储介质和电子设备，该方法包括：确定其余星敏感器相对于第一星敏感器的姿态关系矩阵和地面主光轴夹角；提取星敏感器在各自测量坐标系中的主光轴矢量；将其余星敏感器各自测量坐标系中的主光轴矢量转换至第一星敏感器测量坐标系，并确定主光轴矢量夹角；根据主光轴矢量夹角、地面主光轴夹角、主光轴矢量，确定热变形修正实时矩阵；其余星敏感器通过相应的姿态关系矩阵和热变形修正实时矩阵将各自测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，进行三轴姿态解算。本发明可实现共基准面安装平台的热变形的在轨测量，并具有较高的普适性。

Description

多星敏感器在轨热变形修正方法以及存储介质和电子设备
技术领域
本发明涉及航天器姿态测量技术领域，具体涉及一种多星敏感器在轨热变形修正方法以及存储介质和电子设备。
背景技术
随着遥感卫星的技术发展和遥感图像分辨率的提高，遥感图像数据的应用领域正在扩展，各领域对遥感数据的要求也在不断升高。对于星敏感器与遥感载荷共基准面安装的微纳卫星而言，由于受制于体积重量的约束，热控手段、设备、效果均无法与大型卫星平台相比，因此安装平台的热变形与大卫星相较程度更大，对成像载荷的影响也更严重。为保证遥感卫星成像链路中成像模型精度提高，需要微纳卫星平台能够脱离载荷测量值，独立获取安装结构热变形的参数。现有的卫星在轨标定多是关注星敏感器与载荷相机之间的关联角变化，无法直接反映安装平台的热变形程度；或是需要惯性敏感器配合星敏感器进行系统常值误差校正，但无法将热变形的影响从整体误差中剥离出。
发明内容
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。为此，本发明的第一个目的在于提供一种多星敏感器在轨热变形修正方法，可实现共基准面安装平台的热变形的在轨测量，并可辅助载荷成像模型的精确建立，进而提高遥感图像品质，提高定位精度，并在应用上具有较高的普适性。
本发明的第二个目的在于提供一种计算机可读存储介质。
本发明的第三个目的在于提供一种电子设备。
为达到上述目的，本发明通过以下技术方案实现：
一种多星敏感器在轨热变形修正方法，包括：步骤S1：建立三台星敏感器的各自测量坐标系，基于所述测量坐标系，以第一星敏感器为基准，确定第二星敏感器和第三星敏感器分别相对于所述第一星敏感器的姿态关系矩阵和地面主光轴夹角；步骤S2：三台星敏感器在同一时刻对各自视场内的恒星成像，并提取当前时刻下三台星敏感器在各自测量坐标系中的主光轴矢量；步骤S3：将所述第二星敏感器和所述第三星敏感器各自测量坐标系中的主光轴矢量分别转换至所述第一星敏感器测量坐标系；步骤S4：确定当前时刻转换至所述第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量与第一星敏感器测量坐标系下的主光轴矢量之间的主光轴矢量夹角；步骤S5：根据所述主光轴矢量夹角、所述地面主光轴夹角、第一星敏感器测量坐标系下的所述主光轴矢量和转换至所述第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量，确定当前时刻其余星敏感器相对于所述第一星敏感器的热变形修正实时矩阵；步骤S6：获取各个星敏感器各自测量坐标系下的恒星星光矢量，其余星敏感器通过相对于所述第一星敏感器的姿态关系矩阵和当前时刻下的热变形修正实时矩阵将各自的恒星星光矢量投影至所述第一星敏感器测量坐标系，以进行三轴姿态解算。
可选的，所述步骤S3中，所述第二星敏感器和所述第三星敏感器各自测量坐标系中的主光轴矢量分别通过相应的所述姿态关系矩阵转换至所述第一星敏感器测量坐标系。
可选的，所述步骤S4中，主光轴矢量夹角采用如下公式表示：
Figure BDA0003347200180000021
Figure BDA0003347200180000022
其中，α为第二星敏感器转换后的主光轴矢量夹角，β为第三星敏感器转换后的主光轴矢量夹角，
Figure BDA0003347200180000023
为第一星敏感器测量坐标系下的主光轴矢量，
Figure BDA0003347200180000024
为第二星敏感器测量坐标系下的主光轴矢量，
Figure BDA0003347200180000025
为第三星敏感器测量坐标系下的主光轴矢量，
Figure BDA0003347200180000026
和
Figure BDA0003347200180000027
分别为第二星敏感器和第三星敏感器相对于第一星敏感器的姿态关系矩阵。
可选的，所述步骤S5中，当前时刻其余星敏感器相对于所述第一星敏感器的热变形修正实时矩阵采用如下公式表示：
Figure BDA0003347200180000028
Figure BDA0003347200180000029
其中，αi和βi分别为第i时刻第二星敏感器和第三星敏感器的主光轴矢量夹角，α标和β标分别为第二星敏感器和第三星敏感器相对于第一星敏感器的地面主光轴夹角，
Figure BDA0003347200180000031
和
Figure BDA0003347200180000032
分别为第i时刻第二星敏感器和第三星敏感器相对于第一星敏感器的热变形修正实时矩阵，
Figure BDA0003347200180000033
为转换至第一星敏感器测量坐标系下的第二星敏感器的主光轴矢量，
Figure BDA0003347200180000034
为转换至第一星敏感器测量坐标系下的第三星敏感器的主光轴矢量。
可选的，所述方法还包括：重复步骤S2～S5获取若干个时刻第二星敏感器相对于第一星敏感器的若干个热变形修正实时矩阵，以及第三星敏感器相对于第一星敏感器的若干个热变形修正实时矩阵。
可选的，所述方法还包括：获取各个星敏感器各自测量坐标系下的恒星星光矢量，其余星敏感器通过相对于所述第一星敏感器的姿态关系矩阵和若干个时刻下的若干个热变形修正实时矩阵将各自的恒星星光矢量投影至所述第一星敏感器测量坐标系，并从投影后的若干个时刻下的恒星星光矢量中选取空间优越的恒星星光矢量进行三轴姿态解算。
为达到上述目的，本发明第二方面提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，实现上述所述的多星敏感器在轨热变形修正方法。
为达到上述目的，本发明第三方面提供了一种电子设备，包括处理器和存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现上述所述的多星敏感器在轨热变形修正方法。
本发明至少具有以下技术效果：
通过以第一星敏感器为基准，确定第二星敏感器和第三星敏感器分别相对于第一星敏感器的姿态关系矩阵和地面主光轴夹角，然后提取当前时刻下各个星敏感器在各自测量坐标系中的主光轴矢量，并将其余星敏感器各自测量坐标系中的主光轴矢量分别转换至第一星敏感器测量坐标系，再确定当前时刻转换至第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量与第一星敏感器测量坐标系下的主光轴矢量之间的主光轴矢量夹角，以根据主光轴矢量夹角、地面主光轴夹角、第一星敏感器测量坐标系下的主光轴矢量和转换至第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量，确定当前时刻其余星敏感器相对于第一星敏感器的热变形修正实时矩阵，从而可根据当前时刻热变形修正实时矩阵和姿态关系矩阵将其余星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，以进行三轴姿态解算，完成在轨热变形修正。
进一步的，本发明还可获取多个时刻下的其余星敏感器相对于第一星敏感器的多个热变形修正实时矩阵，然后通过姿态关系矩和多个时刻下的多个热变形修正实时矩阵将其余星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，并从中选取空间优越的恒星星光矢量进行三轴姿态解算，从而可实现共基准面安装平台的热变形的在轨精确测量，并可辅助载荷成像模型的精确建立，进而提高遥感图像品质，提高定位精度，并在应用上具有较高的普适性。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
图1为本发明一实施例提供的多星敏感器在轨热变形修正方法的流程图。
具体实施方式
下面详细描述本实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
目前，现有的针对星敏感器的热变形修正，有的采用在地面进行模拟空间热环境试验，获得星敏感器的热变形参数，通过热变形参数拟合出修正矩阵进行修正，有的采用改进星敏感器的材料，控制热变形，有的基于在轨遥测数据，对星敏感器在轨结构热变形的真实情况进行评估，并采用傅里叶级数建立热变形模型并进行参数辨识，然后根据辨识得到的模型参数对星敏感器进行在轨修正。而本发明实施例的多星敏感器在轨热变形修正方法采用的是有效利用惯性空间中的星光矢量信息，使用多台星敏感器对星空成像观测完成安装结构热变形的在轨量测，辅助载荷成像模型的精确建立，进而使遥感图像品质得到提高，获得更高的定位精度。并且，该方法能够直接反映共基准面安装平台的变形程度，而现有的卫星在轨标定多是关注星敏感器与载荷相机之间的关联角变化，无法直接反映安装平台的热变形程度，或是需要惯性敏感器配合星敏感器进行系统常值误差校正，无法将热变形的影响从整体误差中剥离出。由于，该方法能够直接反映共基准面安装平台的变形程度，所以该方法可适用于多星敏感器与载荷共基准面安装的微纳卫星平台。因为，结合当前微纳卫星平台的轻小型设计要求，多数微纳卫星的星敏感器与载荷相机均为共基准面安装，因此该方法具有很高的普遍适用性。
下面参考附图描述本实施例的多星敏感器在轨热变形修正方法以及存储介质和电子设备。
图1为本发明实施例提供的多星敏感器在轨热变形修正方法的流程图。
如图1所示，该方法包括：
步骤S1：建立三台星敏感器的各自测量坐标系，基于测量坐标系，以第一星敏感器为基准，确定第二星敏感器和第三星敏感器分别相对于第一星敏感器的姿态关系矩阵和地面主光轴夹角。
具体的，可在地面定义星敏感器测量坐标系CS，并以三台星敏感器中的一台星敏感器如A1为基准，获得其它两台星敏感器相对于星敏感器A1的姿态关系矩阵。
例如，在恒温洁净光学暗室内，将三台星敏感器如A1、A2和A3安装后，可对三台星敏感器分别进行地面标定，各自完成测量校准。然后，建立三台星敏感器的各自测量坐标系CA1、CA2、CA3，其中，定义星敏感器的主光轴为坐标系的X轴。进一步的，可基于各自的测量坐标系，以第一星敏感器A1为基准，确定第二星敏感器A2和第三星敏感器A3分别相对于第一星敏感器A1的姿态关系矩阵
Figure BDA0003347200180000051
并且，以第一星敏感器A1的主光轴为基准，确定第二星敏感器A2和第三星敏感器A3的主光轴分别相对于第一星敏感器A1的主光轴的夹角即地面主光轴夹角，分别为α标、β标。
步骤S2：三台星敏感器在同一时刻对各自视场内的恒星成像，并提取当前时刻下三台星敏感器在各自测量坐标系中的主光轴矢量。
例如，第一星敏感器测量坐标系下的主光轴矢量可表示为
Figure BDA0003347200180000052
第二星敏感器测量坐标系下的主光轴矢量可表示为
Figure BDA0003347200180000053
第三星敏感器测量坐标系下的主光轴矢量可表示为
Figure BDA0003347200180000054
步骤S3：将第二星敏感器和第三星敏感器各自测量坐标系中的主光轴矢量分别转换至第一星敏感器测量坐标系。
本实施例中，可将三台星敏感器各自测量坐标系下的三个主光轴矢量
Figure BDA0003347200180000061
和
Figure BDA0003347200180000062
统一至第一星敏感器A1的测量坐标系下，即将第二星敏感器A2和第三星敏感器A3各自测量坐标系中的主光轴矢量分别转换至第一星敏感器A1的测量坐标系。具体而言，第二星敏感器A2和第三星敏感器A3各自测量坐标系中的主光轴矢量可分别通过相应的姿态关系矩阵转换至第一星敏感器测量坐标系。三台星敏感器统一后的主光轴矢量分别为
Figure BDA0003347200180000063
Figure BDA0003347200180000064
步骤S4：确定当前时刻转换至第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量与第一星敏感器测量坐标系下的主光轴矢量之间的主光轴矢量夹角。
具体的，可以第一星敏感器A1测量坐标系下的主光轴矢量为基准，同时提取另外两个星敏感器主光轴矢量与其的夹角。
其中，可根据如下公式确定主光轴矢量夹角：
Figure BDA0003347200180000065
其中，α为第二星敏感器转换后的主光轴矢量夹角，β为第三星敏感器转换后的主光轴矢量夹角，
Figure BDA0003347200180000066
为第一星敏感器测量坐标系下的主光轴矢量，
Figure BDA0003347200180000067
为第二星敏感器测量坐标系下的主光轴矢量，
Figure BDA0003347200180000068
为第三星敏感器测量坐标系下的主光轴矢量，
Figure BDA0003347200180000069
和
Figure BDA00033472001800000610
分别为第二星敏感器和第三星敏感器相对于第一星敏感器的姿态关系矩阵。
步骤S5：根据主光轴矢量夹角、地面主光轴夹角、第一星敏感器测量坐标系下的主光轴矢量和转换至第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量，确定当前时刻其余星敏感器相对于第一星敏感器的热变形修正实时矩阵。
具体的，可根据如下公式确定当前时刻其余星敏感器相对于第一星敏感器的热变形修正实时矩阵：
Figure BDA0003347200180000071
其中，αi和βi分别为第i时刻第二星敏感器和第三星敏感器的主光轴矢量夹角，α标和β标分别为第二星敏感器和第三星敏感器相对于第一星敏感器的地面主光轴夹角，
Figure BDA0003347200180000072
和
Figure BDA0003347200180000073
分别为第i时刻第二星敏感器和第三星敏感器相对于第一星敏感器的热变形修正实时矩阵，
Figure BDA0003347200180000074
为转换至第一星敏感器测量坐标系下的第二星敏感器的主光轴矢量，
Figure BDA0003347200180000075
为转换至第一星敏感器测量坐标系下的第三星敏感器的主光轴矢量。
具体而言，可将当前时刻如第i时刻第二星敏感器的主光轴矢量夹角αi与步骤S1中确定的第二星敏感器相对于第一星敏感器的地面主光轴夹角α标进行对比，如进行作差得到第二星敏感器的主光轴矢量偏角。然后，根据该主光轴矢量偏角确定当前时刻第二星敏感器相对于第一星敏感器的热变形修正实时矩阵
Figure BDA0003347200180000076
类似的，可将当前时刻如第i时刻第三星敏感器的主光轴矢量夹角βi与步骤S1中确定的第三星敏感器相对于第一星敏感器的地面主光轴夹角β标进行对比，如进行作差得到第三星敏感器的主光轴矢量偏角。然后，根据该主光轴矢量偏角确定当前时刻第三星敏感器相对于第一星敏感器的热变形修正实时矩阵
Figure BDA0003347200180000077
步骤S6：获取各个星敏感器各自测量坐标系下的恒星星光矢量，其余星敏感器通过相对于第一星敏感器的姿态关系矩阵和当前时刻下的热变形修正实时矩阵将各自的恒星星光矢量投影至第一星敏感器测量坐标系，以进行三轴姿态解算。
具体的，可采用三台星敏感器实时对星空成像，然后分别获取各个星敏感器各自测量坐标系下的恒星星光矢量。第二星敏感器通过相对于第一星敏感器的姿态关系矩阵和当前时刻下的相应的热变形修正实时矩阵将第二星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，第三星敏感器通过相对于第一星敏感器的姿态关系矩阵和当前时刻下的相应的热变形修正实时矩阵将第三星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，以进行三轴姿态解算。
进一步的，本发明实施例的多星敏感器在轨热变形修正方法还包括：重复步骤S2～S5获取若干个时刻第二星敏感器相对于第一星敏感器的若干个热变形修正实时矩阵，以及第三星敏感器相对于第一星敏感器的若干个热变形修正实时矩阵。
具体的，可提取多时刻的三个主光轴矢量，并进行主光轴矢量夹角提取，以提取角变化的周期性规律，并将提取出来的周期性规律体现在热变形对星敏感器测量矩阵的影响中。
可选的，除了通过单一时刻下的热变形修正实时矩阵进行在轨热变形修正之外，还可通过多个时刻下的热变形修正实时矩阵进行在轨热变形修正。即言，除了获取当前时刻下的三台星敏感器在各自测量坐标系中的主光轴矢量之外，还可获取多个时刻如t1、t2、t3、t4......下的三台星敏感器在各自测量坐标系中的主光轴矢量，并将多个时刻下的第二星敏感器和第三星敏感器各自测量坐标系中的主光轴矢量分别转换至第一星敏感器测量坐标系，进行主光轴矢量夹角计算，得到多个时刻下第二星敏感器和第三星敏感器的主光轴矢量夹角如(α1、β1)、(α2、β2)、(α3、β3)、(α4、β4)、......，然后将提取出来的周期性规律体现在热变形对星敏感器测量矩阵的影响中。
具体而言，可分别将各个时刻下第二星敏感器和第三星敏感器的主光轴矢量夹角与地面主光轴夹角(α标、β标)进行对比，获得各个时刻的角度差值为(α1-α标、β1-β标)、(α2-α标、β2-β标)、(α3-α标、β3-β标)、(α4-α标、β4-β标)……，然后确定第二星敏感器相对于第一星敏感器的热变形修正实时矩阵
Figure BDA0003347200180000081
第三星敏感器相对于第一星敏感器的热变形修正实时矩阵
Figure BDA0003347200180000082
在本发明的一个实施例中，所述方法还包括：获取各个星敏感器各自测量坐标系下的恒星星光矢量，其余星敏感器通过相对于第一星敏感器的姿态关系矩阵和若干个时刻下的若干个热变形修正实时矩阵将各自的恒星星光矢量投影至第一星敏感器测量坐标系，并从投影后的若干个时刻下的恒星星光矢量中选取空间优越的恒星星光矢量进行三轴姿态解算。
其中，可采用三台星敏感器实时对星空成像，分别获取多个恒星星光矢量在星敏感器中的恒星星光矢量，并且通过坐标转换矩阵
Figure BDA0003347200180000083
以及各个时刻的热变形修正矩阵
Figure BDA0003347200180000084
以及
Figure BDA0003347200180000091
将恒星星光矢量投影至第一星敏感器的测量坐标系下，并选取空间分布优越的恒星星光矢量，进行三轴姿态解算。
具体而言，可采用三台星敏感器实时对星空成像，然后分别获取各个星敏感器各自测量坐标系下的恒星星光矢量。第二星敏感器通过相对于第一星敏感器的姿态关系矩阵和各个时刻下的热变形修正实时矩阵
Figure BDA0003347200180000092
将第二星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，第三星敏感器通过相对于第一星敏感器的姿态关系矩阵和各个时刻下的热变形修正实时矩阵
Figure BDA0003347200180000093
将第三星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，然后从投影后的各个时刻下的恒星星光矢量中选取空间优越的恒星星光矢量进行三轴姿态解算。
本发明实施例的多星敏感器在轨热变形修正方法通过以第一星敏感器为基准，确定第二星敏感器和第三星敏感器分别相对于第一星敏感器的姿态关系矩阵和地面主光轴夹角，然后提取当前时刻下各个星敏感器在各自测量坐标系中的主光轴矢量，并将其余星敏感器各自测量坐标系中的主光轴矢量分别转换至第一星敏感器测量坐标系，再确定当前时刻转换至第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量与第一星敏感器测量坐标系下的主光轴矢量之间的主光轴矢量夹角，以根据主光轴矢量夹角、地面主光轴夹角、第一星敏感器测量坐标系下的主光轴矢量和转换至第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量，确定当前时刻其余星敏感器相对于第一星敏感器的热变形修正实时矩阵，从而可根据当前时刻热变形修正实时矩阵和姿态关系矩阵将其余星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，以进行三轴姿态解算，完成在轨热变形修正。
进一步的，本发明还可获取多个时刻下的其余星敏感器相对于第一星敏感器的多个热变形修正实时矩阵，然后通过姿态关系矩和多个时刻下的多个热变形修正实时矩阵将其余星敏感器测量坐标系下的恒星星光矢量投影至第一星敏感器测量坐标系，并从中选取空间优越的恒星星光矢量进行三轴姿态解算，从而可有效提高星敏感器的姿态确定精度，并且针对星敏感器载荷共基准面安装的卫星平台，无需载荷参与测量过程，可获得共基准面的热变形量，为载荷的高品质应用奠定基础，且具有较高的普适性。
进一步地，本实施例还可提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时，可实现上述的多星敏感器在轨热变形修正方法。
具体地，存储介质可以是非临时性计算机可读存储介质，例如，非临时性计算机可读存储介质可以是只读存储器、随机存取存储器、只读光盘、磁带、软盘和光数据存储设备等。
进一步地，本实施例还可提供一种电子设备，其包括处理器和存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，可实现上述的多星敏感器在轨热变形修正方法。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
尽管本发明的内容已经通过上述优选实施例作了详细介绍，但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后，对于本发明的多种修改和替代都将是显而易见的。因此，本发明的保护范围应由所附的权利要求来限定。

Claims (8)

1.一种多星敏感器在轨热变形修正方法，其特征在于，包括：
步骤S1：建立三台星敏感器的各自测量坐标系，基于所述测量坐标系，以第一星敏感器为基准，确定第二星敏感器和第三星敏感器分别相对于所述第一星敏感器的姿态关系矩阵和地面主光轴夹角；
步骤S2：三台星敏感器在同一时刻对各自视场内的恒星成像，并提取当前时刻下三台星敏感器在各自测量坐标系中的主光轴矢量；
步骤S3：将所述第二星敏感器和所述第三星敏感器各自测量坐标系中的主光轴矢量分别转换至所述第一星敏感器测量坐标系；
步骤S4：确定当前时刻转换至所述第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量与第一星敏感器测量坐标系下的主光轴矢量之间的主光轴矢量夹角；
步骤S5：根据所述主光轴矢量夹角、所述地面主光轴夹角、第一星敏感器测量坐标系下的所述主光轴矢量和转换至所述第一星敏感器测量坐标系下的其余星敏感器的主光轴矢量，确定当前时刻其余星敏感器相对于所述第一星敏感器的热变形修正实时矩阵；
步骤S6：获取各个星敏感器各自测量坐标系下的恒星星光矢量，其余星敏感器通过相对于所述第一星敏感器的姿态关系矩阵和当前时刻下的热变形修正实时矩阵将各自的恒星星光矢量投影至所述第一星敏感器测量坐标系，以进行三轴姿态解算。
2.如权利要求1所述的多星敏感器在轨热变形修正方法，其特征在于，所述步骤S3中，所述第二星敏感器和所述第三星敏感器各自测量坐标系中的主光轴矢量分别通过相应的所述姿态关系矩阵转换至所述第一星敏感器测量坐标系。
3.如权利要求2所述的多星敏感器在轨热变形修正方法，其特征在于，所述步骤S4中，主光轴矢量夹角采用如下公式表示：
Figure FDA0003347200170000021
Figure FDA0003347200170000022
其中，α为第二星敏感器转换后的主光轴矢量夹角，β为第三星敏感器转换后的主光轴矢量夹角，
Figure FDA0003347200170000023
为第一星敏感器测量坐标系下的主光轴矢量，
Figure FDA0003347200170000024
为第二星敏感器测量坐标系下的主光轴矢量，
Figure FDA0003347200170000025
为第三星敏感器测量坐标系下的主光轴矢量，
Figure FDA0003347200170000026
和
Figure FDA0003347200170000027
分别为第二星敏感器和第三星敏感器相对于第一星敏感器的姿态关系矩阵。
4.如权利要求3所述的多星敏感器在轨热变形修正方法，其特征在于，所述步骤S5中，当前时刻其余星敏感器相对于所述第一星敏感器的热变形修正实时矩阵采用如下公式表示：
Figure FDA0003347200170000028
Figure FDA0003347200170000029
其中，αi和βi分别为第i时刻第二星敏感器和第三星敏感器的主光轴矢量夹角，α标和β标分别为第二星敏感器和第三星敏感器相对于第一星敏感器的地面主光轴夹角，
Figure FDA00033472001700000210
和
Figure FDA00033472001700000211
分别为第i时刻第二星敏感器和第三星敏感器相对于第一星敏感器的热变形修正实时矩阵，
Figure FDA00033472001700000212
为转换至第一星敏感器测量坐标系下的第二星敏感器的主光轴矢量，
Figure FDA00033472001700000213
为转换至第一星敏感器测量坐标系下的第三星敏感器的主光轴矢量。
5.如权利要求1所述的多星敏感器在轨热变形修正方法，其特征在于，还包括：
重复步骤S2～S5获取若干个时刻第二星敏感器相对于第一星敏感器的若干个热变形修正实时矩阵，以及第三星敏感器相对于第一星敏感器的若干个热变形修正实时矩阵。
6.如权利要求5所述的多星敏感器在轨热变形修正方法，其特征在于，还包括：
获取各个星敏感器各自测量坐标系下的恒星星光矢量，其余星敏感器通过相对于所述第一星敏感器的姿态关系矩阵和若干个时刻下的若干个热变形修正实时矩阵将各自的恒星星光矢量投影至所述第一星敏感器测量坐标系，并从投影后的若干个时刻下的恒星星光矢量中选取空间优越的恒星星光矢量进行三轴姿态解算。
7.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现如权利要求1-6中任一项所述的多星敏感器在轨热变形修正方法。
8.一种电子设备，其特征在于，包括处理器和存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现如权利要求1-6中任一项所述的多星敏感器在轨热变形修正方法。
CN202111326115.6A 2021-11-10 2021-11-10 多星敏感器在轨热变形修正方法以及存储介质和电子设备 Active CN114234962B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111326115.6A CN114234962B (zh) 2021-11-10 2021-11-10 多星敏感器在轨热变形修正方法以及存储介质和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111326115.6A CN114234962B (zh) 2021-11-10 2021-11-10 多星敏感器在轨热变形修正方法以及存储介质和电子设备

Publications (2)

Publication Number Publication Date
CN114234962A true CN114234962A (zh) 2022-03-25
CN114234962B CN114234962B (zh) 2023-09-12

Family

ID=80748914

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111326115.6A Active CN114234962B (zh) 2021-11-10 2021-11-10 多星敏感器在轨热变形修正方法以及存储介质和电子设备

Country Status (1)

Country Link
CN (1) CN114234962B (zh)

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0795806A2 (de) * 1995-08-11 1997-09-17 Daimler-Benz Aerospace Aktiengesellschaft Vorrichtung zur Lageregelung und Stabilisierung eines Raumfahrzeuges mit Entkopplung der Roll- und Gierdynamik, und Verfahren dazu
US6470270B1 (en) * 2000-11-07 2002-10-22 The Boeing Company Lost-in-space stellar attitude acquisition using multiple star trackers
FR2833364A1 (fr) * 2001-12-12 2003-06-13 Centre Nat Etd Spatiales Procede d'estimation d'attitude de dispositif spatial a precision amelioree
US6691033B1 (en) * 2000-07-26 2004-02-10 Hughes Electronics Corporation System and method for calibrating inter-star-tracker misalignments in a stellar inertial attitude determination system
US20060149474A1 (en) * 2005-01-03 2006-07-06 Needelman David D Real-time refinement method of spacecraft star tracker alignment estimates
CN102252673A (zh) * 2011-06-03 2011-11-23 哈尔滨工业大学 一种星敏感器在轨光行差的修正方法
CN106338296A (zh) * 2016-11-04 2017-01-18 上海航天控制技术研究所 一种双星敏感器在轨实时交互的修正方法
CN106441283A (zh) * 2016-10-12 2017-02-22 上海航天控制技术研究所 一种自主标较的星矢量级多视场星敏感器数据融合方法
CN106568436A (zh) * 2016-10-26 2017-04-19 上海航天控制技术研究所 一种利用陀螺修正星敏感器轨道周期系统误差的方法
CN106940196A (zh) * 2017-03-30 2017-07-11 上海航天控制技术研究所 一种星敏感器安装热变形修正方法
CN107024228A (zh) * 2017-04-12 2017-08-08 上海航天控制技术研究所 一种星敏感器非高频误差在轨修正方法
CN107389098A (zh) * 2017-08-22 2017-11-24 长光卫星技术有限公司 一种航天器星敏感器安装矩阵在轨标定方法
CN109696179A (zh) * 2018-11-15 2019-04-30 上海航天控制技术研究所 一种遥感卫星星敏感器热弹性误差估计方法
CN110553667A (zh) * 2019-09-11 2019-12-10 上海航天控制技术研究所 一种对星敏感器进行热变形补偿的方法
CN111380567A (zh) * 2020-03-24 2020-07-07 上海卫星工程研究所 一种卫星星敏感器指向在轨热变形修正系统
CN112061425A (zh) * 2020-09-08 2020-12-11 上海航天控制技术研究所 一种敏捷小卫星星敏感器受地气光干扰的规避方法
CN112414398A (zh) * 2020-10-22 2021-02-26 中国西安卫星测控中心 一种在轨卫星标定星敏感器测量精度的方法
CN112781616A (zh) * 2020-12-24 2021-05-11 中国人民解放军国防科技大学 星敏感器在轨测量低频误差分析方法、装置和存储介质
CN113252075A (zh) * 2021-05-12 2021-08-13 中国科学院微小卫星创新研究院 一种基于轨道周期特性的星敏感器精度修正方法
CN113405567A (zh) * 2021-05-31 2021-09-17 中国人民解放军61540部队 一种重力卫星星敏感器安装矩阵在轨定标方法及系统

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0795806A2 (de) * 1995-08-11 1997-09-17 Daimler-Benz Aerospace Aktiengesellschaft Vorrichtung zur Lageregelung und Stabilisierung eines Raumfahrzeuges mit Entkopplung der Roll- und Gierdynamik, und Verfahren dazu
US6691033B1 (en) * 2000-07-26 2004-02-10 Hughes Electronics Corporation System and method for calibrating inter-star-tracker misalignments in a stellar inertial attitude determination system
US6470270B1 (en) * 2000-11-07 2002-10-22 The Boeing Company Lost-in-space stellar attitude acquisition using multiple star trackers
FR2833364A1 (fr) * 2001-12-12 2003-06-13 Centre Nat Etd Spatiales Procede d'estimation d'attitude de dispositif spatial a precision amelioree
US20060149474A1 (en) * 2005-01-03 2006-07-06 Needelman David D Real-time refinement method of spacecraft star tracker alignment estimates
CN102252673A (zh) * 2011-06-03 2011-11-23 哈尔滨工业大学 一种星敏感器在轨光行差的修正方法
CN106441283A (zh) * 2016-10-12 2017-02-22 上海航天控制技术研究所 一种自主标较的星矢量级多视场星敏感器数据融合方法
CN106568436A (zh) * 2016-10-26 2017-04-19 上海航天控制技术研究所 一种利用陀螺修正星敏感器轨道周期系统误差的方法
CN106338296A (zh) * 2016-11-04 2017-01-18 上海航天控制技术研究所 一种双星敏感器在轨实时交互的修正方法
CN106940196A (zh) * 2017-03-30 2017-07-11 上海航天控制技术研究所 一种星敏感器安装热变形修正方法
CN107024228A (zh) * 2017-04-12 2017-08-08 上海航天控制技术研究所 一种星敏感器非高频误差在轨修正方法
CN107389098A (zh) * 2017-08-22 2017-11-24 长光卫星技术有限公司 一种航天器星敏感器安装矩阵在轨标定方法
CN109696179A (zh) * 2018-11-15 2019-04-30 上海航天控制技术研究所 一种遥感卫星星敏感器热弹性误差估计方法
CN110553667A (zh) * 2019-09-11 2019-12-10 上海航天控制技术研究所 一种对星敏感器进行热变形补偿的方法
CN111380567A (zh) * 2020-03-24 2020-07-07 上海卫星工程研究所 一种卫星星敏感器指向在轨热变形修正系统
CN112061425A (zh) * 2020-09-08 2020-12-11 上海航天控制技术研究所 一种敏捷小卫星星敏感器受地气光干扰的规避方法
CN112414398A (zh) * 2020-10-22 2021-02-26 中国西安卫星测控中心 一种在轨卫星标定星敏感器测量精度的方法
CN112781616A (zh) * 2020-12-24 2021-05-11 中国人民解放军国防科技大学 星敏感器在轨测量低频误差分析方法、装置和存储介质
CN113252075A (zh) * 2021-05-12 2021-08-13 中国科学院微小卫星创新研究院 一种基于轨道周期特性的星敏感器精度修正方法
CN113405567A (zh) * 2021-05-31 2021-09-17 中国人民解放军61540部队 一种重力卫星星敏感器安装矩阵在轨定标方法及系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JØRGENSEN P S: "Denver T and Betto M. 2005b. In-flight quality and accuracy of attitude measurements from the CHAMP advanced stellar compass", 《ACTA ASTRONAUTICA》, vol. 56, no. 1, pages 181 - 186 *
MEI, S., JIANG, R., JI, J., SUN, J., PENG, Y., & ZHANG, Y: "Invariant feature extraction for image classification via multi-channel convolutional neural network", 《 2017 INTERNATIONAL SYMPOSIUM ON INTELLIGENT SIGNAL PROCESSING AND COMMUNICATION SYSTEMS (ISPACS)》, pages 491 - 495 *
胡雄超;毛晓楠;吴永康;闫晓军;余路伟;王兆龙;: "基于亚像元坐标的像素频率误差补偿方法", 《红外与激光工程》, vol. 46, no. 07, pages 1 - 4 *
薛喜平等: "深空探测天文自主导航技术综述", 《天文研宄与技术》, vol. 14, no. 3, pages 382 - 391 *
钟金凤等: "双星敏感器在轨相对热变形分析及修正", 《遥感学报》, vol. 22, no. 2, pages 198 - 202 *

Also Published As

Publication number Publication date
CN114234962B (zh) 2023-09-12

Similar Documents

Publication Publication Date Title
CN106403902B (zh) 一种星地协同的光学卫星在轨实时几何定位方法及系统
US6735348B2 (en) Apparatuses and methods for mapping image coordinates to ground coordinates
CN102741706B (zh) 地理参照图像区域的方法
JP3428539B2 (ja) 人工衛星の姿勢角センサ較正装置
CN107192376B (zh) 基于帧间连续性的无人机多帧图像目标定位校正方法
CN109708649B (zh) 一种遥感卫星的姿态确定方法及系统
CN107728182B (zh) 基于相机辅助的柔性多基线测量方法和装置
US7821535B2 (en) Information processing method and apparatus
CN107144293A (zh) 一种视频卫星面阵相机的几何定标方法
CN106885585B (zh) 一种基于光束法平差的星载摄影测量系统一体化检校方法
CN110986998B (zh) 一种基于有理函数模型的卫星视频相机在轨几何定标方法
CN112050806B (zh) 一种移动车辆的定位方法及装置
CN111912430B (zh) 高轨光学卫星的在轨几何定标方法、装置、设备及介质
CN111473802A (zh) 一种基于线阵推扫的光学传感器内方位元素定标方法
CN112581544B (zh) 基于参数优化的无公共视场相机标定方法
CN110793542A (zh) 一种基于广义探元指向角的面阵光学遥感卫星在轨几何定标方法
Zhou et al. Novel autonomous on-orbit calibration method for star sensors
CN108154535B (zh) 基于平行光管的摄像机标定方法
CN114838740A (zh) 一种考虑不同经纬度区域的卫星图像几何定标方法
Pi et al. On-orbit geometric calibration using a cross-image pair for the linear sensor aboard the agile optical satellite
CN114754798A (zh) 一种陀螺误差特性参数在轨辨识与标定方法
KR102265081B1 (ko) 벡터내적과 3차원 좌표변환을 이용한 카메라의 위치 및 자세를 결정하는 시스템
CN114234962B (zh) 多星敏感器在轨热变形修正方法以及存储介质和电子设备
CN115578448B (zh) 一种基于ccd观测数据批处理的天文定位方法及系统
Woodbury et al. Stellar positioning system (Part II): Improving accuracy during implementation

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