留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

L2范式距离的星图识别方法

王泽文 王广君 刘小波 佘锦华

王泽文, 王广君, 刘小波, 佘锦华. L2范式距离的星图识别方法[J]. 红外与激光工程, 2020, 49(10): 20200040. doi: 10.3788/IRLA20200040
引用本文: 王泽文, 王广君, 刘小波, 佘锦华. L2范式距离的星图识别方法[J]. 红外与激光工程, 2020, 49(10): 20200040. doi: 10.3788/IRLA20200040
Wang Zewen, Wang Guangjun, Liu Xiaobo, She Jinhua. Star map recognition method of L2 normal distance[J]. Infrared and Laser Engineering, 2020, 49(10): 20200040. doi: 10.3788/IRLA20200040
Citation: Wang Zewen, Wang Guangjun, Liu Xiaobo, She Jinhua. Star map recognition method of L2 normal distance[J]. Infrared and Laser Engineering, 2020, 49(10): 20200040. doi: 10.3788/IRLA20200040

L2范式距离的星图识别方法

doi: 10.3788/IRLA20200040
基金项目: 中国地质大学(武汉)中央高校基本科研业务费资助项目(G1323519204)
详细信息
    作者简介:

    王泽文(1994-),男,硕士生,主要从事控制理论与应用研究。Email:wangzewen@cug.edu.cn

    通讯作者: 王广君(1963-),男,教授,博士,主要从事天文导航、遥感图像处理、电子信息技术和智能仪器方面的教学和研究工作。Email:gjwang@cug.edu.cn

Star map recognition method of L2 normal distance

  • 摘要: 在Hausdorff距离的基础上,提出一种不依赖于星敏感器的旋转方向和焦距等因素的星图识别方法。在构造Hausdorff距离的数据点集合时,采用基于L2范数对应的相对欧氏距离作为集合元素,解决星敏感器滚动对星图识别的影响;另一方面,由于受星敏感器焦距的影响,星敏感器图像与标准参考图也会存在误差。在构造标准数据点元素时,考虑到如果一个数据点集包含另一个数据点集,在这两个数据点集之间至少有两个数据点之间的L2范式距离是相同的。对L2范式Hausdorff距离进行比例化处理,每个集合中的相对空间距离除以本集合中最小的相对空间距离,构成一种新的数据点集。这种方法不需要对星敏感器图像由于焦距不同进行标定,避免了星敏感器焦距对星图识别的影响。给出了距离的计算公式和实现步骤,并给出了实验结果。结果表明:在星敏感器转动、尺度变换等情况下,该算法可以正确得到星图识别结果,从而获得星敏感器的姿态信息。
  • 图  1  模拟星敏感器图像

    Figure  1.  Image of analog star sensor

    图  2  标准星库下基于L2范式距离的识别结果

    Figure  2.  Recognition results based on L2 normal form distance under standard star Library

    图  3  L1范式下的最小Hausdorff距离和

    Figure  3.  Sum of minimum Hausdorff distance based on L1 normal form

    图  4  L2范式下的最小Hausdorff距离和

    Figure  4.  Sum of minimum Hausdorff distance based on L2 normal form

    图  5  L1范式下滚动角分别为0.5°、1.0°、1.5°和2.0°时的最小HD距离

    Figure  5.  Minimum HD distance at rolling angle is 0.5°, 1.0°, 1.5°and 2.0° in L1 normal form respectively

    图  6  L1范式下可识别性随滚动角的变换

    Figure  6.  Transformation of identifiability with rolling angle under normal form L1

    图  7  转动角0.5°、1.0°、1.5°和2.0°时,L2范式下的最小HD距离

    Figure  7.  Minimum HD distance under normal form L2 when rotation angle is 0.5°, 1.0°, 1.5° and 2.0° respectively

    图  8  L2范式下HD最小距离与次最小距离之间的差

    Figure  8.  Difference between the HD minimum distance and the secondary minimum distance in L2 normal form

    图  9  星敏感器图像膨胀1、1.05、1.10、1.15倍时,L2范式下的最小HD距离

    Figure  9.  Minimum HD distance under normal form L2 when the image of star sensor expands 1,1.05, 1.10, 1.15 times respectively

    图  10  星敏感器图像膨胀1、2、3、4倍时,范式L2范式下的最小HD距离

    Figure  10.  Minimum HD distance under L2 normal form when the image of star sensor expands 1, 2, 3, 4 times respectively

  • [1] Wang G, Li J, Wei X. Star identification based on Hash Map [J]. IEEE Sensors Journal, 2018, 18(4): 1591-1599. doi:  10.1109/JSEN.2017.2772272
    [2] Pham M D, Low K S, Chen S S. An autonomous star recognition algorithm with optimized database [J]. IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(3): 1467-1475. doi:  10.1109/TAES.2013.6557999
    [3] Li J, Wei X, Zhang G J. Iterative algorithm for autonomous star identification [J]. IEEE Transactions on Aerospace and Electronic Systems, 2015, 51(1): 536-547. doi:  10.1109/TAES.2014.130729
    [4] 朱海龙, 梁斌, 张涛. 基于星点有序集的全天球自主星图识别算法[J]. 光学 精密工程, 2017, 25(6): 1577-1586. doi:  10.3788/OPE.20172506.1577

    Zhu Hailong, Liang Bin, Zhang Tao. All-sky automation star identification algorithm based on ordered star points set [J]. Optics and Precision Engineering, 2017, 25(6): 1577-1586. (in Chinese) doi:  10.3788/OPE.20172506.1577
    [5] Zhao Y, Wei X G, Li J, et al. Star identification algorithm based on k-l transformation and star walk formation [J]. IEEE Sensors Journal, 2016, 16(13): 5202-5210. doi:  10.1109/JSEN.2016.2553245
    [6] 王军, 何昕, 魏仲慧, 等. 基于多特征匹配的快速星图识别[J]. 光学 精密工程, 2019, 27(8): 1870-1879. doi:  10.3788/OPE.20192708.1870

    Wang Jun, He Xin, Wei Zhonghui, et al. Fast atar identification algorithm based on multi-feature matching [J]. Optics and Precision Engineering, 2019, 27(8): 1870-1879. (in Chinese) doi:  10.3788/OPE.20192708.1870
    [7] Zhao Y, Wei X G, Li J, et al. A fast triangular star map recognition algorithm [J]. Journal of Applied Optics, 2018, 39(1): 71-75. doi:  10.5768/JAO201839.0102006
    [8] 刘烟, 席红霞, 曹珺, 等. 结合自组织映射网络及三角形算法的星图识别方法[J]. 中国空间科学与技术, 2018, 38(4): 1-10.

    Liu Yan, Xi Hongxia, Cao Jun, et al. A star pattern recognition method baesd on self-organizing mapnetwork and triangle algorithm [J]. Chinese Space Science and Technology, 2018, 38(4): 1-10. (in Chinese)
    [9] Li X L, Yang J H, Zhang L, et al. Further optimization of selected observation triangles technique before star pattern matching [J]. Journal of Astronautics, 2015, 36(1): 76-81.
    [10] 踪华, 刘嬿, 高晓颖, 等. 星点坐标辅助的全天区三角形星图识别算法[J]. 航空学报, 2018, 39(9): 322045.

    Zong Hua, Liu Yan, Gao Xiaoying, et al. All-sky triangle algorithm of stars identification aided by star coordinates [J]. Acta Aeronautica et Astronautica Sinica, 2018, 39(9): 322045. (in Chinese)
    [11] 高扬, 赵金宇, 陈涛, 等. 添加补偿码的快速径向伴星特征星图识别[J]. 光学 精密工程, 2017, 25(6): 1627-1634. doi:  10.3788/OPE.20172506.1627

    Gao Yang, Zhao Jinyu, Chen Tao, et al. Radial neighbor feature with compensate code star pattern recognition algorithm [J]. Optics and Precision Engineering, 2017, 25(6): 1627-1634. (in Chinese) doi:  10.3788/OPE.20172506.1627
    [12] Li B H, Sun Q, Zhang T S. A star pattern recognition algorithm for the double-FOV star sensor [J]. IEEE Aerospace and Electronic Systems Magazine, 2015, 30(8): 24-31. doi:  10.1109/MAES.2015.140208
    [13] 苟斌, 程咏梅, 赵明艳, 等. 三视场星敏感器的多级星图识别算法[J]. 西北工业大学学报, 2019, 37(3): 541-546. doi:  10.1051/jnwpu/20193730541

    Gou Bin, Cheng Yongmei, Zhao Mingyan, et al. Multi-stage star image identification method of three field-of-view star sensor [J]. Journal of Northwestern Polytechnical University, 2019, 37(3): 541-546. (in Chinese) doi:  10.1051/jnwpu/20193730541
    [14] Sun L, Jiang J, Zhang G J, et al. A discrete HMM-based feature sequence model approach for star identification [J]. IEEE Sensors Journal, 2016, 16(4): 931-940. doi:  10.1109/JSEN.2015.2496662
    [15] Yoon H, Paek S W, Lim Y, et al. New star pattern identification with vector pattern matching for attitude determination [J]. IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(2): 1108-1118. doi:  10.1109/TAES.2013.6494402
    [16] 王广君, 房建成. 基于Hausdorff距离的星图识别方法[J]. 北京航空航天大学学报, 2005, 31(5): 508-510. doi:  10.3969/j.issn.1001-5965.2005.05.006

    Wang Guangjun, Fang Jiancheng. New star pattern recognition approach based on Hausdorff distance [J]. Journal of Beijing University of Aeronautics and Astronautics, 2005, 31(5): 508-510. (in Chinese) doi:  10.3969/j.issn.1001-5965.2005.05.006
    [17] Quan W, Xu L, Fang J C. A new star identification algorithm based on improved Hausdorff distance for star sensors [J]. IEEE Transactions on Aerospace and Electronic Systems, 2013, 49(3): 2101-2109. doi:  10.1109/TAES.2013.6558045
  • [1] 王彦, 程东升, 蒋超, 葛子阳, 金萍.  基于FBG柔性传感器的滑觉信号特性识别 . 红外与激光工程, 2023, 52(3): 20220587-1-20220587-11. doi: 10.3788/IRLA20220587
    [2] 徐晓冰, 左涛涛, 孙百顺, 李奇越, 吴刚.  基于热释电红外传感器的人体动作识别方法 . 红外与激光工程, 2022, 51(4): 20210188-1-20210188-8. doi: 10.3788/IRLA20210188
    [3] 冯婕, 李豫东, 文林, 郭旗.  CMOS图像传感器辐射损伤导致星敏感器性能退化机理 . 红外与激光工程, 2020, 49(5): 20190555-20190555-7. doi: 10.3788/IRLA20190555
    [4] 郑循江, 何峰, 姜丽辉, 曹卫卫, 肖东东.  火星探测环绕段自主光学导航算法研究(特约) . 红外与激光工程, 2020, 49(5): 20190458-20190458-9. doi: 10.3788/IRLA20190458
    [5] 苏本跃, 郑丹丹, 汤庆丰, 盛敏.  单传感器数据驱动的人体日常短时行为识别方法 . 红外与激光工程, 2019, 48(2): 226003-0226003(9). doi: 10.3788/IRLA201948.0226003
    [6] 陈津津, 苏君红, 金宁, 浦恩昌, 张皓, 苏俊波, 周立钢, 明景谦, 徐曼, 杨开宇, 宋治航.  基于中波制冷型碲镉汞探测器的远距离探测/识别连续变焦热像仪 . 红外与激光工程, 2018, 47(4): 404004-0404004(8). doi: 10.3788/IRLA201847.0404004
    [7] 张宇飞, 贺岩, 刘梦庚, 陈卫标.  基于伪随机码调制的测距通信一体化激光雷达 . 红外与激光工程, 2018, 47(9): 930003-0930003(6). doi: 10.3788/IRLA201847.0930003
    [8] 刘先一, 周召发, 张志利, 刘殿剑, 朱文勇.  基于数字天顶仪的视场角分析 . 红外与激光工程, 2016, 45(6): 617001-0617001(6). doi: 10.3788/IRLA201645.0617001
    [9] 李志军, 刘松林, 牛照东, 陈曾平.  基于梯度相位和显著性约束的Hausdorff 距离模板匹配方法 . 红外与激光工程, 2015, 44(2): 775-780.
    [10] 王昊京, 王建立, 吴量, 姚凯男.  三视场定位定向设备的视场确定 . 红外与激光工程, 2015, 44(6): 1890-1897.
    [11] 孙剑明.  剖分特征集星识别法在天文导航中的应用 . 红外与激光工程, 2015, 44(11): 3330-3335.
    [12] 褚永辉, 李茂登, 黄翔宇, 王大轶.  基于陆标和脉冲星测量的自主导航方法 . 红外与激光工程, 2015, 44(S1): 143-146.
    [13] 廖育富, 钟建勇.  一种新的方法用于星敏感器近红外星图星点提取 . 红外与激光工程, 2014, 43(5): 1667-1671.
    [14] 张磊, 何昕, 魏仲慧, 梁国龙.  星图分布对星敏感器最小二乘姿态精度的影响 . 红外与激光工程, 2014, 43(6): 1836-1841.
    [15] 魏新国, 李延鹏, 李健, 江洁.  多视场星敏感器近地轨道自主定位导航方法 . 红外与激光工程, 2014, 43(6): 1812-1817.
    [16] 梁冬生, 刘朝晖, 刘文, 袁辉, 刘夫成.  航空飞行器天文自主导航定位技术 . 红外与激光工程, 2014, 43(9): 3020-3025.
    [17] 毛海岑, 刘爱东, 王亮.  采用混合粒子群算法的星图识别方法 . 红外与激光工程, 2014, 43(11): 3762-3766.
    [18] 蔡义, 汪红熳, 亓波.  星地激光通信无波前传感器优化算法仿真分析 . 红外与激光工程, 2013, 42(4): 1063-1068.
    [19] 陈聪, 王宏力, 崔祥祥, 陆敬辉, 姜伟.  基于预测跟踪星表提高星敏感器实时性的方法 . 红外与激光工程, 2013, 42(8): 2190-2196.
    [20] 李葆华, 陈希军, 郑循江, 泮宏梁.  星敏感器高动态下自主星跟踪算法 . 红外与激光工程, 2012, 41(1): 190-195.
  • 加载中
图(10)
计量
  • 文章访问数:  381
  • HTML全文浏览量:  133
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-01-12
  • 修回日期:  2020-02-04
  • 网络出版日期:  2020-11-03
  • 刊出日期:  2020-11-03

L2范式距离的星图识别方法

doi: 10.3788/IRLA20200040
    作者简介:

    王泽文(1994-),男,硕士生,主要从事控制理论与应用研究。Email:wangzewen@cug.edu.cn

    通讯作者: 王广君(1963-),男,教授,博士,主要从事天文导航、遥感图像处理、电子信息技术和智能仪器方面的教学和研究工作。Email:gjwang@cug.edu.cn
基金项目:  中国地质大学(武汉)中央高校基本科研业务费资助项目(G1323519204)

摘要: 在Hausdorff距离的基础上,提出一种不依赖于星敏感器的旋转方向和焦距等因素的星图识别方法。在构造Hausdorff距离的数据点集合时,采用基于L2范数对应的相对欧氏距离作为集合元素,解决星敏感器滚动对星图识别的影响;另一方面,由于受星敏感器焦距的影响,星敏感器图像与标准参考图也会存在误差。在构造标准数据点元素时,考虑到如果一个数据点集包含另一个数据点集,在这两个数据点集之间至少有两个数据点之间的L2范式距离是相同的。对L2范式Hausdorff距离进行比例化处理,每个集合中的相对空间距离除以本集合中最小的相对空间距离,构成一种新的数据点集。这种方法不需要对星敏感器图像由于焦距不同进行标定,避免了星敏感器焦距对星图识别的影响。给出了距离的计算公式和实现步骤,并给出了实验结果。结果表明:在星敏感器转动、尺度变换等情况下,该算法可以正确得到星图识别结果,从而获得星敏感器的姿态信息。

English Abstract

    • 天文导航由于不存在累计误差常被用于飞行器的自主导航或组合导航中。飞行器自主天文导航的重要技术就是星图的快速识别方法,通过实时星图识别方法获取飞行器的姿态。星图识别算法很多,有基于Hash Map的星图识别算法[1]、自主星图识别算法[2-4]、基于KL变换算法[5]、多特征匹配算法[6]、三角形匹配算法[7-10]以及其他算法[11-15]。目前研究最多的是三角形匹配算法。但是,三角形匹配算法受星图提取位置精度的影响容易造成误识别。特别是在大视场情况下计算速度慢,实时性差。近年来基于Hausdorff距离的星图识别算法由于其抗干扰能力强被广泛关注[16-17]。 在对基于Hausdorff距离的星图识别方法研究中,主要是由星敏感器获得图像和参考星库中星图中恒星的空间位置构成两个不同的数据点集,然后通过计算两个点集之间的Hausdorff距离进行识别。在理想情况下,当一个数据点集包含另一个数据点集时,它们之间的Hausdorff距离为零。但是,由于星敏感器的成像不仅与方向有关,而且与星敏感器的滚动角和光学焦距都有关系。当星敏感器发生转动时,恒星在星敏感器图像中的位置会发生变化;当星敏感器的焦距不同时,图像的大小也不同。参考文献[16-17]是利用L1范数距离进行局部Hausdorff识别,但是,L1范数对应的距离是有方向的矢量距离,尽管识别准确性高,但是,需要有先验知识,适合于与惯导结合的组合导航中。如果没有先验知识,利用L1距离范数进行的HD距离计算方法需要考虑星敏感器旋转问题,其计算速度大大降低。为了使HD距离具有更好的鲁棒性,文中介绍了一种利用L2范数的HD距离星图识别方法,这种方法不依赖于星敏感器的旋转和星敏感器焦距,可满足任意状态下自主导航对星敏感器的要求。

    • Hausdorff距离是两个点集之间相似度的一种描述方法,假设两个点集分别为A={a1a2,…aN}和B={b1b2,…, bM},AB之间的双向Hausdorff距离描述为:

      $$H(A,B) = \max \left\{ {h(A,B),h(B,A)} \right\}$$ (1)

      式中:h(A,B) 称为AB的单向Hausdorff距离,h(B,A)称作BA的单向Hausdorff距离,它们分别表示为:

      $$h(A,B) = \mathop {\max }\limits_{a \in A} \mathop {\min }\limits_{b \in B} \left\| {a - b} \right\|$$ (2)
      $$h(B,A) = \mathop {\max }\limits_{b \in B} \mathop {\min }\limits_{a \in A} \left\| {b - a} \right\|$$ (3)

      这里$\left\| { \cdot } \right\|$是一种距离范式,参考文献[16-17]中采用L1范式距离作为度量。L1范式是指向量中各个元素绝对值之和,由于向量对方向非常敏感,其应用受到限制。当采用L2范式时,对应的是欧氏距离。欧氏距离没有方向性,对星敏感器旋转不敏感。从公式(1)可以看出,双向$H(A,B)$距离取单向$h(A,B)$距离和$h(B,A)$距离的最大值,这个值越大表示AB的不相似性越大。反过来,这个值越小,表示AB的相似性越大。参考文献[16-17]中数据点集是利用图像的空间直角坐标进行构造,即:

      $$A = \left\{ {\left( {{x_i} - {x_1},{y_i} - {y_1}} \right){\rm{|}}b{r_{ai}}} \right\},\;(i = 1,2,...,N)$$ (4)

      即点集A是由第一颗星为原点的相对星图坐标构建,其中$b{r_{ai}}$为图像灰度,它与星库中的星等相对应。在实际应用中,星图的空间坐标依赖于星敏感器的转动方向。在惯性与天文组合导航或已知先验知识的情况下,这种方法具有很强的星图识别能力。但是,在纯自主天文导航中,飞行器状态一旦丢失,它需要很长时间对所有星敏感器滚动角进行计算才能重新获得星敏感器的状态,自主恢复能力较差。如果采用范数L2对应的欧氏距离构造数据点集A,则可以消除滚动角对星图识别的影响。由此,公式(4)改为:

      $$A=\left\{ {{d}_{{{a}_{i}}}}|b{{r}_{i}} \right\}(i=1,2,\cdot \cdot \cdot ,N)$$ (5)

      式中:${d_{{a_{i1}}}} = \sqrt {{{\left( {{x_i} - {x_1}} \right)}^2} + {{\left( {{y_i} - {y_1}} \right)}^2}} $表示星图中任意一恒星点到第一颗恒星点之间的欧氏距离。

    • 在星图识别中,参考星图是由标准星库构成的。标准星库中的恒星一般是由经纬度坐标和星等构成。为了构造与星敏感器星图一样的点集,需要把星库中天球坐标系下的经纬度坐标转换成星敏感器的平面坐标。由于恒星到地球的距离可以认为是无穷远处,对恒星方位角来说,从地心赤道坐标系到星敏感器坐标系的转换无需顾及坐标平移所带来的误差。设原点O指向地心,Ox轴指向春分点,Oz轴指向北极,Oy为右手准则拇指方向,从天球坐标系到平面坐标系转换的基本公式为[16]

      $$\left\{ \begin{aligned} x =\; & \frac{{{N_x}}}{2} \times \frac{1}{{\tan \frac{{{{FO}}{{{V}}_x}}}{{\rm{2}}}}} \times \\ &\frac{{\cos {\delta _i}\sin \left( {{\alpha _i} - {\alpha _0}} \right)}}{{\sin {\delta _i}\sin {\delta _0} + \cos {\delta _i}\cos {\delta _0}{\rm{cos}}({\alpha _{\rm{i}}} - {\alpha _0})}} \\ y =\; & \frac{{{N_y}}}{2} \times \frac{1}{{\tan \frac{{{{FO}}{{{V}}_y}}}{{\rm{2}}}}} \times \\ &\frac{{\sin {\delta _i}\cos {\delta _0} - \cos {\delta _i}\sin {\delta _0}\cos \left( {{\alpha _i} - {\alpha _0}} \right)}}{{\sin {\delta _i}\sin {\delta _0} + \cos {\delta _i}\cos {\delta _0}\cos ({\alpha _{\rm{i}}} - {\alpha _0})}} \end{aligned} \right.$$ (6)

      其中

      $$x \in \left( - \frac{{Nx}}{2},\frac{{Nx}}{2}\right),y \in \left( - \frac{{Ny}}{2},\frac{{Ny}}{2}\right)$$ (7)

      式中:(${{N}_{x}},{{N}_{y}}$)分别是在xy方向上星库中恒星的经纬度坐标转换成星敏感器平面坐标系后的最大像素数,它与星敏感器的焦距有关;(${{FO}}{{{V}}_{x}},{{FO}}{{{V}}_{y}}$)分别是在xy方向上的视场角;($\delta ,\alpha $)为恒星的经纬度坐标;(${{\delta }_{0}},{{\alpha }_{0}}$)为作为参考方向的经纬度坐标。用上面同样的方法可以构造相对于参考星点K的L2范式距离点集B

      $${B_k} = \left\{ {d_{bjk}^{}|{br}'_j} \right\},\left( {j = 1,2, \cdot \cdot \cdot ,M} \right),\left( {1 \leqslant k \leqslant M} \right)$$ (8)

      式中:$d_{bjk}^{} = \sqrt {{{\left( {x_{bj}^{} - x_{bk}^{}} \right)}^2} + {{\left( {y_{bj}^{} - y_{bk}^{}} \right)}^2}} $表示任意一点到参考点K之间的欧氏距离;$br_{bj}^{}$为与星等对应的图像灰度,需要与传感器图像进行定标,文中在星图识别中忽略了星等的影响。参考点K不同,点集B也不同,遍及星库中所有恒星作为参考点,可得得到M个不同的点集B。在实际计算中,为了提高计算速度,根据星敏感器视场角大小,选取参考星周围锥形范围内的星即可。

    • 上面算法尽管对星敏感器的滚动角没有要求,但是,从公式(6)中可以看出,构造参考星图时与参数(${{N}_{x}},{{N}_{y}}$)有关,也即与星敏感器的焦距有关。在进行图像匹配时需要对不同的星敏感器进行单独标定,既麻烦又带来一定的标定误差。为此,文中在上述基础上采用相对距离方法构建数据点集AB

      假设在星敏感器中恒星距离构成的点集A中,距离$\left\{ {{d}_{{{a}_{i}}}} \right\},(i=1,2,\cdot \cdot \cdot ,N)$是由小到大排序的,星敏感器A的点集可以改写为:

      $$\begin{split} & A = \left\{ {\frac{{{d_{{a_i}}}}}{{{d_{a2}}}}|b{r_{ai}}} \right\} = \left\{ {\left( {0,1,\frac{{{d_{a3}}}}{{{d_{a2}}}}, \cdot \cdot \cdot ,\frac{{{d_{aN}}}}{{{d_{a2}}}}} \right){\rm{|}}b{r_{ai}}} \right\},\\ & \quad (i = 1,2, \cdot \cdot \cdot ,N) \end{split}$$ (9)

      因为集合中A中的元素是所有星敏感器中对应所有星的位置相对于第一颗星之间的欧氏距离,所以,第一个值为零。

      同样,假设星敏感器的图像不变,参考星图在星敏感器坐标系中表示为:

      $$\left\{ \begin{aligned} x =\; & a \times \frac{{{N_x}}}{2} \times \frac{1}{{\tan \frac{{{{FO}}{{{V}}_x}}}{{\rm{2}}}}} \times \\ & \frac{{\cos {\delta _i}\sin \left( {{\alpha _i} - {\alpha _0}} \right)}}{{\sin {\delta _i}\sin {\delta _0} + \cos {\delta _i}\cos {\delta _0}{\rm{cos}}({\alpha _{\rm{i}}} - {\alpha _0})}} = \\ \; & a \times {x_0} \\ y =\; & a \times \frac{{{N_y}}}{2} \times \frac{1}{{\tan \frac{{{{FO}}{{{V}}_y}}}{{\rm{2}}}}} \times \\ &\frac{{\sin {\delta _i}\cos {\delta _0} - \cos {\delta _i}\sin {\delta _0}\cos \left( {{\alpha _i} - {\alpha _0}} \right)}}{{\sin {\delta _i}\sin {\delta _0} + \cos {\delta _i}\cos {\delta _0}{\rm{cos}}({\alpha _{\rm{i}}} - {\alpha _0})}} =\\ \; & a \times {y_0} \end{aligned} \right.$$ (10)

      则,点集合B中点数据元素可以表示为:

      $${d_{jk}} = a \times {d_0}_{jk}$$ (11)

      同样,对标准点集B中相对k的L2范式距离排序,并改写B为:

      $$\begin{split} {B_k}(j) = \; &\left\{ {\frac{{{d_{{b_{jk}}}}}}{{{d_{b{}_{{\rm{2}}k}}}}}|b{r_j}} \right\} = \left\{ {\left( {0,\frac{{{d_{{b_{2k}}}}}}{{{d_{{b_{2k}}}}}},\frac{{{d_{{b_{3k}}}}}}{{{d_{{b_{2k}}}}}}, \cdot \cdot \cdot ,\frac{{{d_{bMk}}}}{{{d_{{b_{2k}}}}}}} \right){\rm{|}}b{r_j}} \right\}, \\ & \left( {{{j}} = 1,2, \cdot \cdot \cdot ,M} \right) \end{split} $$ (12)

      从公式(10)可以看出,相对距离${B_k}\left( i \right)$中的元素不依赖于a${B_k}\left( i \right)$表示以参考星k为方向,所有相对于k的欧氏距离排序,然后除以序列i的欧氏距离。

    • 对于范数L2点集A中任意一数据点到点集B之间的最小Hausdorff距离定义为:

      $${D_{\min (ik)}}\left( {A,B} \right) = \min \left( {\left\{ {{w_1}{\rm{|}}{d_{ai1}} - d_{bjk}^{}{{\rm{|}}^m} + w_2^{}{\rm{|}}b{r_{ai}} - br_{bj}^{}{|^m}} \right\}} \right)$$ (13)

      式中:W1W2为权值;m为一个距离范数,当m=2时对应的距离就是欧氏距离;${D_{\min (ik)}}$表示集合A中的任意点到参考点为kB集合的最小Hausdorff距离集合。如果集合B包含集合A,理论上,A中的任何点到B中的最小距离都为零。星敏感器不同指向对应的星图中的恒星个数不同,取最小Hausdorff距离的平均和作为评判标准,表示为:

      $${{H}}{{{D}}_k} = \frac{1}{N}\sum\limits_{i = 1}^N {{D_{\min (ik)}}} ,\;\;(k = 1,2, \cdot \cdot \cdot ,M)$$ (14)

      理论上当星图匹配时,${{H}}{{{D}}_k}$应为零。考虑到噪声影响,取${{H}}{{{D}}_k}$式中最小值对应的k即为识别结果。

      $$H{D_{\min }} = \min \left\{ {H{D_k}} \right\}$$ (15)
    • 根据上述方法,在没有任何先验知识的情况下,根据上面方法进行了实验。实验过程中,没有考虑星等的影响,相当于w2为零。仿真实验中,星敏感器视场角为6°×6°。为了提高计算速度,选取大于二倍以上星敏感器视场区域作为参考区域进行计算。识别过程中,随机选取星敏感器的方向,并对识别结果进行统计分析,其具体计算步骤如下:

      (1)在星库中,任意选取一颗星作为星敏感器的方向,并选取该星6°×6°区域内的所有恒星;

      (2) 利用公式(6)产生星敏感器平面图像$\left\{ {{x_{ai}},{y_{ai}}} \right\}$,然后把该图像旋转任意角度,并对每个像素增加0.5像素的随机噪声,生成仿真用的星星敏感器图像;

      (3) 利用公式(9)生成星敏感器星图相对距离集A

      (4) 在星库中,任意选取第k颗星作为标准星图的参考方向,并选取该星30°×30°区域内的所有恒星作为参考星;

      (5)利用公式(12)产生参考星相对距离集B

      (6)利用公式(13)计算相对于k的最小Hausdorff距离$H{D_k}\left( i \right)$

      (7)重复步骤(4)~(6),直到计算完所有参考星为止。

      (8)利用公式(14)求最小距离和;

      (9)利用公式(15),求Hausdorff最小距离和的最小值,此最小值对应的参考星即为星敏感器对应的恒星。

      理想情况下,星敏感器图像信息与星库中某一位置信息完全一致时,Hausdorff最小距离为零,所有星的Hausdorff最小距离和也为零。图1为在星库中,经纬度坐标为(30,5)、视场大小为6°×6°时利用公式(6)转换得到的星敏感器坐标系下的平面图像,此图像模拟作为星敏感器图像。图2为利用L2范式距离计算得到的在所有标准星库中获得的最小Hausdorff距离和,纵坐标为最小Hausdorff距离和,横坐标为恒星序列数。标准星库中星等从−1.09~6.95之间的恒星个数为14 581颗。从图中可以看出:当星图匹配时,最小HD距离和为零,其他位置的HD距离和非常大,从计算结果中很容易获得识别结果。但是,由于标准星库中恒星个数太多,计算速度非常慢。为了提高计算速度,验证该算法的有效性,在标准星库中任意选取了30°×30°的局部区域所为标准参考区域进行实验。为了对L1范式距离和L2范式距离两种方法的识别结果进行比较,在相同条件下分别对两种方法进行了仿真实验。图3图4分别为利用L1范式距离和L2范式距离在局部星库的识别结果。从图中可以看出,如果没有旋转,无论是利用L1范式距离方法或是利用L2范式距离方法,其对应的最小Hausdorff距离和几乎都为零,而其他位置的最小Hausdorff距离的和非常大,并且具有随机性,说明在没有旋转情况下,两种方法都可以获得正确的识别结果。图5为在旋转角分别为0.5°、1.0°、1.5°、2.0°情况下利用L1范式距离对应的最小Hausdorff距离和。最小Hausdorff距离和随着旋转角的增加而增加,这说明随着旋转角的增加,可识别性降低。图6为不同旋转角下利用L1范式距离计算的正确识别位置Hausdorff距离和与次最小Hausdorff距离和之间的差值,此差值越大表明可识别性越好(纵坐标为正确识别位置的最小HD减次最小HD距离)。从图中可以看出:当旋转角度超过1.5°时,采用L1范式距离计算的最小Hausdorff距离识别方法的可识别性大大下降,当旋转角超过2°时几乎无法识别。说明利用L1范式距离进行星图识别时,星敏感器不能随意转动。一旦发生较大转动,通过L1范式距离的Hausdorff识别方法将无法获得正确的识别结果。图7为在L2范式距离下星敏感器转动对最小Hausdorff距离识别结果的影响。在旋转角分别为0.5°、1.0°、1.5°和2.0°情况下,最小Hausdorff距离和的位置没有任何改变。图8为在0~180°旋转情况下,利用L2范式距离计算的最小Hausdorff距离和与次最小HD距离和之间的差值,这个差值几乎是一条水平线。这个差值说明,利用L2范式进行的最小Hausdorff距离识别方法不受旋转角的影响。沿星敏感器轴向无论如何转动,都可以获得正确的识别结果。所以,在利用L2范式距离对星图识别时,不用考虑前一时刻的状态。星敏感器焦距不同获得的星敏感器图像大小也不同。图910分别显示了星敏感器图像大小不一致时两种范式距离的识别结果。图9是星敏感器图像大小分别为原始图像的1、1.05、1.1、1.15倍时利用L1范式距离得到的相对最小Hausdorff距离,从图中可以看出,当图像大小发生变换时,图像的匹配特征快速消失,可识别性变差。图10是星敏感器图像大小分别为原始图像的1、2、3、4倍时利用L2范式距离得到的相对最小Hausdorff距离,从图中可以看出,无论图像大小如何变换,图像的匹配特征基本不变。这说明L2范式相对距离方法不需要对星敏感器的焦距进行校正。

      图  1  模拟星敏感器图像

      Figure 1.  Image of analog star sensor

      图  2  标准星库下基于L2范式距离的识别结果

      Figure 2.  Recognition results based on L2 normal form distance under standard star Library

      图  3  L1范式下的最小Hausdorff距离和

      Figure 3.  Sum of minimum Hausdorff distance based on L1 normal form

      图  4  L2范式下的最小Hausdorff距离和

      Figure 4.  Sum of minimum Hausdorff distance based on L2 normal form

      图  5  L1范式下滚动角分别为0.5°、1.0°、1.5°和2.0°时的最小HD距离

      Figure 5.  Minimum HD distance at rolling angle is 0.5°, 1.0°, 1.5°and 2.0° in L1 normal form respectively

      图  6  L1范式下可识别性随滚动角的变换

      Figure 6.  Transformation of identifiability with rolling angle under normal form L1

      图  7  转动角0.5°、1.0°、1.5°和2.0°时,L2范式下的最小HD距离

      Figure 7.  Minimum HD distance under normal form L2 when rotation angle is 0.5°, 1.0°, 1.5° and 2.0° respectively

      图  8  L2范式下HD最小距离与次最小距离之间的差

      Figure 8.  Difference between the HD minimum distance and the secondary minimum distance in L2 normal form

      图  9  星敏感器图像膨胀1、1.05、1.10、1.15倍时,L2范式下的最小HD距离

      Figure 9.  Minimum HD distance under normal form L2 when the image of star sensor expands 1,1.05, 1.10, 1.15 times respectively

      图  10  星敏感器图像膨胀1、2、3、4倍时,范式L2范式下的最小HD距离

      Figure 10.  Minimum HD distance under L2 normal form when the image of star sensor expands 1, 2, 3, 4 times respectively

    • 文中研究的基于L2范式Hausdorff星图识别方法与基于L1范式的Hausdorff距离识别方法不完全相同,基于L1范式的Hausdorff距离识别方法对星敏感器的初始状态有较高要求,适合于与惯导的组合导航中。文中研究的基于L2范式的Hausdorff距离识别方法既可以不考虑星敏感器光学焦距的影响,也不用考虑星敏感器滚动带来的影响,具有完全自主的星图识别能力。

参考文献 (17)

目录

    /

    返回文章
    返回