CMOS探测器星敏感器的跟踪窗口星图处理方法.pdf

第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）CMOS探测器星敏感器的跟踪窗口星图处理方法王燕清，郑循江，钟金凤，杜伟峰，张磊（上海航天控制技术研究所，上海 201109）摘要:近年来，商业航天微小卫星已经实现了轻小型化、低成本、快速研制等目标，但是在卫星姿态控制方面受制于星敏感器等测量敏感器的精度水平，无法实现高精度的测量与控制。为满足微小卫星在轨飞行中高精度姿态测量和大角速度运动的要求，本文提出了适用于 CMOS 探测器星敏感器的跟踪窗口星图处理算法方案，设计了基于姿态信息和导航星库，通过坐标系旋转和透视投影变换，预测探测星点成像位置，设置 CMOS探测器跟踪窗口的方法，根据跟踪窗口设置，完成预测恒星所在图像子区域的像素灰度数据采集和背景灰度值的计算。采用这种方法后，恒星成像质心精度得到了提高，加快了星点提取质心的图像处理速度。通过试验室动态光星模和外场观星试验验证。试验结果表明：该算法提升了星敏感器的精度（2/16）和动态性能（3（）/s），同时提高了全天识别的概率和跟踪数据有效率。关键词:星敏感器；星跟踪；星图处理；跟踪窗口中图分类号:TN 241 文献标志码:A DOI:10.19328/ki.20968655.2023.04.016A Star Image Processing Algorithm of Tracking Window for CMOS Detector Star SensorsWANG Yanqing，ZHENG Xunjiang，ZHONG Jinfeng，DU Weifeng，ZHANG Lei（Shanghai Institute of Spaceflight Control Technology，Shanghai 201109，China）Abstract:In recent years，commercial aerospace micro-satellites have achieved the goals such as light weight，small size，low cost，and rapid development.However，in the aspect of satellite attitude control，owing to the accuracy level of measurement sensors such as star sensors，it is hard to achieve high-precision measurement and control.In order to meet the requirements of high-precision attitude measurement and large angular velocity motion for micro-satellites in orbit，in this paper，a star image processing algorithm of tracking window for CMOS detector star sensors is proposed.Based on the attitude information and navigation star library，the imaging positions of the detecting star points can be predicted through coordinate system rotation and perspective projection transformation.The CMOS detector tracking window is set according to the tracking window setting，by which the pixel grayscale data are collected and the background grayscale values for predicting the sub regions of stars in the image are calculated.With the proposed algorithm，the accuracy of star centroid is improved，and the image processing speed of star point centroid extraction is accelerated.Laboratory dynamic light star mode and field observation experiments are carried out for verification.The results show that the proposed algorithm improves the accuracy（2/16）and dynamic performance （3（）/s）of the star sensor，while enhances the probability of all-day recognition and tracking data efficiency.Key words:star sensor；star tracking；star image processing；tracking window收稿日期：20230312；修回日期：20230406基金项目：国家自然科学基金（62005158）作者简介：王燕清（1985），男，硕士，高级工程师，主要研究方向为空间光电探测等。119第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）0引言 姿态测量系统是航天飞行器的重要组成部分，它是控制航天器稳定或平稳运动的前提，而姿态测量系统中的首要关键部件就是姿态敏感器。随着航天技术的不断发展，对航天飞行器姿态的测量精度也提出了更高的要求。尤其商业微小卫星受体积、质量、功率、成本等条件约束，需要配置低功耗、轻小型、低成本的姿态敏感器。星敏感器是一种高度自主的姿态敏感器，通过探测器对星空成像，测量恒星矢量在星敏感器坐标系中的分量，并利用已知恒星的精确位置来确定卫星相对于惯性坐标系的三轴姿态。由于星敏感器是迄今为止最精密且漂移最小的姿态测量部件，已经成为卫星上必不可少的姿态敏感器1-3。星敏感器的测量精度和动态性能是衡量星敏感器性能的最重要指标。分析影响测量精度和动态性能的因素，提高星敏感器的姿态测量精度和动态性能，是星敏感器研究中需要解决的关键问题。对地遥感测量卫星的分辨率已经达到了亚米级，对于卫星的控制和测量精度的要求已经优于 5（）/s。卫星初始入轨、快速机动、大角度调姿等条件对星敏感器动态性能提出了很高要求4-5。传统的星敏感器在精度和动态性方面无法满足现有需求6-9。星敏感器开始工作时为初始姿态捕获模式，进行全天的星图匹配和识别，所以搜索和匹配花的时间比较多，一般识别过程需要秒级的时间，这对于实时姿态的输出是不利的。经过全天星图识别获得初始姿态后，星敏感器自动进入星跟踪模式。全天识别初始捕获时间一般小于 5 s，随后的大部分时间星敏感器处于星跟踪模式，而跟踪模式不需要进行全天搜索和识别，所以跟踪模式处理时间短。初始姿态捕获只有在起始或跟踪丢失时才进行，如果跟踪模式非常稳定，则经过初始姿态捕获，星敏感器就一直处于实时跟踪的状态。因此实时跟踪是星敏感器的主要工作模式之一，直接影响星敏感器的整体性能。星跟踪模式下星图成像质量和快速处理的能力的提高，对星敏感器整体的精度和动态性能的影响非常明显，因此研究此星图处理算法是非常有意义的10-13。目前国内外主流星敏感器公司主要有 2 类产品，其中，1 类产品是高精度高可靠星敏感器，其主要特点为精度高、寿命长、质量体积大和价格高；另1 类产品是微小型星敏感器，其主要特点为质量尺寸小、价格便宜，精度、数据更新率等性能一般。根据目前商业卫星对控制和测量精度的要求，国内外的星敏感器在精度和动态性能指标上都不能满足新的航天飞行器的要求。如德国 Jena 的 ASTRO CL 星 敏 感 器 的 静 态 测 量 精 度 为 6/48，法 国Sodern 公司的 AURIGA SA 星敏感器的静态测量精度为 6/40。现有的星跟踪算法主要有 3 类，第 1 类为窗口法，第 2类为质心跟踪法，第 3类为滤波算法。文献14 提出的星跟踪算法需要先对图像根据阈值粗筛检测到星点像素，然后在其周围开窗口求取精确的质心坐标。该算法需要先求取阈值，并对全帧图像数据进行筛选，处理时间长。文献 15-16 提出的星跟踪算法主要根据之前连续两帧的姿态数据来计算角速度和预报当前的姿态。该算法提供了预报的星像坐标，并没有给出具体的指导星敏感器在该像素坐标开窗、探测器的驱动和采集等详细方法。文献 17 提出的主要是基于基于卡尔曼预测的跟踪算法，该算法需要借助先验信息对星点轨迹的建模，在计算机平台上完成了星跟踪仿真实验，并未给出该算法应用到具体星敏感器硬件平台的具体实施方法和测试效果。文献 18 提出的星跟踪算法，利用星矢量与角速率间关系进行星点预测，但是在大动态角速度下，角速率估计会误差随之增大。文中的星图处理算法，深入分析了影响星敏感器精度和动态性能的成因，并根据星敏感器的电路框架以及 CMOS探测器特点，给出了一种星跟踪模式下的星图处理算法，根据之前帧的姿态数据预测当前的姿态，并利用预测姿态信息指导星敏感器下一帧的开窗区域，给出了探测器的图像驱动和采集具体方法，应用到实际的星敏感器硬件平台上，并通过动态光星模和外场观星试验验证，数据分析结果证明该算法可有效提高星敏感器整机的精度和动态性能。1星图处理算法 1.1星跟踪原理星敏感器在全天识别成功后，星敏感器即进入星跟踪模式，在星跟踪模式下的图像处理算法由以120第 40 卷 2023 年第 4 期王燕清，等：CMOS 探测器星敏感器的跟踪窗口星图处理方法下步骤组成：1）跟踪准备，角速率估计；2）预测光轴指向；3）预测视场内的导航星；4）预测的导航星与实测星点匹配；5）定姿星筛选；6）基于 ESOQ2算法进行姿态解算；7）Kalman滤波；8）更新姿态四元数和角速率。星跟踪的目的是：在具备了姿态先验信息后，再快速建立数字星图中已提取出的星点与导航星库中存储恒星的对应关系。全天识别成功后，已算得全天捕获成功后的四元数 Q=q0，q1，q2，q3，（q0为标量，下同），星敏感器即自主转入星跟踪模式。为了方便描述该星跟踪算法，后面论述中假设星跟踪过程中的第 k-1 帧星图的识别结果已知，对当前第 k帧的观测星图（星敏感器采集到的星图）进行跟踪识别，星跟踪算法原理框图如图 1所示，图中对每个步骤进行详细描述。图中 ESOQ2（Second EStimator of the Optimal Quaternion）为姿态确定算法，Kalman 滤波器是一种利用线性系统状态方程，通过系统输入输出观测数据，对系统状态进行最优估计的算法。步骤 1）描述星跟踪的准备工作主要根据前几帧星图的姿态四元数估计出当前帧的星敏感器的角速率。步骤 2）描述预测光轴指向是基于第 k-1帧算得的姿态四元数Qk-1，可以计算对应的星敏感器光轴指向r k-1。步骤 3）根据星表中保存有导航星在天球坐标系下的坐标（赤经、赤纬）、亮度（星等）、星号等信息，得到在光轴指向r k-1下某一定天区范围（星敏感器视场）的导航星坐标集合；步骤 4）经过星敏感器透视投影变换，即将此范围内导航星的坐标从天球坐标系向图像坐标系的变换（将星表中星体的赤经、赤纬信息转化为图像传感器感探测器面上的位置坐标信息），可以得到第 k 帧参考星图，此参考星图中所有的导航星信息均已知，建立跟踪星列表，与探测星列表进行匹配。以各导航星坐标为中心，在半径为 2R（R 为星跟踪半径，取值与星敏感器运动的角速度有关）的圆形范围内，查看第 k帧星点质心数据包中是否存在相应的实测星点，存在且唯一，则完成匹配；若预测的导航星与实测星点匹配失败，则退回至全天识别模式，若匹配成功执行步骤 5）步骤 8），随后进行下一拍星跟踪。利用第 k帧观测星图的识别结果产生第 k+1帧参考星图，用以实现对第 k+1帧观测星图的识别。如此循环，可以实现星体的跟踪。1.2跟踪窗口设置星敏感器进行全天球识别时，因无当前天区的先验姿态位置信息，所以需要驱动星敏感器的探测器处理模块对拍摄的整幅星图进行处理，并提取星点，即采用全帧式图像处理。全天识别成功并积累多拍姿态信息后星敏进入稳态星跟踪模式，在星跟踪模式下星敏感器根据先验姿态和角速度信息，驱图 1星敏感器星跟踪算法原理框Fig.1Principle block diagram of the star sensor tracking algorithm121第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）动探测器模块只对部分图像进行星点提取处理，即采用局部式图像处理。本文利用星跟踪模式下已知的姿态和角速度先验信息，将当前天区视场内星库中的星点（导航星）转换到探测器平面内，并根据该导航星理论位置对该区域像素进行处理和星点提取。以大小为 2 0482 048 像素的探测器面阵为例，该窗口跟踪设置方式具体描述如下：1）根据恒星在星敏感器探测器像平面上成高斯弥散斑的大小特征19，星敏感器系统的 80%的能量会集中在 44像素的区域内，一般恒星的能量分布范围都在 88像素的区域内，因此将探测器面阵以 88像素为单元，均等划分为 65 536个网格。2）根据探测器图像数据逐行依次读出的物理特性，行方向对齐的 256个网格组成 1个图像条带，整个探测器面阵被分割成 256个条带，如图 2所示3）根据当前星敏姿态和角速度信息预测星敏视场天区内导航星位置，并记录其所在的网格编号作为待处理的窗口波门图。当导航星点坐标与网格边界的距离较近时，相邻网格也被记录，因此一个波门图像可能由 1个、2个或 4个网格组成。波门采集如图 3 所示，设某颗恒星的预测位置由黑色十字标示，该恒星所在网格由阴影标示，提出波门采集要求时，应在图像条带中选择包含预测星点弥散斑的波门，将波门对应的行坐标与列坐标记录，作为每个网格的位置（x，y）坐标。4）因为将所有导航星对应的网格图像进行存储、处理会增大存储器成本和 CPU 处理的运算量，通常会对处理的网格总数或条带总数进行限制。当导航星点对应的波门网格数超过限制时，按星等优先级进行控制。5）根据预测的导航星位置，对上述波门图像进行处理并提取星点，完成波门窗口设置以及波门内星点提取。以星敏感器 CMOS 探测器 CMV4000 为例，从存储空间资源、图像读出时间和质心定位精度等几个方面，来分析传统图像处理方式和基于跟踪窗口设置的图像处理方式之间的差异。传统的星图处理算法主要是分为 2 个部分，首先，处理器读取图像数据并对其进行阈值分割、连通性分析和内插细分定位后，得到星点的有效坐标信息并保存，完成对图像的星点提取工作。其次，处理器根据星点信息采用三角形算法，使星点构成的三角形模式与导航星库中的三角形进行匹配。这种传统的星图处理算法由于需要对整个图像进行读取和处理，需要的时间较长，而且容易引入探图 2探测器面阵分割Fig.2Division diagram of the detector array图 3预测星点波门Fig.3Predicted star gate diagram122第 40 卷 2023 年第 4 期王燕清，等：CMOS 探测器星敏感器的跟踪窗口星图处理方法测器的噪声和环境干扰，产生假星点影响整体精度，同时计算速度慢影响动态性能26-27。传统图像处理方式需要将整个探测器图像进行存储，探测器有 2 048 行，2 048 列像素，每个像素按 照 8 bit 存 储，需 要 2 0482 0488 bit=33 554 432 bit存储空间。跟踪窗口设置图像处理方式，根据实际星敏感器姿态算法需要的探测星数要求 30颗（上限）为例，并且按照星点处于 4 个网格的交界处最差情况，最多需要的存储空间为 304888=61 440 bit存储空间。从图像输出时间分析，根据图像传感器CMV4000的数据手册可以查询到，读出时间为T=（129Tm16Noutput）Nline（1）式中：T 为图像读出时间；Tm为探测器输入主时钟周 期；Noutput为 图 像 输 出 通 道 数；Nline为 图 像 输 出行数。传统图像处理方式需要驱动探测器输出整张图像，采集后进行处理，探测器输入主时钟频率为 5 MHz，时钟周期 Tm为 200 ns，按照 4 通道输出，Noutput为 4，整图输出行数 Nline为 2 048。因此可以计算得到，图像传感器输出整张图像时间 T为 211.4 ms。跟踪窗口设置图像处理方式由于只要输出预测星点附近的图像数据，根据实际星敏感器姿态算法需要的探测星数要求 30颗（上限）为例，并且按照星点处于 4 个网格的交界处最差情况，最多需要的输出的图像行数为 3028=480 行，图像输出时间 T1为 49.5 ms。星敏感器是姿态测量敏感器，在相同的姿态识别算法条件下，姿态精度取决于测量星点坐标的准确度，星点坐标的计算中阈值起到关键作用，阈值用来分割图像中的星点和背景，其取值的合理会影响到星点质心的精度。跟踪窗口设置图像处理算法是以每个导航星 4个网格的形式（1616 pix）存储波门图像，以预测的导航星坐标为中心，处理 1010 pix 范围内的图像数据，在提取星点前，首先确定星点提取阈值，阈值在计算得到的灰度均值基础上，累加灰度偏移量得到，计算方法如下：1）计算 1616 pix.范围内的灰度均值 m；2）星点提取阈值T=m+Offset（2）式中：m 为该像素所在子区域的灰度均值；Offset为灰度偏移量。在传统的星图处理算法中，Offset的值根据灰度均值自适应调整，考虑到整个图像阵列内未知的噪声和坏像元等因素，为了避免错误信息的干扰，一般取灰度均值 1/2，传统的星图处理算法中的 m值的选取是整张星图的灰度平均值，由于星空背景大多情况下是黑暗的，因此这个灰度均值相对于恒星星点附近能量的均值存在误差，通常是过低，均值过低会导致星点提取阈值偏低，从而会受到噪声干扰，提取出假星点，影响星图识别的姿态精度。在跟踪窗口设置图像处理算法中，由于星点信息是之前四元数姿态信息预测计算得到，因此只要考虑星点附近窗口内的像素，受到像素随机噪声、坏像元和环境干扰的因素降低，在求得子区域的灰度均值后，灰度偏移量的值可以取小一些。因此，跟踪窗口图像处理方法既避免了全局图像灰度均值带来的噪声、坏像元等环境因素的干扰，同时在窗口范围内计算得到更加精确的均值，再通过降低阈值偏移量，使得星点弥散斑的效应充分得到利用，筛选出的有效像素数量提升 4 倍。在后续应用质心法计算星点中心坐标时，质心定位精度可以从传统的 1/20 像素提升到 1/50像素。星点质心精度的提升，可以加快三角形角距匹配速度，减少误匹配的概率，为最后姿态识别和精度提升提供保证，传统星图处理算法和跟踪窗口设置图像处理算法的对比见表 1。相对于传统方式，将整个探测器图像进行存储并从整个探测器图像中通过逐行扫描提取星点方式，基于跟踪窗口设置的图像处理方式大幅降低了图像存储空间和处理时间，在降低存储器成本的同时提高姿态数据更新频率，从而实现星敏感器动态性能和的提升。同时预测姿态信息选择合适的窗表 1传统星图处理算法和跟踪窗口设置图像处理算法的对比Tab.1Comparison between the traditional star image processing algorithm and the proposed tracking window algorithm图像处理方式传统图像处理方式跟踪窗口设置图像处理方式存储空间/B33 554 43261 440图像读出时间/ms211.449.5数据更新率/Hz510质心定位精度/像素1/201/50姿态精度/()s-1)4.52.1123第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）口范围，选取合适的阈值，充分利用星点弥散效应，提高星点质心定位精度，从而实现提高星敏感器姿态精度的效果。2试验结果 2.1试验方法星敏感器主要性能指标很大程度上与星图处理算法设计相关。基于外场观星试验对本文所述算法进行验证，验证内容如下：1）静态精度；2）动态性能；3）数据有效率。精度测试试验方法：在地速条件下，垂直对天，利用地速的稳定性。分析星敏感器垂直对天时的观星实测数据，画出数据曲线。以多项式拟合四元数样本点，得到参考四元数曲线，将实测四元数与参考四元数做差，并转化为三轴姿态角，可得到姿态误差曲线20-25。动态性能测试试验方法：控制二维观星转台绕产品 X 或 Y 轴以一定角速度旋转运动，观察不同速度下遥测数据中的数据有效标志以及星点数。2.2实际测试外场观星对比试验中，应用传统星图算法的星敏感器可稳定跟踪的最高角速度为 0.8（）/s，采用跟踪窗口设置图像处理算法的星敏感器可稳定跟踪的最高角速度为 3（）/s。测试方法：控制二维观星转台绕产品 X 或 Y 轴以一定角速度旋转。控制转台，使星敏感器按照 0.6、0.7、0.8、0.9、1.0、2.3、2.8、3.0（）/s的速度运动，观察不同速度下遥测数据中的数据有效标志以及星点数，有效标志以及星点数曲线如图 4所示。图 4传统图像处理算法星敏感器在不同角速度下数据有效标志以及星点数曲线Fig.4Valid data marks and star points of the star sensor obtained by the traditional star image processing algorithm under different angular velocities124第 40 卷 2023 年第 4 期王燕清，等：CMOS 探测器星敏感器的跟踪窗口星图处理方法由图 4 可得星敏感器在以角速度 0.6、0.7（）/s机动过程中，可以稳定跟踪，数据有效率为 100%，星敏感器在以角速度 0.8（）/s机动过程中，可以稳定跟踪，数据有效率为 99.2%，在星敏感器 0.9（）/s机动过程中，数据有效率为 73.3%，不能稳定跟踪。用跟踪窗口处理算法星敏感器三轴角速度曲线如图 5 所示，在不同角速度下数据有效标志以及星点数曲线如图 6所示。图 5采用跟踪窗口处理算法星敏感器角速度曲线Fig.5Angular velocity curves of the star sensor processed by the tracking window algorithm图 6跟踪窗口处理算法星敏感器在不同角速度下数据有效标志以及星点数曲线Fig.6Valid data marks and star points of the star sensor obtained by the tracking window processing algorithm algorithm under different angular velocities125第 40 卷 2023 年第 4 期上海航天（中英文）AEROSPACE SHANGHAI（CHINESE&ENGLISH）由图 6 可得，星敏感器以角速度 1.0、2.3、2.8、3.0（）/s 转动时，数据有效率 100%，可以稳定跟踪。传统星图处理算法的星敏感器三轴姿态精度如图 7所示。如图 7 所得，星敏感器使用传统星图算法，垂直对天观测数据进行分析，计算三轴姿态精度，得到 三 轴 姿 态 精 度 分 别 为 Z 轴 4.41（）/s，Y 轴 4.64（）/s，X 轴 27.23（）/s。应用跟踪窗口设置星 图 处 理 算 法 的 星 敏 感 器 三 轴 姿 态 精 度 如 图 8所示。如图 8 所得，星敏感器使用跟踪窗口星图处理算法，垂直对天观测数据进行分析，计算三轴姿态精度，得到三轴姿态精度分别为 Z 轴 2.05（）/s，Y轴 2.17（）/s，X轴 15.91（）/s。星敏感器采用不同星图处理算法的主要性能参数对比见表 2。图 7传统星图处理算法的星敏感器姿态精度Fig.7Attitude accuracy of the star sensor based on the traditional star image processing algorithm图 8应用跟踪窗口设置星图处理算法的星敏感器姿态精度Fig.8Attitude accuracy of the star sensor obtained by the tracking window algorithm表 2星敏感器不同星图处理算法主要性能参数对比表Tab.2Main performance parameters of the star sensor obtained by different algorithms测试项目静态精度（3）动态性能的跟踪数据有效率/%Z轴Y轴X轴0.8()/s0.9()/s1.0()/s2.3()/s传统星图处理算法4.414.6427.2399.277.250.0无法全天识别跟踪窗口设置星图处理算法2.052.1715.91100100100100126第 40 卷 2023 年第 4 期王燕清，等：CMOS 探测器星敏感器的跟踪窗口星图处理方法3结束语 本文提出的这种星图处理算法，在星跟踪环节，基于姿态信息和导航星库得到预测导航星，并对其升序排列，通过坐标系旋转和透视投影变换，得到预测星点坐标。根据星敏感器选用 CMV4000探测器特点，对预测星点相邻区域进行合并得到探测器开窗的起始和结束位置，星图读出时间减少 1/3，提高数据更新率。对预测星图与实测星图进行匹配，求取新的姿态信息。该星图处理算法通过外场观星试验数据的分析，星敏感器的主要性能指标精度从4.5（3）提升到 2.1（3），动态性能从 0.8（）/s提升到 3.0（）/s，证明在精度和动态性能方面均有显著效果。给出了一种星敏感器的星图处理方法和试验验证结果，为提高星敏感器的姿态精度、动态跟踪以及数据有效率等核心性能提供了参考。参考文献1 李新鹏，任平川，高原，等.多探头星敏感器分布式视场融合方法 J.飞控与探测，2018（1）：071-076.2 程会艳，郑然，武延鹏，等.一种提高多探头星敏感器姿态测量精度的方法J.空间控制技术与应用，2022（4）：78-85.3 王燕清，吴永康，沈杰，等.国产元器件在星敏感器的应用设计和验证测试J.导弹与航天运载技术（中英文），2023（1）：59-65.4 练达，周琦，余路伟，等.高动态条件下星点像斑建模与补偿 J.飞控与探测，2020，3（3）：086-094.5 陈纾，张广军，郑循江，等.小型化星敏感器技术 J.上海航天，2013（4）：69-78.6 郑循江.轻小型高动态星敏感器技术研究 D.上海：上海交通大学，2012：15-17.7 王军，何昕，魏仲慧，等.基于多特征匹配的快速星图识别 J.光学精密工程，2019（8）：1870-1879.8 金荷，翟正一，杜伟峰，等.星敏感器热稳定性的试验分析方法 J.中国激光，2020（4）：204-209.9 隋杰，程会艳，余成武，等.星敏感器光轴热稳定性仿真分析方法 J.空间控制技术与应用，2017（4）：37-41.10 江洁，李霄，张广军，等.一种星敏感器中快速星跟踪方法 J.宇航学报，2006（27）：953-955.11 LIEBE C C.Accuracy performance of star tracker atutorial J.IEEE Trans on AES，2002，38（2）：587-599.12 余路伟，毛晓楠，金荷，等.星敏感器抗杂光背景滤波图像处理方法研究 J.上海航天，2016（4）：26-31.13 余路伟，毛晓楠，周琦，等.采用最大背景估计的星敏感器图像处理方法 J.激光与红外，2017（7）：889-895.14 管宇，张迎春，王常虹.一种适用于星敏感器的大角速度下星跟踪算法J.中国惯性技术学报.2012（4）：440-443.15 李葆华，陈希军，郑循江，等.星敏感器高动态下自主星跟踪算法 J.红外与激光工程，2012，41（1）：191-19516 金雁，江洁，张广军.高动态星体目标提取方法 J.红外与激光工程，2011，40（11）：2281-2285.17 金荷，毛晓楠，孙少勇，等.基于恒星矢量差分的星敏感器星跟踪算法研究 J.上海航天，2017（2）：99-105.18 郑循江，张广军，毛晓楠.一种甚高精度星敏感器精度测试方法 J.红外与激光工程，2015（5）：1605-1609.19 刘太阳，王仕成，李庆会，等.基于高斯弥散参数拟合的伪星检测 J.激光与红外，2011（6）：678-682.20 郭敬明，杨轻云，贾建禄.高精度星敏感器动态测量精度检测方法 J.科技创新导报，2015（26）：75-76.21 任平川，周琦，金荷，等.星敏感器动态精度评价方法J.红外与激光工程，2022（3）：349-356.22 龚德铸，武延鹏，卢欣.一种提高星敏感器动态性能的方法 J.空间控制技术与应用，2009（6）：19-23.23 向明江，信思博，姜通，等.高分五号卫星星敏感器在轨测量精度评估 J.上海航天，2019（36）：51-55.24 胡雄超，吴永康，闫晓军，等.一种星敏感器标定误差修正方法研究 J.上海航天，2016（6）：88-92.25 徐樱，吴德安，汪礼成，等.星敏感器慢变误差校准方法研究 J.上海航天，2016（4）：63-69.26 李杰.星敏感器系统的研究与实现 M.长春：吉林人 民出版社，2008.27 李德良，阮锦.一种适用于星敏感器的星点提取方法 J.激光与红外，2009（12）：1348-1350.续表 2星敏感器不同星图处理算法主要性能参数对比表Continue tab.2Main performance parameters of the star sensor obtained by different algorithms测试项目2.8()/s3()/s传统星图处理算法无法全天识别无法全天识别跟踪窗口设置星图处理算法100100127