The Satellite’s On-Orbit Attitude System Error Compensation Technique Based on Stereo Models

According to the problem that the low measurement accuracy of TH-1 satellite star sensor, the low frequency and “slow drift” error which cannot be ig-nored in the attitude determination system, resulting in obvious random error in the horizontal position and elevation direction, and the change of the error with time and latitude, cannot be calibrated by the ground field of the real problem. In this paper, a low frequency detection model is established by using the principle of relative orientation, and the low frequency error is obtained by parallax elimination. Finally, the satellite attitude is compensated and the more accurate exterior orientation elements are obtained, thus im-proving the positioning accuracy and stability. The experimental results show that: the proposed methods are feasible, and by using the model to dynamically calibrate the exterior orientation angle elements on orbit, the plane and elevation errors of the ground points can be basically eliminated. The global uncontrollable positioning accuracy and stability of the photogrammetry satellite are improved.


Introduction
It is an important problem that high-precision positioning accuracy without GCPs, which is plaguing the acquisition of high-precision geographic information products around the world [1] [2] [3]. Research status and engineering practices show that the high positioning accuracy of satellite imagery has little to do with image resolution [1]. The key is the acquisition of high-precision attitude data. If the external orientation elements of the satellite on-orbit detection are highly accurate and reliable, then it is simple to realize the global photogrammetry problem without GCPs [4] [5]. The geometric calibration technique is used to calibrate the error of the external orientation elements in the calibration field [6] [7]. In this way, the requirements for positioning without GCPs can be achieved. However, due to the low accuracy of satellite star sensors such as TH-1 and ZY-3, the low frequency errors and "slow drift" errors are not negligible, which will result in significant random errors in the horizontal and altitude directions [8] [9]. Different from the system error, the error of the attitude measurement system changes with time and latitude, resulting in its short-period systematization in the horizontal position direction [10]. The static calibration in orbit cannot give a stable solution. It causes great difficulties in positioning without ground control points (GCPs) [11].
The TH-1 satellite adopts the unique EFP (Equivalent Frame Photo) bundle adjustment, and the low-frequency error compensation technology is proposed on this basis, which realizes the global uncontrolled high-precision positioning, achieving the positioning accuracy of 3.7 m on plane and 2.4 m on elevation which meets the mapping accuracy requirement of 1:25000 geographic information product [12]. EFP bundle adjustment is based on the special three-line-array LMCCD aerospace surveying camera system; the low-frequency error compensation with this method is not reported in other satellite projects. In addition, in the domestic study on the positioning accuracy of ZY-3, a few GCPs are needed to fit the system error and improve the positioning accuracy [13].
In this paper, our work is spirited by above method, we use the vertical parallax between the stereo image pair to eliminate the system error of image's pitch angle and yaw angle, and then the relative orientation of two images is restored. The positioning accuracy of satellite images is improved by compensating low frequency error of the pitch angle and yaw angle. The proposed model is simple to be constructed and its transplantation to the general aerospace photogrammetry satellite positioning method is convenient.
The low-frequency error compensation model and compensation process are also given in this paper. The proposed method is tested by the actual data of TH-1 satellite, which proves its correctness and feasibility [14].

Low-Frequency Error Compensation Model for Exterior
Orientation Angle Elements

Fundamental Theory
By using the principle of relative orientation, the low frequency detection model is established, and the low frequency error of attitude is obtained by eliminating parallax. The satellite attitude is compensated to obtain more accurate foreign position elements and improve the accuracy of satellite positioning. As shown in Figure 1, 1 a and 2 a are a pair of homonymy light, homonymy light pair is located in the same epipolar plane as the photographic baseline B , that is, 1 a , 2 a and B three lines are coplanar. From the knowledge of spatial analytic geometry, if three lines are coplanar, the mixed product of their corresponding vectors is zero, shown as (1) [15].

The Basic Calculation Process
The coordinates of the three vectors in the image space assisted coordinate system are respectively ( ) , , , , , , X Y Z ，The coplanar condition equation can be represented by coordinates as (2).
Whether the coplanar condition equation is established is the standard to complete the relative orientation. The analytic relative orientation is to solve elements of relative orientation according to the coplanar condition equation. When relative directional parsing is calculated, the photographic baseline B is usually rewritten to b, which is called projection baseline. Here, where m is the photographic scale denominator; ( ) , ,   (4), µ is the deflection angles of the baseline (that is, the yaw corner of low frequency error) and ν is the inclination (that is, the pitch angle of low frequency error). Substituting the above formula into the coplanar condition equation.
F will be expanded by the Taylor series, holding the first order terms, and the linear formula of the coplanar equation is the derivation is finally available: the upper formula is converted into an error equation:

Provisions for Changes in the Inner Orientation Elements of Camera
Taking the TH-1 satellite as an example, the three-line-array CCD space photogrammetry camera adopts the form of three-lens three-line array, and the three-lens

The Correct Model of dµ, dv
In the error Equation (8), the vertical parallax of homonymy points originates not only from the system value produced by the low frequency error of the attitude measurement system, but also from the accidental value generated by the random error of the attitude measurement system. In the actual low frequency error compensation, using vertical parallax of all points participating in the adjustment as the basis, the random error can be weakened by the adjustment method (about 0.6 factors) and its size can be ignored, that is, dϕ , dω and dκ are about zero.
At the same time, centering on the central axis of the flight direction of strip image, match two rows of homonymy points on the upper and lower sides of the same strip image for constraining. Taking one pair of homonymy points in the upper and lower rows as an example, the vertical parallax can be converted to In the equation the symbol U and D respectively represents the upper and lower sides of the same strip image for constraining. Because of the upper and lower symmetry of the homonymy points, there is bring Equation (11) into the Equation (12) we can get 2 The vertical parallax of the homonymy points is further transformed into 1) The special correction of d µ available by Formula (14) (Figure 3) the correction value of the obtained d µ is 2) The special correction of dν available by Formula (16) (Figure 4): the correction value of the obtained dν is:

Low Frequency Error Compensation Process
The process of low frequency error compensation is as follows and the for details can be seen in Figure 5.
Step 1: Fully automatic match the three-line-array image, and obtain the image points coordinates of the upper and lower rows on the forward-looking, nadir-looking and backward-looking images.
Step 2: Coordinate transformation of the original pose data, and convert attitude data between CGCS2000 and earth observation camera to the local and camera coordinate system for the vertical parallax calculation.
Step 3: According to the relative orientation principle, calculate the parallax of the upper and lower rows points, and the correction amount corresponding to d µ of each row, and use the sieve strategy to eliminate the outlier to ensure the uniformity of the upper and lower parallax distribution of the upper and lower rows points. Thus, the average values of d µ and dν are obtained.
Step 4: Correct d µ and dν , and calculate the parallax again, check whether it meets the threshold.
Step 5: After the correction is completed, the pose data is converted to the data between CGCS2000 and the camera.
Step 6: Use the corrected attitude data for bundle adjustment and assess the positioning accuracy.

Sensitivity Analysis of Low Frequency Error Compensation
In this paper, we used five test sites of TH01-01 and TH01-02's image to validate our methods. After static on-orbit geometry inspection of five test sites of TH01-01 star and TH01-02 star, the system error is basically eliminated, that is, when each test site star sensor data is corrected by its own calibration parameters, the low frequency measurement error should be eliminated. Under the condition of tilt error compensation of three-line-array image plane, the low frequency error of five test sites should theoretically be 0. The low frequency error detection of TH01-01 star and TH01-02 star is obtained by means of dynamic calibration technique of exterior orientation angle elements, and the low-frequency error correction value is given as shown in Table 1. 1) The pitch angle and yaw angle error of each test site fluctuates around 0, showing randomness, through the error calculation, the sensitivity of low-frequency error compensation is nearly 1″, which is a high value.
2) The sensitivity of yaw angle is higher than the pitch angle, and its value is basically zero. Through the position method of the satellite images, the system error of the yaw angle will cause elevation accuracy error [14]. Therefore, the reduction of elevation accuracy caused by low-frequency measurement error can be basically eliminated by dynamic calibration.

Correction of Tilt Error of Third Line Array Image Plane
After static on-orbit geometry calibration, the real pitch angle and yaw angle correction values of the same test sites are obtained by the dynamic calibration technique of exterior orientation angle elements, which is shown in Table 2.
Y. Zhou et al. 1) The pitch angle error of TH01-01 star and TH01-02 star is systematic, the error of TH01-01 star's test sites fluctuates around the mean value 2.51", and the error of TH01-01 star's test sites revolves around the mean value 10.60".
2) The partial system error of the pitch angle of TH01-01 star and TH01-02 star is the tilt error of the three-line-array image plane.
3) After correcting the tilt error of the three-line-array image plane of TH01-01 star and TH01-02 star, the real low-frequency measurement error can be obtained by dynamic calibration, as shown in Table 2.

Low Frequency Compensation Performance Verification
After dynamic calibration of the test site data, the measurement error of the star sensor is corrected, and the positioning without GCPs is carried out. The positioning accuracy results of TH01-01 star and TH01-02 star are shown as Tables 3-5.
Comparing the experimental results above, we can summarize that:   2) Correct the low-frequency measurement error by dynamically calibration, the positioning accuracy of plane and elevation is improved, especially the elevation accuracy can be improved to the accuracy of each test site corrected by the calibration parameters of each test site.
3) Correct low-frequency measurement error by dynamic calibration, and the plane positioning accuracy is still lower than the accuracy of each test sites corrected by the calibration parameters of each test site.
Y. Zhou et al.

Conclusion
Through the technique of exterior orientation angle elements dynamic calibration, the low-frequency error detection of TH01-01 star and TH01-02 star can basically eliminate the plane and elevation error of the ground points, and the sensitivity of the low-frequency error compensation is about 1", which is a high sensitivity. Under the condition of low precision of domestic satellite attitude measurement system, it is of great significance to the improvement of global positioning accuracy without GCPs.