WO2021135161A1 - 基于指向自动测定的空间碎片实时天文定位和测光方法 - Google Patents
基于指向自动测定的空间碎片实时天文定位和测光方法 Download PDFInfo
- Publication number
- WO2021135161A1 WO2021135161A1 PCT/CN2020/102323 CN2020102323W WO2021135161A1 WO 2021135161 A1 WO2021135161 A1 WO 2021135161A1 CN 2020102323 W CN2020102323 W CN 2020102323W WO 2021135161 A1 WO2021135161 A1 WO 2021135161A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- star
- stars
- theoretical
- space debris
- image
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/02—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means
- G01C21/025—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by astronomical means with the use of startrackers
Definitions
- the invention relates to the technical field of space debris positioning, in particular to a method for real-time astronomical positioning and photometry of space debris based on automatic pointing determination.
- Absolute positioning is the use of the telescope shaft system to achieve space debris measurement, which is affected by the processing accuracy of the telescope shaft system, the accuracy of atmospheric refraction correction, and temperature deformation. , It does not rely on background stars.
- Relative positioning is based on the relative position of space debris and background stars to achieve space debris measurement.
- the telescope pointing accuracy does not directly affect the measurement results, but when the telescope pointing and image plane installation errors are large, it will cause the theoretical coordinates of the stars on the image.
- absolute measurement there are also two methods for precise photometric measurement of space debris: absolute measurement and relative measurement.
- multiple standard photometric calibration stars in different fields of view are used to obtain a photometric model, which is called absolute measurement model, also called absolute photometry; using multiple photometric calibration stars in the same field of view,
- the obtained metering model is called a relative pair measurement model, also known as relative metering.
- the light metering method it is based on the successful matching of the actual gray scale of the star and the theoretical magnitude of the star, which also requires more accurate station temperature, relative humidity, and atmospheric pressure. Even so, in areas with low elevation angles, the accuracy of the correction of the telescope's pointing error is low due to the correction accuracy of atmospheric refraction.
- the present invention provides a real-time astronomical positioning method for space targets based on automatic pointing determination, which can automatically optimize the film model according to the size of the observation field of view, and according to the given time on the image and the orientation of the image center, There is no need to input pointing information, and the image center orientation and image plane rotation angle are automatically determined, so as to realize the automatic matching of the theoretical coordinates of stars and the actual measured coordinates, thereby realizing real-time astronomical positioning and photometry of space debris.
- the purpose of the present invention is to provide a real-time astronomical positioning and photometric method for space debris based on automatic orientation determination, which can automatically select the film model according to the size of the observation field of view, and automatically determine the image center orientation and image plane rotation angle according to a given time on the image , To realize the automatic matching of the theoretical coordinates of the stars and the actual measured coordinates, to realize the automatic matching of the gray value of the star on the image and its theoretical magnitude, so as to realize the real-time astronomical positioning and photometry of the space debris.
- this method reduces the requirements for telescope shafting accuracy, reduces the requirements for telescope field installation and commissioning, reduces the requirements for inputting environmental temperature parameters at the station, and reduces the need for pointing calibration before observation. Requirements.
- this method can automatically determine the zero point difference of the telescope's two axes without precise astronomical latitude and longitude, and can also realize real-time astronomical positioning and relative photometry of space debris. More importantly, this method can achieve high-precision astronomical positioning and relative photometry on mobile sites without astronomical latitude and longitude (only geographic latitude and longitude).
- the present invention proposes a real-time astronomical positioning and photometric method for space debris based on automatic orientation determination.
- the real-time astronomical positioning and photometric method includes the following steps:
- S1 Generate the astronomical positioning star library and the first index data used to express the information of all the stars contained in the astronomical positioning star library; based on the astronomical positioning star library, generate a theoretical star map of the whole sky and use it for presentation The second index data of the angular distance information between the stars contained in the theoretical star map of the whole sky;
- S2 Receive at least one frame of image including space debris and background stars, and obtain star information of stars and space debris within the preset detection threshold on the image; based on the obtained star information, calculate the angle between any two stars Distance, generate the star map of the star;
- S3 Determine the upper and lower bounds of the measured star map of stars according to the theoretical star map and the second index data of the whole sky area, and obtain the maximum angular distance Minimum angular distance Intermediate angular distance Combined maximum angular distance Minimum angular distance Intermediate angular distance According to the preset matching rules, several stellar images included in the actual star map that match the theoretical star map of the entire sky are calculated;
- the metering model is:
- G i is the gray value of the i-th matched star after subtracting the background
- i 1, 2,..., N 3
- N 3 is the total number of matched star images
- a and B are calculated by the least square method Relative metering model coefficients.
- step S2 the process of generating an actual star map of stars includes the following steps:
- N 1 star images from the image and define it as the first candidate star.
- f is the focal length of the telescope
- (x j , y j ) are the two-dimensional plane coordinates of the j-th first candidate star
- (x k , y k ) are the two-dimensional plane coordinates of the k-th first candidate star.
- step S1 the generation of the astronomical positioning stellar library and the first index data used to express the self-information of all stars contained in the astronomical positioning stellar library refer to:
- step S1 the astronomically positioned star database is used to generate a theoretical star map of the entire sky region, and a second index for expressing the angular distance information between the stars contained in the theoretical star map of the entire sky region Data refers to:
- N 2 stars in the whole sky area and define them as the second candidate star.
- select three second candidate stars to form a triangular star map generate a theoretical star map of the whole sky area, and sort according to the angular distance of each triangle, and generate the corresponding index data.
- step S1 the following formula is used to calculate the angular distance between any two second candidate stars:
- ( ⁇ u , ⁇ u ) are the right ascension and declination of the u-th second candidate star
- ( ⁇ v , ⁇ v ) are the right ascension and declination of the v-th second candidate star
- step S3 according to the preset matching rules, calculating and obtaining several stellar image information of the stars included in the actual measured constellation map that matches the theoretical constellation map of the entire sky region includes the following steps:
- ⁇ 1 and ⁇ 2 are both preset angular distance thresholds
- step S4 the process of calculating the center pointing deviation and the image plane rotation angle based on the information of a number of successfully matched stars includes the following steps:
- step S4 the process of obtaining the film constant model includes the following steps:
- the related information of the star meeting the given magnitude threshold includes its corresponding two-dimensional plane coordinate theoretical values (X, Y), ascension and declination theoretical values ( ⁇ s , ⁇ s ), ideal coordinates
- the shooting information corresponding to the image includes the shooting time of the image, pointing information, station latitude and longitude, station altitude, station temperature, station humidity, atmospheric pressure, and given field of view size;
- the constant model calculation is performed, and the film constant model is automatically selected according to the positioning accuracy of the star, and the selected film constant model is automatically stored.
- the calculation of the constant model includes:
- the six-constant model, the twelve-constant model, and the fourteen-constant model are respectively selected for constant model calculations, among which:
- the six-constant model corresponds to at least 3 calibration stars:
- the twelve-constant model corresponds to at least 6 calibration stars:
- the fourteen constant model corresponds to at least 7 calibration stars:
- the real-time astronomical positioning method further includes:
- ( ⁇ T , ⁇ T ) are the ideal coordinates of the space debris, which are obtained by substituting (x T , y T ) into the six-constant, twelve-constant or fourteen-constant model.
- the film model can be automatically selected according to the size of the observation field of view, and the image center direction and the image plane rotation angle can be automatically determined according to the given time and image center direction on the image, so as to realize the automatic matching of the theoretical and measured coordinates of the star, thus realizing the space Real-time astronomical positioning and photometry of debris.
- this method reduces the requirements for the machining accuracy of the telescope shafting system, reduces the requirements for the installation and commissioning of the telescope field, reduces the input requirements for the ambient temperature parameters of the station, and reduces the need for pre-observation Point to the requirements of the standard school.
- this method can also realize real-time astronomical positioning and relative photometry of space debris when the telescope pointing calibration cannot be achieved without precise astronomical latitude and longitude. Therefore, this method is a very good real-time astronomical positioning and photometric method for space debris.
- This method can achieve high-precision astronomical positioning and relative photometry on mobile sites without astronomical latitude and longitude (only geographic latitude and longitude), and has good actual processing effects, and can be widely used in scientific research and engineering fields.
- the computer system can provide real-time space debris astronomical positioning results, stellar astronomical positioning results, pointing image rotation measurement results, and star retrieval results on the image. These results can be used in a wide range of applications. For example, they can be displayed by the display system and stored in the storage medium of the computer system. They can also be used for orbit determination and precise orbit determination for space debris cataloging, and they can also be used to modify the space based on the results of the pointing measurement. The predicted location of the debris is conducive to improving the success rate of space debris capture and tracking.
- Fig. 1 is a flowchart of the real-time astronomical positioning and photometric method of space debris based on automatic pointing determination of the present invention.
- the real-time astronomical positioning and photometric method includes the following steps:
- S1 Generate the astronomical positioning star library and the first index data used to express the information of all the stars contained in the astronomical positioning star library; based on the astronomical positioning star library, generate a theoretical star map of the whole sky and use it for presentation The second index data of the angular distance information between the stars contained in the all-sky theoretical star map.
- S2 Receive at least one frame of image including space debris and background stars, and obtain star information of stars and space debris within the preset detection threshold on the image; based on the obtained star information, calculate the angle between any two stars Distance to generate the measured star map of the star.
- S3 Determine the upper and lower bounds of the measured star map of stars according to the theoretical star map and the second index data of the whole sky area, and obtain the maximum angular distance Minimum angular distance Intermediate angular distance Combined maximum angular distance Minimum angular distance Intermediate angular distance According to the preset matching rules, a number of stellar images contained in the actual star map that match the theoretical star map of the entire sky area are calculated.
- the metering model is:
- G i is the gray value of the i-th matched star after subtracting the background
- i 1, 2,..., N 3
- N 3 is the total number of matched star images
- a and B are calculated by the least square method Relative metering model coefficients.
- the stars in the entire sky area of a given magnitude are stored in regions in the order of increasing right ascension and increasing declination, and indexing is formed to generate astronomical positioning star database and first index data for use in star retrieval.
- indexing is formed to generate astronomical positioning star database and first index data for use in star retrieval.
- select N 2 stars in the whole sky area and use the following formula to calculate the angular distance between any two stars according to the star’s right ascension and declination information:
- ( ⁇ u , ⁇ u ) are the right ascension and declination of the u-th second candidate star
- ( ⁇ v , ⁇ v ) are the right ascension and declination of the v-th second candidate star
- the given threshold such as the minimum and maximum angular distance
- choose three stars to form a triangular star map generate a theoretical star map of the whole sky, and sort according to the angular distance of each triangle.
- the star information of the stars and space debris within the detection threshold on the image is obtained, including its two-dimensional plane coordinates (x, y), the number of pixels, the gray scale sum, and the image is compared in the order of decreasing the number of pixels. put in order.
- any space debris acquisition method in the prior art can be used to obtain the star information of stars and space debris within the detection threshold on the image.
- the right side of the image is the increasing direction of the x axis
- the lower side of the image is the increasing direction of the y axis
- x is the distance between the position of the star in the image and the origin of the coordinate on the x axis
- Distance, y is the distance between the position of the star in the image and the origin of the coordinate in the y-axis direction.
- N 1 stellar constellations According to the given threshold, select N 1 stellar constellations. According to the two-dimensional plane coordinates (x, y) of the stellar constellations, the focal length of the telescope is f.
- the following formula is used to calculate the angular distance between any two stars, and three stars are selected to form a triangular star chart to generate the actual star chart of the star. Use the following formula to calculate the angular distance between any two first candidate stars:
- (x j , y j ) are the two-dimensional plane coordinates of the j-th second candidate star
- (x k , y k ) are the two-dimensional plane coordinates of the k-th second candidate star.
- the order of angular distance is According to the index of the theoretical star map, the upper and lower bounds n of the candidate star map can be quickly realized.
- the order of the corresponding angular distances of the triangles in the upper and lower bounds is Perform matching judgments.
- Star information for a given magnitude threshold including theoretical values of two-dimensional plane coordinates (X, Y), theoretical values of right ascension and declination ( ⁇ s , ⁇ s ), theoretical values of ideal coordinates ( ⁇ s , ⁇ s ), theoretical values
- the magnitude is M and sorted in the order of increasing theoretical magnitude.
- the ideal coordinates ( ⁇ s , ⁇ s ) satisfy the following formula:
- the film constant model is automatically selected according to the positioning accuracy of the star (only need to be selected once, and the optimization result is automatically stored).
- ( ⁇ T , ⁇ T ) are the ideal coordinates of the space debris, which are obtained by substituting (x T , y T ) into the six-constant, twelve-constant or fourteen-constant model.
- the computer system Based on the above input data, the computer system provides real-time space debris astronomical positioning and photometric results, stellar astronomical positioning and photometric results, pointing image rotation measurement results, and star retrieval results on the image in real time. These results can be displayed by the display system and stored in the storage medium of the computer system. They can be used for orbit determination of space debris cataloging and precise orbit determination. The predicted position of space debris can be corrected according to the pointing measurement results, which is beneficial to improve The success rate of space debris capture and tracking; it can be used for space debris identification and evaluation of the working status and rotation status of space targets with attitude control.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Telescopes (AREA)
- Image Processing (AREA)
Abstract
一种基于指向自动测定的空间碎片实时天文定位和测光方法,该方法包括:理论星图生成;获取恒星及空间碎片星象;实测星图生成;理论星图及实测星图匹配;指向及像面旋转测定;恒星检索;底片模型优选;测光模型计算;空间碎片天文定位和测光。该方法能够根据观测视场大小自动优选底片模型,按照图像上给定时间,自动测定图像中心指向及像面旋转角,实现恒星理论坐标和实测坐标的自动匹配,实现图像上恒星灰度值和其理论星等的自动匹配,从而实现空间碎片的实时天文定位及测光。
Description
本发明涉及空间碎片定位技术领域,具体而言涉及一种基于指向自动测定的空间碎片实时天文定位和测光方法。
在科研、军事等许多领域,都需要对空间碎片进行监视,一方面测定空间碎片的每一个观测时刻在天空中的位置及其变化,确定空间碎片的运行轨道,从而获取空间碎片精确的信息。基于此需求,空间碎片的精确测量是非常重要的基础环节,没有空间碎片的精确测量,空间碎片轨道识别,编目定轨,及精密定轨都无法实现。
空间碎片的位置精确测量有两种方法:绝对定位和相对定位,其中绝对定位就是利用望远镜的轴系实现空间碎片测量,它受到望远镜轴系加工精度、大气折射修正精度、温度变形等因素的影响,它不依赖背景恒星。相对定位是根据空间碎片和背景恒星的相对位置实现空间碎片的测量,望远镜指向精度不直接影响测量结果,但是当望远镜指向及像面安装误差大的情况下,就会造成恒星在图像上理论坐标和恒星在图像上实测坐标相差较大,尤其对于有像面误差的图像,边缘部分误差更大,无法满足给定匹配门限,因此造成恒星理论星图和实测星图匹配失败,无法实现相对定位。
空间碎片的光度精确测量也有两种方法:绝对测量和相对测量。其中利用不同视场中的多个标准测光定标星,得到的测光模型,称之为绝对测量模型,也称为绝对测光;利用同一视场中的多个测光定标星,得到的测光模型,称之为相对对测量模型,也称为相对测光。无论是哪一种测光方式,它都是建立在恒星实测灰度和与恒星理论星等匹配成功的基础上,而这同样需要较为精确的测站温度、相对湿度、及大气压强。即使如此,在仰角较低的天区,受大气折射改正精度影响,望远镜指向误差修正精度偏低。
基于现有天文定位方法的不足,本发明给出一种基于指向自动测定的空间目标实时天文定位方法,它能够根据观测视场大小自动优选底片模型,按照图像上给定时间和图像中心指向,无需指向信息输入,自动测定图像中心指向及像面旋转角,实现恒星理论坐标和实测坐标的自动匹配,从而实现空间碎片的实时天文定位及测光。
发明内容
本发明目的在于提供一种基于指向自动测定的空间碎片实时天文定位和测光方法,能够根据观测视场大小自动优选底片模型,按照图像上给定时间,自动测定图像中心指向及像面旋转角,实现恒星理论坐标和实测坐标的自动匹配,实现图像上恒星灰度值和其理论星等的自动匹配,从而实现空间碎片的实时天文定位及测光。对于固定站址(有精密的天文经纬度)的望远镜,该方法降低了望远镜轴系加工精度要求,降低望远镜外场安装调试要求,减少了测站环境温度参输入要求,减少了观测前需要指向标校的要求。对于可移动望远镜,该方法在无精密天文经纬度不能实现望远镜指向标校的情况下,能够自动测定望远镜两轴零点差, 也能实现空间碎片实时天文定位及相对测光。更重要的是该方法能够在没有天文经纬度的移动站址上(只有地理经纬度)上实现高精度天文定位及相对测光。
为达成上述目的,结合图1,本发明提出一种基于指向自动测定的空间碎片实时天文定位和测光方法,所述实时天文定位和测光方法包括以下步骤:
S1:生成天文定位恒星星库、和用于表述天文定位恒星星库所包含的所有恒星自身信息的第一索引数据;基于天文定位恒星星库,生成全天区理论星图、和用于表述全天区理论星图所包含的恒星之间角距信息的第二索引数据;
S2:接收至少一帧包括空间碎片和背景恒星的图像,获得图像上在预设检测门限内的恒星和空间碎片的星象信息;基于获取的恒星星象信息,计算得到任意两颗恒星之间的角距,生成恒星实测星图;
S3:根据全天区理论星图和第二索引数据,确定恒星实测星图的上下界限,获取最大角距
最小角距
中间角距
结合最大角距
最小角距
中间角距
按照预设的匹配规则,计算得到恒星实测星图中所包含的与全天区理论星图相匹配的若干个恒星星象;
S4:以匹配成功的若干个恒星星象信息为基础,计算得到中心指向偏差、像面旋转角、底片常数模型、相对测光模型;
其中,测光模型为:
其中,G
i是第i颗匹配成功的恒星星象扣除背景后的灰度值,
是第i颗匹配成功的恒星星象对应的理论星等,i=1,2,...,N
3,N
3是匹配成功的恒星星象的总数,A和B是采用最小二乘法计算得到的相对测光模型系数。
进一步的实施例中,步骤S2中,所述生成恒星实测星图的过程包括以下步骤:
按照给定门限,从图像中选择N
1颗恒星星象,定义成第一候选恒星,结合第一候选恒星在图像上的二维平面坐标和望远镜的焦距,计算得到任意两颗第一候选恒星之间的角距,选择三颗第一候选恒星组成三角形星图,生成恒星实测星图;
其中,采用下述公式计算任意两颗第一候选恒星之间的角距:
其中,f是望远镜的焦距,(x
j,y
j)是第j颗第一候选恒星的二维平面坐标,(x
k,y
k)是第k颗第一候选恒星的二维平面坐标。
进一步的实施例中,步骤S1中,所述生成天文定位恒星星库、和用于表述天文定位恒星星库所包含的所有恒星自身信息的第一索引数据是指:
将给定星等的全天区恒星按照赤经增加及赤纬增加的顺序分区存放,并形成索引,生成天文定位恒星星库及索引数据。
进一步的实施例中,步骤S1中,所述基于天文定位恒星星库,生成全天区理论星图、和用于表述全天区理论星图所包含的恒星之间角距信息的第二索引数据是指:
按照给定星等门限,选择全天区N
2颗恒星星象,定义成第二候选恒星,结合第二候选恒星的赤经和赤纬,计算得到任意两颗第二候选恒星之间的角距,按照给定角距门限,任选三颗第二候选恒星组成三角形星图,生成全天区理论星图,并按照每个三角形的角距大小进行排序,生成相应的索引数据。
进一步的实施例中,步骤S1中,采用下式计算任意两颗第二候选恒星之间的角距:
进一步的实施例中,步骤S3中,所述按照预设的匹配规则,计算得到恒星实测星图中所包含的与全天区理论星图相匹配的若干个恒星星象信息包括以下步骤:
第j颗定标星和第k颗定标星之间满足下式:
任意三颗i,j,k定标星之间满足下式:
其中,ε
1和ε
2均为预设的角距门限;
进一步的实施例中,步骤S4中,所述以匹配成功的若干个恒星星象信息为基础,计算得到中心指向偏差、像面旋转角的过程包括以下步骤:
设匹配成功的N
3颗恒星在图像上的二维平面坐标为(x
i,y
i),i=1,2,...N
3,对应的理论二维平面坐标坐标为(X
i,Y
i),i=1,2,...N
3;
利用下式,采用最小二乘方法,计算出系数a,b,c,d,e,f,从而得到中心指向偏差及像面旋转角:
进一步的实施例中,步骤S4中,所述底片常数模型的获取过程包括以下步骤:
结合天文定位恒星星库和第一索引数据,根据图像对应的拍摄信息和全天区星图指向测定结果(α
p,δ
p),检索出视场中满足给定星等门限的所有恒星的相关信息,所述满足给定星等门限的恒星的相关信息包括其所对应的二维平面坐标理论值(X,Y)、赤经和赤纬理论值(α
s,δ
s)、理想坐标理论值(ξ
s,ζ
s)、理论星等M,按照理论星等由小到大的顺序对检索出的恒星进行排序;其中,所述理想坐标(ξ
s,ζ
s)满足以下公式:
其中,所述图像对应的拍摄信息包括图像的拍摄时间、指向信息、测站经纬度、测站海拔高度、测站温度、测站湿度、大气压强、给定视场大小;
根据图像对应视场大小,结合匹配成功的N
3颗恒星在图像上的二维平面坐标(x
i,y
i),及理想坐标
i=1,2,...N
3,进行常数模型计算,根据恒星的定位精度,自动优选底片常数模型,并且自动存储优选出的底片常数模型。
进一步的实施例中,所述进行常数模型计算包括,
结合匹配成功的定标星数量,分别选用六常数模型、十二常数模型、十四常数模型进行常数模型计算,其中:
所述六常数模型对应至少3个以上的定标星:
所述十二常数模型对应至少6个以上的定标星:
所述十四常数模型对应至少7个以上的定标星:
进一步的实施例中,所述实时天文定位方法还包括:
S5:根据空间碎片的二维平面坐标实测值(x
T,y
T),采用以下公式获得空间碎片的赤经和赤纬(α
T,δ
T):
其中,(ξ
T,ζ
T)为空间碎片的理想坐标,由(x
T,y
T)代入六常数、十二常数或者十四常数模型获得。
以上本发明的技术方案,与现有相比,其显著的有益效果在于,
(1)能够根据观测视场大小自动优选底片模型,按照图像上给定时间和图像中心指向,自动测定图像中心指向及像面旋转角,实现恒星理论坐标和实测坐标的自动匹配,从而实现空间碎片的实时天文定位及测光。
(2)对于固定站址(有精密的天文经纬度)的望远镜,该方法降低了望远镜轴系加工精度要求,降低望远镜外场安装调试要求,减少了测站环境温度参输入要求,减少了观测前需要指向标校的要求。对于可移动望远镜,该方法在无精密天文经纬度不能实现望远镜指向标校的情况下,也能实现空间碎片实时天文定位和相对测光。因此该方法是一种的非常好的空间碎片实时天文定位和测光方法。
(3)该方法能够在没有天文经纬度的移动站址上(只有地理经纬度)上实现高精度天文定位和相对测光,实际处理效果好,能够广泛地应用到科研、及工程领域中。
(4)计算机系统能够实时给出图像上的空间碎片天文定位结果、恒星天文定位结果、指向像面旋转测定结果、恒星检索结果。这些结果应用场景广泛,例如,可以提供给可以通过显示系统显示出来,以及存储在计算机系统的存储介质中,也可以供空间碎片编目定轨及精密定轨使用,还可以根据指向测定结果修正空间碎片的预报位置,有利于提高空间碎片的捕获和跟踪成功率。
应当理解,前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外,所要求保护的主题的所有组合都被视为本公开的发明主题的一部分。
结合附图从下面的描述中可以更加全面地理解本发明教导的前述和其他方面、实施例和特征。本发明的其他附加方面例如示例性实施方式的特征和/或有益效果将在下面的描述中显见,或通过根据本发明教导的具体实施方式的实践中得知。
附图不意在按比例绘制。在附图中,在各个图中示出的每个相同或近似相同的组成部分可以用相同的标号表示。为了清晰起见,在每个图中,并非每个组成部分均被标记。现在,将通过例子并参考附图来描述本发明的各个方面的实施例,其中:
图1是本发明的基于指向自动测定的空间碎片实时天文定位和测光方法的流程图。
为了更了解本发明的技术内容,特举具体实施例并配合所附图式说明如下。
结合图1,本发明提及一种基于指向自动测定的空间碎片实时天文定位和测光方法,所述实时天文定位和测光方法包括以下步骤:
S1:生成天文定位恒星星库、和用于表述天文定位恒星星库所包含的所有恒星自身信息的第一索引数据;基于天文定位恒星星库,生成全天区理论星图、和用于表述全天区理论星图所包含的恒星之间角距信息的第二索引数据。
S2:接收至少一帧包括空间碎片和背景恒星的图像,获得图像上在预设检测门限内的恒 星和空间碎片的星象信息;基于获取的恒星星象信息,计算得到任意两颗恒星之间的角距,生成恒星实测星图。
S3:根据全天区理论星图和第二索引数据,确定恒星实测星图的上下界限,获取最大角距
最小角距
中间角距
结合最大角距
最小角距
中间角距
按照预设的匹配规则,计算得到恒星实测星图中所包含的与全天区理论星图相匹配的若干个恒星星象。
S4:以匹配成功的若干个恒星星象信息为基础,计算得到中心指向偏差、像面旋转角、底片常数模型、相对测光模型。
其中,测光模型为:
其中,G
i是第i颗匹配成功的恒星星象扣除背景后的灰度值,
是第i颗匹配成功的恒星星象对应的理论星等,i=1,2,...,N
3,N
3是匹配成功的恒星星象的总数,A和B是采用最小二乘法计算得到的相对测光模型系数。
简单来说,本发明的技术方案包括以下几个工作步骤:
(1)理论星图生成。
(2)恒星及空间碎片星象。
(3)实测星图生成。
(4)理论星图及实测星图匹配。
(5)指向及像面旋转测定。
(6)恒星检索。
(7)底片模型优选。
(8)测光模型计算。
(9)空间碎片天文定位。
在实际应用中,针对采集到连续的观测图像,通过空间目标检测获得每一帧图像上恒星及空间碎片的星象信息之后,依次采用上述9个步骤,获得空间碎片的天文定位及测光数据。更优化和更具体描述以上各步骤如下:
一、理论星图生成
按照给定星等的全天区恒星按照赤经增加及赤纬增加的顺序,分区存放,并形成索引,生成天文定位恒星星库及第一索引数据,用于恒星检索使用。按照给定星等门限,选择全天区N
2颗恒星星象,根据恒星的赤经和赤纬信息,采用下式计算任意两颗恒星之间的角距:
其中,(α
u,δ
u)是第u颗第二候选恒星的赤经和赤纬,(α
v,δ
v)是第v颗第二候选恒星的赤经和赤纬
按照给定门限,比如最小及最大角距,任选三颗恒星组成三角形星图,生成全天区理论星图,并按照每个三角形的角距大小进行排序。生成全天区理论星图库及相应的第二索引数 据,供星图检索使用。
二、恒星及空间碎片星象
按照空间碎片检测方法,获得图像上在检测门限内恒星及空间碎片的星象信息,包括其二维平面坐标(x,y)、像素个数、灰度和,按照像素个数减少的顺序对图像进行排序。优选的,可以采用现有技术中的任意一种空间碎片获取方法来获得图像上在检测门限内恒星及空间碎片的星象信息。其中,设图像的左上角为坐标原点(0,0),图像右侧为x轴增加方向,图像下侧为y轴增加方向,x为图像中星象所在位置在x轴方向上与坐标原点的距离,y为图像中星象所在位置在y轴方向上与坐标原点的距离。
三、实测星图生成
按照给定门限,选择N
1颗恒星星象,根据恒星星象的二维平面坐标(x,y),望远镜的焦距为f。采用下式计算任意两颗恒星之间的角距,人选三颗恒星组成三角形星图,生成恒星实测星图。采用下述公式计算任意两颗第一候选恒星之间的角距:
其中,(x
j,y
j)是第j颗第二候选恒星的二维平面坐标,(x
k,y
k)是第k颗第二候选恒星的二维平面坐标。
四、理论星图及实测星图匹配
由于f可能不是太准,但是第j颗定标星和第k颗定标星之间依然能够满足下式:
五、指向及像面旋转测定。
假设匹配成功了N
3颗恒星,图像上星象二维平面坐标(x
i,y
i),i=1,2,...N
3,星象的理论二维平面坐标坐标(X
i,Y
i),i=1,2,...N
3。
利用下式,采用最小二乘方法,计算出系数a,b,c,d,e,f,从而得到中心指向偏差及像面旋转角:
六、恒星检索。
根据图像对应的时间和全天区星图指向测定结果(α
p,δ
p),测站经纬度、海拔高度,测站温度、湿度、大气压强,给定视场大小,检索出视场中满足给定星等门限的恒星信息,包括二维平面坐标理论值(X,Y)、赤经和赤纬理论值(α
s,δ
s)、理想坐标理论值(ξ
s,ζ
s)、理论星等M,并按照理论星等增加的顺序进行排序。所述理想坐标(ξ
s,ζ
s)满足以下公式:
七、底片模型优选。
根据图像对应视场大小,假设匹配成功了N
3颗恒星,图像上星象二维平面坐标(x
i,y
i),及理想坐标
i=1,2,...N
3。采用以下六常数、十二常数、十四常数模型计算,根据恒星的定位精度,自动优选底片常数模型(只需优选一次,自动存储优选结果)。
六常数模型(需要3个以上的定标星)
十二常数模型(需要6个以上的定标星)
十四常数模型(需要7个以上的定标星)
八、测光模型计算
九、空间碎片天文定位。
根据空间碎片的二维平面坐标实测值(x
T,y
T),采用以下公式获得空间碎片的赤经和赤纬(α
T,δ
T):
其中,(ξ
T,ζ
T)为空间碎片的理想坐标,由(x
T,y
T)代入六常数、十二常数或者十四常数模型获得。
计算机系统根据上述输入数据,实时给出了图像上的空间碎片天文定位及测光结果、恒星天文定位及测光结果、指向像面旋转测定结果、恒星检索结果。这些结果可以提供给可以通过显示系统显示出来,以及存储在计算机系统的存储介质中,可以供空间碎片编目定轨及 精密定轨使用,可以根据指向测定结果修正空间碎片的预报位置,有利于提高空间碎片的捕获和跟踪成功率;可以供空间碎片识别使用,评估有姿控空间目标的工作状态及旋转状态等。
在本公开中参照附图来描述本发明的各方面,附图中示出了许多说明的实施例。本公开的实施例不必定义在包括本发明的所有方面。应当理解,上面介绍的多种构思和实施例,以及下面更加详细地描述的那些构思和实施方式可以以很多方式中任意一种来实施,这是因为本发明所公开的构思和实施例并不限于任何实施方式。另外,本发明公开的一些方面可以单独使用,或者与本发明公开的其他方面的任何适当组合来使用。
虽然本发明已以较佳实施例揭露如上,然其并非用以限定本发明。本发明所属技术领域中具有通常知识者,在不脱离本发明的精神和范围内,当可作各种的更动与润饰。因此,本发明的保护范围当视权利要求书所界定者为准。
Claims (10)
- 一种基于指向自动测定的空间碎片实时天文定位和测光方法,其特征在于,所述实时天文定位和测光方法包括以下步骤:S1:生成天文定位恒星星库、和用于表述天文定位恒星星库所包含的所有恒星自身信息的第一索引数据;基于天文定位恒星星库,生成全天区理论星图、和用于表述全天区理论星图所包含的恒星之间角距信息的第二索引数据;S2:接收至少一帧包括空间碎片和背景恒星的图像,获得图像上在预设检测门限内的恒星和空间碎片的星象信息;基于获取的恒星星象信息,计算得到任意两颗恒星之间的角距,生成恒星实测星图;S3:根据全天区理论星图和第二索引数据,确定恒星实测星图的上下界限,获取最大角距 最小角距 中间角距 结合最大角距 最小角距 中间角距 按照预设的匹配规则,计算得到恒星实测星图中所包含的与全天区理论星图相匹配的若干个恒星星象;S4:以匹配成功的若干个恒星星象信息为基础,计算得到中心指向偏差、像面旋转角、底片常数模型、相对测光模型;其中,测光模型为:
- 根据权利要求1所述的基于指向自动测定的空间碎片实时天文定位和测光方法,其特征在于,步骤S1中,所述生成天文定位恒星星库、和用于表述天文定位恒星星库所包含的所有恒星自身信息的第一索引数据是指,将给定星等的全天区恒星按照赤经增加及赤纬增加的顺序分区存放,并形成索引,生成天文定位恒星星库及索引数据。
- 根据权利要求1所述的基于指向自动测定的空间碎片实时天文定位和测光方法,其特征在于,步骤S1中,所述基于天文定位恒星星库,生成全天区理论星图、和用于表述全天区 理论星图所包含的恒星之间角距信息的第二索引数据是指,按照给定星等门限,选择全天区N 2颗恒星星象,定义成第二候选恒星,结合第二候选恒星的赤经和赤纬,计算得到任意两颗第二候选恒星之间的角距,按照给定角距门限,任选三颗第二候选恒星组成三角形星图,生成全天区理论星图,并按照每个三角形的角距大小进行排序,生成相应的索引数据。
- 根据权利要求1所述的基于指向自动测定的空间碎片实时天文定位和测光方法,其特征在于,步骤S3中,所述按照预设的匹配规则,计算得到恒星实测星图中所包含的与全天区理论星图相匹配的若干个恒星星象信息包括以下步骤:第j颗定标星和第k颗定标星之间满足下式:任意三颗i,j,k定标星之间满足下式:其中,ε 1和ε 2均为预设的角距门限。
- 根据权利要求1所述的基于指向自动测定的空间碎片实时天文定位和测光方法,其特征在于,步骤S4中,所述底片常数模型的获取过程包括以下步骤:结合天文定位恒星星库和第一索引数据,根据图像对应的拍摄信息和全天区星图指向测定结果(α p,δ p),检索出视场中满足给定星等门限的所有恒星的相关信息,所述满足给定星等门限的恒星的相关信息包括其所对应的二维平面坐标理论值(X,Y)、赤经和赤纬理论值(α s,δ s)、理想坐标理论值(ξ s,ζ s)、理论星等M,按照理论星等由小到大的顺序对检索出的恒星进行排序;其中,所述理想坐标(ξ s,ζ s)满足以下公式:其中,所述图像对应的拍摄信息包括图像的拍摄时间、指向信息、测站经纬度、测站海拔高度、测站温度、测站湿度、大气压强、给定视场大小;
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911416733.2 | 2019-12-31 | ||
CN201911416733.2A CN111156990B (zh) | 2019-12-31 | 2019-12-31 | 基于指向自动测定的空间碎片实时天文定位和测光方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2021135161A1 true WO2021135161A1 (zh) | 2021-07-08 |
Family
ID=70560257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2020/102323 WO2021135161A1 (zh) | 2019-12-31 | 2020-07-16 | 基于指向自动测定的空间碎片实时天文定位和测光方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN111156990B (zh) |
WO (1) | WO2021135161A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111156990B (zh) * | 2019-12-31 | 2021-03-16 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位和测光方法 |
CN113552648A (zh) * | 2021-07-09 | 2021-10-26 | 中国测绘科学研究院 | 空间目标的天文定位及自动跟踪方法、系统及电子设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7218384B2 (en) * | 2002-02-27 | 2007-05-15 | Pentax Corporation | Surveying system |
CN102374856A (zh) * | 2011-11-15 | 2012-03-14 | 中国科学院紫金山天文台 | 天文望远镜实时精确指向的测定方法 |
CN102506829A (zh) * | 2011-11-15 | 2012-06-20 | 中国科学院紫金山天文台 | 快速星图匹配方法 |
CN103925922A (zh) * | 2013-08-05 | 2014-07-16 | 中国航天科工集团第三研究院第八三五八研究所 | 一种适用于高动态条件下iccd星图的恒星识别方法 |
CN104331888A (zh) * | 2014-10-30 | 2015-02-04 | 西北工业大学 | 基于最小二乘拟合的星等标定方法 |
CN111156990A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位和测光方法 |
CN111156989A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位方法 |
CN111156991A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向误差自动测定的空间碎片实时天文定位方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2929395B1 (fr) * | 2008-03-27 | 2010-05-21 | Centre Nat Etd Spatiales | Methode d'estimation d'attitude d'un senseur stellaire |
CN102981262A (zh) * | 2012-11-20 | 2013-03-20 | 中国科学院紫金山天文台 | 大口径透射式望远镜系统 |
CN104776848B (zh) * | 2015-04-20 | 2018-05-22 | 李智 | 一种空间目标识别、定位、跟踪方法 |
-
2019
- 2019-12-31 CN CN201911416733.2A patent/CN111156990B/zh active Active
-
2020
- 2020-07-16 WO PCT/CN2020/102323 patent/WO2021135161A1/zh active Application Filing
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7218384B2 (en) * | 2002-02-27 | 2007-05-15 | Pentax Corporation | Surveying system |
CN102374856A (zh) * | 2011-11-15 | 2012-03-14 | 中国科学院紫金山天文台 | 天文望远镜实时精确指向的测定方法 |
CN102506829A (zh) * | 2011-11-15 | 2012-06-20 | 中国科学院紫金山天文台 | 快速星图匹配方法 |
CN103925922A (zh) * | 2013-08-05 | 2014-07-16 | 中国航天科工集团第三研究院第八三五八研究所 | 一种适用于高动态条件下iccd星图的恒星识别方法 |
CN104331888A (zh) * | 2014-10-30 | 2015-02-04 | 西北工业大学 | 基于最小二乘拟合的星等标定方法 |
CN111156990A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位和测光方法 |
CN111156989A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向自动测定的空间碎片实时天文定位方法 |
CN111156991A (zh) * | 2019-12-31 | 2020-05-15 | 中国科学院紫金山天文台 | 基于指向误差自动测定的空间碎片实时天文定位方法 |
Non-Patent Citations (1)
Title |
---|
SUN RONGYU, LU YAO;ZHAO CHANG-YIN: "A Method for Correcting Telescope Pointing Error in Optical Space Debris Surveys", ACTA ASTRONOMICA SINICA, vol. 56, no. 3, 1 May 2015 (2015-05-01), pages 253 - 263, XP055828229, ISSN: 0001-5245, DOI: 10.15940/j.cnki.0001-5245.2015.03.006 * |
Also Published As
Publication number | Publication date |
---|---|
CN111156990A (zh) | 2020-05-15 |
CN111156990B (zh) | 2021-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111156988B (zh) | 基于指向误差自动测定的空间碎片天文定位和测光方法 | |
CN103047985B (zh) | 一种空间目标的快速定位方法 | |
CN111156989B (zh) | 基于指向自动测定的空间碎片实时天文定位方法 | |
US11635481B2 (en) | Stellar atmospheric refraction measurement correction method based on collinearity of refraction surfaces | |
CN106871927B (zh) | 一种无人机光电吊舱安装误差标校方法 | |
CN109708649B (zh) | 一种遥感卫星的姿态确定方法及系统 | |
CN101246590B (zh) | 星载相机空间畸变图像几何校正方法 | |
CN110940354B (zh) | 一种光电跟踪系统捷联惯导安装姿态的标定方法 | |
CN101539397B (zh) | 物体三维姿态的精密光学测量方法 | |
CN104913780A (zh) | 集成gnss和ccd天顶筒的高精度垂线偏差快速测量方法 | |
WO2021135161A1 (zh) | 基于指向自动测定的空间碎片实时天文定位和测光方法 | |
CN114001756B (zh) | 一种小视场星敏感器外场地面寻星方法 | |
CN115079728A (zh) | 小视场相机快速恒星瞄准跟踪装置、方法及探测装置 | |
CN108398123A (zh) | 一种全站仪及其度盘标定方法 | |
CN112489091B (zh) | 一种基于直瞄模板的全捷联图像导引头目标跟踪方法 | |
CN108225294A (zh) | 一种组合式船用平台罗经系泊态航向标定方法 | |
CN113218577A (zh) | 一种星敏感器星点质心位置精度的外场测量方法 | |
CN111156991B (zh) | 基于指向误差自动测定的空间碎片实时天文定位方法 | |
CN108154535A (zh) | 基于平行光管的摄像机标定方法 | |
CN103852079B (zh) | 一种基于双星顶点剖分弧度集合模糊匹配的船舶天文导航方法 | |
CN110887474B (zh) | 一种用于精密跟踪型望远镜的星图识别方法 | |
CN113916217A (zh) | 一种基于分区平流层大气折射模型的星光定位方法 | |
CN104458653B (zh) | 一种测量大天顶距处的大气折射值的方法和系统 | |
CN115077560A (zh) | 一种船载可见光及中波红外系统光轴平行度动态检测方法 | |
CN114998423B (zh) | 一种望远镜自主定向方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20908890 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20908890 Country of ref document: EP Kind code of ref document: A1 |