CN111156991A - 基于指向误差自动测定的空间碎片实时天文定位方法 - Google Patents

Abstract

本发明公开了一种基于指向误差自动测定的空间碎片实时天文定位方法，包括：恒星及空间碎片星象；恒星检索；理论星图生成；实测星图生成；理论星图及实测星图匹配；指向及像面旋转测定；底片模型优选；空间碎片天文定位。本发明能够根据观测视场大小自动优选底片模型，按照图像上给定时间和图像中心指向，自动测定图像中心指向及像面旋转角，实现恒星理论坐标和实测坐标的自动匹配，从而实现空间碎片的实时天文定位。

Description

基于指向误差自动测定的空间碎片实时天文定位方法

技术领域

本发明涉及空间碎片定位技术领域，具体而言涉及一种基于指向误差自动测定的空间碎片实时天文定位方法。

背景技术

在科研、军事等许多领域，都需要对空间碎片进行监视，一方面测定空间碎片的每一个观测时刻在天空中的位置及其变化，确定空间碎片的运行轨道，从而获取空间碎片精确的信息。基于此需求，空间碎片的精确测量是非常重要的基础环节，没有空间碎片的精确测量，空间碎片轨道识别，编目定轨，及精密定轨都无法实现。

空间碎片的位置精确测量有两种方法：绝对定位和相对定位。其中绝对测定位就是利用望远镜的轴系实现空间碎片测量，它受到望远镜轴系加工精度、大气折射修正精度、温度变形等因素的影响，它不依赖背景恒星的位置。相对定位是根据空间碎片和背景恒星的相对位置实现空间碎片的测量，望远镜指向精度不直接影响测量结果，但是当望远镜指向及像面安装误差大的情况下，就会造成恒星在图像上理论坐标和恒星在图像上实测坐标相差较大，尤其对于有像面误差的图像，边缘部分误差更大，无法满足给定匹配门限，因此造成恒星理论星图和实测星图匹配失败，无法实现相对定位。通常为了实现空间碎片天文定位，必须实事先完成以下三件事：

(1)人工测量获得像面旋转角。

(2)通过观测几十颗恒星获得望远镜指向修正模型。

(3)优选底片常数模型。

上述的第二件事需要经常做，而且需要较为精确的测站温度、相对湿度、及大气压强。即使如此，在仰角较低的天区，受大气折射改正精度影响，望远镜指向误差修正精度偏低。

在此基础上，现有空间碎片天文定位方法的计算步骤如下：

(1)根据望远镜指向修正模型，获得图像中心指向误差。

(2)通过时间和修正后的指向，通过恒星检索获得恒星在图像上理论坐标。

(3)按照给定匹配门限实现理论坐标和实测坐标的匹配。

(4)按照优选的底片常数模型实现空间碎片实时天文定位。

基于现有天文定位方法的不足，本发明给出一种基于指向误差自动测定的空间目标实时天文定位方法，它能够根据观测视场大小自动优选底片模型，按照图像上给定时间和图像中心指向，自动测定图像中心指向及像面旋转角，实现恒星理论坐标和实测坐标的自动匹配，从而实现空间碎片的实时天文定位。

发明内容

本发明目的在于提供一种基于指向误差自动测定的空间碎片实时天文定位方法，能够根据观测视场大小自动优选底片模型，按照图像上给定时间和图像中心指向，自动测定图像中心指向及像面旋转角，实现恒星理论坐标和实测坐标的自动匹配，从而实现空间碎片的实时天文定位。对于固定站址(有精密的天文经纬度)的望远镜，该方法降低了望远镜轴系加工精度要求，降低望远镜外场安装调试要求，减少了测站环境温度参输入要求，减少了观测前需要指向标校的要求，降低了图像时间标识和指向标识因不严格对应造成的指向误差。对于可移动望远镜，该方法在无精密天文经纬度不能实现望远镜指向标校的情况下，也能实现空间碎片实时天文定位。

为达成上述目的，结合图1，本发明提出一种基于指向误差自动测定的空间碎片实时天文定位方法，所述实时天文定位方法包括以下步骤：

S1：接收至少一帧包括空间碎片和背景恒星的图像，获得图像上在预设检测门限内的恒星和空间碎片的星象信息，所述星象信息包括每个恒星和空间碎片的二维平面坐标(x,y)、像素个数、灰度和，按照像素个数减少的顺序对所有星象进行排序；其中，设图像的左上角为坐标原点(0，0)，图像右侧为x轴增加方向，图像下侧为y轴增加方向，x为图像中星象所在位置在x轴方向上与坐标原点的距离，y为图像中星象所在位置在y轴方向上与坐标原点的距离；

S2：结合天文定位恒星星库，根据图像对应的拍摄信息和图像中心含有指向误差的指向信息(α_p,δ_p)，检索出视场中满足给定星等门限的所有恒星的相关信息，所述满足给定星等门限的恒星的相关信息包括其所对应的二维平面坐标理论值(X,Y)、赤经和赤纬理论值(α^s,δ^s)、理想坐标理论值(ξ^s,ζ^s)、理论星等M，按照理论星等由小到大的顺序对检索出的恒星进行排序；

S3：按照给定门限，从检索出的恒星中选择N₁颗恒星星象，定义成第一候选恒星，结合第一候选恒星在图像上的二维平面坐标和望远镜的焦距，计算得到任意两颗第一候选恒星之间的角距，选择三颗第一候选恒星组成三角形星图，生成恒星实测星图；

S4：按照给定门限，从检索出的恒星中选择N₂颗恒星星象，定义成第二候选恒星，结合第二候选恒星的二维平面理论坐标、赤经和纬经，以及望远镜的焦距，计算得到任意两颗第二候选恒星之间的角距，选择三颗第二候选恒星组成三角形星图，生成恒星理论星图；

S5：结合实际测得的二维平面坐标，计算得到恒星实测星图上任意三颗定标星i,j,k组成的三角形的边长为

角距大小顺序为

按照预设顺序对计算结果进行排序以生成定标星的实测星图特征库；

结合二维平面理论坐标，计算得到恒星理论星图上任意三颗定标星l,m,n组成的三角形的边长为

角距大小顺序为

按照预设顺序对计算结果进行排序生成定标星的理论星图特征库；

根据下述匹配条件，结合定标星的实测星图特征库和理论星图特征库，对理论星图和实测星图中的恒星进行匹配，直至匹配成功N₃颗恒星：

第j颗定标星和第k颗定标星之间满足下式：

任意三颗i,j,k定标星之间满足下式：

其中，ε₁和ε₂均为预设的角距门限；

S6：设匹配成功的N₃颗恒星在图像上的二维平面坐标为(x_i,y_i)，i＝1,2,…N₃，对应的理论二维平面坐标坐标为(X_i,Y_i)，i＝1,2,…N₃；

利用下式，采用最小二乘方法，计算出系数a,b,c,d,e,f，从而得到中心指向偏差及像面旋转角：

S7：根据图像对应视场大小，结合匹配成功的N₃颗恒星在图像上的二维平面坐标(x_i,y_i)，及理想坐标(ξ_i,ζ_i)，i＝1,2,…N₃，进行常数模型计算，根据恒星的定位精度，自动优选底片常数模型。

进一步的实施例中，步骤S2中，所述理想坐标(ξ^s,ζ^s)满足以下公式：

进一步的实施例中，步骤S3中，采用下述公式计算任意两颗第一候选恒星之间的角距：

其中，f是望远镜的焦距，(x_u,y_u)是第u颗第一候选恒星的二维平面坐标，(x_v,y_v)是第v颗第一候选恒星的二维平面坐标。

进一步的实施例中，步骤S4中，采用下述公式计算任意两颗第二候选恒星之间的角距：

其中，f是望远镜的焦距，(X_u,Y_u)是第u颗第二候选恒星的二维平面理论坐标，(Y_v,Y_v)是第v颗第二候选恒星的二维平面理论坐标。

进一步的实施例中，步骤S4中，采用下述公式计算任意两颗第二候选恒星之间的角距：

其中，f是望远镜的焦距，(α_u,δ_u)是第u颗第二候选恒星的赤经和纬经，(α_v,δ_v)是第v颗第二候选恒星的赤经和赤纬。

进一步的实施例中，步骤S2中，所述图像对应的拍摄信息包括图像的拍摄时间、指向信息、测站经纬度、测站海拔高度、测站温度、测站湿度、大气压强、给定视场大小。

进一步的实施例中，所述天文定位恒星星库的生成方法包括：

将给定星等的全天区恒星按照赤经增加及赤纬增加的顺序，分区存放，并形成索引，生成天文定位恒星星库及索引数据。

进一步的实施例中，步骤S7中，结合匹配成功的定标星数量，分别选用六常数模型、十二常数模型、十四常数模型进行常数模型计算，其中：

所述六常数模型对应至少3个以上的定标星：

所述十二常数模型对应至少6个以上的定标星：

所述十四常数模型对应至少7个以上的定标星：

进一步的实施例中，步骤S7还包括：

自动存储优选底片常数模型。

进一步的实施例中，所述实时天文定位方法还包括：

S8：根据空间碎片的二维平面坐标实测值(x^T,y^T)，采用以下公式获得空间碎片的赤经和赤纬(α^T,δ^T)：

其中，(ξ^T,ζ^T)为空间碎片的理想坐标，由(x^T,y^T)代入六常数、十二常数或者十四常数模型获得。

以上本发明的技术方案，与现有相比，其显著的有益效果在于：

能够根据观测视场大小自动优选底片模型，按照图像上给定时间和图像中心指向，自动测定图像中心指向及像面旋转角，实现恒星理论坐标和实测坐标的自动匹配，从而实现空间碎片的实时天文定位。

对于固定站址(有精密的天文经纬度)的望远镜，该方法降低了望远镜轴系加工精度要求，降低望远镜外场安装调试要求，减少了测站环境温度参输入要求，减少了观测前需要指向标校的要求，降低了图像时间标识和指向标识因不严格对应造成的指向误差。

对于可移动望远镜，该方法在无精密天文经纬度不能实现望远镜指向标校的情况下，也能实现空间碎片实时天文定位。

应当理解，前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外，所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。

结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见，或通过根据本发明教导的具体实施方式的实践中得知。

附图说明

附图不意在按比例绘制。在附图中，在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见，在每个图中，并非每个组成部分均被标记。现在，将通过例子并参考附图来描述本发明的各个方面的实施例，其中：

图1是本发明的基于指向误差自动测定的空间碎片实时天文定位方法的流程图。

具体实施方式

为了更了解本发明的技术内容，特举具体实施例并配合所附图式说明如下。

结合图1，本发明提及一种基于指向误差自动测定的空间碎片实时天文定位方法，所述实时天文定位方法包括以下步骤：

S1：接收至少一帧包括空间碎片和背景恒星的图像，获得图像上在预设检测门限内的恒星和空间碎片的星象信息，所述星象信息包括每个恒星和空间碎片的二维平面坐标(x,y)、像素个数、灰度和，按照像素个数减少的顺序对所有星象进行排序；其中，设图像的左上角为坐标原点(0，0)，图像右侧为x轴增加方向，图像下侧为y轴增加方向，x为图像中星象所在位置在x轴方向上与坐标原点的距离，y为图像中星象所在位置在y轴方向上与坐标原点的距离。

S2：结合天文定位恒星星库，根据图像对应的拍摄信息和图像中心含有指向误差的指向信息(α_p,δ_p)，检索出视场中满足给定星等门限的所有恒星的相关信息，所述满足给定星等门限的恒星的相关信息包括其所对应的二维平面坐标理论值(X,Y)、赤经和赤纬理论值(α^s,δ^s)、理想坐标理论值(ξ^s,ζ^s)、理论星等M，按照理论星等由小到大的顺序对检索出的恒星进行排序。

S3：按照给定门限，从检索出的恒星中选择N₁颗恒星星象，定义成第一候选恒星，结合第一候选恒星在图像上的二维平面坐标和望远镜的焦距，计算得到任意两颗第一候选恒星之间的角距，选择三颗第一候选恒星组成三角形星图，生成恒星实测星图。

S4：按照给定门限，从检索出的恒星中选择N₂颗恒星星象，定义成第二候选恒星，结合第二候选恒星的二维平面理论坐标、赤经和纬经，以及望远镜的焦距，计算得到任意两颗第二候选恒星之间的角距，选择三颗第二候选恒星组成三角形星图，生成恒星理论星图。

S5：结合实际测得的二维平面坐标，计算得到恒星实测星图上任意三颗定标星i,j,k组成的三角形的边长为

角距大小顺序为

按照预设顺序对计算结果进行排序以生成定标星的实测星图特征库。

结合二维平面理论坐标，计算得到恒星理论星图上任意三颗定标星l,m,n组成的三角形的边长为

角距大小顺序为

按照预设顺序对计算结果进行排序生成定标星的理论星图特征库。

根据下述匹配条件，结合定标星的实测星图特征库和理论星图特征库，对理论星图和实测星图中的恒星进行匹配，直至匹配成功N₃颗恒星：

第j颗定标星和第k颗定标星之间满足下式：

任意三颗i,j,k定标星之间满足下式：

其中，ε₁和ε₂均为预设的角距门限。

S6：设匹配成功的N₃颗恒星在图像上的二维平面坐标为(x_i,y_i)，i＝1,2,…N₃，对应的理论二维平面坐标坐标为(X_i,Y_i)，i＝1,2,…N₃。

利用下式，采用最小二乘方法，计算出系数a,b,c,d,e,f，从而得到中心指向偏差及像面旋转角：

S7：根据图像对应视场大小，结合匹配成功的N₃颗恒星在图像上的二维平面坐标(x_i,y_i)，及理想坐标(ξ_i,ζ_i)，i＝1,2,…N₃，进行常数模型计算，根据恒星的定位精度，自动优选底片常数模型。

简单来说，本发明的技术方案包括以下几个工作步骤：

(1)恒星及空间碎片星象。

(2)恒星检索。

(3)理论星图生成。

(4)实测星图生成。

(5)理论星图及实测星图匹配。

(6)指向及像面旋转测定。

(7)底片模型优选。

(8)空间碎片天文定位。

在实际应用中，针对采集到连续的观测图像，通过空间目标检测获得每一帧图像上恒星及空间碎片的星象信息之后，依次采用上述8个步骤，获得空间碎片的天文定位及测光数据。更优化和更具体描述以上各步骤如下：

一、恒星及空间碎片星象

按照空间碎片检测方法，获得图像上在检测门限内恒星及空间碎片的星象信息，包括其二维平面坐标(x,y)、像素个数、灰度和，按照像素个数减少的顺序对图像进行排序。优选的，可以采用现有技术中的任意一种空间碎片获取方法来获得图像上在检测门限内恒星及空间碎片的星象信息。

二、恒星检索

可以事先生成天文定位恒星星库。其中一种天文定位恒星星库及索引数据的生成方式为：按照给定星等的全天区恒星按照赤经增加及赤纬增加的顺序，分区存放，并形成索引，生成天文定位恒星星库及索引数据，用于恒星检索使用。

根据图像对应的时间和指向信息，测站经纬度、海拔高度，测站温度、湿度、大气压强，给定视场大小，根据加载的恒星星库及索引，快速检索出视场中满足给定星等门限的恒星信息，包括二维平面坐标(x,y)、赤经和赤纬(α,δ)、理想坐标(ξ,ζ)、理论星等，按照理论星等增加的顺序对检索出的恒星进行排序。

所述理想坐标(ξ^s,ζ^s)满足以下公式：

三、实测星图生成

按照给定门限，选择N₁颗恒星星象，根据恒星星象的二维平面坐标(x,y)，望远镜的焦距为f。采用下式计算任意两颗恒星之间的角距，人选三颗恒星组成三角形星图，生成恒星实测星图。

其中，f是望远镜的焦距，(x_u,y_u)是第u颗恒星的二维平面坐标，(x_v,y_v)是第v颗恒星的二维平面坐标。

四、理论星图生成

按照给定门限，选择N₂颗恒星星象，根据恒星星象的二维平面坐标(X,Y)，(α,δ)，望远镜的焦距为f。采用下式计算任意两颗恒星之间的角距，任选三颗恒星组成三角形星图，生成恒星实测星图：

其中，f是望远镜的焦距，(X_u,Y_u)是第u颗第二候选恒星的二维平面理论坐标，(Y_v,Y_v)是第v颗第二候选恒星的二维平面理论坐标。

或者，采用下式计算任意两颗恒星之间的角距，人选三颗恒星组成三角形星图，生成恒星实测星图：

其中，f是望远镜的焦距，(α_u,δ_u)是第u颗第二候选恒星的赤经和纬经，(α_v,δ_v)是第v颗第二候选恒星的赤经和赤纬。

五、理论星图及实测星图匹配

对于恒星实测星图上任意三颗定标星i,j,k组成的三角形的边长为

角距大小顺序为

按照一定的顺序进行排列生成定标星的实测星图特征库。

对于恒星理论星图上任意三颗定标星l,m,n组成的三角形的边长为

角距大小顺序为

按照一定的顺序进行排列生成定标星的理论星图特征库。

由于f可能不是太准，本发明选择下述匹配条件进行恒星实测星图和恒星理论星图的匹配：

第j颗定标星和第k颗定标星之间满足下式：

任意三颗i,j,k定标星之间满足下式：

其中，ε₁和ε₂均为预设的角距门限。

六、指向及像面旋转测定

假设匹配成功了N₃颗恒星，图像上星象二维平面坐标(x_i,y_i)，i＝1,2,…N₃，星象的理论二维平面坐标坐标(X_i,Y_i)，i＝1,2,…N₃。

利用下式，采用最小二乘方法，计算出系数a,b,c,d,e,f，从而得到中心指向偏差及像面旋转角。

七、底片模型优选

根据图像对应视场大小，假设匹配成功了N₃颗恒星，图像上星象二维平面坐标(x_i,y_i)，及理想坐标(ξ_i,ζ_i)，i＝1,2,…N₃。采用以下六常数、十二常数、十四常数模型计算，根据恒星的定位精度，自动优选底片常数模型(只需优选一次，自动存储优选结果)。

(1)六常数模型(需要3个以上的定标星)

(2)十二常数模型(需要6个以上的定标星)

(3)十四常数模型(需要7个以上的定标星)

八、空间碎片天文定位；

根据空间碎片的二维平面坐标实测值(x^T,y^T)，采用以下公式获得空间碎片的赤经和赤纬(α^T,δ^T)：

其中，(ξ^T,ζ^T)为空间碎片的理想坐标，由(x^T,y^T)代入六常数、十二常数或者十四常数模型获得。

计算机系统根据上述输入数据，实时给出了图像上的空间碎片天文定位结果、恒星天文定位结果、指向像面旋转测定结果、恒星检索结果。这些结果可以提供给可以通过显示系统显示出来，以及存储在计算机系统的存储介质中，可以供空间碎片编目定轨及精密定轨使用，可以根据指向测定结果修正空间碎片的预报位置，有利于提高空间碎片的捕获和跟踪成功率。

在本公开中参照附图来描述本发明的各方面，附图中示出了许多说明的实施例。本公开的实施例不必定义在包括本发明的所有方面。应当理解，上面介绍的多种构思和实施例，以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施，这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外，本发明公开的一些方面可以单独使用，或者与本发明公开的其他方面的任何适当组合来使用。

虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。

Claims (10)

1.一种基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，所述实时天文定位方法包括以下步骤：

S1：接收至少一帧包括空间碎片和背景恒星的图像，获得图像上在预设检测门限内的恒星和空间碎片的星象信息，所述星象信息包括每个恒星和空间碎片的二维平面坐标(x,y)、像素个数、灰度和，按照像素个数减少的顺序对所有星象进行排序；其中，设图像的左上角为坐标原点(0，0)，图像右侧为x轴增加方向，图像下侧为y轴增加方向，x为图像中星象所在位置在x轴方向上与坐标原点的距离，y为图像中星象所在位置在y轴方向上与坐标原点的距离；

S2：结合天文定位恒星星库，根据图像对应的拍摄信息和图像中心含有指向误差的指向信息(α_p,δ_p)，检索出视场中满足给定星等门限的所有恒星的相关信息，所述满足给定星等门限的恒星的相关信息包括其所对应的二维平面坐标理论值(X,Y)、赤经和赤纬理论值(α^s,δ^s)、理想坐标理论值(ξ^s,ζ^s)、理论星等M，按照理论星等由小到大的顺序对检索出的恒星进行排序；

S3：按照给定门限，从检索出的恒星中选择N₁颗恒星星象，定义成第一候选恒星，结合第一候选恒星在图像上的二维平面坐标和望远镜的焦距，计算得到任意两颗第一候选恒星之间的角距，选择三颗第一候选恒星组成三角形星图，生成恒星实测星图；

S4：按照给定门限，从检索出的恒星中选择N₂颗恒星星象，定义成第二候选恒星，结合第二候选恒星的二维平面理论坐标、赤经和纬经，以及望远镜的焦距，计算得到任意两颗第二候选恒星之间的角距，选择三颗第二候选恒星组成三角形星图，生成恒星理论星图；

S5：结合实际测得的二维平面坐标，计算得到恒星实测星图上任意三颗定标星i,j,k组成的三角形的边长为

角距大小顺序为

按照预设顺序对计算结果进行排序以生成定标星的实测星图特征库；

结合二维平面理论坐标，计算得到恒星理论星图上任意三颗定标星l,m,n组成的三角形的边长为

角距大小顺序为

按照预设顺序对计算结果进行排序生成定标星的理论星图特征库；

根据下述匹配条件，结合定标星的实测星图特征库和理论星图特征库，对理论星图和实测星图中的恒星进行匹配，直至匹配成功N₃颗恒星：

第j颗定标星和第k颗定标星之间满足下式：

任意三颗i,j,k定标星之间满足下式：

其中，ε₁和ε₂均为预设的角距门限；

S6：设匹配成功的N₃颗恒星在图像上的二维平面坐标为(x_i,y_i)，i＝1,2,…N₃，对应的理论二维平面坐标坐标为(X_i,Y_i)，i＝1,2,…N₃；

利用下式，采用最小二乘方法，计算出系数a,b,c,d,e,f，从而得到中心指向偏差及像面旋转角：

S7：根据图像对应视场大小，结合匹配成功的N₃颗恒星在图像上的二维平面坐标(x_i,y_i)，及理想坐标

进行常数模型计算，根据恒星的定位精度，自动优选底片常数模型。

2.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S2中，所述理想坐标(ξ^s,ζ^s)满足以下公式：

3.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S3中，采用下述公式计算任意两颗第一候选恒星之间的角距：

其中，f是望远镜的焦距，(x_u,y_u)是第u颗第一候选恒星的二维平面坐标，(x_v,y_v)是第v颗第一候选恒星的二维平面坐标。

4.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S4中，采用下述公式计算任意两颗第二候选恒星之间的角距：

其中，f是望远镜的焦距，(X_u,Y_u)是第u颗第二候选恒星的二维平面理论坐标，(Y_v,Y_v)是第v颗第二候选恒星的二维平面理论坐标。

5.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S4中，采用下述公式计算任意两颗第二候选恒星之间的角距：

其中，f是望远镜的焦距，(α_u,δ_u)是第u颗第二候选恒星的赤经和纬经，(α_v,δ_v)是第v颗第二候选恒星的赤经和纬经。

6.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S2中，所述图像对应的拍摄信息包括图像的拍摄时间、指向信息、测站经纬度、测站海拔高度、测站温度、测站湿度、大气压强、给定视场大小。

7.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，所述天文定位恒星星库的生成方法包括：

将给定星等的全天区恒星按照赤经增加及赤纬增加的顺序，分区存放，并形成索引，生成天文定位恒星星库及索引数据。

8.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S7中，结合匹配成功的定标星数量，分别选用六常数模型、十二常数模型、十四常数模型进行常数模型计算，其中：

所述六常数模型对应至少3个以上的定标星：

所述十二常数模型对应至少6个以上的定标星：

所述十四常数模型对应至少7个以上的定标星：

9.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，步骤S7还包括：

自动存储优选底片常数模型。

10.根据权利要求1所述的基于指向误差自动测定的空间碎片实时天文定位方法，其特征在于，所述实时天文定位方法还包括：

S8：根据空间碎片的二维平面坐标实测值(x^T,y^T)，采用以下公式获得空间碎片的赤经和赤纬(α^T,δ^T)：

其中，(ξ^T,ζ^T)为空间碎片的理想坐标，由(x^T,y^T)代入六常数、十二常数或者十四常数模型获得。

Images