Distance Elevation and Azimuth Calculation » History » Version 9

LANVIN, Jean-baptiste, 12/15/2015 11:16 AM

1 3 LANVIN, Jean-baptiste
2 3 LANVIN, Jean-baptiste
{{toc}}
3 3 LANVIN, Jean-baptiste
4 1 LANVIN, Jean-baptiste
h1. Distance Elevation and Azimuth Calculation
5 1 LANVIN, Jean-baptiste
6 1 LANVIN, Jean-baptiste
h3.    ECI coordinates
7 1 LANVIN, Jean-baptiste
8 1 LANVIN, Jean-baptiste
The coordinates of the satellite given by the propagator are given in the Earth-Centered Inertial (ECI) coordinate system . This system is a cartesian coordinate system whose origin is located at the center of the earth (at the center of mass to be precise). The z axis is orthogonal to the equatorial plane pointing north, the x axis is pointing towards the vernal equinox, and the y axis is such that the system remains a direct orthogonal system The x and y axis are located in the equatorial plane as shown in the following figure.
9 1 LANVIN, Jean-baptiste
10 1 LANVIN, Jean-baptiste
11 5 LANVIN, Jean-baptiste
p=. !ECI.jpg!
12 8 LANVIN, Jean-baptiste
Figure 1 : ECI coordinates
13 1 LANVIN, Jean-baptiste
14 1 LANVIN, Jean-baptiste
This system is convenient to represent the positions and velocities of space objects rotating around the earth considering firstly that the origin of the system is the center of mass of the earth, and that the system does not rotate with the earth. Indeed, "inertial" means that the coordinate system is not accelerating, therefore not rotating: considering the way the three axis are defined, the system is fixed in space regarding the stars.
15 1 LANVIN, Jean-baptiste
The problem here is that our ground is station is located on the surface of the earth, so its coordinates are given in the geodetic coordinate system, which is obviously a system that is rotating with the earth.
16 1 LANVIN, Jean-baptiste
We had therefore two options: either calculate the coordinates of the ground station in the ECI coordinate system, or calculate the coordinates of the satellite in the geodetic coordinate system. We chose the first option .
17 1 LANVIN, Jean-baptiste
18 1 LANVIN, Jean-baptiste
h3. θ(τ) calculation
19 1 LANVIN, Jean-baptiste
20 1 LANVIN, Jean-baptiste
The main problem to go from the geodetic system to the ECI system is to calculate the angle between the ground station meridian (longitude) and the vernal equinox direction at the time of interest τ. This angle, also called the local sidereal time is a measure of time that depends on the stars, and not on the sun and that is why it is a bit touchy to calculate. For our implementation and tests, we relied on the algorithms and example given in the magazine _Satellite Times_ ,in the column Orbital Coordinate Systems, Part II. θ(τ) is actually defined as the sum between the ground station east longitude and the greenwich sidereal time GST, which is the angle between the greenwich meridian and the vernal equinox that we will note  θg(τ).
21 1 LANVIN, Jean-baptiste
θg(τ)can be calculated using the formula:
22 1 LANVIN, Jean-baptiste
23 1 LANVIN, Jean-baptiste
<pre>
24 1 LANVIN, Jean-baptiste
θg(τ) = θg(0h) + ωe·Δτ (1)
25 1 LANVIN, Jean-baptiste
</pre>
26 1 LANVIN, Jean-baptiste
27 1 LANVIN, Jean-baptiste
28 1 LANVIN, Jean-baptiste
where Δτ is the number of seconds elapsed since 0h at the time of the calculation τ and ωe = 7.29211510 × 10-5 radians/second is the Earth's rotation rate. θg(0h) is the GST calculated at 0h that day and it is given as:
29 1 LANVIN, Jean-baptiste
30 1 LANVIN, Jean-baptiste
<pre>
31 1 LANVIN, Jean-baptiste
θg(0h) = 24110.54841 + 864018.812866 Tu + 0.093104 Tu2 - 6.2 × 10-6 Tu3 (2)
32 1 LANVIN, Jean-baptiste
</pre>
33 1 LANVIN, Jean-baptiste
34 1 LANVIN, Jean-baptiste
where Tu = du/36525 and du is the number of days of Universal Time elapsed since the julian date 2451545.0 (2000 January 1, 12h UT1). Therefore, to calculate θg(τ), the first thing we need is the julian date of the day, which we can deduce from the julian date of the year as follow:
35 1 LANVIN, Jean-baptiste
36 1 LANVIN, Jean-baptiste
<pre>
37 1 LANVIN, Jean-baptiste
JD = Julian_Date_of_Year() + Number of day since the first of January
38 1 LANVIN, Jean-baptiste
</pre>
39 1 LANVIN, Jean-baptiste
40 1 LANVIN, Jean-baptiste
41 1 LANVIN, Jean-baptiste
This calculation is done in the VI julian date.vi where the calculation of the Julian date of the year is done using the Meeus' approach.
42 1 LANVIN, Jean-baptiste
43 1 LANVIN, Jean-baptiste
Once we had the Julian date, we could calculate du and therefore Tu, and that way we calculated θg(0h) using (2). From there, we were able to calculate θg(τ) using (1) and since we know the east longitude of the antenna, we obtained θ(τ) by a simple sum. Those calculations are done in the VI tetag.vi where θ(τ) is calculated in radians.
44 1 LANVIN, Jean-baptiste
45 1 LANVIN, Jean-baptiste
46 1 LANVIN, Jean-baptiste
47 1 LANVIN, Jean-baptiste
48 1 LANVIN, Jean-baptiste
49 1 LANVIN, Jean-baptiste
h3. Geodetic to ECI conversion
50 1 LANVIN, Jean-baptiste
51 1 LANVIN, Jean-baptiste
Once we had the local sidereal time, we were able to make the change of system using the following formulas:
52 1 LANVIN, Jean-baptiste
53 1 LANVIN, Jean-baptiste
<pre>
54 1 LANVIN, Jean-baptiste
x=Rcos θ
55 1 LANVIN, Jean-baptiste
y=Rsin θ
56 1 LANVIN, Jean-baptiste
z=Re sin φ
57 1 LANVIN, Jean-baptiste
</pre>
58 1 LANVIN, Jean-baptiste
59 1 LANVIN, Jean-baptiste
where Re = 6378.135 km, φ is the latitude of the antenna, θ is the local sidereal time and R = Re cos φ
60 1 LANVIN, Jean-baptiste
Those calculations are done in the VI _Antenna coordinates ECI.vi_ and x,y,and z are in kilometers.
61 1 LANVIN, Jean-baptiste
62 1 LANVIN, Jean-baptiste
h3. Distance
63 1 LANVIN, Jean-baptiste
64 1 LANVIN, Jean-baptiste
Now that we had the coordinates of the ground station in the ECI coordinate system, and since the propagator gives us the coordinates of the satellite, the distance was pretty trivial to calculate. Indeed, as we said before, the ECI coordinate system is a cartesian system, and therefore the distance between two points A(x,y,z) and B(x',y',z') is defined by:
65 1 LANVIN, Jean-baptiste
66 1 LANVIN, Jean-baptiste
<pre>
67 1 LANVIN, Jean-baptiste
d= sqrt((x-x')²+(y-y')²+(z-z')²)
68 1 LANVIN, Jean-baptiste
</pre>
69 1 LANVIN, Jean-baptiste
70 1 LANVIN, Jean-baptiste
We used this formula to calculate the distance satellite to ground station in the VI distance.vi, which is very important to know for the link budget (for the free space loss) as we will see after.
71 1 LANVIN, Jean-baptiste
   
72 1 LANVIN, Jean-baptiste
h3. Elevation and azimuth
73 1 LANVIN, Jean-baptiste
74 2 LANVIN, Jean-baptiste
The elevation and the azimuth are very important to know as well in order to analyse the satellite link, because if the elevation is too low (below 10 degrees) 
75 1 LANVIN, Jean-baptiste
it is not even worth calculating the link budget because the communication cannot be settled.
76 2 LANVIN, Jean-baptiste
The first step for this calculation is to calculate the range vector r, which is defined by:
77 1 LANVIN, Jean-baptiste
78 1 LANVIN, Jean-baptiste
<pre>
79 1 LANVIN, Jean-baptiste
 r = [rx, ry, rz] = [xs - xa, ys - ya, zs - za].
80 1 LANVIN, Jean-baptiste
(xa,ya,za) being the antenna ECI coordinates and (xs,ys,zs) the satellite ECI coordinates
81 1 LANVIN, Jean-baptiste
</pre>
82 1 LANVIN, Jean-baptiste
83 1 LANVIN, Jean-baptiste
But r is in the ECI coordinate system which is not adapted for the elevation and the azimuth calculation. Therefore, we had to make another change of coordinates from the ECI to the Topocentric-Horizon Coordinate System that is defined in the figure below. 
84 1 LANVIN, Jean-baptiste
85 1 LANVIN, Jean-baptiste
p=. !Topocentric-Horizon.jpg!
86 8 LANVIN, Jean-baptiste
Figure 2 : Topocentric-Horizon Coordinate System
87 1 LANVIN, Jean-baptiste
88 1 LANVIN, Jean-baptiste
In the figure, θ is the local sidereal time, and φ is the latitude of the antenna. The new coordinates of the range vector r(Rs,Re,Rz) are therefore defined by  
89 1 LANVIN, Jean-baptiste
90 1 LANVIN, Jean-baptiste
<pre>
91 1 LANVIN, Jean-baptiste
92 1 LANVIN, Jean-baptiste
rS = sin φ cos θ rx + sin φ sin θ ry - cos φ rz
93 1 LANVIN, Jean-baptiste
94 1 LANVIN, Jean-baptiste
rE = -sin θ rx + cos θ ry
95 1 LANVIN, Jean-baptiste
96 1 LANVIN, Jean-baptiste
rZ = cos φ cos θ rx + cos φ sin θ ry + sin φ rz
97 1 LANVIN, Jean-baptiste
98 1 LANVIN, Jean-baptiste
</pre>
99 1 LANVIN, Jean-baptiste
100 1 LANVIN, Jean-baptiste
The range to the satellite is defined by:
101 1 LANVIN, Jean-baptiste
102 1 LANVIN, Jean-baptiste
<pre>
103 1 LANVIN, Jean-baptiste
r = √ [rS2 + rE2 + rZ2]
104 1 LANVIN, Jean-baptiste
</pre>
105 1 LANVIN, Jean-baptiste
106 1 LANVIN, Jean-baptiste
The elevation and the azimuth are then given by
107 1 LANVIN, Jean-baptiste
108 1 LANVIN, Jean-baptiste
<pre>
109 1 LANVIN, Jean-baptiste
El = arcsin(rZ / r)
110 1 LANVIN, Jean-baptiste
111 1 LANVIN, Jean-baptiste
Az = arctan(-rE / rS)
112 1 LANVIN, Jean-baptiste
</pre>
113 1 LANVIN, Jean-baptiste
114 1 LANVIN, Jean-baptiste
h3. Tests
115 1 LANVIN, Jean-baptiste
116 1 LANVIN, Jean-baptiste
In order to check the validity of our program we had to test each step of our calculation.
117 1 LANVIN, Jean-baptiste
Considering that the local sidereal time is not a very familiar value, to test the calculation, we used the examples given in Satellite Times to compare our values. We first check the julian date of interest and then the value of θg(τ) and θ(τ).
118 1 LANVIN, Jean-baptiste
Once that was ok, the distance was easier to test because we knew that we were supposed to find something around 850 km.
119 9 LANVIN, Jean-baptiste
120 9 LANVIN, Jean-baptiste
h3. References
121 9 LANVIN, Jean-baptiste
122 9 LANVIN, Jean-baptiste
ECI coordinates, Azimuth, Elevation, Distance & Julian Date Calculation  : htps://celestrak.com/columns/v02n02/
123 9 LANVIN, Jean-baptiste
htps://celestrak.com/columns/v02n01/