中国科技核心期刊      中国指挥与控制学会会刊     军事装备类重点期刊
多源信息融合

星敏感器星图数据质量的研究与应用

  • 沈宇腾 ,
  • 张弛
展开
  • 上海航天控制技术研究所, 上海 201108

沈宇腾(1990—),男,工程师,研究方向为星敏感测试。

张 弛(1998—),男,助理工程师。

Copy editor: 胡前进

收稿日期: 2024-11-03

  修回日期: 2024-12-22

  网络出版日期: 2025-07-28

Research and apply of the quality of star image for star tracker

  • SHEN Yuteng ,
  • ZHANG Chi
Expand
  • Shanghai Institute of Spaceflight Control Technology, Shanghai 201108, China

Received date: 2024-11-03

  Revised date: 2024-12-22

  Online published: 2025-07-28

摘要

近年来,随着商业航天技术的快速发展,对商业卫星用星敏感器的生产效能和指标一致性等研产质量的要求越来越高。利用了星敏感器生产中积累的海量星图数据,通过数据建模等理论研究,开发了一套量化评价星敏感器光机装调中星点图像的“7参数图像量化评价标准”。利用这一标准,在皮型星敏感器的光机装调中进行应用,并用外场观星测试验证了产品的最终测量精度,最终试验结果表明该方法能有效指导星敏感器光机组件装调,并为提高后续产品研产的测量精度一致性夯实基础。

本文引用格式

沈宇腾 , 张弛 . 星敏感器星图数据质量的研究与应用[J]. 指挥控制与仿真, 2025 , 47(4) : 80 -86 . DOI: 10.3969/j.issn.1673-3819.2025.04.012

Abstract

In recent years,with the development of commercial spaceflight, the requirements of the production quality for star tracker in commercial satellite become higher and higher. In this paper,“Seven Parameter Image Quantitative Evaluation Standard”is developed, which makes use of the massive star image data accumulated in the production of star tracker and through data modeling theory research. The adjustment of optical machine for Pi star tracker use this standard, and the final measurement accuracy of the product for Pi star tracker was verified by stargazing test, which prove the standard is effectively and can improve the consistency of measurement accuracy of the product for Pi star tracker.

星敏感器是航天器实现姿态感知、自主控制的重要仪器,具有精度高、可靠性高等特点,星敏感器的核心构成包含用于星点成像的光学系统、用于定点提取的图像传感器电路、星图识别与姿态解算数据处理电路[1-3]。传统的总装质量检测包括光机组件装调标定、整机合盖前后的光星模测试等,其中光机组件装调标定是星敏感器研产过程中装配操作最复杂、装配精度最高、安全风险最大的操作环节,同时也是不同工序设备使用冲突最严重的环节。
在当前商业航天加速发展的背景下,星敏感器在横向市场的需求量突增,为了达成高质高效产出,批量化服务商业卫星应用[4-5],本文开发一套基于星点图像数据的“7参数图像量化评价标准”,通过对星点图像数据应用分析预估光机组件标定精度情况,为提高大批量研产星敏感器测量精度一致性夯实基础,继而为商业卫星发展提供强有力的支撑[6]

1 星图数据的研究

在“天体-像照-电信号-位姿信息”的转换中,天体成像是核心环节,因此星图测试是评估星敏感器产品质量性能的关键手段之一,在总装生产、外场观星等过程中多次依靠星图数据来判别星敏感器产品的质量性能[7]。利用星图测试方法,能诊断恒星辐射能量波动、光学系统像差、光学系统结构热稳定性、视场杂散光干扰、星图视场能量分布极度不均、大视场成像角分辨率低等故障问题;外场观星星图数据基本“反映”产品的使用性能,如图1所示,暗室的模拟观星星图数据是对星敏感器产品性能的全面检验。星图数据与总装生产过程的工艺参数密切相关,其隐含着丰富的产品质量信息和生产质量信息。
图1 星图数据示意图

Fig.1 The schematic diagram of star image data

尽管星敏感器产品装调中的星图数据蕴含着丰富的产品质量信息和生产质量信息,但之前缺乏对星图数据的系统研究,没有形成量化的评价产品光机组件装调生产质量的方法,只能依赖操作人员的装调经验对星图图像的对称性、清晰度进行判读,其人为评价误差导致光机组件装调反复操作。而星敏感器产品测量精度一致性只能依靠后续外场观星测试环节确定结果,在外场观星阶段发现精度不达标往往会导致耗费大量人力物力进行光机组件的重复装调与标定,其评判结果一致性较差。因此,急需对星图数据进行研究,得出评价星图质量的量化指标。

1.1 星点图像建模

星点图像是恒星辐射信息经光学系统传递,在像面输出灰度信息的能量转换过程[8]。因此星点图上的灰度信息既可以表征星点分布特征,又可以反映星点能量特征。考虑星点图的实际特征,构建的特征模型如下:
(1)像斑大小(pxs_size):像斑大小即为星点图中每一个星点所占的像元数目。
(2)平均灰度(average_gray):平均灰度表示星点的灰度平均值,即星点占据像元的灰度值之和与星点占据的像元数的比值。
(3)灰度集中度(gray_con):灰度集中度表示星点内某一确定大小的同心像元区域的灰度和占星点灰度总和的百分比。
(4)信噪比(snr):信噪比为星图中星点目标的信号峰值与均方根误差的比值。
基于灰度特征分析的星图数字特征建模过程如图2所示,利用像斑大小、平均灰度值、灰度集中度、信噪比这4个参数描述一个像图,再对图1的产品星图测试结果进行产品质量特征建模,中心轴对称型星敏感器产品的质量特征参数主要包括:
(1)中心点灰度极值;
(2)中心点灰度集中值;
(3)周对称系数场。
图2 星点图像灰度特征提取

Fig.2 The gray feature extraction of star image

1.2 星点特征数据观测

作者分别以三个指标为轴绘制散点图,如图3所示,能够看出各产品不同视场的星点图指标数据具有一定的聚集性,相比标定数据也有更好的分辨力。
图3 散点图

Fig.3 The scatter diagram

1.3 视场中心星点图灰度峰值与理论中心点的偏差建模

图4所示,根据视场中心星点图的灰度峰值与理论中心点的偏差建模,其中偏差值计算公式为:a= ( x - x 0 ) 2 + ( y - y 0 ) 2图像偏差值越小,则星图质量越好。
图4 视场中心星点图

Fig.4 The center of field in star image

1.4 左、右星图与上、下星图对称性建模

图5所示,基于左、右星图与上、下星图对称性建模,因为4项指标特征表征星点分布特征和能量特征,与具体的位置无关,则对称性比较可以转换为相似性比较。
图5 星图对称性

Fig.5 The symmetry of star image

左、上指标向量:a1=(像斑大小,平均灰度,灰度集中度,信噪比)=(x1,y1,z1,u1)
右、下指标向量:a2=(像斑大小,平均灰度,灰度集中度,信噪比)=(x2,y2,z2,u2)
左、右(上、下)相似性(向量内积):θ=arccos a 1 · a 2 a 1 * a 2
星图的向量夹角越趋于零,则越相似,具体产品应用结果如图6图7图8所示。
图6 偏差值-左右相似性-上下相似性散点图

Fig.6 The scatter diagram of deviation value, left-right similarity and up-down similarity

图7 向量内积相似性计算的指标值

Fig.7 The index value of inner product similarity in vectors

图8 结构相似性计算的指标值

Fig.8 The index value of structural similarity

2 星图数据的应用

2.1 图像量化评价标准确立

根据星点图像建模与产品装调质量参数的大数据对比,以及综合考虑星点能量特征的量化标准[9-10],本文提出全视场点目标像斑图像的灰度和的均方差graystd、全视场点目标像斑图像灰度值的能量一致性graybias、全视场点目标像斑图像集中度grayconcent、灰度均值graymean、最大灰度值graymax、最小灰度值graymin和对称性graybalance共计7个指标参数作为评价星敏感器光机组件调试质量的判据。7个指标参数的具体含义、计算方法如下:
(1)全视场点目标像斑图像的灰度和的均方差graystd,表示图像M1~M4M6~M9中星点能量的方差,计算方法为
$\ { graystd }=\frac{\sum_{i=1}^{i=4}\left(g_{i}-\ { graymean }\right)^{2}+\sum_{i=6}^{i=9}\left(g_{i}-\ { graymean }\right)^{2}}{8}$
(2)全视场点目标像斑图像灰度值的能量一致性graybias,表示图像M1~M9中星点能量分别与灰度均值graymean的偏差,图像Mi中星点的能量一致性graybiasi计算方法为
$\ { graybias }_{i}=\frac{g_{i}-\ { graymean }}{\ { graymean }}$
(3)全视场点目标像斑图像集中度grayconcent,表示星点能量的聚集程度。理论的星点大小应为3*3,因此以提取到的质心Si为中心,开一个3*3大小的窗口,计算3*3窗口内的灰度和,为表征图像Mi中理论星点的能量。能量集中度是分别计算图像M1~M9中理论星点的能量占实际星点能量的比值,该值越接近1,说明星点能量越聚集。图像Mi中星点的能量集中度grayconcent计算方法为
$\text { grayconcent }_{i}=\frac{g_{-} n e w_{i}}{g_{i}}$
(4)灰度均值graymean表示图像M1~M4M6~M9中星点能量的灰度均值,计算方法为
$\ { graymean }=\frac{\sum_{i=1}^{i=9} g_{i}-g_{5}}{8}$
(5)最大灰度值graymax表示图像M1~M4M6~M9中星点能量的最大值。
(6)最小灰度值graymin表示图像M1~M4M6~M9中星点能量的最小值。
(7)对称性graybalance表示九宫图中,左-右,上-下,左上-右下,右上-左下这4组图像中星点能量的对称性,计算方法为每组图像中2个星点能量的比值。
本文研究基于以往大批量星敏感器装调数据分析和星图数据建模研究的基础上,提出了“7参数图像量化评价标准”的指标约束条件,其具体的约束条件如表1所示,若星点图像同时满足下述指标约束条件,则认为该产品光机组件标定结果满足产品最终测量精度要求。
表1 “7参数图像量化评价标准”的指标约束条件

Tab.1 The index constraint of “Seven Parameter Image Quantitative Evaluation Standard”

评价指标 约束条件
灰度均值graymean ≥400
方差graystd ≤100
最大灰度值graymax ≤750
最小灰度值graymin ≥350
能量一致性graybias -0.25~+0.25
能量集中度grayconcent ≥0.7
对称性graybalance 0.75~1.25

2.2 图像量化评价标准在皮型星敏感器产品的应用

本文选择了商业卫星型谱产品-皮型星敏感器(以下简称皮星敏)的产品1和产品2使用“7参数图像量化评价标准”进行光机组件调试和标定,产品的量化指标如表2所示。
表2 皮星敏产品“7参数图像量化评价标准”指标

Tab.2 The index constraint of “Seven Parameter Image Quantitative Evaluation Standard”

本文对皮星敏产品1和产品2进行外场观星测试,外场观星测试如图9所示,测试所得产品精度如图10图11所示。
图9 商业卫星用皮星敏批量外场观星测试照片

Fig.9 The outfield stargazing photos of Pi Star Tracker in commercial satellite

图10 皮星敏产品1精度图

Fig.10 The precision of Pi Star Tracker 1

图11 皮星敏产品2精度图

Fig.11 The precision of Pi Star Tracker 2

图10图11可知,皮星敏产品1三轴总误差分别为x=4.84,y=4.50,z=39.81;产品2三轴总误差分别为x=4.63,y=4.79,z=37.02。该2台产品均满足测量精度要求X轴、Y轴5,Z轴50的指标要求。

3 结束语

为了提升星敏感器生产制造中测量精度一致性,提高星敏感器研产水平,本文利用星敏感器研产中星点图像的大数据分析和建模手段,提出了“7参数图像量化评价标准”。作者按照这个标准对皮星敏产品光机组件装调进行应用,通过后续外场观星测试验证满足产品的测量精度要求,为提高后续产品研产的测量精度一致性夯实基础。
[1]
郑循江. 轻小型高动态星敏感器技术研究[D]. 上海: 上海交通大学, 2012.

ZHENG X J. Research on light and small high dynamic star sensor technology[D]. Shanghai: Shanghai Jiao Tong University, 2012.

[2]
潘迪, 温兆伦, 任平川, 等. 高动态下的星体目标质心提取方法[J]. 飞控与探测, 2022, 5(5): 19-25.

PAN D, WEN Z L, REN P C, et al. Star spot centroid extraction method based on difference hash[J]. Flight Control & Detection, 2022, 5(5): 19-25.

[3]
练达, 周琦, 余路伟, 等. 高动态条件下星点像斑建模与补偿[J]. 飞控与探测, 2020, 3(3): 86-94.

LIAN D, ZHOU Q, YU L W, et al. Modeling and compensation of star spot in high dynamic condition[J]. Flight Control & Detection, 2020, 3(3): 86-94.

[4]
曾思康, 赵汝进, 马跃博, 等. 基于事件的高动态星敏感器星点提取方法[J]. 光子学报, 2022, 51(9): 0912003.

ZENG S K, ZHAO R J, MA Y B, et al. An event-based method for extracting star points from high dynamic star sensors[J]. Acta Photonica Sinica, 2022, 51(9): 0912003.

[5]
侯娅贤, 赵汝进, 马跃博, 等. 一种基于自适应滤波的高动态星敏感器在轨校正方法[J]. 光子学报, 2021, 50(2): 163-173.

HOU Y X, ZHAO R J, MA Y B, et al. An on-orbit correction method for high dynamic APS star tracker based on adaptive filtering[J]. Acta Photonica Sinica, 2021, 50(2): 163-173.

[6]
李小燕. 高精度星敏感器光学系统装调工艺与检测方法研究[D]. 西安: 西安电子科技大学, 2018.

LI X Y. Research on assembly and adjustment technology and detection method of high-precision star sensor optical system[D]. Xi'an: Xidian University, 2018.

[7]
魏新国, 张广军, 樊巧云, 等. 利用仿真星图的星敏感器地面功能测试方法[J]. 红外与激光工程, 2008, 37(6): 1 087-1 091.

WEI X G, ZHANG G J, FAN Q Y, et al. Ground function test method of star sensor using simulated sky image[J]. Infrared and Laser Engineering, 2008, 37(6): 1 087-1 091.

[8]
陈新华, 季轶群, 沈为民. 基于星点图像的小像差复原[J]. 光学精密工程, 2012, 20(4): 706-711.

CHEN X H, JI Y Q, SHEN W M. Small-aberration retrieval based on spot images[J]. Optics and Precision Engineering, 2012, 20(4): 706-711.

[9]
张新宇, 郝云彩. 星敏感器星点能量分布数学模型及其对精度的影响分析[J]. 空间控制技术与应用, 2013, 39(3): 14-18, 23.

ZHANG X Y, HAO Y C. Analysis of the star image energy distribution mathematical model and its effect on the accuracy for a startracker[J]. Aerospace Control and Application, 2013, 39(3): 14-18, 23.

[10]
吴庭翱. 基于能量分布的星点提取方法[J]. 信息技术, 2014, 38(11): 178-182.

WU T A. A method of star extraction based on the energy distribution of the starpoint[J]. Information Technology, 2014, 38(11): 178-182.

文章导航

/