^{1}

^{1}

^{2}

^{3}

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 ignored 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 improving 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.

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 [

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 [

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 [

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

B ⋅ ( a 1 × a 2 ) = 0 (1)

The coordinates of the three vectors in the image space assisted coordinate system are respectively ( B x , B y , B z ) , ( X 1 , Y 1 , Z 1 ) and ( X 2 , Y 2 , Z 2 ) ，The coplanar condition equation can be represented by coordinates as (2).

F = | B x B y B z X 1 Y 1 Z 1 X 2 Y 2 Z 2 | = 0 (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,

B = m ⋅ b (3)

where m is the photographic scale denominator; ( b x , b y , b z ) are the components of the projection baseline. From the geometric relationship in

{ b y = b x ⋅ tan μ ≈ b x ⋅ μ b z = b x cos μ ⋅ tan ν ≈ b x ⋅ ν (4)

in (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 = | b x b x μ b x ν X 1 Y 1 Z 1 X 2 Y 2 Z 2 | = b x | 1 μ ν X 1 Y 1 Z 1 X 2 Y 2 Z 2 | = 0 (5)

F will be expanded by the Taylor series, holding the first order terms, and the linear formula of the coplanar equation is

F = F 0 + ∂ F ∂ μ d μ + ∂ F ∂ ν d ν + ∂ F ∂ ϕ d ϕ + ∂ F ∂ ω d ω + ∂ F ∂ κ d κ = 0 (6)

the derivation is finally available:

Q = b x d μ − Y 2 Z 2 b x d ν − X 2 Y 2 Z 2 N 2 d ϕ − ( Z 2 + Y 2 2 Z 2 ) N 2 d ω + X 2 N 2 d κ (7)

the upper formula is converted into an error equation:

V Q = b x d μ − Y 2 Z 2 b x d v − X 2 Y 2 Z 2 N 2 d φ − ( Z 2 + Y 2 2 Z 2 ) N 2 d ω + X 2 N 2 d κ − Q (8)

where X 2 , Y 2 , Z 2 is the coordinate of homonymy points in right image in the photogrammetry coordinate system, b x represents the component of baseline b in the Axis X, N 2 represents the right image projection coefficient, Q is the vertical parallax. μ , υ represent the low-frequency detection error values. The correction amount of attitude in pitch and yaw direction is calculated by global matching connection point, and the error compensation is carried out.

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 three-line array constitutes an image plane equivalent to the same frame, with three CCD linear array spaces parallel and horizontally placed. That is, the on-orbit calibration based on the change values of the inner orientation elements is under the assumption that the image plane is horizontal.

However, in the long flight, the vibration and change of temperature of satellite platform will affect the three-line-array CCD camera’s change, it causes the virtual image plane is not horizontal, showing the change that the tilt up and down.

As shown in

It can be seen from

camera, referred to as the pitch angle error) is all in the same expression form. Therefore, when the angle between star sensor and earth observation cameras geometrically corrected, the tilt error of the three-line-array image plane and low-frequency measurement error will be confused with pitch angle error. When the low-frequency error of the exterior orientation angle elements is compensated, the tilt error of the three-line image plane will be confused to low-frequency measurement error.

In fact, the tilt error of the three-line-array image plane, the pitch angle error and the low-frequency measurement error are different in mechanism, and the error properties are different.

Repeated detection or leak detection without reasonable discrimination will result in lower ground positioning accuracy. Research shows that the tilt error of the three-line image plane is up to tens of seconds, and its influence on the positioning accuracy cannot be ignored.

Therefore, when there is the tilt error of three-line-array image plane, in order to avoid repeated detection or leak detection, it is necessary to re-specify the inner orientation elements in the form of

{ f l = F l ⋅ cos ( α l + α v ) x l = f l ⋅ tan ( α l + α v ) f v = F v x v = 0 f r = F r ⋅ cos ( α r − α v ) x r = − f l ⋅ tan ( α l − α v ) (9)

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

{ Q U = b x d μ U − Y 2 U Z 2 U b x d ν U Q D = b x d μ D − Y 2 D Z 2 D b x d ν U (10)

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

{ x U = x D ， y U = − y D d μ U = d μ D ， d ν U = d ν D (11)

bring Equation(11) into the Equation (12)

[ X 2 Y 2 Z 2 ] = m R [ x y − f ] (12)

we can get X 2 U = X 2 D ,Y 2 U = − Y 2 D ， Z 2 U = Z 2 D .

The vertical parallax of the homonymy points is further transformed into

{ Q U = b x d μ U − Y 2 U Z 2 U b x d ν U = b x d μ + Y 2 Z 2 b x d ν Q D = b x d μ D − Y 2 D Z 2 D b x d ν D = b x d μ - Y 2 Z 2 b x d ν (13)

1) The special correction of d μ available by Formula (14) (

Q U + Q D = 2 b x d μ (14)

the correction value of the obtained d μ is

d μ = Q U + Q D 2 b x (15)

2) The special correction of d ν available by Formula (16) (

Q U − Q D = 2 Y 2 Z 2 b d υ (16)

the correction value of the obtained d ν is:

d υ = Q U − Q D 2 Y 2 Z 2 (17)

The process of low frequency error compensation is as follows and the for details can be seen in

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.

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

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 [

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

Test sites | Dynamic calibration value of TH01 | Dynamic calibration value of TH-02 | ||
---|---|---|---|---|

dv | dµ | dv | dµ | |

Xinjiang | 0.45 | −0.59 | −0.69 | 0.12 |

Beijing | −0.77 | 0.13 | 1.4 | −0.05 |

Heilongjiang | 0.91 | −0.09 | −0.12 | 0.38 |

Anhui | 0.15 | 0.54 | 1.13 | 0.27 |

Heilongjiang and Jilin | 1.42 | −0.11 | −1.98 | −0.42 |

RMS | 0.86 | 0.37 | 1.24 | 0.29 |

Test sites | Dynamic calibration value of TH01 | Dynamic calibration value of TH-02 | ||
---|---|---|---|---|

dv | dµ | dv | dµ | |

Xinjiang | 2.36 | −0.58 | 9.77 | 0.02 |

Beijing | 1.22 | 0.12 | 12.21 | −0.05 |

Heilongjiang | 2.92 | −0.13 | 10.42 | 0.38 |

Anhui | 2.18 | 0.48 | 11.65 | 0.19 |

Heilongjiang and Jilin | 3.34 | −0.21 | 8.54 | −0.38 |

RMS | 2.51 | 0.36 | 10.60 | 0.26 |

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

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:

1) Because of the existence of low-frequency measurement error, the positioning accuracy without GCPs of each test site corrected by the test parameters of Beijing calibration field is lower than the accuracy of each test site corrected by the calibration parameters of each test site.

Test sites | TH01-01 star | TH01-02 star | ||||||
---|---|---|---|---|---|---|---|---|

X | Y | XY | H | X | Y | XY | H | |

Xinjiang | 14.06 | 11.31 | 18.04 | 6.85 | 8.35 | 10.03 | 13.05 | 10.28 |

Beijing | 6.58 | 5.28 | 8.44 | 5.88 | 6.42 | 8.12 | 10.35 | 10.56 |

Heilongjiang | 11.83 | 11.24 | 16.32 | 7.48 | 6.77 | 9.22 | 11.44 | 8.75 |

Anhui | 7.42 | 6.28 | 9.72 | 8.45 | 9.12 | 9.98 | 13.52 | 11.12 |

Heilongjiang and Jilin | 13.33 | 12.14 | 18.10 | 8.34 | 4.23 | 9.78 | 10.66 | 10.35 |

RMS | 11.08 | 9.68 | 14.71 | 7.46 | 7.18 | 9.45 | 11.87 | 10.24 |

Test sites | TH01-01 star | TH01-02 star | ||||||
---|---|---|---|---|---|---|---|---|

X | Y | XY | H | X | Y | XY | H | |

Xinjiang | 4.28 | 5.79 | 7.20 | 6.89 | 6.96 | 5.89 | 9.12 | 5.12 |

Beijing | 3.12 | 5.38 | 6.22 | 5.63 | 7.28 | 6.12 | 9.51 | 10.87 |

Heilongjiang | 3.45 | 4.12 | 5.37 | 5.12 | 3.55 | 3.39 | 4.91 | 3.87 |

Anhui | 2.98 | 3.09 | 4.29 | 7.41 | 3.38 | 5.12 | 6.14 | 4.28 |

Heilongjiang and Jilin | 5.97 | 4.88 | 7.71 | 6.12 | 4.50 | 4.67 | 6.49 | 2.89 |

RMS | 4.11 | 4.75 | 6.28 | 6.29 | 5.40 | 5.13 | 7.45 | 6.10 |

Test sites | TH01-01 star | TH01-02 star | ||||||
---|---|---|---|---|---|---|---|---|

X | Y | XY | H | X | Y | XY | H | |

Xinjiang | 6.12 | 6.09 | 8.63 | 5.12 | 7.12 | 6.12 | 9.39 | 5.12 |

Beijing | 3.22 | 5.80 | 6.63 | 6.02 | 10.02 | 5.48 | 11.42 | 10.33 |

Heilongjiang | 4.89 | 7.31 | 8.79 | 4.83 | 4.12 | 6.23 | 7.47 | 4.08 |

Anhui | 3.12 | 6.97 | 7.64 | 5.21 | 3.09 | 6.04 | 6.78 | 3.99 |

Heilongjiang and Jilin | 5.89 | 4.30 | 7.29 | 6.89 | 6.89 | 5.87 | 9.05 | 2.78 |

RMS | 4.82 | 6.18 | 7.84 | 5.66 | 6.71 | 5.95 | 8.97 | 5.89 |

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.

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.

The authors declare no conflicts of interest regarding the publication of this paper.

Zhou, Y., Hu, X., Du, P. and Wang, Y. (2020) The Satellite’s On-Orbit Attitude System Error Compensation Technique Based on Stereo Models. Journal of Computer and Communications, 8, 42-53. https://doi.org/10.4236/jcc.2020.86005