一种基于光学测角数据的星链卫星初轨确定方法和装置与流程

1.本发明涉及天文测量的技术领域，尤其是涉及一种基于光学测角数据的星链卫星初轨确定方法和装置。


背景技术：

2.当前针对低轨卫星的光学数据的初轨确定常用方法，主要是对“方位俯仰”测角数据进行的初轨确定，针对“赤经赤纬”数据定初轨的方法相对较少，而且初轨确定精度与观测时间相关；在实际应用的低轨卫星（轨道高度在200~600公里）短弧段光学数据定初轨方法中，计算得到的误差比较大，甚至无法获取轨道根数，而根据两个点定初轨的方法更是凤毛麟角，因为测量数据太少，往往无法定出结果。
3.另外，目前还没有专门针对星链卫星的光学数据进行初轨确定的方法。
4.星链，是美国太空探索技术公司的一个项目，太空探索技术公司计划在2019年至2024年间在太空搭建由约1.2万颗卫星组成的“星链”网络提供互联网服务，其中1584颗将部署在地球上空550千米处的近地轨道，并从2020年开始工作。
5.针对上述问题，还未提出有效的解决方案。


技术实现要素：

6.有鉴于此，本发明的目的在于提供一种基于光学测角数据的星链卫星初轨确定方法和装置，以缓解了现有的初轨确定方法的精度和效率较低的技术问题。
7.第一方面，本发明实施例提供了一种基于光学测角数据的星链卫星初轨确定方法，包括： 获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道；基于所述最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨。
8.进一步地，所述光学测角数据包括：测量时间，赤经和赤纬；所述测站的位置信息包括：经度，纬度和高度。
9.进一步地，基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距，包括：第一计算步骤，基于所述预设观测距离粗值，以及初始时刻的赤经和赤纬，计算出所述星链卫星在j2000惯性坐标系中的初始位置信息；第一确定步骤，基于所述初始时刻的测站的位置信息，确定出所述测站在地固坐标系中的位置信息，并将所述测站在地固坐标系中的位置信息转换为所述测站在j2000惯性坐标
系中的位置信息；第二计算步骤，基于所述星链卫星在j2000惯性坐标系中的初始位置信息和所述测站在j2000惯性坐标系中的位置信息，计算出所述目标轨道的初始地心距；若所述初始地心距不处于预设范围内，则调整所述预设观测距离粗值，并重复执行所述第一计算步骤，第一所述确定步骤和所述第二计算步骤，直至所述初始地心距处于所述预设范围内，并将处于所述预设范围内的初始地心距确定为所述目标轨道的地心距。
10.进一步地，所述标准卫星轨道的轨道参数包括：半长轴，倾角，偏心率，升交点赤经，近地点幅角，平近点角和第一初始历元，其中，所述半长轴为所述目标轨道的地心距，所述倾角为53
°
，所述偏心率、所述升交点赤经、所述近地点幅角和所述平近点角均为0，所述第一初始历元为初始时刻；所述预设轨道的位置速度包括：所述测站在j2000惯性坐标系中的位置，所述星链卫星在所述标准卫星轨道中初始时刻的速度和所述第一初始历元。
11.进一步地，基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道 ，包括：第一修正步骤，基于所述标准卫星轨道在初始时刻的星下点位置信息和所述预设轨道在初始时刻的星下点位置信息，计算出第一平近点角，并基于所述第一平近点角对所述标准卫星轨道进行修正，得到第一中间标准卫星轨道；第二修正步骤，基于所述第一中间标准卫星轨道中初始时刻时所述星链卫星的经度和所述预设轨道中初始时刻时所述星链卫星的经度，计算出第一升交点赤经，并基于所述第一升交点赤经对所述第一中间标准卫星轨道进行修正，得到所述目标标准卫星轨道。
12.进一步地，基于最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨，包括：第二确定步骤，确定出所述目标标准卫星轨道中最后时刻时所述星链卫星的目标赤经和目标赤纬；执行步骤，若所述目标赤经和所述目标赤纬与所述光学测角数据中最后时刻对应的赤经和赤纬之间的差值大于门限值，则计算所述目标标准卫星轨道在最后时刻的可见性距离，并基于所述可见性距离和所述目标赤经和所述目标赤纬，计算出所述星链卫星在最后时刻的位置信息和星下点，以及基于所述星链卫星在最后时刻的位置信息和星下点，调整所述目标标准卫星轨道的半长轴，得到最终标准卫星轨道；第三修正步骤，基于所述最终标准卫星轨道在最后时刻的星下点位置信息和所述预设轨道在最后时刻的星下点位置信息，计算出第二平近点角，并基于所述第二平近点角对所述最终标准卫星轨道进行修正，得到第二中间标准卫星轨道；第四修正步骤，基于所述第二中间标准卫星轨道中最后时刻时所述星链卫星的经度和所述预设轨道中初始时刻时所述星链卫星的经度，计算出第二升交点赤经，并基于所述第一升交点赤经对所述第二中间标准卫星轨道进行修正，得到所述第三中间标准卫星轨道；第三确定步骤，确定出所述第三中间标准卫星轨道中初始时刻时所述星链卫星的最终赤经和最终赤纬；若所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值大于门限值，则计算所述第三中间标准卫星轨道在初始时刻的可见性距离，并基于所述初始时刻的可见性距离和所述初始时刻对应的赤经和赤纬，确定出目标星下点，将所述第三中间标准卫星轨道确定为所述标准卫星轨道，以及将所述目标星下点确定为所述预设轨道的星下点，重复执行所述第一修正步骤，所述第二修正步骤，所述第二确定步骤，所述执行步骤，所述第三修正步骤，所述第四修正步骤和所述第三确定步骤，直至所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值小于或等于门限值，并
将所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值小于或等于门限值对应的第三中间标准卫星轨道确定为所述星链卫星的初轨。
13.第二方面，本发明实施例还提供了一种基于光学测角数据的星链卫星初轨确定装置，包括：获取单元，确定单元，构建单元，修正单元和优化单元，其中，所述获取单元，用于获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；所述确定单元，用于基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；所述构建单元，用于基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；所述修正单元，用于基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道；所述优化单元，用于基于所述最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨。
14.进一步地，所述光学测角数据包括：测量时间，赤经和赤纬；所述测站的位置信息包括：经度，纬度和高度。
15.第三方面，本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述第一方面中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
16.第四方面，本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序。
17.在本发明实施例中，通过获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道；基于最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨，达到了利用光学测角数据确定星链卫星初轨的目的，进而解决了现有的初轨确定方法的精度和效率较低的技术问题，从而实现了提高确定星链卫星初轨的精度和效率的技术效果。
18.本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
19.为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
20.为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的
附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
21.图1为本发明实施例提供的一种基于光学测角数据的星链卫星初轨确定方法的流程图；图2为本发明实施例提供的一种基于光学测角数据的星链卫星初轨确定装置的示意图；图3为本发明实施例提供的一种电子设备的示意图。
具体实施方式
22.为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
23.实施例一：根据本发明实施例，提供了一种基于光学测角数据的星链卫星初轨确定方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
24.图1是根据本发明实施例的一种基于光学测角数据的星链卫星初轨确定方法的流程图，如图1所示，该方法包括如下步骤：步骤s102，获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；需要说明的是，上述的光学测角数据包括：测量时间，赤经和赤纬；上述的测站的位置信息包括：经度，纬度和高度，其中，一个测量时间对应一个赤经和一个赤纬。
25.步骤s104，基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；步骤s106，基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；步骤s108，基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道；步骤s110，基于所述最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨。
26.在本发明实施例中，通过获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨
道进行修正，得到目标标准卫星轨道；基于最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨，达到了利用光学测角数据确定星链卫星初轨的目的，进而解决了现有的初轨确定方法的精度和效率较低的技术问题，从而实现了提高确定星链卫星初轨的精度和效率的技术效果。
27.在本发明实施例中，步骤s104包括如下步骤：第一计算步骤，基于所述预设观测距离粗值，以及初始时刻的赤经和赤纬，计算出所述星链卫星在j2000惯性坐标系中的初始位置信息；第一确定步骤，基于所述初始时刻的测站的位置信息，确定出所述测站在地固坐标系中的位置信息，并将所述测站在地固坐标系中的位置信息转换为所述测站在j2000惯性坐标系中的位置信息；第二计算步骤，基于所述星链卫星在j2000惯性坐标系中的初始位置信息和所述测站在j2000惯性坐标系中的位置信息，计算出所述目标轨道的初始地心距；若所述初始地心距不处于预设范围内，则调整所述预设观测距离粗值，并重复执行所述第一计算步骤，第一所述确定步骤和所述第二计算步骤，直至所述初始地心距处于所述预设范围内，并将处于所述预设范围内的初始地心距确定为所述目标轨道的地心距。
28.在本发明实施例中，首先给定一个预测观测距离粗值，需要说明的是观测距离用于表征测站与星链卫星之间的距离，根据预测观测距离粗值、初始时刻赤经（α）、赤纬（δ），星链卫星在j2000惯性坐标系中的初始位置信息，具体计算公式如下：; 。
29.接着，基于初始时刻的测站的位置信息，确定出测站在地固坐标系中的位置信息，并将测站在地固坐标系中的位置信息转换为测站在j2000惯性坐标系中的位置信息()。
30.然后，基于星链卫星在j2000惯性坐标系中的初始位置信息和测站在j2000惯性坐标系中的位置信息，计算出目标轨道的初始地心距，其中，初始地心距r=，()= + ()。
31.最后，确定初始地心距不处于预设范围内，则调整预设观测距离粗值，并重复执行第一计算步骤，第一确定步骤和第二计算步骤，直至初始地心距处于预设范围内，并将处于预设范围内的初始地心距确定为目标轨道的地心距。
32.需要说明的是预设范围的上限值为地球赤道半径+轨道高度最大值，预设范围的下限值为地球赤道半径-轨道高度最小值。
33.在本发明实施例中，在获取到测站的地心距之后，根据测站的地心距和目标数据构建标准卫星轨道和预设轨道。
34.标准卫星轨道的轨道参数包括：半长轴，倾角，偏心率，升交点赤经，近地点幅角，平近点角和第一初始历元，其中，所述半长轴为所述目标轨道的地心距，所述倾角为53
°
，所述偏心率、所述升交点赤经、所述近地点幅角和所述平近点角均为0，所述第一初始历元为所述初始时刻。
35.预设轨道的轨道参数包括：所述测站在j2000惯性坐标系中的位置，所述星链卫星在所述标准卫星轨道中初始时刻的速度和所述第一初始历元。
36.在本发明实施例中，步骤s108包括如下步骤：第一修正步骤，基于所述标准卫星轨道在初始时刻的星下点位置信息和所述预设轨道在初始时刻的星下点位置信息，计算出第一平近点角，并基于所述第一平近点角对所述标准卫星轨道进行修正，得到第一中间标准卫星轨道；第二修正步骤，基于所述第一中间标准卫星轨道中初始时刻时所述星链卫星的经度和所述预设轨道中初始时刻时所述星链卫星的经度，计算出第一升交点赤经，并基于所述第一升交点赤经对所述第一中间标准卫星轨道进行修正，得到所述目标标准卫星轨道。
37.在本发明实施例中，首先，将标准卫星轨道和预设轨道由j2000惯性系转换为地固坐标系，然后，将地固坐标系下的准卫星轨道和预设轨道的位置转换为星下点（即，用于经纬度表示）。
38.接着，根据星下点，寻找标准卫星轨道中纬度等于预设轨道初始时刻的纬度的时刻，计算j2000位置，结合初始时刻j2000星链卫星的位置，计算第一平近点角，修正标准卫星轨道的得到第一中间标准卫星轨道，其中，第一平近点角m1计算过程如下：由于，原始数据太短，因此需要先对标准卫星轨道进行二体预报，生成j2000位置速度；接着，重新将标准卫星轨道j2000位置速度转换为星下点。
39.然后，寻找标准卫星轨道的经纬度报表中纬度与预设轨道的初始时刻纬度相同的时刻，计算该时刻标准轨道的j2000位置速度；根据位置速度计算平近点角m。
40.其中，为标准卫星轨道在初始时刻的j2000位置；为标准卫星轨道在纬度与预设轨道初始时刻纬度相同时刻的j2000位置。
41.在得到第一中间标准卫星轨道之后，利用第一中间标准卫星轨道中初始时刻时星链卫星的经度和预设轨道中初始时刻时星链卫星的经度，计算出第一升交点赤经，并基于第一升交点赤经对第一中间标准卫星轨道进行修正，得到目标标准卫星轨道，其中，第一升交点赤经raan=预设轨道的初始时刻经度-第一中间标准卫星轨道的初始时刻经度。
42.需要说明的是，目标标准卫星轨道为利用光学测角数据中的初始时刻的赤经和赤纬进行星链卫星的定初轨之后得到的卫星轨道。
43.在本发明实施例中，步骤s110包括如下步骤：第二确定步骤，确定出所述目标标准卫星轨道中最后时刻时所述星链卫星的目标赤经和目标赤纬；执行步骤，若所述目标赤经和所述目标赤纬与所述光学测角数据中最后时刻对应的赤经和赤纬之间的差值大于门限值，则计算所述目标标准卫星轨道在最后时刻的可见性距离，并基于所述可见性距离和所述目标赤经和所述目标赤纬，计算出所述星链卫星在最后时刻的位置信息和星下点，以及基于所述星链卫星在最后时刻的位置信息和星下点，调整所述目标标准卫星轨道的半长轴，得到最终标准卫星轨道；第三修正步骤，基于所述最终标准卫星轨道在最后时刻的星下点位置信息和所述预设轨道在最后时刻的星下点位置信息，计算出第二平近点角，并基于所述第二平近点角
对所述最终标准卫星轨道进行修正，得到第二中间标准卫星轨道；第四修正步骤，基于所述第二中间标准卫星轨道中最后时刻时所述星链卫星的经度和所述预设轨道中初始时刻时所述星链卫星的经度，计算出第二升交点赤经，并基于所述第一升交点赤经对所述第二中间标准卫星轨道进行修正，得到所述第三中间标准卫星轨道；第三确定步骤，确定出所述第三中间标准卫星轨道中初始时刻时所述星链卫星的最终赤经和最终赤纬；若所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值大于门限值，则计算所述第三中间标准卫星轨道在初始时刻的可见性距离，并基于所述初始时刻的可见性距离和所述初始时刻对应的赤经和赤纬，确定出目标星下点，将所述第三中间标准卫星轨道确定为所述标准卫星轨道，以及将所述目标星下点确定为所述预设轨道的星下点，重复执行所述第一修正步骤，所述第二修正步骤，所述第二确定步骤，所述执行步骤，所述第三修正步骤，所述第四修正步骤和所述第三确定步骤，直至所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值小于或等于门限值，并将所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值小于或等于门限值对应的第三中间标准卫星轨道确定为所述星链卫星的初轨。
44.在本发明实施例中，在得到目标标准卫星轨道之后，对目标标准卫星轨道行轨道预报，计算目标标准卫星轨道最后时刻的目标赤经和目标赤纬，与光学测角数据中最后时刻的赤经和赤纬对比，若所述目标赤经和所述目标赤纬与所述光学测角数据中最后时刻对应的赤经和赤纬之间的差值小于或等于门限值，则星链卫星的初轨确定结束，目标标准卫星轨道为星链卫星的初轨。
45.需要说明的是，目标赤经和目标赤纬的计算过程如下：首先，计算最后时刻测站的在j2000惯性坐标系中的位置；然后，计算目标标准卫星轨道相对于测站的j2000位置速度（x，y，z；v
x
，vy，vz）= 目标标准卫星轨道在最后时刻j2000惯性坐标系中的位置
ꢀ‑
最后时刻测站在j2000惯性坐标系中的位置;计算目标赤经和目标赤纬：； ()。
46.若所述最终赤经和所述最终赤纬与所述光学测角数据中初始时刻对应的赤经和赤纬之间的差值大于门限值，则对目标标准卫星轨道进行迭代优化，具体过程如下：使用测站对星链卫星进行二体预报，计算目标标准卫星轨道的可见性，得到最后时刻的可见性距离；根据可见性距离和目标赤经，目标赤纬，计算目标标准卫星轨道在最后时刻的星下点（经纬度）：1）先计算目标标准卫星轨道在j2000惯性坐标下的位置：; 2）再计算目标标准卫星轨道相对于测站的j2000惯性坐标系中的位置(
)=轨道j2000下的位置+测站j2000位置;3）将目标标准卫星轨道相对于测站的j2000惯性坐标系中的位置转换为地固坐标系中的位置；4）根据地固坐标系中的位置转换位星下点（经纬度）。
47.根据光学测量数据中最后时刻赤经赤纬，计算位置和半长轴过程：使用计算的最后时刻的可见性距离，光学测角数据中最后时刻赤经赤纬，计算星下点（经纬度）。
48.1）先计算目标标准卫星轨道j2000惯性坐标系下的位置：; 2）再计算目标标准卫星轨道相对于测站的j2000惯性坐标系中的位置()= 目标标准卫星轨道在j2000惯性坐标系中的位置+测站在j2000惯性坐标系中的位置;3）将目标标准卫星轨道相对于测站在j2000惯性坐标系中的位置转换为地固坐标系中的位置；4）根据地固坐标系中的位置转换为星下点（经纬度）。
49.最后，计算调整的半长轴，并对目标标准卫星轨道进行调整，得到最终标准卫星轨道；。
50.然后，寻找最终标准卫星轨道下点中纬度与光学测角数据中最后时刻纬度相同的点，计算第二升交点赤经，修正最终标准卫星轨道，得到第二中间标准卫星轨道；调整的第二升交点赤经计算过程：第二升交点赤经=第一升交点赤经+（光学测角数据中最后时刻的经度-目标标准卫星轨道的最后时刻经度）。
51.根据第二中间标准卫星轨道最后时刻在j2000惯性坐标系中的位置和光学测角数据中最后时刻星链卫星在j2000惯性坐标系中的位置，计算第二平近点角，修正第二中间标准卫星轨道，得到第三中间标准卫星轨道。
52.第二平近点角m计算过程：方法：计算轨道的最后时刻与原始数据最后时刻对应的位置夹角:cosm=《m0,m1》其中，为第二中间标准卫星轨道在最后时刻的j2000位置；为光学测角数据中最后时刻赤经赤纬转换的j2000位置。
53.第三中间标准卫星轨道为利用最后时刻光学测角数据优化后得到的初轨。
54.接着，将第三中间标准卫星轨道转换为初始历元的轨道根数，计算赤经赤纬，与光学测角数据中初始时刻赤经赤纬对比，判断是否小于门限值，若小于，则定初轨结束，即，第三中间标准卫星轨道为星链卫星的初轨。
55.具体的，使用二体轨道预报器计算三中间标准卫星轨道的初始时刻的轨道根数；
将初始时刻第三中间标准卫星轨道六根数转为三中间标准卫星轨道j2000位置速度，并计算测站在j2000惯性坐标系中的位置；计算第三中间标准卫星轨道相对于测站的在j2000关系坐标系中的位置速度=轨道j2000位置速度-测站j2000位置；计算最终赤经和最终赤纬：最终；最终()。
56.若仍大于，重新将第三中间标准卫星轨道设为标准轨道，计算第三中间标准卫星轨道在初始时刻的可见性，生成新的可见性距离，结合光学测角数据中初始时刻赤经和赤纬计算生成的目标星下点，并将目标星下点作为新的预设轨道的星下点，重复执行第一修正步骤第二修正步骤，第二确定步骤，执行步骤，第三修正步骤，第四修正步骤和第三确定步骤，直至最终赤经和最终赤纬与光学测角数据中初始时刻对应的赤经和赤纬之间的差值小于或等于门限值，并将最终赤经和最终赤纬与光学测角数据中初始时刻对应的赤经和赤纬之间的差值小于或等于门限值对应的第三中间标准卫星轨道确定为星链卫星的初轨。
57.实施例二：本发明实施例还提供了一种基于光学测角数据的星链卫星初轨确定装置，该基于光学测角数据的星链卫星初轨确定装置用于执行本发明实施例上述内容所提供的基于光学测角数据的星链卫星初轨确定方法，以下是本发明实施例提供的装置的具体介绍。
58.如图2所示，图2为上述基于光学测角数据的星链卫星初轨确定装置的示意图，该基于光学测角数据的星链卫星初轨确定装置包括：获取单元10，确定单元20，构建单元30，修正单元40和优化单元50。
59.所述获取单元，用于获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；所述确定单元，用于基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；所述构建单元，用于基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；所述修正单元，用于基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道；所述优化单元，用于基于最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨。
60.在本发明实施例中，通过获取星链卫星对应的目标数据，其中，所述目标数据包括：初始时刻的光学测角数据，最终时刻的光学测角数据和测站的位置信息；基于所述初始时刻的光学测角数据，所述测站的位置信息和预设观测距离粗值，确定出所述目标轨道的地心距；基于所述地心距和所述目标数据，构建标准卫星轨道和所述星链卫星的预设轨道，并计算出所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息；基于所述标准卫星轨道的星下点位置信息和所述预设轨道的星下点位置信息，对所述标准卫星轨道进行修正，得到目标标准卫星轨道；基于最终时刻的光学测角数据对所述目标标准卫星轨道进行迭代优化，得到所述星链卫星的初轨，达到了利用光学测角数据确定星链卫星初轨的目的，进而解决了现有的初轨确定方法的精度和效率较低的技术问题，从而实现了提
高确定星链卫星初轨的精度和效率的技术效果。
61.优选的，所述光学测角数据包括：测量时间，赤经和赤纬；所述测站的位置信息包括：经度，纬度和高度。
62.实施例三：本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述实施例一中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
63.参见图3，本发明实施例还提供一种电子设备100，包括：处理器60，存储器61，总线62和通信接口63，所述处理器60、通信接口63和存储器61通过总线62连接；处理器60用于执行存储器61中存储的可执行模块，例如计算机程序。
64.其中，存储器61可能包含高速随机存取存储器（ram，random access memory），也可能还包括非不稳定的存储器（non-volatile memory），例如至少一个磁盘存储器。通过至少一个通信接口63（可以是有线或者无线）实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。
65.总线62可以是isa总线、pci总线或eisa总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图3中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
66.其中，存储器61用于存储程序，所述处理器60在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器60中，或者由处理器60实现。
67.处理器60可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器60中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器60可以是通用处理器，包括中央处理器(central processing unit，简称cpu)、网络处理器(network processor，简称np)等；还可以是数字信号处理器(digital signal processing，简称dsp)、专用集成电路(application specific integrated circuit，简称asic)、现成可编程门阵列(field-programmable gate array，简称fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器61，处理器60读取存储器61中的信息，结合其硬件完成上述方法的步骤。
68.实施例四：本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器运行时执行上述实施例一中所述方法的步骤。
69.另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是
两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
70.在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
71.在本技术所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
72.所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
73.另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
74.最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。