The Scientific World Journal

Volume 2013 (2013), Article ID 237189, 7 pages

http://dx.doi.org/10.1155/2013/237189

Research Article

An Autonomous Navigation Algorithm for High Orbit Satellite Using Star Sensor and Ultraviolet Earth Sensor

Li Baohua,¹ Lai Wenjie,¹ Chen Yun,¹ and Liu Zongming²

¹Space Control and Inertial Technology Research Center, Harbin Institute of Technology, Harbin 150001, China

²Shanghai Aerospace Control and Engineering Institute, Shanghai 200233, China

Received 31 July 2013; Accepted 10 September 2013

Academic Editors: H. Li and C. Zhang

Copyright © 2013 Li Baohua et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

An autonomous navigation algorithm using the sensor that integrated the star sensor (FOV1) and ultraviolet earth sensor (FOV2) is presented. The star images are sampled by FOV1, and the ultraviolet earth images are sampled by the FOV2. The star identification algorithm and star tracking algorithm are executed at FOV1. Then, the optical axis direction of FOV1 at J2000.0 coordinate system is calculated. The ultraviolet image of earth is sampled by FOV2. The center vector of earth at FOV2 coordinate system is calculated with the coordinates of ultraviolet earth. The autonomous navigation data of satellite are calculated by integrated sensor with the optical axis direction of FOV1 and the center vector of earth from FOV2. The position accuracy of the autonomous navigation for satellite is improved from 1000 meters to 300 meters. And the velocity accuracy of the autonomous navigation for satellite is improved from 100 m/s to 20 m/s. At the same time, the period sine errors of the autonomous navigation for satellite are eliminated. The autonomous navigation for satellite with a sensor that integrated ultraviolet earth sensor and star sensor is well robust.

1. Introduction

An autonomous navigation system of the high orbit satellite can be less reliant on, and eventually does not depend, on the ground under the condition of system support [1]. That’s to say the on-orbit real-time position and velocity of satellite are calculated by themselves with the multiform sensors and navigation algorithm. The merits of the satellite autonomous navigation system not only reduce dependence on the ground level but also improve the system ability to survive. In addition, the satellite autonomous navigation system can also effectively reduce the burden of the ground control stations, ground support to reduce costs, thereby reducing the whole space program development [2].

A star sensor is a camera device which measures the direction of a star in the spacecraft coordinate system. Star sensors are widely used in spacecraft attitude determination because they produce higher accuracy attitude estimates than any other existing sensor [3, 4]. But the star sensor can only output attitude angle of aircraft relative to the inertial coordinate system for aircraft control system. However, the position and velocity information must be provided for the satellite from other sensors.

Infrared earth sensor is an optical measurement instrument, which measures the earth and the sky difference of infrared radiation and accesses satellite attitude information. Infrared earth sensors are usually used for measuring the pitch angle of aircraft attitude angle and the rolling attitude angle [5, 6]. The current satellite autonomous navigation system uses star sensor and infrared earth sensor. The starlight angular distance observed information is obtained through infrared earth sensor installation matrix vector under vehicle body coordinate system to geocentric direction of projection, and it is depending on the radius of the earth. Then, the observed variable is starlight angular distance, and the radius of earth, as the combination of satellite orbit dynamics model and filtering technology, can determine the position of the vehicle information. However, due to the fact that earth's infrared radiation characteristic is not stable, the measurement accuracy of earth radius calculated with infrared earth sensor is reduced. So, the satellite navigation information accuracy is reduced because the radius of the earth precision is reduced [7, 8].

Earth edge energy at ultraviolet radiation, no matter of day or at night, is adapted to the detection of ultraviolet. The variety of earth edge energy at ultraviolet radiation extreme height is small, along with the longitude and latitude and season change. So, the edge energy at ultraviolet radiation is of enough accuracy to provide a signal for the establishment of the image sensor. Therefore, the autonomous navigation for satellite is used with earth ultraviolet band and stars visible light band [9].

2. Navigation Principium with the Star Sensor and Ultraviolet Earth Sensor

2.1. The Characteristic of Ultraviolet Radiation of the Earth

The analysis of variation between atmospheric transmittance and ozone transmittance over 0.2 km height above sea level at ultraviolet wavelengths are shown in Figure 1. The energy radiation at ultraviolet wavelengths (200–300 nm) within the scope of atmospheric transmittance is zero. But the energy radiation at ultraviolet wavelengths (200–300 nm) is certain atmospheric transmittance. However, the energy radiation is determinate over 55 km height above sea level at ultraviolet wavelengths (200–300 nm) in Figure 2. So the energy radiation of ultraviolet can be determinate under 80 km height above sea level. The transmittance of ultraviolet wavelengths is almost 100% above 80 km height above sea level from Figure 3. So it’s feasible that the autonomous navigation of satellite is used with ultraviolet earth sensors, which detect the edge of the earth at ultraviolet wavelengths.

Figure 1: The analysis of variation between atmospheric transmittance and O₃ transmittance.

Figure 2: The refractivity of ultraviolet rays in the upper atmosphere.

Figure 3: The analysis of O₃ Concentration variation ultraviolet rays in the atmosphere.

2.2. Navigation Principium

The autonomous navigation system is constitutive of two Field of Views (FOVs). The first FOV (FOV1) is used to observe the stars which is called star sensor, and another FOV (FOV2) is used to observe ultraviolet radiation of the earth which is called ultraviolet earth sensor. The angle between optical axis of star sensor and optical axis of ultraviolet radiation of the earth is 90∘. The star identification algorithm and star tracking algorithm are executed at FOV1. And the quaternion at J2000.0 coordinate system is calculated at FOV1. Then, the optical axis direction of FOV1 at J2000.0 coordinate system is calculated. The ultraviolet image of earth is sampled by FOV2. The coordinates of ultraviolet earth is obtained from ultraviolet image of earth of FOV2. The center vector of earth at FOV2 coordinate system is calculated with the coordinates of ultraviolet earth. The observed vector is obtained according to the geometric relationships among the earth, the optical axis of FOV1, and satellite. The position vector of satellite and direction vector of satellite can be used alone; also we can use a combination of them. The accuracy of the autonomous navigation for satellite can greatly be improved with the observation combining related to the observation value of ultraviolet earth sensor and observation value of star sensor.

2.3. Obtaining the Coordinates of Ultraviolet Earth Image

Based on the ultraviolet earth properties, ultraviolet energy of the earth's surface is mainly caused by atmospheric reflect of the sun's ultraviolet light. So the ultraviolet image of earth in the FOV2 can never be a complete circular gaussian distribution when the satellites move around the earth. The ultraviolet image may be a circular image, or only images of the arc. So the coordinate of ultraviolet earth image is not directly obtained with centroiding algorithm.

If the shape of ultraviolet earth in the ultraviolet earth image is loop in Figure 5, the arbitrary three points called , , and are obtained from the loop. The line segment is done with and . The vertical line of the called is done. And the vertical line is through the center point of the line segment . In the same way, the line segment is done with and . The vertical line of the called is done. And the vertical line is through the center point of the line segment . The point of intersection between the vertical line and vertical line 23, which is the coordinates of ultraviolet earth image in the FOV2 coordinate system, is obtained. Then, the center vector of earth at FOV2 coordinate system is calculated with the point of intersection .

If the shape of ultraviolet earth in the ultraviolet earth image is arc in Figure 6, the endpoint coordinates of the arc called and are obtained from the arc. And the center coordinates of the arc called are obtained from the arc. The line segment is done with and . The vertical line of the called is done. And the vertical line is through the center point of the line segment . In the same way, the line segment is done with and . The vertical line of the called is done. And the vertical line is through the center point of the line segment . The point of intersection between the vertical line and vertical line 23, which is the coordinates of ultraviolet earth image in the FOV2 coordinate system, is obtained. Then, the center vector of earth at FOV2 coordinate system is calculated with the point of intersection .

3. Autonomous Algorithm

3.1. State Model of the Autonomous Algorithm

Consider the following.

Autonomous navigation algorithm state model of the satellite is the satellite orbit dynamics equation. There is a variety of forms of the satellite orbit dynamics equation. The most commonly the satellite orbit dynamics equation is used with celestial navigation system of rectangular coordinate expression of perturbation motion equation and Newton equations of motion among them.

The satellite orbit dynamic equation of the rectangular coordinate system is chosen at epoch (J2000.0) geocentric equatorial coordinate system. The autonomous navigation algorithm state model of the satellite is chosen as formula (1) [10].

That’s to say

where, is the state vector,   is the position scalar quantity of axis at J2000.0 coordinate system, is the position scalar quantity of axis at J2000.0 coordinate system, is the position scalar quantity of axis at J2000.0 coordinate system, is the velocity scalar quantity of axis at J2000.0 coordinate system, is the velocity scalar quantity of axis at J2000.0 coordinate system, is the velocity scalar quantity of axis at J2000.0 coordinate system, is the earth’ gravitational constant, is the position vector at J2000.0 coordinate system, is the coefficient of earth’s gravity, are the disturbance factor from the moon, the pressure of the sun, and atmosphere of the earth, and is the radius of the earth.

3.2. Observed Equation

The star light angle between star vector observed from the satellite and the position vector at J2000.0 is used for observed variable. According the geometric relationships from Figure 4, the star light angle is calculated as follows:

So the observed equation is as follows:

where the obtained from at FOV2 is the satellite position vector at J2000.0 coordinate system, the vector obtained from at FOV1 is the optical axis direction of FOV1 at J2000.0 coordinate system, and is the observed noise.

Figure 4: The autonomous theory with star light and ultraviolet earth.

Figure 5: The ring image of ultraviolet earth.

Figure 6: The arc image of ultraviolet earth.

3.3. Filter Equation

Formula (1) and formula (4) must be discrete and linearized because they are nonlinear equation, and the Extended Kalman Filter (EKF) is used for the autonomous navigation algorithm of the satellite. The discrete equation of formula (1) is as follows:

The second order Taylor series of the formula (5) at called linearization is as follows:

where is the sampling time.

The discrete equation and the linearization of formula (4) at are as follows:

That’s to say

The observed variable is the angle between the optical axis direction of FOV1 and optical axis of ultraviolet radiation of the earth FOV2. The observed matrix is as follows:

The EKF filter equation based error is as follows:

where , , and .

So the most merit filter equation is as follows:

4. Experiment Result and Analysis

In order to validate the proposed autonomous navigation algorithm of high orbit satellite using star sensor and ultraviolet earth sensor, the navigation accuracy is tested in comparison with the star sensor and Infrared earth sensor. All the experiments are done with star sensor whose FOV is 14∘ × 14∘, whose size of CCD is 1024 × 1024 pixels of 13.2 × 13.2 μm pixel size, and whose update rate up is 10 Hz and the ultraviolet earth sensor whose FOV is 20∘ × 20∘, whose size is 2048 × 2048 pixels of 5.5 × 5.5 μm pixel size, and whose update rate up is 8 Hz (Figure 7). The autonomous navigation algorithm of high orbit satellite using star sensor and ultraviolet earth sensor is used for the earth synchronous orbit satellite. The orbit period of the earth synchronous orbit satellite is 24 hours. The feasibility, robustness, and navigation accuracy of the autonomous navigation algorithm of high orbit satellite using star sensor and ultraviolet earth sensor are validated under the sun in different positions. In order to reduce the amount of data sent to the PC, the sampling period of the data is 15 seconds. However, the control period of the autonomous navigation algorithm of high orbit satellite using star sensor and ultraviolet earth sensor is 200 milliseconds.

Figure 7: The experiment method of integrated sensor.

A sensor that integrated ultraviolet earth sensor and star sensor is laid on the test bench. The star simulator is laid in front of the lens of star sensor. The star simulator is adjusted, which results in the angle between the optical axis vector of star simulator, and the optical axis vector of star sensor is . The ultraviolet simulator is laid in front of the lens of ultraviolet earth sensor. The ultraviolet simulator is adjusted, which results in the angle between the optical axis vector of ultraviolet simulator, and the optical axis vector of ultraviolet earth sensor is . Let the sensor that integrated ultraviolet earth sensor and star sensor put on. The RS422 is connected between the integrated ultraviolet earth sensor and star sensor and the PC. Let the star simulator and ultraviolet simulator put on. The real-time satellite orbit parameters are sent to star simulator and ultraviolet simulator with Ethernet network from orbit dynamics simulator. The star image is sampled by the star sensor of integrated sensor, and the ultraviolet image is sampled by ultraviolet earth sensor of integrated sensor. The autonomous navigation data of satellite are calculated by integrated sensor with the star image and ultraviolet image. Then, the autonomous navigation data of satellite are sent to the PC with RS422 from the integrated sensor. At the same time, the ideal navigation data of satellite to the PC are from orbit dynamics simulator. Then the errors of navigation data of satellite are calculated and saved between autonomous navigation data and ideal navigation data of satellite in the PC. The errors curves of position and velocity of satellite are showed with offline errors of navigation data in Figure 8.

(a) The errors of satellite position

(b) The errors of satellite velocity

Figure 8: The autonomous result with a sensor that integrated ultraviolet earth sensor and star sensor.

In order to further compare the errors accuracy of navigation using a sensor that integrated ultraviolet earth sensor and star sensor with the errors accuracy of navigation using a sensor that integrated infrared earth sensor and star sensor. The errors of navigation data of satellite are tested with a sensor that integrated infrared earth sensor and star sensor. The errors curves of position and velocity of satellite are showed with offline errors of navigation data in Figure 9.

(a) The errors of satellite position

(b) The errors of satellite velocity

Figure 9: The autonomous result with a sensor that integrated infrared earth sensor and star sensor.

The autonomous navigation for satellite is calculated with a sensor that integrated infrared earth sensor and star sensor. Firstly, the infrared earth image is sampled with infrared earth sensor. Most infrared earth sensor sample the infrared earth image with CO₂ of wavelengths 14–16 μm. But the energy radiation of the infrared wavelengths is variational with day or night, longitude and latitude, and season. So the accuracy of earth radius vector is reduced with infrared earth sensor. As a result, the accuracy of autonomous navigation for satellite is reduced. The position errors of autonomous navigation for satellite are about 1000 meters, and the velocity errors of autonomous navigation for satellite are about 100 meters/second with a sensor that integrated infrared earth sensor and star sensor in Figure 9. However, earth edge energy at ultraviolet radiation, no matter at day or at night, is adapted to the detection of ultraviolet. The variety of earth edge energy at ultraviolet radiation extreme height is small, along with the longitude and latitude and season change. So the position errors of autonomous navigation for satellite are less than 300 meters, and the velocity errors of autonomous navigation for satellite are less than 20 meters/second with a sensor that integrated ultraviolet earth sensor and star sensor in Figure 8.

5. Summary

Earth edge energy at ultraviolet radiation, no matter at day or at night, is adapted to the detection of ultraviolet. The variety of earth edge energy at ultraviolet radiation extreme height is small, along with the longitude and latitude and season change. So the edge energy at ultraviolet radiation is of enough accuracy to provide a signal for the establishment of the image sensor. Therefore, an autonomous navigation algorithm for high orbit satellite using star sensor and ultraviolet earth sensor is presented. A sensor that integrated ultraviolet earth sensor (FOV2) and star sensor (FOV1) is developed. The star images are sampled by FOV1 of the integrated sensor, and the ultraviolet earth images are sampled by the FOV2 of the integrated sensor. The autonomous navigation data of satellite are calculated by integrated sensor with the star image and ultraviolet image. The star identification algorithm and star tracking algorithm are executed at FOV1, and the quaternion at J2000.0 coordinate system is calculated at FOV1. Then, the optical axis direction of FOV1 at J2000.0 coordinate system is calculated. The ultraviolet image of earth is sampled by FOV2. The coordinates of ultraviolet earth is obtained from ultraviolet image of earth of FOV2. The center vector of earth at FOV2 coordinate system is calculated with the coordinates of ultraviolet earth. The observed vector is obtained according to the geometric relationships among the earth, the optical axis of FOV1, and satellite. The position vector of satellite and direction vector of satellite can be used only; also we can use a combination them. The accuracy of the autonomous navigation for satellite is doubled. And the period errors of the autonomous navigation for satellite are eliminated. So the autonomous navigation for satellite with a sensor that integrated ultraviolet earth sensor and star sensor is well robust.

References

M. A. Chory, D. P. Hoffman, and J. L. LeMay, “Satellite autonomous navigation status and history,” IEEE Postion Location and Navigation SymposiumNV, vol. 47, pp. 110–121, 1986.

G. Shorshi and I. Y. Bar-Itzhack, “Satellite autonomous navigation based on magnetic field measurements,” Journal of Guidance, Control, and Dynamics, vol. 18, no. 4, pp. 843–850, 1995.

C. C. Liebe, E. W. Dennison, B. Hancock, R. C. Stirbl, and B. Pain, “Active pixel sensor (APS) based star tracker,” in Proceedings of the IEEE Aerospace Conference, pp. 119–127, March 1998.

C. C. Liebe, “Star trackers for attitude determination,” IEEE Aerospace and Electronic Systems Magazine, vol. 10, no. 6, pp. 10–16, 1995.

P. C. Calhoun, R. Adhami, and L. L. Joiner, “Application of horizon sensor measurements in Kalman filtering techniques for on-board orbit navigation during aerobraking,” Circuits, Systems, and Signal Processing, vol. 21, no. 4, pp. 399–413, 2002.

D. Pledger, J. B. Ross, and W. Saylor, “Development of honeywell's earth reference attitude determination system,” in Proceeding of the 7th Annual AIAA/USU Conference on Small Satellites (ERADS '93), 1993.

R. Gounley, R. White, and E. Gai, “Autonomous satellite navigation by stellar refraction,” Journal of Guidance, Control, and Dynamics, vol. 7, no. 2, pp. 129–134, 1984.

D. Pledger and J. B. Ross, “Development of Honeywell' Earth Reference Attitude Determination System,” in Proceedings of the of 7th Annual AIAA/USU Conference on Small Satellites, 1993.

S. J. Janz, “Ultraviolet atmospheric and space remote sensing: methods and Instrumentation,” Proceeding Rayleigh Scattering Attitude Sensor, vol. 2831, pp. 146–153, 1996.

R. L. White, S. W. Thurman, and F. A. Barnes, “Autonomous satellite navigation using observation of starlight atmospheric refraction,” Navigation, Journal of the Institute of Navigation, vol. 32, no. 4, pp. 317–333, 1985.