Rapid star pattern matching for precisely tracking telescopes

精密跟踪型望远镜适配的快速星图匹配 Rapid star pattern matching for precisely tracking telescopes corresp first-author E-mail： liudl@cho.ac.cn E-mail： liudl@cho.ac.cn 刘德龙（1988-），男，黑龙江哈尔滨人，博士，中科院特别研究助理，2011年于哈尔滨工业大学获得学士学位，2015年于中国科学院大学获得硕士学位，2021年于法国图卢兹第三大学获得博士学位，主要从事目标天文定位及光谱特征分析的研究工作。E-mail： liudl@cho.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226267&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226277&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226272&type= 22.01333237 32.00399780 刘德龙 （1988-），男，黑龙江哈尔滨人，博士，中科院特别研究助理，2011年于哈尔滨工业大学获得学士学位，2015年于中国科学院大学获得硕士学位，2021年于法国图卢兹第三大学获得博士学位，主要从事目标天文定位及光谱特征分析的研究工作。E-mail： liudl@cho.ac.cn liudl@cho.ac.cn 杨文波（1980-），男，吉林省吉林市人，博士，副研究员，2014年于中国科学院长春光学精密机械与物理研究所获得博士学位，主要从事天体测量及自动控制方面的研究。E-mail： yangwb@cho.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226282&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226287&type= 22.01332855 32.00399780 杨文波 （1980-），男，吉林省吉林市人，博士，副研究员，2014年于中国科学院长春光学精密机械与物理研究所获得博士学位，主要从事天体测量及自动控制方面的研究。E-mail： yangwb@cho.ac.cn yangwb@cho.ac.cn 中国科学院 国家天文台长春人造卫星观测站，吉林 长春 130117 中国科学院 国家天文台长春人造卫星观测站，吉林 长春 130117 Changchun Observatory， National Astronomical Observatories， Chinese Academy of Sciences， Changchun 130117， China Changchun Observatory， National Astronomical Observatories， Chinese Academy of Sciences， Changchun 130117， China 空间目标的精密定轨需要高精度的天文定位技术，而星图匹配是天文定位的基础。本文针对精密跟踪型望远镜提出一种快速星图匹配算法，它包含优化三角形匹配和序列图像修正两部分。论文主要叙述了该算法原理、流程与实现，并对其速度与精度等方面进行了研究。首先，利用编码器轴系定位划取指向天区的星表数据，经筛选和归算列为导航星表。接着，应用降维查表的方法加速三角形匹配并通过理想坐标底片常数之间的关系校验得出首帧匹配结果。然后，修正算法结合后续各帧中望远镜指向的变化量计算导航星的理想坐标，并应用上一帧的底片常数，匹配观测星和导航星。最后，对底片常数计算得出的导航星和目标定位结果进行统计对比。经实验，计算39组星对时采用快速三角形匹配可以将时间缩短至近1/300；对后续图像（每帧含大约100组匹配星对）使用序列图像修正算法，均可以在0.04 s以内完成匹配；采用快速星图匹配算法获取的匹配星对用于对中高轨激光星定位，其平均误差在0.5″左右。由此可见，快速星图匹配算法充分满足精密跟踪型望远镜天文定位精度高、速度快的要求。 Star pattern matching is critical when using high-precision astronomical positioning technique， which， in turn， is required for the precise determination of space objects’ orbits. This paper proposes a rapid star-pattern-matching method for precisely tracking telescopes， including its principles， program flow and implementation. First， the catalog of pointed sky region was filtered out with the shafting positioning from the encoder； then， after selection and reduction， the stars in the area were compiled into a navigation stars list. After which， the optimized triangle matching method， based on dimension reduction and table look-up methods， was employed in the first frame， where the relation between plate constants verifies the success of matching. The changes of the telescope pointing were considered to calculate the standard coordinates of navigation stars， and the observed stars were matched through the last plate constants. Finally， the coordinates of navigation stars calculated were compared by the plate constants with their reduced values from the catalog， and the position of one space object was calculated through celestial positioning. Experimental results indicated that 39 star pairs could be matched by applying the optimized triangle matching proposed， however， with only approximately 1/300 of the time taken compared to the conventional method. When calculating the following sequential frames （approximately 100 star pairs per-frame）， the matching could be finished in less than 0.04 s. The matched star pairs obtained by the rapid star pattern matching algorithm were used to locate an MEO laser satellite with an average error of approximately 0.5″. As a result， the method proposed fully satisfies the requirements of high accuracy and speed in astronomical positioning for precisely tracking telescopes. 星图匹配 天文定位 三角形匹配 图像序列 底片常数 star pattern matching astronomical positioning triangle matching image sequence Plate constants 1 引 言 1 近年来，大气层外围的卫星、空间碎片、火箭箭体和载荷等人造物体不断增加，空间环境不断恶化，严重威胁了航天器的正常运行，为此需要对目标进行精密定轨、编目、预报及特征分析。基于地基望远镜的光电观测目前仍是最有效的空间目标研究手段之一，它具有观测轨道高度广、组网成本低、无接触及免干扰等优势。其中，精密跟踪型望远镜通过较长弧段观测，取得已知轨道目标的定位数据，为其精密定轨服务 ［ 1 1 ］ 。天文定位基于背景恒星的绝对位置克服了轴系定位中机械误差、脱靶量测量误差和坐标系转换等带来的影响，可显著提高空间目标的定位精度。 星图匹配是天文定位中不可或缺的一环，它将观测图像中恒星的摄影坐标和标准星表中的天球坐标关联起来，以此解算观测系统的底片模型。星图匹配主要依靠子图同构 ［ 2 2 - 3 3 2-3 ］ 和模式识别 ［ 4 4 - 5 5 4-5 ］ 。子图同构中的三角形相似算法仍是至今使用最广泛和有效的星图匹配算法，然而传统算法由于特征维数较低，并且需遍历所有星象，导致匹配效率低，导致既耗时又浪费存储空间 ［ 6 6 ］ 。随着望远镜光学口径和视场的不断增大，视场中的星象呈指数级增加。同时天文相机技术也不断进步，越来越多望远镜配备了高帧频、高分辨率的科学级CMOS相机。这导致在一个观测周期内图像序列的帧数和每帧数据急剧增加，进一步凸显了传统三角形匹配算法的不足。近年来见诸报道的一些改进算法包括：2019年王军等人 ［ 7 7 ］ 采用三角形外接圆和内切圆半径乘积作为匹配特征，初步筛选后再进行三角形三边匹配，2020年刘先一 ［ 8 8 ］ 采用径向特征量缩小了待匹配导航星数量，增强了匹配的针对性；2021年李变等人 ［ 9 9 ］ 公开的发明专利提出基础三角形的概念，预筛选观测三角形，降低了匹配量；2021年徐伟等人 ［ 10 10 ］ 选择观测星中星等最小的三颗组成特征三角形，后采用三边组合特征值的方式进行导航星表的遍历匹配；2021年熊琨等人 ［ 11 11 ］ 的专利也提出缩减导航星，并建立导航星角距库的思想。这些算法增加了星表或观测星预处理的复杂度，提高了计算成本；前期预选还有可能造成匹配失败，因此不满足序列图像的首帧星图匹配中容易实现和成功率高的要求。 精密跟踪（简称精跟）型望远镜具有高精度的控制伺服和较大的相机靶面。其图像序列的特点首先是目标和背景恒星的相对运动，由于望远镜跟踪目标，恒星星象会出现拉长。因此需使用短曝光或采用星象模糊还原算法，前者增加了图像噪声，后者使得提取位置发生改变。这均需要匹配算法对于不同天区、轨道和类型的目标具有一定鲁棒性，能够自动选取误差容限进行运算；第二，码盘可在拍摄图像的同时获取望远镜指向，即图像中心的轴系定位。该指向较为粗略，但可用于定位星表。同一观测圈次内序列图像帧间的指向改变量非常精准，可用于快速匹配；第三，望远镜指向连续变化，并且帧间变化较小。因此相邻帧的底片常数变化很小，可用上帧常数进行匹配，之后通过匹配更新底片常数；第四，相对于固定方向观测，精跟型望远镜在观测周期内动作范围较大，因此分区指向模型可能会造成结果误差的不连续。基于以上特点本文提出一种专门适配的星图匹配算法。 2 星图匹配模型 1 望远镜电控系统根据引导数据驱动跟踪，同时采集图像。后经图像处理和目标识别给出目标和背景恒星的位置。对于精跟型望远镜较小的视场（<2.5°），六参数模型即可完成匹配。 由于拍摄图像和实际星图只存在比例、平移、旋转或镜像关系 ［ 12 12 - 13 13 12-13 ］ ，只需要3个对应点，就可以确定底片常数。首先采用三点拟合计算图像中心对应的天球坐标 <math id="M1"><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226299&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226294&type= 11.17599964 4.48733330 ，利用该点的天球坐标计算得到导航星的理想坐标，即 <math id="M2"><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226310&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226304&type= 10.41399956 4.48733330 到 <math id="M3"><mo stretchy="false">(</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226322&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226315&type= 9.90600014 4.82600021 ： <math display="block" id="M4"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>ζ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226333&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226327&type= 67.81800079 20.82799911 . 理想坐标实际上是空间星象的地心投影（Gnomonic Projection），其本质是针孔相机的成像关系，即带入较小畸变误差地模拟长焦望远镜。理想坐标将天球系统转化为笛卡尔平面坐标，得出了导航星在各帧不同的位置坐标。目标定位时，将其理想坐标 <math id="M5"><mo stretchy="false">(</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>,</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226341&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226337&type= 10.07533360 4.82600021 转换为天球坐标 <math id="M6"><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>,</mo><msub><mrow><mi>δ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226352&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226346&type= 10.58333302 4.48733330 的变换如 式（2） 式（2） ： <math display="block" id="M7"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mrow><mi>ζ</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>s</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>s</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>s</mi></mrow></msub><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226365&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226356&type= 63.07666779 20.82799911 . 仅依靠3个点可能出现误匹配的情况，此外三点拟合精度也较低，因此必须对考察列表进行校验，并进行多点拟合。对于正确匹配，当图像 <math id="M8"><mi>x</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226376&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226371&type= 2.20133328 2.79399991 和 <math id="M9"><mi>y</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226386&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226381&type= 2.11666679 3.72533321 方向的比例尺相同时，恒星的理想坐标 <math id="M10"><mo stretchy="false">(</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226472&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226469&type= 9.99066734 4.82600021 和其度量坐标 <math id="M11"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226407&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226402&type= 10.75266743 4.91066647 应满足底片常数的对应关系，如 公式（3） 公式（3） （当拍摄图像为镜像时，转换矩阵元素应为相应的下组符号）。通过该组关系可以判定匹配是否正确（ <math id="M12"><mi mathvariant="bold-italic">T</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226417&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226412&type= 2.87866688 2.79399991 是否为酉矩阵 ［ 14 14 ］ ），并且将6底片常数减为3底片常数（ <math id="M13"><mi>a</mi><mo>,</mo><mi>b</mi><mo>,</mo><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226428&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226422&type= 9.22866726 3.21733332 ）： <math display="block" id="M14"><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>ζ</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mi>a</mi></mtd></mtr><mtr><mtd><mi>b</mi></mtd></mtr></mtable></mrow></mfenced><mo>+</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mo>±</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi></mtd><mtd><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi></mtd></mtr><mtr><mtd><mo>∓</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext> </mtext><mi>θ</mi></mtd><mtd><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>θ</mi></mtd></mtr></mtable></mrow></mfenced><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced><mo>=</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><mi>a</mi></mtd></mtr><mtr><mtd><mi>b</mi></mtd></mtr></mtable></mrow></mfenced><mo>+</mo><mi mathvariant="bold-italic">T</mi><mo>·</mo><mfenced open="[" close="]" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226438&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226433&type= 73.66000366 10.24466610 ， 式中： <math id="M15"><mi>a</mi><mo>,</mo><mi>b</mi><mtext> </mtext></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226446&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226442&type= 6.34999990 3.21733332 表示理想坐标系和度量坐标系间的平移量， <math id="M16"><mi>θ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226456&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226451&type= 1.86266661 2.79399991 表示两坐标系间的夹角。当已知一组 <math id="M17"><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226465&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226461&type= 10.75266743 4.91066647 和其对应的 <math id="M18"><mo stretchy="false">(</mo><msub><mrow><mi>ζ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>,</mo><msub><mrow><mi>η</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226472&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226469&type= 9.99066734 4.82600021 时，底片常数便可由最小二乘法得出。 3 算法构成 1 3.1　图像获取与预处理 2 本研究工作所有实际观测都依托于中国科学院国家天文台长春人造卫星观测站运维的1.2米口径地平式光电望远镜的主焦点 ［ 15 15 ］ 。该望远镜位于吉林省吉林市（ 图1 图1 ）。望远镜动态峰值误差为方位轴：2.555″，高度轴：2.170″，为典型的精跟型望远镜。望远镜搭载的相机为Andor iKon-XL 230，4 096×4 108 pixel，像素尺寸为15 μm，视场约为1.5°。望远镜的焦距约为2.2 m。 拍摄得到的图像如 图2 图2 ，应用SExtractor （Source-Extractor）软件 ［ 16 16 ］ 进行星象质心位置的提取。之后采用静态特征和运动特征差异来初选疑似目标并通过连续帧关联确定目标。最后根据望远镜指向筛选Tycho-2星表，剔除双星并选出视星等介于8~10、自行小于10 mas/yr的恒星进行匹配。 3.2　星表数据的资料归算 2 本研究采用Tycho-2恒星星表，该星表是建立于ICRS下的，历元为J2000.0。因此需将星表中的恒星坐标进行自行、周年视差、周年光行差和光线偏转的归算，解算到恒星在观测时刻的GCRS坐标进行星图匹配。本文应用ANCI C版本的SOFA（Standards Of Fundamental Astronomy， from IAU）软件包 ［ 17 17 ］ 中的程序做以上的换算。 几点说明：基于望远镜拍摄图像的天文定位给出的是目标的站心空固坐标，由于恒星站心和地心间的周日视差可忽略，可采用恒星的GCRS坐标进行星图匹配。恒星由岁差—章动、极移引起视位置差别很小，不影响理论与观测星图的匹配。大气折射同时影响背景恒星和目标，精跟型观测中该较差可被忽略。 3.3　优化三角形匹配算法 2 本程序包括两部分算法：优化三角形匹配和序列图像修正。优化三角形匹配算法通过降维查表的方法提高运算速度，得到首帧的底片常数和较为精确的中心位置；该中心位置叠加码盘获得的望远镜指向改变量获得相应帧底片中心位置，并结合上一帧的底片常数进行新一轮的匹配。流程可见 图3 图3 。底片常数可用来解算各像素对应的天球坐标。 本研究中三角形匹配算法利用“边-边-边”模式匹配观测星和导航星（如 图4 图4 所示）。但是传统算法特征维数较低并且组构三角形过多，造成计算大量的冗余匹配，浪费计算时间和存储空间。 优化三角形匹配算法活动图如 图5 图5 所示，观测端编号并计算每两颗星象之间的“边”——角距 ［ 18 18 ］ ，按照自然编号的顺序存储。经归算后星表中的导航星，编号并计算两两之间的角距，按照角距大小排列存储，称为排序导航星表。具有相同观测星的两个星对（例如观测星对： <math id="M19"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226578&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226575&type= 1.10066664 2.79399991 ， <math id="M20"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226571&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226564&type= 1.60866666 3.64066648 和 <math id="M21"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226578&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226575&type= 1.10066664 2.79399991 ， <math id="M22"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226593&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226589&type= 1.69333339 2.79399991 ），即为待测观测星对，分别放缩一定的误差σ后划取排序导航星表中的两个片段。判断这两个片段内是否存在共享导航星的两个星对，如果存在（例如导航星对：A，B和A，C），即为待选导航星对，计算这两个导航星对中不同导航星（即导航星B，C）间的角距。哈希查表（哈希函数如 式（4） 式（4） ，其中 <math id="M23"><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226562&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226559&type= 3.30200005 2.87866688 为筛选星表中导航星数目）得到观测星对两相异星角距（即 <math id="M24"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226571&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226564&type= 1.60866666 3.64066648 ， k 之间角距），放缩该误差后所划定的导航星表范围，判断该范围是否包括对应相异导航星（即导航星B，C）的角距。是就将这些观测星和导航星按照对应关系放入待考察列表中（即 <math id="M25"><mi>i</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226578&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226575&type= 1.10066664 2.79399991 ， <math id="M26"><mi>j</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226586&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226583&type= 1.60866666 3.64066648 ， <math id="M27"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226593&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226589&type= 1.69333339 2.79399991 对应A，B，C）： <math display="block" id="M28"><mi>h</mi><mo>=</mo><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mi>N</mi><mo>-</mo><mn mathvariant="normal">1</mn><mo>-</mo><mi>j</mi></mrow></mfenced><mfrac><mrow><mi>j</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo>+</mo><mi>k</mi><mo>-</mo><mi>j</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226612&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226610&type= 49.52999878 7.87400007 . 完成一组匹配后即得到粗略的定位关系，可求得图像中心对应的天球位置。该位置用来求得导航星的理想坐标，并将这组匹配进行参数拟合得到底片常数。用 式（2） 式（2） 验证该底片常数，满足校验关系后进行观测星和导航星间的匹配，当有足够多匹配时，该底片常数方为正确。当无法得到正确匹配时，需要放大误差再次进行匹配，直到得到正确匹配或误差不可接受为止。放大误差后，导航星表范围增加，需屏蔽已考察的列表片段避免重复匹配，提高匹配效率。 3.4　序列图像修正算法 2 地心投影通过各帧图像中心将导航星绝对的天球坐标投影为本帧理想坐标；对于较高帧频的精跟型望远镜，帧间的底片常数变化很小，因此本帧理想坐标和上一帧参数可完成匹配。但是为了不引入累积误差，完成匹配后还需再一次的多点拟合修正底片常数。序列图像修正算法在不损失精度的前提下点对点地完成匹配，极大地节约存储空间，提高计算速度，这对于实时性解析目标位置有着重要意义。除此之外，对于大视场的望远系统，较多的背景恒星会造成三角形匹配中内存消耗过大、效率过低，对此可以根据望远镜指向筛选导航星并划取图像中心区域的观测星进行优化三角形匹配，再将得出的底片常数应用序列图像修正算法推广到整帧图像，进行底片常数计算。需要注意的是，这种情况下，常采用十二参数模型进行匹配，可以在一定程度上克服相机图像的畸变。序列图像修正算法流程如 图6 图6 所示。 望远镜码盘给出了其中心指向，但一些原因，如自然条件、热胀冷缩、配重变化、重力变形、结构改变（本例中望远镜可以调节主焦点和卡塞格林焦点）等会造成该指向不准；然而在同一圈次的观测中，指向的相对精度较高，因此配合首帧的修正中心就可以大幅提高底片的中心位置精度。该中心位置对恒星理想坐标运算和目标位置解算有着重要意义，同时还可对望远镜指向偏差提供参考。 4 实验结果 1 本研究采用Windows系统的计算平台，CPU：AMD Ryzen 5 @ 3.4 GHz，内存：8 GB DDR4，编译环境为VS 2019。 将优化算法和传统的三角形匹配算法作比较，对于39颗观测星对应39颗导航星的实验，传统算法耗时2.105 s，而优化算法耗时0.044 s。 将算法应用于精跟测试中轨激光星的事后分析中，观测时间为2019年4月，指向范围：（1 h30 m16.2 s，53°0′50.2″）~（2 h15 m13.0 s， 50°33′40.0″）。首帧识别到326颗观测星，筛选出66颗导航星，首先对所有涉及的背景恒星进行了归算，共消耗时间0.005 s，后采用优化的三角形匹配算法对首帧数据进行匹配。经4次放大，理想坐标匹配误差为0.000 04时，可匹配60对星象，共需时间0.75 s。经最小二乘法拟合出底片和地平取向基本成镜像。匹配结果如 图7 图7 ，其中红圈表示成功匹配的星象（彩图见期刊电子版）。 1,2; 对后续帧采用序列图像修正算法，将计算结果列于 表1 表1 ，可见该算法耗时均在0.04 s以内，在识别到充分多观测星时，可实现绝大多数导航星的匹配。 No. Observed stars Navigation stars Matched pairs Time/s No. Observed stars Navigation star Matched pairs Time/s 1 326 66 60 0.75 16 589 96 92 0.018 2 396 66 63 0.024 17 591 87 86 0.016 3 381 72 69 0.024 18 526 91 90 0.029 4 389 68 63 0.013 19 545 86 85 0.021 5 391 70 68 0.018 20 566 79 77 0.01 6 401 65 62 0.037 21 531 73 70 0.017 7 429 72 69 0.017 22 514 84 81 0.025 8 472 87 84 0.032 23 577 95 94 0.026 9 510 104 104 0.015 24 597 100 97 0.026 10 474 114 110 0.019 25 611 102 96 0.028 11 529 103 99 0.019 26 633 88 82 0.012 12 528 87 85 0.015 27 675 96 93 0.018 13 533 87 86 0.024 28 659 104 104 0.025 14 545 102 101 0.023 29 622 105 102 0.021 15 573 101 97 0.021 30 630 79 76 0.037 匹配星对用于拟合，得出各帧的底片常数。底片常数的缓变是序列星图匹配算法的主要依据，30帧的常数（ 式（3） 式（3） 中 <math id="M29"><mi>a</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226670&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226666&type= 1.86266661 2.79399991 和 <math id="M30"><mi>b</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226679&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226675&type= 1.86266661 2.79399991 ）变化如 图8 图8 所示，可见常数变化较为均匀，可用线性回归来表示，其皮尔森相关系数均达到99.6%以上；并且经30帧，两参数均减小了不超过10%。 由拟合的底片常数计算得出的观测星的天球位置和实际导航星之间的误差可由 图9 图9 表示（彩图见期刊电子版），某一典型帧赤纬方向的误差由绿色球表示，赤经方向由红色方块表示。可见除个别边缘点外，绝大多数的星象误差在6角秒以内。由于精跟型拍摄目标多位于视场中心，其误差会进一步缩小，对于观测的中轨激光星，30帧的目标统计误差如 图10 图10 所示。可以看出，除第一帧以外，绝大部分目标的误差在2″以内，所有帧的平均误差为 <math id="M35"><mover><mrow><mo>∆</mo><mi>α</mi><mo>⋅</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>δ</mi></mrow><mo>¯</mo></mover><mo>=</mo><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">53</mn></mrow><mrow><mi>'</mi><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226722&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226718&type= 28.02466965 3.72533321 ， <math id="M36"><mover><mrow><mo>∆</mo><mi>δ</mi></mrow><mo>¯</mo></mover><mo>=</mo><mo>-</mo><mn mathvariant="normal">0.34</mn><mi>'</mi><mi>'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226731&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226726&type= 22.35199928 3.72533321 。由此可证明本文所述方法的有效性。 5 结 论 1 本文所述的星图匹配程序包括基于降维查表优化的三角形匹配和基于上一帧底片常数的序列图像修正，该程序特别适用于精密跟踪型的大视场望远镜。其中优化三角形匹配算法相较于原始算法大幅提高了匹配速度（39颗观测星对应39颗导航星，时间从2.105 s降到0.044 s），而序列图像底片常数修正进一步加速，使得100颗左右的匹配星对在0.04 s以内完成，这说明了算法的高效性；本方法对中轨激光星定位的结果显示其平均误差为 <math id="M37"><mover><mrow><mo>∆</mo><mi>α</mi><mo>⋅</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mi>δ</mi></mrow><mo>¯</mo></mover><mo>=</mo><mn mathvariant="normal">0</mn><mo>.</mo><msup><mrow><mn mathvariant="normal">53</mn></mrow><mrow><mi>'</mi><mi>'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226762&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226759&type= 28.44799995 3.72533321 ， <math id="M38"><mover><mrow><mo>∆</mo><mi>δ</mi></mrow><mo>¯</mo></mover><mo>=</mo><mo>-</mo><mn mathvariant="normal">0.34</mn><mi>'</mi><mi>'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226770&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37226766&type= 22.52133369 3.72533321 ，完全满足天文定位的精度要求。 参考文献 吴连大 . 人造卫星与空间碎片的轨道和探测 ［M］. 2版 . 北京 ： 中国科学技术出版社 ， 2012 . 人造卫星与空间碎片的轨道和探测 WU L D . The Orbit and Detection of Artificial Satellites and Space Debris ［M］. Beijing ： Science and technology of China Press ， 2012 . The Orbit and Detection of Artificial Satellites and Space Debris MORTARI D ， SAMAAN M A ， BRUCCOLERI C ， et al . The pyramid star identification technique ［J］. Navigation ， 2004 ， 51 （ 3 ）： 171 - 183 . doi: 10.1002/j.2161-4296.2004.tb00349.x http://dx.doi.org/10.1002/j.2161-4296.2004.tb00349.x The pyramid star identification technique LIU H ， WEI X G ， LI J ， et al . A star identification algorithm based on simplest general subgraph ［J］. Acta Astronautica ， 2021 ， 183 ： 11 - 22 . doi: 10.1016/j.actaastro.2021.02.026 http://dx.doi.org/10.1016/j.actaastro.2021.02.026 A star identification algorithm based on simplest general subgraph 翟优 ， 曾峦 ， 熊伟 . 基于不变特征描述符实现星点匹配 ［J］. 光学 精密工程 ， 2012 ， 20 （ 11 ）： 2531 - 2539 . doi: 10.3788/ope.20122011.2531 http://dx.doi.org/10.3788/ope.20122011.2531 基于不变特征描述符实现星点匹配 ZHAI Y ， ZENG L ， XIONG W . Star matching based on invariant feature descriptor ［J］. Opt. Precision Eng. ， 2012 ， 20 （ 11 ）： 2531 - 2539 . （in Chinese） . doi: 10.3788/ope.20122011.2531 http://dx.doi.org/10.3788/ope.20122011.2531 Star matching based on invariant feature descriptor LU W T ， HUA G X ， ZHAO Y F ， et al . Design and optimization of star recognition algorithm based on hierarchical CNN ［J］. Journal of Physics： Conference Series ， 2021 ， 2132 （ 1 ）： 012009 . doi: 10.1088/1742-6596/2132/1/012009 http://dx.doi.org/10.1088/1742-6596/2132/1/012009 Design and optimization of star recognition algorithm based on hierarchical CNN 时圣革 ， 雷肖剑 ， 于长海 . 星图识别三角形算法综述 ［J］. 光电技术应用 ， 2014 ， 29 （ 5 ）： 1 - 6 . doi: 10.3969/j.issn.1673-1255.2014.05.001 http://dx.doi.org/10.3969/j.issn.1673-1255.2014.05.001 星图识别三角形算法综述 SHI SH G ， LEI X J ， YU CH H . Triangle algorithm overview of star pattern recognition ［J］. Electro-Optic Technology Application ， 2014 ， 29 （ 5 ）： 1 - 6 . （in Chinese） . doi: 10.3969/j.issn.1673-1255.2014.05.001 http://dx.doi.org/10.3969/j.issn.1673-1255.2014.05.001 Triangle algorithm overview of star pattern recognition 王军 ， 何昕 ， 魏仲慧 ， 等 . 基于多特征匹配的快速星图识别 ［J］. 光学 精密工程 ， 2019 ， 27 （ 8 ）： 1870 - 1879 . doi: 10.3788/ope.20192708.1870 http://dx.doi.org/10.3788/ope.20192708.1870 基于多特征匹配的快速星图识别 WANG J ， HE X ， WEI Z H ， et al . Fast star identification algorithm based on multi-feature matching ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 8 ）： 1870 - 1879 . （in Chinese） . doi: 10.3788/ope.20192708.1870 http://dx.doi.org/10.3788/ope.20192708.1870 Fast star identification algorithm based on multi-feature matching 刘先一 ， 张志利 ， 周召发 ， 等 . 径向特征下的改进三角形星图识别算法 ［J］. 国防科技大学学报 ， 2020 ， 42 （ 6 ）： 42 - 46 . doi: 10.11887/j.cn.202006005 http://dx.doi.org/10.11887/j.cn.202006005 径向特征下的改进三角形星图识别算法 LIU X Y ， ZHANG ZH L ， ZHOU ZH F ， et al . Modified algorithm of triangle star identification based on radial feature ［J］. Journal of National University of Defense Technology ， 2020 ， 42 （ 6 ）： 42 - 46 . （in Chinese） . doi: 10.11887/j.cn.202006005 http://dx.doi.org/10.11887/j.cn.202006005 Modified algorithm of triangle star identification based on radial feature 李变 ， 屈俐俐 ， 高玉平 ， 等 . 基于冗余参考星的三角形星图匹配方法 ： CN112880666A ［P］. 2021-06-01 . 基于冗余参考星的三角形星图匹配方法 LI B ， QU L L ， GAO Y P ， et al . Triangular star map matching method based on redundant reference star ： CN112880666A ［P］. 2021-06-01 . （in Chinese） Triangular star map matching method based on redundant reference star 徐伟 ， 魏铭 ， 杨秀彬 ， et al . Rapid star map identification method based on triangle matching ： CN110879063B ［P］. 2021-02-26 . （in Chinese） Rapid star map identification method based on triangle matching 徐伟 ， 魏铭 ， 杨秀彬 ， 等 . 基于三角形匹配的快速星图识别方法 ： CN110879063B ［P］. 2021-02-26 . 基于三角形匹配的快速星图识别方法 熊琨 ， 王春喜 ， 姜云翔 ， 等 . 一种三角形匹配式星敏感器星库的构建方法 ： CN107883946B ［P］. 2021-01-05 . 一种三角形匹配式星敏感器星库的构建方法 XIONG K ， WANG CH X ， JIANG Y X ， et al . Construction method of star database of triangle matching-type star sensor ： CN107883946B ［P］. 2021-01-05 . （in Chinese） Construction method of star database of triangle matching-type star sensor PÁL A ， BAKOS G Á . Astrometry in wide-field surveys ［J］. Publications of the Astronomical Society of the Pacific ， 2006 ， 118 （ 848 ）： 1474 - 1483 . doi: 10.1086/508573 http://dx.doi.org/10.1086/508573 Astrometry in wide-field surveys 廖明哲 ， 吴谨 ， 朱磊 . 基于ResNet和RF-Net的遥感影像匹配 ［J］. 液晶与显示 ， 2020 ， 35 （ 9 ）： 972 - 980 . doi: 10.37188/YJYXS20203509.0972 http://dx.doi.org/10.37188/YJYXS20203509.0972 基于ResNet和RF-Net的遥感影像匹配 LIAO M ZH ， WU J ， ZHU L . Remote sensing images matching based on ResNet and RF-Net ［J］. Chinese Journal of Liquid Crystals and Displays ， 2020 ， 35 （ 9 ）： 972 - 980 . （in Chinese） . doi: 10.37188/YJYXS20203509.0972 http://dx.doi.org/10.37188/YJYXS20203509.0972 Remote sensing images matching based on ResNet and RF-Net 兀颖 ， 葛亮 ， 卢晓猛 ， 等 . 基于GPU的小行星光学观测图像实时处理 ［J］. 光学 精密工程 ， 2020 ， 28 （ 11 ）： 2527 - 2534 . doi: 10.37188/OPE.20202811.2527 http://dx.doi.org/10.37188/OPE.20202811.2527 基于GPU的小行星光学观测图像实时处理 WU Y ， GE L ， LU X M ， et al . A method of real-time data reduction of asteroid optical observation image based on GPU ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 11 ）： 2527 - 2534 . （in Chinese） . doi: 10.37188/OPE.20202811.2527 http://dx.doi.org/10.37188/OPE.20202811.2527 A method of real-time data reduction of asteroid optical observation image based on GPU 邓诗宇 . 空间目标光谱测量技术与识别方法应用研究 ［D］. 北京 ： 中国科学院大学 ， 2021 . doi: 10.3788/LOP202158.2230001 http://dx.doi.org/10.3788/LOP202158.2230001 空间目标光谱测量技术与识别方法应用研究 DENG SH Y . Research on Application of Spectral Measurement Technology and Recognition Method of Space Target ［D］. Beijing ： University of the Chinese Academy of Sciences ， 2021 . （in Chinese） . doi: 10.3788/LOP202158.2230001 http://dx.doi.org/10.3788/LOP202158.2230001 Research on Application of Spectral Measurement Technology and Recognition Method of Space Target BERTIN E ， ARNOUTS S . SExtractor： Software for source extraction ［J］. Astronomy and Astrophysics Supplement Series ， 1996 ， 117 （ 2 ）： 393 - 404 . doi: 10.1051/aas:1996164 http://dx.doi.org/10.1051/aas:1996164 SExtractor： Software for source extraction Software Routines from the IAU SOFA Collection were used ［EB/OL］. Copyright © International Astronomical Union Standards of Fundamental Astronomy ［ 2021-05-12 ］. http：//www.iausofa.org http://www.iausofa.org . Software Routines from the IAU SOFA Collection were used 李振伟 ， 张涛 ， 张楠 ， 等 . 暗弱空间目标的高精度定位 ［J］. 光学 精密工程 ， 2015 ， 23 （ 9 ）： 2627 - 2634 . doi: 10.3788/ope.20152309.2627 http://dx.doi.org/10.3788/ope.20152309.2627 暗弱空间目标的高精度定位 LI ZH W ， ZHANG T ， ZHANG N ， et al . High precision orientation of faint space objects ［J］. Opt. Precision Eng. ， 2015 ， 23 （ 9 ）： 2627 - 2634 . （in Chinese） . doi: 10.3788/ope.20152309.2627 http://dx.doi.org/10.3788/ope.20152309.2627 High precision orientation of faint space objects 10.37188/OPE.20223022.2952 P123.2 P123.3 A 1004-924X（2022）22-2952-10 刘德龙,杨文波,柳鸣等.精密跟踪型望远镜适配的快速星图匹配[J].光学精密工程,2022,30(22):2952-2961. 刘德龙,杨文波,柳鸣等.精密跟踪型望远镜适配的快速星图匹配[J].光学精密工程,2022,30(22):2952-2961. LIU Delong,YANG Wenbo,LIU Ming,et al.Rapid star pattern matching for precisely tracking telescopes[J].Optics and Precision Engineering,2022,30(22):2952-2961. LIU Delong,YANG Wenbo,LIU Ming,et al.Rapid star pattern matching for precisely tracking telescopes[J].Optics and Precision Engineering,2022,30(22):2952-2961. 国家自然科学基金项目资助（No.12003052） 2022-05-17 2022-07-08 2022-11-25 2022-11-24 光学精密工程 22 30