METHOD FOR JOINT ESTIMATION OF STELLAR ATMOSPHERIC REFRACTION AND STAR TRACKER ATITUDE

Abstract

Disclosed is a method for joint estimation of stellar atmospheric refraction and star tracker attitude, including: capturing, by a star tracker, an image of stars, and recognizing the image to obtain a matching relationship between an observed star and a reference star; transforming the reference star to a terrestrial reference frame based on time and position of observation to obtain a true zenith distance of the reference star; calculating an estimated stellar atmospheric refraction based on a simplified atmospheric refraction model, the true zenith distance, and an initial atmospheric refractive coefficient, and subjecting the reference star to stellar atmospheric refraction compensation; calculating the star tracker attitude, and re-projecting the observed star to the terrestrial reference frame to calculate the observed stellar atmospheric refraction and the stellar atmospheric refraction error; if the observed stellar atmospheric refraction is misaligned with the estimated stellar atmospheric refraction, adjusting the atmospheric refractive coefficient based on the stellar atmospheric refraction error, and then recompensating the reference star and calculating the attitude, till the observed stellar atmospheric refraction is aligned with the estimated stellar atmospheric refraction, whereby joint estimation results of the stellar atmospheric refraction and the star tracker attitude are obtained. The disclosure realizes real-time, autonomous estimation and cancelation of stellar atmospheric refraction without being limited by external sensors or empirical formulae.

Claims

1. A method for joint estimation of stellar atmospheric refraction and star tracker attitude, comprising: capturing, by a star tracker, an image of stars, and recognizing the image to obtain a matching relationship between an observed star and a reference star; transforming the reference star to a terrestrial reference frame based on time and position of observation to obtain a true zenith distance of the reference star; obtaining a stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame; obtaining an attitude cosine matrix from a star tracker coordinate system to the terrestrial reference frame; obtaining a stellar atmospheric refraction error based on the attitude cosine matrix; subjecting the stellar atmospheric refraction error to linear least squares fitting to obtain a fit slope and a fit intercept, wherein if the absolute value of the fit slope is greater than a cutoff condition, an estimated atmospheric refractive coefficient is superimposed with a feedback gain that is positively correlated to the fit slope, and the method returns to obtaining a stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame; if the absolute value of the fit slope is smaller than the cutoff condition, iteration ends, whereby an optimal atmospheric refractive coefficient and an atmospheric refraction error compensated star tracker attitude are obtained.

2. The method of claim 1, wherein the capturing, by a star tracker, an image of stars, and recognizing the image to obtain a matching relationship between an observed star and a reference star comprises: obtaining, by the star tracker, the image of the stars to acquire coordinates of the stars on an imaging plane, whereby all observed stellar vectors in the star tracker coordinate system are obtained; and performing image recognition to obtain the matching relationship between the observed star and the reference star on a navigational star catalogue.

3. The method of claim 2, wherein the transforming the reference star to a terrestrial reference frame based on time and position of observation to obtain a true zenith distance of the reference star comprises: transforming the reference star from a celestial reference frame adopted in the navigational star catalogue to the terrestrial reference frame at the position of observation of the star tracker to obtain a reference stellar vector free of stellar atmospheric refraction in the terrestrial reference frame, whereby a true zenith distance of the reference star is further obtained.

4. The method of claim 3, wherein the obtaining a stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame comprises: calculating an estimated stellar atmospheric refraction based on a simplified atmospheric refraction model, the true zenith distance, and an initial atmospheric refractive coefficient, and subjecting the reference star to stellar atmospheric refraction compensation, whereby a stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame is obtained.

5. The method of claim 4, wherein the obtaining an attitude cosine matrix from a star tracker coordinate system to the terrestrial reference frame comprises: calculating an attitude quaternion from the star tracker coordinate system to the terrestrial reference frame by substituting the observed stellar vector in the star tracker coordinate system and the stellar atmospheric refraction compensated reference stellar vector in a QUEST algorithm, whereby the attitude cosine matrix from the star tracker coordinate system to the terrestrial coordinate system is further obtained.

6. The method of claim 4, wherein the subjecting the reference star to stellar atmospheric refraction compensation refers to subtracting corresponding estimated stellar atmospheric refraction from the true zenith distance of the reference star to obtain a new reference stellar vector.

7. The method of claim 5, wherein the obtaining a stellar atmospheric refraction error based on the attitude cosine matrix comprises: transforming the observed stellar vector in the star tracker coordinate system to the terrestrial coordinate system via the attitude cosine matrix, comparing the observed stellar vector with the reference star in the terrestrial coordinate system to calculate a magnitude of an observed stellar atmospheric refraction, and finding a difference between the observed stellar atmospheric refraction and the estimated stellar atmospheric refraction, whereby the stellar atmospheric refraction error is obtained.

8. The method of claim 7, wherein the stellar atmospheric refraction error is subjected to linear least squares fitting to obtain a fit slope and a fit intercept.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0026] FIG. 1 is a flow diagram of a method for joint estimation of stellar atmospheric refraction and star tracker attitude according to the present disclosure;

[0027] FIG. 2 is a diagram of a simplified starlight atmospheric refraction model according to the present disclosure;

[0028] FIG. 3 is a schematic diagram of an atmospheric refractive coefficient feedback mechanism according to the present disclosure;

[0029] FIG. 4 is a diagram of test measurements of 1000 continuous frames of star images captured by the star tracker at a turntable tilt angle of 20° according to the present disclosure;

[0030] FIG. 5 is an effect diagram of stellar atmospheric refraction estimation for single-frame measurement at a turntable tilt angle of 20° according to the present disclosure.

DETAILED DESCRIPTION

[0031] To make the objectives, technical solutions, and advantages of the present disclosure more apparent, the present disclosure will be described in further detail through specific implementations with reference to the drawings. It should be understood that the depictions are only exemplary, not intended for limiting the scope of the present disclosure. In addition, depictions about known structures and technologies is omitted in the description below so as to avoid unnecessary confusion to the concepts of the present disclosure.

[0032] The present disclosure provides a method for joint estimation of stellar atmospheric refraction and star tracker attitude, specifically comprising:

[0033] step S100: capturing, by a star tracker, an image of stars, and recognizing the image to obtain a matching relationship between an observed star and a reference star.

[0034] Specifically, the star tracker captures an image of stars to acquire coordinates of the stars on the imaging plane to obtain all observed stellar vectors in the star tracker coordinate system, and then performs star image recognition to obtain a matching relationship between the observed star and the reference star on a navigational stellar catalogue.

[0035] Specifically, all observed stellar vectors in the star tracker coordinate system are obtained based on internal parameters of the star tracker, a distortion model, and a camera imaging model.

[0036] Specifically, the star tracker captures an image of the observed star to obtain a series of observed vectors W.sub.s in the star tracker coordinate system, where W.sub.s = [W.sub.s,1, W.sub.s,2, W.sub.s,3, ..., W.sub.s,n].

[0037] The star image may be recognized using algorithms such as a subgraph isomorphism algorithm (e.g., triangulation algorithm, pyramid algorithm, etc.), a pattern recognition algorithm (e.g., grid algorithm, recognition algorithm based on the radial companion-circumferential feature, etc.); and other types of recognition algorithms (e.g., neural network-based star image recognition algorithm) are also allowed. The star image recognition algorithm is not unique, and an appropriate algorithm is selected dependent on actual conditions. Irrespective of which algorithm is selected, it needs to ensure correctness of the recognition result. The information about the reference star includes right ascension, declination, proper motion, parallax, and radial velocity.

[0038] Specifically, the star tracker captures an image of observed stars to acquire coordinates [X, Y] of the stars on the imaging plane, where = [X.sub.1,X.sub.2,X.sub.3,•••, X.sub.n].sup.T, Y = (Y.sub.1, Y.sub.2, Y.sub.3, •••, Y.sub.n].sup.T. the observed star vectors W.sub.s = [w.sub.s,.sub.1, w.sub.s,.sub.2, w.sub.s,.sub.3, •••, w.sub.s,.sub.n] in the star tracker coordinate system are obtained based on a lens distortion model, a camera imaging model, and internal parameters of the star tracker (e.g., principle point X.sub.0, Y.sub.0, focal length ƒ, distortion parameters p.sub.1, p.sub.2, q.sub.1, q.sub.2, and pixel sizes d.sub.x, d.sub.y), where:

[00001]$w_{s,i}=\frac{1}{\sqrt{{\widehat{x}}_i{}^2+{\widehat{y}}_i{}^2+f^2}}\left[\begin{array}{c}-{\widehat{x}}_i\\ -{\widehat{y}}_i\\ f\end{array}\right]$

where:

[00002]$\begin{array}{c}{\widehat{x}}_i=x_i\left(1+q_1{r}^2+q_2{r}^4\right)+2p_1{x}_i{y}_i+p_2\left(r^2+2x_i{}^2\right)\\ {\widehat{y}}_i=y_i\left(1+q_1{r}^2+q_2{r}^4\right)+p_1\left(r^2+2y_i{}^2\right)+2p_2{x}_i{y}_i\\ x_i=\left(X_i-X_0\right)dx\\ y_i=\left(Y_i-Y_0\right)dy\\ r=\sqrt{x_i{}^2+y_i{}^2}\end{array}$

[0039] Then, the reference stars V.sub.r = [v.sub.r,1, v.sub.r,2, v.sub.r,3,•••, v.sub.r,n] matched with the observed stars on the stellar catalogue are obtained by performing star image recognition using a modified triangular algorithm, where v.sub.r,i = [RA.sub.i, DE.sub.i, pmRA.sub.i, pmDE.sub.i, Plx.sub.i, rv.sub.i].sup.T, where RA.sub.i denotes the right ascension, DE.sub.i denotes the declination, pmEA.sub.i, denotes the longitude component of proper motion, pmDE.sub.i denotes the latitude component of the proper motion, PLx.sub.i, denotes the parallax, and rv.sub.i denotes the radial velocity.

[0040] Step 200: transforming the reference star to a terrestrial reference frame based on time and position of observation to obtain a true zenith distance of the reference star.

[0041] Specifically, the reference star obtained from step S100 is transformed from the celestial reference system adopted by the navigational stellar catalogue to the terrestrial reference frame at the observation position of star tracker, obtaining reference stellar vectors free of stellar atmospheric refraction in the terrestrial reference frame, whereby the true zenith distance of the reference star is obtained.

[0042] Specifically, the reference star obtained from step S100 is transformed from the celestial reference system adopted by the navigational stellar catalogue to the terrestrial reference frame at the observation position of the star tracker using a high-precision celestial reference system transformation chain based on the imaging time and position of the star tracker and the Earth’s rotation information.

[0043] Specifically, in transformation of the reference star’s coordinate system, the star tracker imaging time refers to the exposure center time, which may be represented by UTC time on the scale of milliseconds; and the position refers to the longitude, latitude, and altitude of the star tracker at the exposure center time, which may be represented by a WGS84 ellipsoid reference coordinate system, precise to a scale of seconds and meters.

[0044] Specifically, the Earth’s rotation information includes the Earth’s polar wander and DUT 1 of the day, which may be provided by IERS (International Earth Rotation Service). To perform reference star coordinate system transformation based on the high-precision coordinate system transformation chain provided by SOFA (Standards Of Fundamental Astronomy), it further requires acquisition of parameters such as the right ascension, declination, proper motion, parallax, and radial velocity of the reference star on the navigation stellar catalogue.

[0045] Specifically, the observation site’s terrestrial reference frame refers to the East-North-Up celestial reference system, wherein the following in the terrestrial reference frame are obtained: reference vector: V.sub.g = [v.sub.g,1, v.sub.g,2, v.sub.g,3, ..., v.sub.g,n], zenith vector U = [0,0,1].sup.T, true zenith distance Z = [Z.sub.1,Z.sub.2,Z.sub.3, •••,Z.sub.n].sup.T, where:

[00003]$Z_i=\text{arccos}\left(U\cdot v_{g,i}\right)$

[0046] Specifically, the UTC time T for the star tracker exposure center time, the longitude-latitude-altitude coordinates P (WGS84 elliptic reference coordinate system) where the star tracker is located, the Earth’s pole wander Pm = [px, py] and dut1 of the day are obtained, among which pm and dut1 may be provided by EIRS (https://iers.org). Then, by transforming the reference star from the celestial reference system adopted by the navigational stellar catalogue to the terrestrial reference frame of the star tracker’s observation site using the high-precision coordinate system transformation chain provided by IAU SOFA (http://www.iausofa.org/), the reference vector V.sub.g = [V.sub.g,1, v.sub.g,2, v.sub.g,3, •••, v.sub.g,n] in the terrestrial reference frame is given as such:

[00004]$V_g=iauCST\left(V_r,T,P,Pm,dut1\right)$

where iau(*) denotes using the iauSOFA function to transform the reference stellar coordinate from the celestial reference system to the terrestrial reference frame, Vr denotes the coordinate in the reference stellar celestial reference system, T denotes the observation time, P denotes the observation position, Pm denotes the Earth’s pole wander of the day, and dut1 denotes the difference constant between UT1 time and UTC time.

[0047] Given the zenith vector U = [0,0,1].sup.T, the true zenith distances Z = [Z.sub.1, Z.sub.2, Z.sub.3, •••, Z.sub.n].sup.T of all reference stars are calculated according to the equation below:

[00005]$Z_i=\text{arccos}\left(U\cdot v_{g,i}\right)$

[0048] Step S300: obtaining a stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame.

[0049] Specifically, a stellar atmospheric refraction calculation formula is obtained based on a simplified starlight atmospheric refraction model, as shown in FIG. 2. The estimated atmospheric refractions of all reference stars are calculated by substituting the estimated atmospheric refractive coefficient and the true zenith distances of all reference stars into the atmospheric refraction calculation formula, and the reference stellar vectors are subjected to stellar atmospheric refraction correction along the elevation angle incremental direction, whereby the stellar atmospheric refraction compensated reference stellar vectors in the terrestrial reference frame are obtained.

[0050] Specifically, the simplified starlight atmospheric refraction model refers to regarding the Earth’s atmospher as being comprised of innumerable concentric spheres, where the atmospheric density and the refractive index increase as the altitude decreases, and the atmosphere in the field of view (FOV) of the star tracker may be regarded as comprising 1~n planar layers from the observation site to the interstellar space irrespective of the curvature of the atmospheric sphere, the incident angle of the starlight outside the atmosphere is referred to as the true zenith distance Z; refraction occurs when the starlight propagates from one layer to the next layer, and the starlight curves convexly toward the vertical normal to the observation site. Supposing the atmospheric refractive index at the outermost layer is n.sub.n, the following is given by the refraction law:

[00006]$sin\left(Z\right)=n_{n}sin\left(Z_n\right)$

where Z.sub.n stands for the true zenith distance of the stellar vector at the outermost layer of atmosphere, and the adjacent next layer is given as follows:

[00007]$n_n\text{sin}\left(Z_n\right)=n_{n-1}\text{sin}\left(Z_{n-1}\right)$

[0051] Till the lowest layer which is given as follows:

[00008]$n_2\text{sin}\left(Z_2\right)=n_1\text{sin}\left(Z_1\right)$

[0052] The above equations can be combined into the following:

[00009]$\text{sin}\left(Z\right)=n_1\text{sin}\left(Z_1\right)$

where Z.sub.1 stands for the true zenith distance when the stellar vector is propagated to the lowest layer; from the perspective of the observation site, it is the apparent zenith distance, represented by Z′, where Z′ < Z. The stellar atmospheric refraction ρ is defined as a difference between the true zenith distance Z and the apparent zenith distance Z′:

[00010]$\rho =Z-Z^{\prime }$

so, the equation (9) may be rewritten as follows:

[00011]$\text{sin}\left(Z\right)=n_1\text{sin}\left(Z-\rho \right)=n_1\sin \left(Z\right)\cos \left(\rho \right)-n_1\cos \left(Z\right)\sin \left(\rho \right)$

since ρ is very small, taking cos(ρ) = 1, sin(ρ) = ρ, then the atmospheric refraction ρ is given as follows:

[00012]$\rho =\frac{n_1-1}{n_1}\text{tan}\left(Z\right)=Ntan\left(Z\right)$

where N denotes the atmospheric refractive coefficient, which is on the scale of arc-second. When N.sub.1 takes the standard atmospheric refractive index 1.000292, N0=60.21″. The initial value of the atmospheric refractive coefficient N may be set to standard atmospheric coefficient N.sub.0 = 60.21. According to the stellar atmospheric refraction calculation equation (12), the estimated stellar atmospheric refractive coefficient N and the true zenith distances Z of all reference stars may be substituted to calculate the estimated atmospheric refractions P = [ρ̂.sub.1, ρ̂.sub.2, ρ̂.sub.3, •••,ρ̂.sub.n].sup.T of all reference stellar vectors, where:

[00013]${\widehat{\rho }}_i=\widehat{N}\tan \left(Z_i\right)$

[0053] Then, the reference stellar vectors are subjected to stellar atmospheric refraction correction along the elevation angle incremental direction, obtaining the stellar atmospheric refraction compensated reference stellar vector V.sub.gc = [v.sub.gc,1, v.sub.gc,2, v.sub.gc,3,•••, v.sub.gc,n] in the terrestrial reference frame, where:

[00014]$v_{gc,i}=\left[cos{\widehat{\rho }}_{i}I+\left(1-cos{\widehat{\rho }}_i\right)n_i{n}_i{}^T+sin{\widehat{\rho }}_i{n}_i{}^{\wedge }\right]v_{g,i}$

where I denotes the unit matrix, n.sub.i denotes the rotation vector when the reference stellar vector v.sub.g,i increments by elevation angle ρ̂.sub.i:

[00015]$n_i=v_{g,i}\times U$

n.sub.i.sup.∧ denotes the antisymmetric matrix composed of vector elements, and n.sub.i(1), n.sub.i(2), and n.sub.i(3) represent the 1.sup.st, 2.sup.nd, and 3.sup.rd elements of the vector n.sub.i:

[00016]$n_i{}^{\wedge }=\left[\begin{array}{ccc}0& -n_i\left(3\right)& n_i\left(2\right)\\ n_i\left(3\right)& 0& -n_i\left(1\right)\\ -n_i\left(2\right)& n_i\left(1\right)& 0\end{array}\right]$

[0054] Step S400: obtaining an attitude cosine matrix from the star tracker coordinate system to the terrestrial reference frame.

[0055] Specifically, the observed stellar vector in the star tracker coordinate system and the atmospheric refraction compensated reference stellar vector in the terrestrial reference frame are substituted into a QUEST algorithm to calculate an attitude quaternion from the star tracker coordinate system to the terrestrial reference frame, whereby the attitude cosine matrix from the star tracker coordinate system to the terrestrial reference frame is obtained.

[0056] Specifically, the observed stellar vector in the star tracker coordinate system and the stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame are substituted into the QUEST algorithm to calculate the attitude quaternion from the star tracker coordinate system to the terrestrial reference frame, where the attitude quaternion is given as follows:

[00017]$Q_{out}=QUEST\left(W_s,V_{gc}\right)$

[0057] where QUEST(*) stands for the QUEST algorithm, W.sub.s denotes the observed stellar vector in the star tracker coordinate system, and V.sub.gc refers to the stellar atmospheric refraction compensated reference stellar vector in the terrestrial reference frame; [0058] in the equation, Q.sub.out = [v, q.sub.4], where v denotes the imaginary part of the quaternion, v = [q.sub.1, q.sub.2, q.sub.3], q.sub.1, q.sub.2, q.sub.3 being elements of the imaginary part of the quaternion, and q.sub.4 denotes the real part of the quaternion; then the star tracker attitude cosine matrix is given as follows:

[0059] Step S500: obtaining a stellar atmospheric refraction error based on the attitude cosine matrix;

[0060] Specifically, the observed stellar vector in the star tracker coordinate system is transformed to the terrestrial reference frame via the attitude cosine matrix and compared with the reference star in the terrestrial reference frame to calculate the magnitude of the observed stellar atmospheric refraction; and then a difference between the observed stellar atmospheric refraction and the estimated stellar atmospheric refraction is solved, whereby the stellar atmospheric refraction error is obtained.

[0061] Specifically, the observed stellar vector in the star tracker coordinate system is transformed to the terrestrial reference frame, obtaining W.sub.g = [w.sub.g,1, w.sub.g,2, w.sub.g,n]:

[00019]$W_g=R_s^g{W}_s$

[0062] Then, the W.sub.g is compared with the reference star in the terrestrial reference frame to calculate the magnitude of the stellar atmospheric refraction P = [ρ́.sub.1, ρ́.sub.2, ρ́.sub.3, •••,ρ́.sub.n].sup.T, where:

[00020]${{\rho }^{\prime }}_i=Z_i-\text{arccos}\left(U\cdot w_{g,i}\right)$

[0063] And then, a difference between the observed stellar atmospheric refraction and the estimated stellar atmospheric refraction is found to obtain the stellar atmospheric refraction error as follows:

[00021]$\Delta \text{P =}{\text{P}}^{\prime }-\widehat{\text{P}}$

[0064] Step S600: subjecting the stellar atmospheric refraction error to linear least squares fitting to obtain a fit slope and a fit intercept; if the absolute value of the fit slope is greater than the cutoff condition, the estimated atmospheric refractive coefficient is superimposed with a feedback gain, the feedback gain being positively correlated to the fit slope, and then the process returns to step S300; if the absolute value of the fit slope is less than the cutoff condition, the iteration is terminated, whereby the optimal atmospheric refractive coefficient and the atmospheric refraction error compensated star tracker attitude are obtained.

[0065] Specifically, the stellar atmospheric refraction error is subjected to linear least squares fitting based on the atmospheric refractive coefficient feedback mechanism to obtain the fit slope and the fit intercept.

[0066] The atmospheric refractive coefficient feedback mechanism goes like this: according to the characteristics of the simplified starlight atmospheric refraction model and the QUEST algorithm least squares, when N< N*(N*is the true value): the elevation angle of the stellar atmospheric refraction compensated reference star is smaller than the true position, with the elevation angle distribution range being greater than that of the observed star, such that the reference star is misaligned with the observed star; after attitude determination by QUEST, the boresight pointing is lower than the true pointing, which indicates an erroneous estimated attitude: however, since the stellar atmospheric refraction error is proportional to the zenith distance and the slope is greater than 0, the more N deviates from N*, the greater the slope of the atmospheric refraction; when N = N*: the elevation angle of the stellar atmospheric refraction compensated reference star is equal to the true position, indicating that the elevation angle distribution range is aligned the observed star; after attitude determination by QUEST, the boresight pointing is aligned with the true pointing, indicating that there is no error in attitude estimation; the stellar atmospheric refraction error is 0, and the slope is 0; when N > N*: the elevation angle of the stellar atmospheric refraction compensated reference star is greater than the true position, the elevation angle distribution range is smaller than that of the observed star, indicating that the reference star is misaligned with the observed star; after attitude determining by QUEST, the boresight pointing is higher than the true pointing, indicating that the attitude estimation has an error; however, since the stellar atmospheric refraction error is inversely proportional to the zenith distance and the slope is less than 0, the more N deviates from N*, the smaller the slope of the stellar atmospheric refraction error. Thus, the atmospheric refractive coefficient feedback mechanism is constructed, and the stellar atmospheric refraction error linear parameter is fitted using the least squares:

[00022]${\left[k,b\right]}^T={\left(M^{T}M\right)}^{-1}{M}^T\Delta \text{P}$

where M = [Z,L],L = [1,1,1,•••,1].sub.1×n.sup.T. If |k| > ε (where ε denotes the termination condition, close to 0), it indicates that N does not reach the optimal and needs to be updated; then gain G(k) is superimposed to N, and the process returns to step S300, where G(k) and k are positively correlated; to achieve rapid convergence of N, a PID gain may be adopted:

[00023]$\text{a}\widehat{N}=\widehat{N}+G\left(k\right)=\widehat{N}+g_p{k}_i+g_I{\displaystyle {.Math.}_{i=1}^n{k}_i+g_D}\left(k_i-k_{i-1}\right)$

where gP, gI,and gD are the proportional, integral, and derivative coefficients, respectively; .sub.ki is the slope of current stellar atmospheric refraction error, k.sub.i-1 is the previous value. If |k| < ε, it may be believed that N reaches the optimal; now, the optimal atmospheric refractive coefficient is obtained, and the atmospheric refraction error compensated star tracker attitude is also obtained.

Implementation Effect

[0067] To verify the effectiveness of the method for estimation of star tracker attitude provided by the present disclosure, a star observation field experiment was conducted. Before the experiment, the star tracker was mounted on a turntable, where the initial tilt angle of the turntable was 0°, and the installation plane normal pointed approximately towards the zenithal direction; the turntable was turned 6 times, about 5° for each turn; after the turntable became stable, data acquisition started, where seven groups of continuous star images were acquired in total, each including 1000 frames. The method provided by the present disclosure was adopted for attitude estimation, and the mean angular distance error (reprojection error, ARE) between the observed vector projection and the reference vector in the terrestrial reference frame was adopted as the evaluation criteron:

[00024]$\Delta d=\frac{1}{\text{b}}{\displaystyle \underset{i=1}{\overset{\text{b}}{.Math.}}arccos\left(w_{g,i}\cdot v_{g,i}\right)}$

[00025]$\Delta D=\frac{1}{\text{B}}{\displaystyle \underset{j=1}{\overset{\text{B}}{.Math.}}\Delta d_j}$

[0068] b represents the number of all observed stars in one frame of star image; and B represents the total number of star images.

[0069] Results of comparison between attitude measurements before and after correction are shown in Table 1, which sets forth the attitude estimation effects of respective groups when the turntable tilt angle was at 0°, 5°, ..., 30°. FIG. 4 shows the comparison results of 1000 test frames when the turntable tilt angle was at 20°, and FIG. 5 shows the stellar atmospheric refraction estimations of one frame thereof; Table 2 shows an iteration process of single-frame measurement in terms of atmospheric refractive coefficient of the star tracker when the turntable tilt angle was at 20°.

TABLE-US-00001 Turntable tilt angle (°) Attitude Euler angle (°) ARE (”) Standard deviation of ARE (”) Before correction After correction Before correction After correction Before correction After correction 0 86.6256 86.6258 22.7287 4.6154 1.3695 0.8484 -5.4339 -5.4324 0.7364 0.7348 5 86.6239 86.6243 28.0878 4.7634 1.4643 0.7284 -5.4609 -5.4578 5.7518 5.7458 10 86.6326 86.6331 27.5967 5.7002 1.5795 0.704 -5.4846 -5.4803 10.6949 10.6861 15 86.6322 86.6325 27.7093 5.0529 1.3457 0.6372 -5.5128 -5.5087 15.856 15.8421 20 86.6338 86.6351 30.8921 6.0615 0.7465 0.5925 -5.5365 -5.5314 20.6404 20.6232 25 86.6309 86.635 31.113 6.6874 1.1355 0.7698 -5.563 -5.556 25.7288 25.705 30 86.6261 86.6307 28.5306 6.3751 1.0973 0.5871 -5.5854 -5.5815 30.7538 30.7239

TABLE-US-00002 i  k b 1 60.21 2.405 -44.6921 2 180.4603 0.12823 -1.7799 3 186.8719 0.0068465 0.50786 4 187.2142 0.00036557 0.63001 5 187.2325 1.952e-05 0.63653

[0070] In view of the above, the present disclosure provides a method for joint estimation of stellar atmospheric refraction and star tracker attitude, comprising: capturing, by a star tracker, an image of stars, recognizing the image to obtain a matching relationship between an observed star and a reference star; transforming the reference star to a terrestrial reference frame based on time and position of observation to obtain a true zenith distance of the reference star; calculating an estimated stellar atmospheric refraction using a simplified stellar atmospheric refraction model, the true zenith distance, and an initial stellar atmospheric refractive coefficient; subjecting the reference star to stellar atmospheric refraction compensation for the reference star; calculating the star tracker attitude, re-projecting the observed star to the terrestrial reference frame, and calculating an error between the observed stellar atmospheric refraction and the estimated stellar atmospheric refraction; adjusting the atmospheric refractive coefficient based on the stellar atmospheric refraction error if the observed stellar atmospheric refraction is misaligned with the estimated stellar atmospheric refraction, so as to re-compensate the reference star and calculate the attitude till the observed stellar atmospheric refraction is aligned with the estimated stellar atmospheric refraction, whereby results of the joint estimation of the stellar atmospheric refraction and the star tracker attitude are obtained. The disclosure realizes real-time, autonomous estimation and elimination of star tracker atmospheric refraction, avoids the limitations of empirical formulas, eliminates the reliance on additional sensors, such as thermometers, barometers, hygrometers, and attitude sensors, and reduces the difficulty and system overheads of a terrestrial star tracker in cancelling the atmospheric refraction; the present disclosure adopts a simplified starlight atmospheric refraction model to construct an atmospheric refractive coefficient feedback mechanism based on the star image recognition algorithm and the celestial reference system transformation chain in consideration of the characteristics of QUEST algorithm for attitude determination and the relationship between stellar atmospheric refraction and optimal attitude. The atmospheric refractive coefficient feedback mechanism enables real-time, optimal estimation of the atmospheric refractive coefficient while obtaining the terrestrial star tracker attitude free of stellar atmospheric refraction noise.

[0071] It is understood that the specific implementations described above are only for illustrating or explaining the principle of the disclosure through examples, not intended for constituting limitation to the disclosure. Accordingly, any modification, equivalent substitution, and improvement without departing from the spirits and scope of the disclosure falls into the scope sought for protection in the disclosure. In addition, the appended claims intend to cover the scope and boundary of the claims or all variations and modifications within equivalents of the scope and boundary.