Stellar atmospheric refraction measurement correction method based on collinearity of refraction surfaces

Abstract

A stellar atmospheric refraction measurement correction method based on collinearity of refraction surfaces, comprising: performing star identification on the basis of observed star vectors in a star sensor and the reference star catalog, to obtain matching relationships between observed stars and reference stars; converting reference star vectors corresponding to the observed stars to a geographic coordinate system before entering the atmosphere to obtain zenith distances and azimuth angles of incident stellar; on the basis of a principle of collinearity of refraction surfaces, performing optimal solving according to imaging coordinates of observation stars, to obtain the optimal position coordinates of the zenith direction on an imaging surface of the star sensor; according to the optimal zenith direction, performing atmospheric refraction correction on all the recognized observed stars by means of the trigonometric cosine formula to obtain corrected star coordinates; and performing optimal solving to obtain the attitude of the star sensor in the geographic coordinate system.

Claims

1. A method for measuring and correcting stellar atmospheric refraction based on collinear refraction surfaces; characterized in that it comprises the following steps: Step 1, based on a star identification method, obtaining a matching relationship between an observation star in the star sensor and a reference star in a reference star catalogue according to the position coordinates of the observation star on the imaging plane of the star sensor and the reference star catalogue; Step 2, according to the first star identification result and the observation time and the observation position, converting the direction vector of the reference star corresponding to the observation star in the star sensor from the reference star coordinate system into a geographic coordinate system before entering the atmosphere layer; ie, the East-North-Up coordinate system, to obtain the zenith distance-altitude and the azimuth angle before entering the atmosphere layer, Step 3, based on the stellar atmospheric refraction surface collinear principle According to the position coordinates of the observed star on the imaging plane of the star sensor and the azimuth angle of the stellar obtained in the second step before entering the atmosphere layer, the position coordinates of the top direction of the earth on the imaging surface of the star sensor are solved; Step 4, performing atmospheric refraction correction on the coordinates of all the observation stars in the star sensor on the imaging plane after the star is identified according to the position coordinates of the earth top direction obtained in the third step and the sky top distance before the stellar obtained in the second step enters the atmosphere layer, so as to obtain theoretical position coordinates of all observed stars on the imaging surface of the star sensor after the atmospheric refraction is removed; Step 5, according to the calculated theoretical position coordinates of all the observation stars on the imaging surface of the star sensor after the atmospheric refraction is removed in the fourth step, the matching relationship between the star and the reference star in the first step, and the direction vector of the reference star in the geographic coordinate system in the second step, the attitude of the star sensor in the geographic coordinate system is solved.

2. The stellar atmospheric measuring and correcting method based on collinear refraction surfaces according to claim 1; wherein in the first step, the star recognition implementation steps are as follows: (a1) setting the position coordinates of the ith observation star on the imaging surface of the star sensor as A.sub.i(x.sub.A, y.sub.A), and obtaining an observation star vector of the observation star in the star sensor coordinate system according to the optical parameters of the star sensor, namely the size, the focal length and the main point of the star sensor O.sub.sA.sub.i, wherein OS is the origin of the star sensor coordinate system; (a2) observing the inter-satellite angular distance between two observation star vectors based on the star identification method, and matching them with the reference angular distance formed by the reference star in the reference star catalogue to obtain the matching relationship between the observation star in the star sensor and the reference star in the reference star catalogue.

3. The stellar atmospheric measuring and correcting method based on collinear refraction surfaces according to claim 1, wherein in the second step, the antenna top distance and the azimuth angle of the incident stellar before entering the atmosphere layer are implemented as follows: (b1) obtaining the observation time of the star sensor and the longitude and latitude of the observation position, and converting the observation time to the earth mechanics time; (b2) taking out the position vector of the reference star matched with the observation star in the star sensor from the reference star table according to the result of the first star identification and the related astronomical parameters thereof, wherein the astronomical parameters comprise the red warp, the red weft, the sun and the parallax; (b3) based on the International Astronomical Union (IAU) basic Standards of Fundamental Astronomy (SOFA), observing the observation time and the observation position and substituting the reference star position vector and the related astronomical parameters of the star into an astronomical conversion formula, and calculating to obtain a position vector of the reference star matched with the observation star in the star sensor in a geographic coordinate system, ie, the East-North-Up coordinate system, before entering the atmosphere layer, wherein the vector is represented by the antenna top distance and the azimuth angle.

4. The stellar atmospheric measuring and correcting method based on collinear refraction surfaces according to claim 1, characterized in that: in the third step, the steps of obtaining the position coordinates of the zenith direction of the earth on the imaging plane of the star sensor are as follows: (c1) setting the theoretical position coordinates of the top direction of the earth on the imaging surface of the star sensor as B(x.sub.B, y.sub.B), and obtaining the direction vector of the top direction of the earth in the star sensor coordinate system according to the optical parameters, namely the size, the focal length and the main point of the star sensor O.sub.sB; (c2) the ith observed star vector O.sub.sA.sub.i in the star sensor and the direction vector O.sub.sB of the earth's zenith form the atmospheric refraction surface O.sub.sA.sub.iB, the vector of the refraction surface method is FA.sub.i=O.sub.sA.sub.i×O.sub.sB, the atmospheric refraction surfaces of different stars in the star sensor intersect with the earth's zenith direction vector O.sub.sB, that is, the refraction face collinear; (c3) the angle between the atmospheric refraction surfaces O.sub.sA.sub.iB and O.sub.sA.sub.iB of the i-th observed star and the j-th observed star in the star sensor is equal to the angle between their normal vectors, since the azimuth angle does not change before and after the atmospheric refraction of the stellar, the angle between the normal vectors of the atmospheric refraction surfaces O.sub.sA.sub.iB and O.sub.sA.sub.jB is equal to the angle between the azimuth angles before the observed stars i and j enter the atmosphere, namely: a cos(FAi.Math.FAj/|FAj||FAj|)=θ.sub.Azimuth,Ai−θ.sub.Azimuth,Aj, where θ.sub.Azimuth,Ai and θ.sub.Azimuth,Aj are azimuth angles of the observed stars i and j before they enter the atmosphere obtained in the second step; (c4) any three observation stars in the star sensor form a refraction surface triangle, and the refractive surface triangles form a least square error model, namely, $\underset{i\ne j=1}{\overset{3}{.Math.}}{.Math.\mathrm{acos}\left(\stackrel{\_}{\mathrm{FAi}}.Math.\stackrel{\_}{\mathrm{FAj}}/.Math.\stackrel{\_}{\mathrm{FAj}}.Math..Math.\stackrel{\_}{\mathrm{FAj}}.Math.\right)-\left({\theta }_{\mathrm{Azimuth},\mathrm{Ai}}-{\theta }_{\mathrm{Azimuth},\mathrm{Aj}}\right).Math.}_2^2;$ (c5) all observed stars in the star sensor form C.sub.N.sup.3 groups of refracting triangles, where N is the number of observed stars identified after star recognition, the least squares error model of all the refractive surface triangles in the star sensor is accumulated to obtain a position coordinate solving model of the top direction of the earth on the imaging surface of the star sensor, that is, $\left(\stackrel{\_}{x_B},\stackrel{\_}{y_B}\right)=\arg \min \underset{n=1}{\overset{M}{.Math.}}\left(\underset{i\ne j=1}{\overset{3}{.Math.}}{.Math.\mathrm{acos}\left(\stackrel{\_}{\mathrm{FAi}}.Math.\stackrel{\_}{\mathrm{FAj}}/.Math.\stackrel{\_}{\mathrm{FAj}}.Math..Math.\stackrel{\_}{\mathrm{FAj}}.Math.\right)-\left({\theta }_{\mathrm{Azimuth},\mathrm{Ai}}-{\theta }_{\mathrm{Azimuth},\mathrm{Aj}}\right).Math.}_2^2\right)$ in the formula M=C.sub.n.sup.3; Finally, theoretical position coordinates of the optimal earth top direction of the whole star image on the imaging surface of the star sensor are obtained B(x.sub.b, y.sub.b).

5. The stellar atmospheric measuring and correcting method based on collinear refraction surfaces according to claim 1, characterized in that: in the fourth step, an observation sun in the star sensor is subjected to atmospheric refraction correction to obtain theoretical position coordinates of all observed stars on the imaging surface of the star sensor after the atmospheric refraction is removed, and the steps are as follows: (d1) the azimuth angle of the i-th observation sun on the imaging surface of the star sensor is A.sub.i(x.sub.A, y.sub.A) and the position coordinates after atmospheric refraction correction are the same as the azimuth angle of the A′.sub.i(x′.sub.A, y′.sub.A), on the imaging surface of the star sensor, the position coordinates A.sub.i(x.sub.A, y.sub.A), A′.sub.i(x′.sub.A, y′.sub.A) and the theoretical position coordinates B(x.sub.B, y.sub.B) on the imaging plane of the star sensor are collinear, that is, BA.sub.i.Math.BA′.sub.i=0; (d2) setting the top distance of the i-th observation star on the imaging surface of the star sensor in the second step as θ.sub.Zenith,Ai, and according to the cosine theorem |A.sub.iA′.sub.i|.sup.2=|O.sub.sA.sub.i|.sup.2+|O.sub.sA′.sub.i|.sup.2−2|O.sub.sA.sub.i||O.sub.sA′.sup.i|cos(θ.sub.Zenith,Ai−∠BO.sub.sA.sub.i); (d3) the above-mentioned implementation steps (d1) and (d2) are combined to form a system equation set, and theoretical position coordinates A′.sub.i(x′.sub.A, y′.sub.A) of the observation sun i on the imaging surface of the star sensor after the atmospheric refraction is removed are solved; (d4) atmospheric refraction correction is carried out on other observation stars in the star sensor according to steps (d1), (d2) and (d3), so that theoretical position coordinates after atmospheric refraction are removed.

6. The stellar atmospheric measuring and correcting method based on collinear refraction surfaces according to claim 1, wherein the step of solving the attitude of the star sensor in the geographic coordinate system comprises the following steps: (e1) obtaining the direction vector of the star i under the star sensor coordinate system after atmospheric refraction correction according to the optical parameters of the star sensor, namely the size, the focal length and the main point of the star sensor O.sub.sA′.sub.i; (e2) obtaining the attitude of the star sensor in the geographic coordinate system by using the attitude solving method according to the matching relationship between the observation star in the star sensor in the first step and the reference star in the reference star table and the position vector before the reference star enters the atmosphere in the geographic coordinate system in the second step.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

(1) FIG. 1 is the realization flow chart of the method of the present invention;

(2) FIG. 2 is a schematic diagram of the collinearity of the refraction surface of the stellar atmosphere, wherein (a) is a schematic diagram of the refraction surface of the star, and (h) is a schematic diagram of the normal vector of the refraction surface;

(3) FIG. 3 is a schematic diagram of atmospheric refraction correction;

(4) FIG. 4 shows the coordinates of point B and the optimization results obtained by solving different refracting surface triangles of a frame of star image;

(5) FIG. 5 shows the average angular distance deviation before and after atmospheric refraction correction for 1000 frames of experimental star images, (a) is a schematic diagram of the average angular distance deviation before atmospheric refraction correction, and (b) is a schematic diagram of the average angular distance deviation after atmospheric refraction correction.

DETAILED DESCRIPTION

(6) The present invention will be described in detail below with reference to the accompanying drawings and embodiments.

(7) As shown in FIG. 1, the present invention proposes a real-time atmospheric refraction correction measurement method based on the collinearity of refracting surfaces. The method is based on the coplanarity of atmospheric refraction, that is, atmospheric refraction only changes the zenith angle of incident light without affecting its incident azimuth. Based on this, different observed stars in the star sensor have their own atmospheric refraction surfaces. Without considering other errors, the atmospheric refraction surfaces of all observed stars intersect on a straight line, that is, the zenith direction of the earth at the observation time and observation position of the star sensor.

(8) Although the pressure, temperature, humidity and other factors at various positions of the stellar's trajectory through the atmosphere will affect the atmospheric refraction changes, they are always in a refraction plane. Therefore, the atmospheric refraction correction process based on the collinearity of the refracting surfaces has nothing to do with the propagation trajectory of stellar in the atmosphere, and its refraction correction accuracy will not be affected by the atmosphere's internal environment (temperature, pressure, humidity, etc.) and the refraction trajectory during the actual observation process.

(9) First, it is assumed that the position coordinate of the i-th observed star on the imaging surface of the star sensor is A.sub.i(x.sub.Ai, y.sub.Ai). According to the optical parameters of the star sensor, that is, the pixel size, focal length and principal point, the observed star vector O.sub.sA.sub.i of the observed star in the star sensor coordinate system is obtained:
O.sub.sAi=[x.sub.Ai−x.sub.0,y.sub.Ai−y.sub.0,−f]  (1)

(10) In the formula, (x.sub.0, y.sub.0) is the principal point of the star sensor, f is the focal length of the star sensor, and O.sub.s is the origin of the star sensor coordinate system, as shown in (a) in FIG. 2.

(11) Based on the star recognition method, the angular distance between the observed stars is composed of two pairs of observed star vectors, and they are matched with the reference angular distance composed of the reference stars in the reference star catalog to obtain the angular distance in the star sensor. Matches between observed stars and reference stars in the reference catalog. The observed star in the star sensor corresponds to the reference star in the reference catalog.

(12) Then, based on the International Astronomical Union (IAU) Basic Standards of Fundamental Astronomy (SOFA), the reference star corresponding to the observed star in the star sensor is converted from the reference star catalog coordinate system to the geographic coordinate system of the longitude and latitude of the observation position at the observation time t (East-North-Up). The reference vector before entering the atmosphere in the coordinate system) is:
V.sub.i(θ.sub.Zenith,Ai,θ.sub.Azimuth,Ai)=R.sub.2(90°−λ)R.sub.3(ϕ)W(t)R.sub.3(−β)R.sub.Σf[g[u.sub.B(t.sub.0+(t−t.sub.0){dot over (u)}.sub.B(t.sub.0)−πE.sub.B(t)]]  (2)

(13) In the formula, θ.sub.Zenith,Ai and θ.sub.Azimuth,Ai are the zenith distance and azimuth angle of the reference star before entering the atmosphere in the geographic coordinate system (East-North-Up coordinate system), respectively, t is the geomechanical time at the time of observation, t.sub.0 is the reference epoch time of the Hipparcos Astronomical Catalogue t.sub.0=J2000.0, and u.sub.B (t.sub.0) is the position vector in the reference star v.sub.i the reference star catalog, ū.sub.B(t.sub.0) is the stellar motion velocity (arcsec/year) at the stellar reference epoch t.sub.0, E.sub.B(t) is the position ephemeris of the Earth in the International Celestial Coordinate System (JCRS) at the time t of observation, g[ . . . ] is the gravitational light deflection function, ƒ[ . . . ] is the optical aberration function, R.sub.Σ is the precession-nutation transformation matrix, R.sub.3(−β) is the earth's rotation matrix, W(t) is the earth's polar shift matrix, R.sub.2(90°−λ)R.sub.3(ϕ) is the transformation matrix from the earth coordinate system to the geographic coordinate system, R.sub.1, R.sub.2 and R.sub.3 represent the unit rotation angle matrix of x, y and z axes respectively.

(14) As shown in FIG. 2(a), it is assumed that the theoretical position coordinate of the earth's zenith direction on the imaging plane of the star sensor is B(x.sub.B, y.sub.B). According to the optical parameters of the star sensor, that is, the pixel size, focal length and principal point, the zenith direction vector O.sub.sB of the earth is obtained:
O.sub.sB=[x.sub.B−x.sub.0,y.sub.B−y.sub.0,−f]  (3)

(15) In the star sensor, the observed star vector O.sub.sA and the earth's zenith direction vector O.sub.sB form the atmospheric refraction surface O.sub.sA.sub.iB. The atmospheric refraction surfaces of different stars in the star sensor intersect at the zenith direction vector O.sub.sB of the earth, that is, the refraction surfaces are collinear. Based on the observed star vector O.sub.sA.sub.i and the earth's zenith direction vector O.sub.sB, the normal vector FAi of the atmospheric refraction surface O.sub.sA.sub.iB is:
FAi=O.sub.sAi×O.sub.sB  (4)

(16) The angle between the refraction surfaces of the atmospheres of different observed stars is equal to the angle between their normal vectors. Because atmospheric refraction only changes the zenith angle of incident light without affecting its incident azimuth, the azimuth angle before and after atmospheric refraction of stellar does not change. Therefore, the angle between the normal vectors of the atmospheric refraction surface is equal to the angle between their azimuths before entering the atmosphere.

(17) Suppose there are atmospheric refraction surfaces O.sub.sA.sub.iB and O.sub.sA.sub.jB of the i-th observed star and the j-th observed star in the star sensor, and their plane normal vectors are FAi and FAj, respectively. The angle between the normal vector of the atmospheric refraction surface O.sub.sA.sub.iB and O.sub.sA.sub.jB satisfies the following formula:

(18) $\begin{array}{cc}\mathrm{acos}\left(\frac{\stackrel{\_}{\mathrm{FAi}}.Math.\stackrel{\_}{\mathrm{FAj}}}{.Math.\stackrel{\_}{\mathrm{FAj}}.Math..Math.\stackrel{\_}{\mathrm{FAj}}.Math.}\right)={\theta }_{\mathrm{Azimuth},\mathrm{Ai}}-{\theta }_{\mathrm{Azimuth},\mathrm{Aj}}& \left(5\right)\end{array}$

(19) In the formula, θ.sub.Azimuth,Ai and θ.sub.Azimuth,Aj are the azimuth angles of the observed stars i and j before they enter the atmosphere, respectively.

(20) As shown in (a) and (b) of FIG. 2, according to the angle between the atmospheric refraction surfaces of the three observed stars, the theoretical position coordinates B(x.sub.B, y.sub.B) of the earth's zenith direction on the imaging plane of the star sensor can be solved as:

(21) $\begin{array}{cc}\{\begin{array}{c}\mathrm{acos}\left(\frac{\stackrel{\_}{\mathrm{FA}1}.Math.\stackrel{\_}{\mathrm{FA}2}}{.Math.\stackrel{\_}{\mathrm{FA}1}.Math..Math.\stackrel{\_}{\mathrm{FA}2}.Math.}\right)=\left({\theta }_{\mathrm{Azimuth},A1}-{\theta }_{\mathrm{Azimuth},A2}\right)\\ \mathrm{acos}\left(\frac{\stackrel{\_}{\mathrm{FA}1}.Math.\stackrel{\_}{\mathrm{FA}3}}{.Math.\stackrel{\_}{\mathrm{FA}1}.Math..Math.\stackrel{\_}{\mathrm{FA}3}.Math.}\right)=\left({\theta }_{\mathrm{Azimuth},A1}-{\theta }_{\mathrm{Azimuth},A3}\right)\\ \mathrm{acos}\left(\frac{\stackrel{\_}{\mathrm{FA}2}.Math.\stackrel{\_}{\mathrm{FA}3}}{.Math.\stackrel{\_}{\mathrm{FA}2}.Math..Math.\stackrel{\_}{\mathrm{FA}3}.Math.}\right)=\left({\theta }_{\mathrm{Azimuth},A2}-{\theta }_{\mathrm{Azimuth},A3}\right)\end{array}& \left(6\right)\end{array}$

(22) In the formula, (θ.sub.Azimuth,A1, θ.sub.Azimuth,A2 and θ.sub.Azimuth,A3 are the azimuth angles of the observed stars A1, A2 and A3 before they enter the atmosphere, respectively.

(23) In theory, any three observed stars in the star sensor form a refracting surface triangle, and their calculated results should all be the same. However, the actual imaging process of the star sensor not only includes the atmospheric refraction error, but also includes the random noise error of the star sensor positioning system error and the atmospheric turbulence error. Therefore, it needs to be optimized to solve:

(24) $\begin{array}{cc}\left(x_B,y_B\right)=\arg \min \underset{i\ne j=1}{\overset{3}{.Math.}}{.Math.\mathrm{acos}\left(\frac{\stackrel{\_}{\mathrm{FAi}}.Math.\stackrel{\_}{\mathrm{FAj}}}{.Math.\stackrel{\_}{\mathrm{FAj}}.Math..Math.\stackrel{\_}{\mathrm{FAj}}.Math.}\right)-\left({\theta }_{\mathrm{Azimuth},\mathrm{Ai}}-{\theta }_{\mathrm{Azimuth},\mathrm{Aj}}\right).Math.}_2^2& \left(7\right)\end{array}$

(25) On the other hand, because the observed stars in different refracting surface triangles are affected by the random noise error of the star sensor positioning system error and the atmospheric turbulence error, the calculation results are different. The observed stars identified after the identification of N star images form C.sub.N.sup.3 groups of refracting surface triangles, and C.sub.N.sup.3 groups of calculation results are obtained. In addition to most of the results clustered within a small range of variation, there are also some singular values. However, for a star image, the observed stars are imaged at the same time, so the zenith direction of all observed stars in the entire star image is the same.

(26) Therefore, the optimal solution of the refracting surface triangle composed of all the observed stars in the star image can be obtained, so that the theoretical position coordinates B of the optimal Earth zenith direction of the entire star image on the imaging surface of the star sensor can be obtained:

(27) $\begin{array}{cc}\left(\stackrel{\_}{x_B},\stackrel{\_}{y_B}\right)=\arg \min \underset{n=1}{\overset{M}{.Math.}}\left(\underset{\mathrm{ni}\ne \mathrm{nj}=1}{\overset{3}{.Math.}}.Math.\mathrm{acos}\left(\frac{\stackrel{\_}{{\mathrm{FA}}_{\mathrm{ni}}}.Math.\stackrel{\_}{{\mathrm{FA}}_{\mathrm{nj}}}}{.Math.\stackrel{\_}{{\mathrm{FA}}_{\mathrm{ni}}}.Math..Math.\stackrel{\_}{{\mathrm{FA}}_{\mathrm{nj}}}.Math.}\right)-\left({\theta }_{\mathrm{Azimuth},\mathrm{Ani}}-{\theta }_{\mathrm{Azimuth},\mathrm{Anj}}\right).Math.\right)& \left(8\right)\end{array}$

(28) In the formula, n.sub.i and n.sub.j are the numbers of the observed stars in the refraction surface triangle, and n is the number of the refraction surface triangle formed by the observed stars identified after the star recognition. If there are N observed stars recognized after the star recognition, t, the number of refracting surface triangles M=C.sub.N.sup.3.

(29) According to Equation (7) and Equation (8), FIG. 3 shows the 16 observed stars in a frame of measured star image after star identification, and obtained the zenith direction of the Earth with different refraction surface triangles on the imaging surface of the star sensor. The theoretical position coordinate B x.sub.B, y.sub.B) is denoted as ‘x’, and the theoretical position coordinate B(x.sub.B, y.sub.B) of the optimal Earth zenith direction of the entire star image on the imaging plane of the star sensor is denoted as ‘o’.

(30) Then, the atmospheric refraction corrections are made for all observed stars in the star sensor. Since the azimuth angle of the incident stellar does not change before and after atmospheric refraction, the azimuth angle of the observed star before and after atmospheric refraction correction is the same. As shown in FIG. 4, the position coordinates A.sub.i(x.sub.A, y.sub.A) of the i-th observed star on the star sensor imaging plane, the position coordinates A′.sub.i(x′.sub.A, y′.sub.A) after atmospheric refraction correction and the theoretical position coordinates B(x.sub.B, y.sub.B) of the earth's zenith direction on the star sensor imaging plane are collinear with three points, which is:
BA.sub.i.Math.BA′.sub.i=(x.sub.A−x.sub.B)(y′.sub.A−y.sub.B)−(y.sub.A−y.sub.B)(x′−x.sub.B)=0  (9)

(31) According to formula (2), the zenith distance of the i-th observed star on the imaging surface of the star sensor before entering the atmosphere is θ.sub.Zenith,Ai, and the position A′.sub.i(x′.sub.A, y′.sub.A) after atmospheric refraction correction is obtained by solving the triangular cosine law:
|A.sub.iA′.sub.i|.sup.2=|O.sub.sA.sub.i|.sup.2+|O.sub.sA′.sub.i|.sup.2−2|O.sub.sA′.sub.i|cos(θ.sub.Zenith,Ai−∠BO.sub.sA.sub.i)  (10)

(32) In the formula, A.sub.iA′.sub.i is the distance between the coordinates A.sub.i(x.sub.A, y.sub.A) and A′.sub.i(x′.sub.A, y′.sub.A), O.sub.sA.sub.i and O.sub.sA′.sub.i are the distance from the coordinates A.sub.i(x.sub.A, y.sub.A) and A′.sub.i(x′.sub.A, y′.sub.A) to the origin of the star sensor coordinate system Os, respectively, and ∠BO.sub.sA.sub.i is the angle between the vector O.sub.sB and O.sub.sA.sub.i.

(33) Simultaneous equations (9) and (10) form a system of equations, and the theoretical position coordinates A′.sub.i(x′.sub.A, y′.sub.A) of the observed star i on the imaging plane of the star sensor are obtained after removing atmospheric refraction. Similarly, performing atmospheric refraction correction on other observed stars in the star sensor to obtain the theoretical position coordinates after removing atmospheric refraction. According to the optical parameters of the star sensor, that is, the pixel size, focal length and principal point, the observed star vector O.sub.sA.sub.i of the observed star in the star sensor coordinate system after atmospheric refraction correction is obtained:
O.sub.sA′.sub.i=[x′.sub.Ai−x.sub.0,u′.sub.Ai−y.sub.0,−f]  (11)

(34) By the star identification method, the observed star A.sub.i in the star sensor corresponds to the reference star V.sub.i in the reference star catalog. Therefore, using the observed star vector corrected by atmospheric refraction and its corresponding reference star vector in the reference star catalog, the attitude transformation matrix C.sub.s.sup.m of the star sensor in the geographic coordinate system is calculated:
O.sub.sA′.sub.i=C.sub.s.sup.mV.sub.i(θ.sub.Zenith,Ai,θ.sub.Azimuth,Ai)  (12)

(35) In the formula, V.sub.i(θ.sub.Zenith,Ai, θ.sub.Azimuth,Ai) is the position vector of the reference star corresponding to the observed star in the geographic coordinate system.

(36) 2. Embodiment Effect

(37) The invention is applied to the stellar atmospheric refraction correction of the star sensor in the atmosphere, and realizes the real-time correction of the position coordinates of all identified stars during the observation process of the star sensor, thereby obtaining accurate attitude information of the star sensor.

(38) In order to verify the correctness and effectiveness of the present invention, an outfield night-sky test was carried out at the Xinglong National Astronomical Observatory in Hebei. The specific content of the test is as follows:

(39) Continuously collect 1000 frames of star images, and perform star identification on the centroid positioning results of each frame of star images to obtain the reference stars corresponding to the observed stars in the reference star list. Then, use the atmospheric refraction method proposed in the present invention to perform atmospheric refraction correction on the identified stars. In order to verify the atmospheric refraction effect of the present invention, the average value of the observed star angular distance and the reference star angular distance composed of all stars before and after the atmospheric refraction correction in each star image is compared.

(40) FIG. 5 shows the average angular distance deviation before and after atmospheric refraction correction for 1000 frames of experimental star images. FIG. 5(a) is the average angular distance deviation before atmospheric refraction correction, the standard deviation is 1.2713 arc seconds, and the mean is 22.483 arc seconds. FIG. 5(b) is the average angular distance deviation after atmospheric refraction correction, the standard deviation is 0.4643 arc seconds, and the mean is 1.7346 arc seconds.

(41) The results show that the stellar atmospheric refraction correction method based on the collinear refraction surface proposed in the present invention is effective and has high atmospheric refraction correction accuracy, which can effectively improve the attitude measurement accuracy of the ground observation star sensor.

(42) The above embodiments are provided for the purpose of describing the present invention only, and are not intended to limit the scope of the present invention. The scope of the invention is defined by the appended claims. Various equivalent replacements and modifications made without departing from the spirit and principle of the present invention should be included within the scope of the present invention.