CN105318871B - 双星敏感器载体安装矩阵动态标定方法 - Google Patents
双星敏感器载体安装矩阵动态标定方法 Download PDFInfo
- Publication number
- CN105318871B CN105318871B CN201510754528.2A CN201510754528A CN105318871B CN 105318871 B CN105318871 B CN 105318871B CN 201510754528 A CN201510754528 A CN 201510754528A CN 105318871 B CN105318871 B CN 105318871B
- Authority
- CN
- China
- Prior art keywords
- star sensor
- star
- matrix
- under
- 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
Links
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
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Abstract
本发明涉及一种双星敏感器载体安装矩阵动态标定的方法,方法包括以下步骤:一是标定双星敏感器之间相对安装矩阵,及获取双星敏感器之间的相对位置关系;二是标定双星敏感器与载体安装矩阵,即通过双星敏感器联合求解获取某一星敏感器相对载体的位置关系。本发明解决了两个星敏感器单独标定安装矩阵时存在的沿星敏感器光轴方向旋转角标定精度低的问题;可以在载体运动情况下实行动态标定,尤其适用于船载星敏感器的动态标定。
Description
技术领域
本发明涉及一种双星敏感器安装矩阵动态标定的方法。属于天文导航领域。
背景技术
星敏感器以恒星为参照物,由探测单元在某一时刻对天空完成星图捕获,经过数据处理单元对星图进行恒星质心提取、星图识别、恒星跟踪、姿态计算等一系列处理,最终获得星敏感器相对惯性空间的高精度姿态信息,是一种高精度的姿态测量装置。
目前星敏感器广泛应用于航空航天飞行器姿态测量,随着高精度大视场星敏感器技术的成熟,星敏感器已经能够完成船体姿态测量。
星敏感器在使用之前,需要对其与载体之间的安装矩阵进行标定,标定精度越高,其获取的载体姿态数据精度越高。
关于传统的安装矩阵获取方法,一般是通过静态方法测量星敏感器的安装角,每次对一个星敏感器的安装矩阵进行单独标定。该方法存在两方面问题:一是在沿星敏感器光轴方向旋转角标定精度比较低;二是标定需要载体处于静态条件获取安装矩阵的条件较为严格。
发明内容
本发明所要解决的技术问题是针对上述现有技术提供一种双星敏感器载体安装矩阵动态标定方法,该方法通过获取星敏感器之间相对安装矩阵,在动态条件下完成星敏感器与载体间安装矩阵的标定。
本发明解决上述问题所采用的技术方案为:一种双星敏感器与载体间安装矩阵动态标定方法,所述方法包括以下步骤:
步骤一、标定双星敏感器之间相对安装矩阵
1、获取恒星在星敏感器坐标系下的观测矢量Wsi和恒星在地心惯性坐标系下的参考矢量Vii。
Vii=[cosαicosδi sinαicosδi sinδi]T (3)
f为星敏感器焦距,(x0,y0)为图像传感器的主点位置(像元数),(xui,yvi)为第i颗恒星像中心像素坐标,D为图像传感器像元尺寸,αi和δi为第i颗恒星在地心惯性坐标系(i系)下的赤经、赤纬分别。为地心惯性系到星敏感器坐标系(s系)的转换矩阵。
2、求取双星敏感应器之间相对安装矩阵
对星敏感器A1、A2的观测量与参考量,可列出如下公式:
其中分别为当地地理坐标系到星敏感器A1、A2坐标系的姿态转移矩阵。
令星敏感器A1、A2之间的相对安装矩阵为表示由A2坐标系到A1坐标系的旋转矩阵。对于A2坐标系下观测量Ws2j,其在A1坐标系下的值Ws1j可由下式计算:
对(6)式右边继续展开,并将(5)式代入其中可得:
对A1坐标系下的两个星敏感器的观测向量Ws1i、Ws1j进行点乘可得到
观察(8)式,右边为已知,左边Ws1i为已知,Ws1j可由包含的(6)式计算,多个点构成联立方程组,根据约束最小二乘即可求得
步骤二、标定双星敏感器与载体间安装矩阵
双星敏感器与载体之间的安装矩阵标定就是求或以为例给出计算方法如下:
1、载体配置了捷联惯导,通过三轴陀螺解算得到载体坐标系到当地地理坐标系(g系)下的姿态转移矩阵
2、根据可将星敏感器A1、A2在g系下的参考量Vgi转换到载体坐标系(b系),公式如下:
3、星敏感器A2在s1系下的观测矢量可通过(6)式计算得到,与星敏感器A1在s1系下的观测矢量共同构成双星敏感器在s系下的观测矢量,记为Ws1i。
4、当双星敏感器视场内总识别星数目大于3时,可根据下式求得安装矩阵
与现有技术相比,本发明的优点在于:
一是标定精度高,解决了两个星敏感器单独标定安装矩阵时存在的沿星敏感器光轴方向旋转角标定精度低的问题;二是标定条件不受限制,可以在载体运动情况下实行动态标定,尤其适用于船载星敏感器的动态标定。
附图说明
图1相对安装角误差与实验次数关系。
图2星敏感器A1安装角误差对比结果。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明涉及一种双星敏感器载体安装矩阵动态标定的方法,方法包括以下步骤:一是标定双星敏感器之间相对安装矩阵,及获取双星敏感器之间的相对位置关系;二是标定双星敏感器与载体安装矩阵,即通过双星敏感器联合求解获取某一星敏感器相对载体的位置关系。
为便于本发明内容描述,在此对发明中涉及的坐标系定义、相关概念与符号含义等进行简要说明。
一、常用坐标系及其转换关系
1、星敏感器固联安装在载体上,根据其安装矩阵可以得到载体坐标系(b系)到星敏感器坐标系(s系)的姿态转移矩阵
2、载体配置了捷联惯导,通过三轴陀螺解算得到载体坐标系到当地地理坐标系(g系)下的姿态转移矩阵
3、地心惯性坐标系(i系)到当地地理坐标系g系的姿态转移矩阵可以由观测时刻和测量点的经纬度等信息计算得到,记为
二、实施过程
步骤一、标定双星敏感器之间相对安装矩阵
1、获取恒星在星敏感器坐标系下的观测矢量Wsi和恒星在地心惯性坐标系下的参考矢量Vii。
Vii=[cosαicosδi sinαicosδi sinδi]T (3)
f为星敏感器焦距,(x0,y0)为图像传感器的主点位置(像元数),(xui,yvi)为第i颗恒星像中心像素坐标,D为图像传感器像元尺寸,αi和δi为第i颗恒星在地心惯性坐标系(i系)下的赤经、赤纬分别。为地心惯性系到星敏感器坐标系(s系)的转换矩阵。
2、求取双星敏感应器之间相对安装矩阵
设星敏感器A1和A2分别识别到4颗恒星,A1恒星编号一次为1、2、3、4,A2恒星编号5、6、7、8。
对星敏感器A1、A2的观测量与参考量,可列出如下公式:
其中分别为当地地理坐标系到星敏感器A1、A2坐标系的姿态转移矩阵。
令星敏感器A1、A2之间的相对安装矩阵为表示由A2坐标系到A1坐标系的旋转矩阵。对于A2坐标系下观测量Ws2j,其在A1坐标系下的值Ws1j可由下式计算:
对(6)式右边继续展开,并将(5)式代入其中可得:
对A1坐标系下的两个星敏感器的观测向量Ws1i、Ws1j进行点乘可得到
观察(8)式,右边为已知,左边Ws1i为已知,Ws1j可由包含的(6)式计算,多个点构成联立方程组,根据约束最小二乘即可求得
步骤二、标定双星敏感器与载体间安装矩阵
双星敏感器与载体之间的安装矩阵标定就是求或以为例给出计算方法如下:
1、载体配置了捷联惯导,通过三轴陀螺解算得到载体坐标系到当地地理坐标系(g系)下的姿态转移矩阵
2、根据可将星敏感器A1、A2在g系下的参考量Vgi转换到b系,公式如下:
3、星敏感器A2在s1系下的观测矢量可通过(6)式计算得到,与星敏感器A1在s1系下的观测矢量共同构成双星敏感器在s1系下的观测矢量,记为Ws1i,i=1...8;
4、当双星敏感器视场内总识别星数目大于3时,可根据下式求得安装矩阵
实施例:
1、仿真试验条件及参数
星敏感器参数如表1所示:
表1星敏感器参数
船体位置姿态信息如表2所示:
表2船体位置姿态信息
恒星赤经、赤纬信息如表3所示,其中恒星1-4在星敏感器A1视场中,恒星5-8在星敏感器A2视场中。
表3恒星位置
观测时刻其它参数如表4所示:
表4观测时刻相关参数
星敏感器安装角如表5所示:
表5星敏感器安装角
2、仿真试验结果及相关数据
3、仿真试验精度分析
根据前述仿真条件,可以计算出8颗恒星在理想情况和实际情况(包含蒙气差影响情况)下在各自星敏感器像面上的成像位置,在仿真实验中,在实际情况的成像点x、y方向分别加入均值为0、标准差为0.1的高斯白噪声,之后根据参考量蒙气差修正方法对每个像点进行蒙气差修正。
根据双星敏感器之间相对安装矩阵动态标定方法对进行计算,随机进行100次实验,对计算得到的相对安装角结果求平均,与仿真中设置的真值比较,角度误差如图1所示。
由图1可以看出,相对安装角误差随试验组数增加而减小,仿真条件下经过50次以后,相对安装角误差可降低到5〞以下,100次以后,相对安装角误差可降低到3〞以下。
双星敏感器之间相对安装矩阵标定以后,可以对双星敏感器与载体之间的安装矩阵进行标定,以星敏感器A1为例,进行100次实验,对计算得到的A1安装角结果求平均,与仿真中设置的真值比较,角度误差如图2所示。
图2中实线表示动态标定计算得到的星敏感器A1安装角误差曲线,虚线表示利用传统的单星敏感器信息计算得到的星敏感器A1安装角误差曲线。由图2可以看出,动态标定结果精度明显高于传统方法,安装角误差随试验组数增加而减小,仿真条件下经过10次以后,安装角误差可降低到1〞以下,20次以后,相对安装角误差可降低到0.5〞以下。
除上述实施例外,本发明还包括有其他实施方式,凡采用等同变换或者等效替换方式形成的技术方案,均应落入本发明权利要求的保护范围之内。
Claims (1)
1.一种双星敏感器与载体间安装矩阵动态标定方法,所述方法包括以下步骤:
步骤一、标定双星敏感器之间相对安装矩阵
1)、获取恒星在星敏感器坐标系下的观测矢量Wsi和恒星在地心惯性坐标系下的参考矢量Vii:
Vii=[cosαicosδi sinαicosδi sinδi]T (3)
f为星敏感器焦距,(x0,y0)为图像传感器的主点位置,(xui,yvi)为第i颗恒星像中心像素坐标,D为图像传感器像元尺寸,αi和δi为第i颗恒星在地心惯性坐标系下的赤经、赤纬分别,为地心惯性系到星敏感器坐标系的转换矩阵;
2)、求取双星敏感应器之间相对安装矩阵
令星敏感器A1、A2之间的相对安装矩阵为表示由A2坐标系到A1坐标系的旋转矩阵,对于A2坐标系下观测量Ws2j,其在A1坐标系下的值Ws1j可由下式计算:
对A1坐标系下的两个星敏感器的观测向量Ws1i、Ws1j进行点乘可得到
其中:Vgi代表星敏感器A1在当地地理坐标系下的参考量;Vgj代表星敏感器A2在当地地理坐标系下的参考量;
由(6)和(8)根据约束最小二乘即可求得
步骤二、标定双星敏感器与载体间安装矩阵
1)载体配置了捷联惯导,通过三轴陀螺解算得到载体坐标系到当地地理坐标系下的姿态转移矩阵
2)根据可将星敏感器A1、A2在g系下的参考量Vgi、Vgj转换到载体坐标系,公式如下:
3)当双星敏感器视场内总识别星数目大于3时,可根据下式求得安装矩阵
其中:Wsi为恒星在星敏感器坐标系下的观测矢量,Vbi为星敏感器在载体坐标系下的参考量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510754528.2A CN105318871B (zh) | 2015-11-09 | 2015-11-09 | 双星敏感器载体安装矩阵动态标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510754528.2A CN105318871B (zh) | 2015-11-09 | 2015-11-09 | 双星敏感器载体安装矩阵动态标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105318871A CN105318871A (zh) | 2016-02-10 |
CN105318871B true CN105318871B (zh) | 2018-03-30 |
Family
ID=55246777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510754528.2A Active CN105318871B (zh) | 2015-11-09 | 2015-11-09 | 双星敏感器载体安装矩阵动态标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105318871B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106441283B (zh) * | 2016-10-12 | 2019-03-26 | 上海航天控制技术研究所 | 一种自主标较的星矢量级多视场星敏感器数据融合方法 |
CN106338296B (zh) * | 2016-11-04 | 2019-03-26 | 上海航天控制技术研究所 | 一种双星敏感器在轨实时交互的修正方法 |
CN108072387B (zh) * | 2016-11-14 | 2021-06-22 | 上海微小卫星工程中心 | 一种低精度敏感器在轨偏差的地面校正方法及系统 |
CN107228683B (zh) * | 2017-06-27 | 2020-04-10 | 上海航天控制技术研究所 | 一种多星敏感器间慢变误差实时在轨修正方法 |
CN109470269B (zh) * | 2018-11-08 | 2020-09-04 | 上海航天控制技术研究所 | 空间目标测量机构的标定方法、标定设备及标定系统 |
CN111637879B (zh) * | 2020-04-22 | 2021-10-01 | 北京控制工程研究所 | 基于多维差异化误差特性分布的双星敏加权姿态确定方法 |
CN111637885B (zh) * | 2020-05-12 | 2021-10-01 | 北京控制工程研究所 | 一种船载白昼星敏感器定位算法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426017A (zh) * | 2011-11-03 | 2012-04-25 | 北京航空航天大学 | 一种基于星敏感器确定载体相对于地理坐标系姿态的方法 |
CN102865866A (zh) * | 2012-10-22 | 2013-01-09 | 哈尔滨工业大学 | 基于双星敏感器的卫星姿态确定方法及定姿误差分析方法 |
CN104406583A (zh) * | 2014-12-05 | 2015-03-11 | 中国人民解放军63680部队 | 双星敏感器联合确定载体姿态方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7451022B1 (en) * | 2006-12-28 | 2008-11-11 | Lockheed Martin Corporation | Calibration of ship attitude reference |
-
2015
- 2015-11-09 CN CN201510754528.2A patent/CN105318871B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426017A (zh) * | 2011-11-03 | 2012-04-25 | 北京航空航天大学 | 一种基于星敏感器确定载体相对于地理坐标系姿态的方法 |
CN102865866A (zh) * | 2012-10-22 | 2013-01-09 | 哈尔滨工业大学 | 基于双星敏感器的卫星姿态确定方法及定姿误差分析方法 |
CN104406583A (zh) * | 2014-12-05 | 2015-03-11 | 中国人民解放军63680部队 | 双星敏感器联合确定载体姿态方法 |
CN104792340A (zh) * | 2015-05-15 | 2015-07-22 | 哈尔滨工业大学 | 一种星敏感器安装误差矩阵与导航系统星地联合标定与校正的方法 |
Non-Patent Citations (3)
Title |
---|
基于双星敏感器的船体姿态测量系统设计;郭敬明等;《光电子技术》;20140331;第34卷(第1期);全文 * |
捷联惯导_星敏感器组合系统的在轨自标定方法研究;杨波等;《航天控制》;20100228;第28卷(第1期);全文 * |
星敏感器自主在轨标定算法;袁彦红等;《上海航天》;20081231(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN105318871A (zh) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105318871B (zh) | 双星敏感器载体安装矩阵动态标定方法 | |
CN105698764B (zh) | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 | |
CN104897175B (zh) | 多相机光学推扫卫星在轨几何定标方法及系统 | |
CN106500731B (zh) | 一种基于恒星模拟系统的船载经纬仪的标校方法 | |
CN111044074B (zh) | 基于外场观星的星敏感器标定装置与观星标定方法 | |
CN107728182B (zh) | 基于相机辅助的柔性多基线测量方法和装置 | |
CN108592945A (zh) | 一种惯性/天文组合系统误差的在线标定方法 | |
CN108845335A (zh) | 一种基于图像和导航信息的无人机地面目标定位方法 | |
CN104406583B (zh) | 双星敏感器联合确定载体姿态方法 | |
CN106708066A (zh) | 基于视觉/惯导的无人机自主着陆方法 | |
CN107806874B (zh) | 一种视觉辅助的捷联惯导极区初始对准方法 | |
CN104698486B (zh) | 一种分布式pos用数据处理计算机系统实时导航方法 | |
CN109459058B (zh) | 一种基于三轴转台的多视场星敏感器的标定方法 | |
CN102636081B (zh) | 一种基于视觉运动建模的传递对准方法及装置 | |
CN107945229B (zh) | 用于静止轨道对地观测卫星面阵仪器的恒星质心提取方法 | |
KR102075028B1 (ko) | 무인기 고속비행 중 정밀위치 영상 획득 장치 및 그를 이용한 정밀위치 획득 방법 | |
CN104864852B (zh) | 一种基于密集控制点的高分辨率卫星姿态颤振探测方法 | |
CN110793542A (zh) | 一种基于广义探元指向角的面阵光学遥感卫星在轨几何定标方法 | |
CN106679695B (zh) | 一种基于星敏感器的焦平面热变形的测试方法 | |
CN114972078A (zh) | 应用sar影像提升国产光学卫星影像无控几何质量方法及系统 | |
CN105424060B (zh) | 一种飞行器星敏感器与捷联惯组安装误差的测量方法 | |
CN102576459A (zh) | 估计机载或空中的推扫式传感器获得的线图像的偏移的方法 | |
CN104977021A (zh) | 一种星光导航环境误差补偿方法 | |
CN105424047B (zh) | 基于路标信息的航天器姿态指向误差辨识方法 | |
CN108106634B (zh) | 一种直接观星的星敏感器内参数标定方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |