一种基于恒星分布的星敏感器导航星库制作方法

叶志龙，孙朔冬，陈纾，吴迪，武斌

（1. 上海航天控制技术研究所，上海 201109；2. 上海市空间智能控制技术重点实验室，上海 201109）

0 引 言

星敏感器作为一种高精度的姿态测量器件，已广泛应用于卫星，深空探测等航天飞行器中，其通过匹配探测恒星信息与导航星库中的恒星信息，并采用QUEST（Quaternion Estimator）算法解算姿态。导航星库是星敏感器进行星图识别和跟踪的关键信息，其合理性和完备性对星敏感器工作的鲁棒性和精度都有着重要的影响。现有的导航星库制作方法主要有星等阈值法、星等加权法、正交网格法、回归选取法、二维精简索引分层构造法等[1-4]。星等阈值法是通过选取星等小于一定阈值的星构成导航星库，该方法虽然简单，也具有良好的匹配性，但是恒星分布的不均匀性会导致星图识别准确性下降。若采用提高星等阈值来改善此问题，则会导致星库容量增大，不利于星敏感器存储使用，也增加星库搜索的时间。星等加权法通过对星等较低的星赋予较大权值来提高导航星库恒星较少天区的星颗数，但依旧存在与星等阈值法类似的问题。正交网格法则是通过映射法将候选导航星映射从球面到平面上，并利用互不重叠的等面积正交网格将其划分，选取每个网格中最亮的星作为导航星。此方法可通过调节网格大小筛选导航，有较强的灵活性。但是未考虑到星敏感器的实际探测能力，导致一些网格中的导航星不可用。回归选取法则是将星等阈值法中的固定阈值变为动态阈值，可有效地筛选出导航星，但动态阈值的自适应性仍不够完善。二维精简索引分层构造法把导航星分为局部亮星和暗星。以局部亮星构造导航三角形，节约了导航星库空间，大大缩短了搜索时间，然而该方法对星敏感器的视场有一定的要求，仅适用于探测星等较低的第二代星敏感器。

由此可见，现有的导航星库多以星等和分布均匀性为依据选取导航星，此类方法虽然可以较为完备地筛选出导航星库，但未考虑星库中恒星的空间分布对于定姿的影响，往往会出现星敏感器视场划过某一天区的时候，其输出姿态的误差有一定幅度增大[5-8]。因此，本文考虑利用恒星的空间分布与定姿精度的关系，并结合星等和恒星分布的均匀性来选取制作导航星库。

1 星库的选择

星库记录了恒星位置、星等、光谱等天文参数。由于编制星库使用仪器、观测条件和处理方法不一样，不同星库中恒星的颗数不一样，恒星的位置精度也会有差异。星库是设计制作导航星库的基础，也是决定星敏感器定姿精度的关键因素之一。表1给出了天文领域实际使用的几种星库。

表1 常见的几种星库Table 1 Several common star datebase

表1中SAO星库和Hipparcos星库较容易获得，其中Hipparcos星库的精度较高，可达到0.000 8″，因此，本文选择Hipparcos星库为基础星库。从表1中可以看出，Hipparcos星库中恒星的数目近12万颗，数目过于庞大，且星敏感器光学系统的敏感能力有限。因此，可根据星敏感器的敏感极限来删减星库。另一方面，星库中存在的双星和变星，对姿态解算有一定影响，可删除双星和变星进一步精简星库。

2 基于姿态精度的选星策略

星敏感器的定姿误差来源主要有星点质心定位误差，光学系统误差以及定姿算法的转换误差等。为便于分析，将星点定位误差和光学系统误差认为是天文观测误差。星敏感器通过获取星点在星敏感器测量坐标系的坐标，以及星点对应于惯性坐标的坐标，通过Quest算法解算得到姿态信息。

2.1 恒星空间几何分布对定姿精度的影响

假设恒星在测量坐标系b下的三维坐标为在惯性坐标系i下的三维坐标为，则

其中：为观测坐标系b到惯性坐标系i的姿态转移矩阵。若记依据上式可知

星敏感器在观测星数不小于3时，可用最小二乘法求解姿态转移矩阵，即

其中：由于存在天文观测误差，恒星的实际观测矢量为其中∆W为天文观测误差。则

若以观测坐标系b为参考，记定姿误差矢量为当定姿误差角为一小量时，可表示为

式中b’为包含误差计算的测量坐标系。依据式（4）和式（5）可得

式中：记M的协方差矩阵为PM，误差矩阵的协方差矩阵为P∆，定姿误差矢量δφ的协方差矩阵为Pφ。

由式（7）和式（8）可知

天文观测误差∆W主要取决于敏感器件及星点提取等误差，与选星无关。因此，假设∆W一定，其方差为由式（3）、式（4）和式（9）可得

由式（9）和式（10）可知

以Quest算法的三星定姿为例，其中为3颗定姿星的惯性坐标系下的单位矢量，矩阵V一定程度上反映了恒星的空间几何分布[9–10]。假设惯性坐标i的原点为O，3颗星对应的单位矢量为Ov1，Ov2，Ov3。3颗恒星与原点构成的四面体体积Tvol为

将上式代入式（10），式（9），令则

由式（13）可知，定姿误差的方差与星矢量构成的体积以及矩阵B的迹相关。为了说明星矢量构成体积和B的迹对定姿精度的影响量级，以（120°，30°）天区为例，在天区探测到的恒星共有30颗，共有4 060种定姿组合，其中由此可见，恒星的空间几何分布是决定定姿误差的主要因素。因此，可考虑恒星的空间几何分布因素作为导航星库选取策略，以提升星敏感器的定姿精度。

2.2 基于恒星分布的导航星库制作方法和步骤

综上所述，本文选取Hipparcos星库作为基础星库来筛选导航星库。

第1步：星库精简，根据星敏感器的敏感极限星等精简星库，删除原始星库中星等大于6.2的恒星，同时删除双星，光学头部的视场为14.5° × 14.5°，探测器尺寸为2 048 × 2 048，若星点成像为3 × 3的高斯分布，则能分辨两颗恒星的最小间距为6个像素，且星跟踪模式下的星跟踪半径最大设为30，则剔除星对角矩小于即删星对角矩小于0.212°的双星。图1为精简前后的星库三维图示。精简后的星库共有4 975颗星。

图1 精简前后的星库三维图Fig. 1 3-D image of star datebase before and after simplification

第2步：随机生成光轴指向，并搜索当前光轴指向所在天区视场内的星，任取3颗星计算其星矢量组成的四面体的体积，分别计算每种组合中恒星矢量组成四面体体积。

第3步：基于恒星几何分布选取导航星。由2.1所述可知，用于定姿的恒星矢量空间构成体积越大，定姿误差越小，定姿精度越高。以（120°，30°）天区为例，在基础星库中该天区共有30颗星，任取3颗导航星，则有4 060种组合。图3给出了姿态误差与四面体体积的变化趋势图，从图中可以看出，四面体体积越大，姿态角的误差越小，即姿态角度的精度越高。图4给出了精简后的星库在（120°，30°）天区的星点分布图，图中红色的3颗星为空间几何分布体积最小的3颗星，可以看出这3星在平面坐标上基本处于一条直线上，其星矢量形成的空间四面体的体积最小，解算出的姿态数据误差也是最大的，可先删除定姿误差较大的3颗星。

第4步：在减少恒星分布对定姿误差影响后，再次遍历全天球，每个视场保留最亮的15颗，以进一步精简星库，精简之后星库共有导航星4 065颗。

图2 姿态误差与体积变化关系Fig. 2 Relationship between attitude error and volume

图3 精简星库恒星分布图Fig. 3 The distribution map of simplified star datebase

图4 视场划分图Fig. 4 Split graph of FOV

第5步：为了保证导航星库在视场中的星分布尽可能均匀，本文将内切圆视场划分成4个面积相等的扇形区域，如图4所示。遍历每个视场，统计每个扇形区域的星颗数，若出现有星颗数为0的扇形区，则将基础星库在该扇形区最亮的星添加到导航星库中。仿真结果表明，共有2 260个天区出现星颗数为0的区域，共添加126颗星，即星库星颗数共为4 191颗。图5～7给出了扇区为0区域加添星的过程。图8为星库的全天球分布图。

图5 删除误差较大的定姿星Fig. 5 Delete large error star for attitude determination

图6 选取视场中最亮的15颗星Fig. 6 Select the most 15 brightest star in FOV

3 试验验证情况

为了验证星库的有效性和完备性，对制作的导航星库进行全天球识别验证，输入15组动力学参数进行全天识别测试全天球的覆盖性，即赤经0°到360°，赤纬分别为–90°到90°的轨道。识别情况如表2所示。为了进一步评价星库的完备性，随机生成1万个光轴指向，统计每个视场内的星颗数情况，如表3所示。

从表2和表3可以看出，在测试的全天球轨道中，本文方法制作的星库可以较好地满足全天球捕获的要求，且颗数大于10颗的天区占到97.64%，覆盖性较好。

为了进一步验证本文方法制作的星库对提升姿态精度的有效性，在全天球随机生成了10万个光轴指向，计算得到理想姿态与实际姿态的误差，如图9所示。

从图中可以看出，利用本文方法制作的星库，删除了定姿误差较大的恒星，提升了星敏的定姿精度，极大地减小了X、Y、Z3个方向的理想姿态与实际误差。与经典的导航星库制作方法相比，利用本文方法制作的导航星库计算出的姿态精度更高，姿态输出也更为稳定，且不会出现如图2所示姿态误差变化较大的情况。

表2 筛选后星库的全天识别情况Table 2 The Recognition of all-sky by navigation star datebase

表3 星库中星数的分布情况Table 3 The distribution of star number in star datebase

4 结 论

图9 理想姿态与实际姿态的误差Fig. 9 Error between ideal attitude and actual attitude

本文分析了恒星分布对定姿精度的影响，进而提出了一种基于恒星几何分布的星敏感器导航星库的制作方法。首先根据星敏感器的敏感极限以及跟踪半径的大小精简星库，计算每个视场内不同恒星分布组合的定姿误差，删除误差较大的星以进一步精简星库。同时对于视场内恒星分布的均匀性进行了优化，添加了星颗数为0的扇区，使得导航星库进一步完善，经全天球识别测试和天区星数统计结果表明，本文方法制作的星库完备性和均匀性较好，同时能够有效提升星敏感器的定姿精度。

参考文献

[1]朱长征，居永忠，杜晓辉. 导航星库制定方法研究[J]. 宇航学报，2010，31（5）：1327-1330.ZHU C Z，JU Y Z，DU X H. Approach to founding a guide star catalogue[J]. Journal of Astronautics，2010，31（5）：1327-1330.

[2]郑万波，刘智，郝志航. 二维精简索引分层导航星库的构造[J]. 吉林大学学报（信息科学版），2003，21（2）：123-127.ZHENG W B，LIU Z，HAO Z H. Construction of delaminating guide star database with 2-D condensed index[J]. Journal of Changchun Post and Telecommunication Institute，2003，21（2）：123-127.

[3]许士文，李葆华，张迎春，等. 用字符匹配进行星图识别的导航星库存储方法[J]. 哈尔滨工业大学学报[J]. 2005，37（6）：819-821.XU S W，LI B H，ZHANG Y C，et al. A store navigation star database method for star map recognition using character match[J]. Journal of Harbin Institute of Technology，2005，37（6）：819-821.

[4]蔡志武，韩春好，陈金平. 一种高精度星载导航星库的构建方法[J].测绘科学技术学报，2006，23（1）：29-32.CAI Z W，HAN C H，CHEN J P. A approach of construction high precision guide star catalogue on board[J]. Journal of Zhengzhou Institute of Surveying and Mapping，2006，23（1）：29-32.

[5]ZHANG H，SANG H. A double-FOV star sensor for high dynamic spacecraft navigation[C]//Proceedings of the SPIE，MIPPR 2009.Yichang，China：[s.n.]，2009.

[6]LIEBE C C. Star tracker for attitude determination[J]. IEEE Aerospace and Electronic Systems Magazine，1995，10（6）：10-16.

[7]SUN Y，XIAO Y，GENG Y. On-obrit calibration of star sensor based on a new lens distortion model[C]// The 32 nd Chinese Control Conference. Xi'an，China：IEEE，2013：4989-4994.

[8]张磊，魏仲慧，何昕. 快速全天自主星图识别[J]. 光学精密工程，2009，17（4）：909-915.ZHANG L. WEI Z H，HE X. Fast all-sky autonmous star identification[J]. Optics and Precision Engineering，2009，17（4）：909-915

[9]王丽娜，赵慧，熊智，等. 基于导航恒星几何分布的天文导航定姿系统误差建模及误差特性分析[J]. 兵工学报，2015，36（10）：1933-1942.WANG L N，ZHAO H，XIONG Z，et al. Modeling and error characteristics analysis of celestial attitude determination error based on geometric configuration of guide stars[J]. Acta Armamentarii，2015，36（10）：1933-1942.

[10]赵慧，熊智，王丽娜，等. 基于恒星几何构型分布的天文定位误差建模及误差特性分析[J]. 兵工学报，2015，36（5）：813-822.ZHAO H，XIONG Z，WANG L N，et al. Modeling of celestial positioning error and analysis of error characteristics based on distribution of guide stars[J]. Acta Armamentarii，2015，36（5）：813-822.