Calculation of Right Ascension and Declination and its conversion to Azimuth and Altitude

Main Page – Real time planet tracking system

The Idea!- Kepler’s Algorithms

Yes! the Real time position of any celestial body can be calculated using some parameters called Orbital Elements (or Osculating Elements or Keplerian Elements). These are the parameters that define an orbit at a particular time.

1. Inclination (i)angle between the plane of the Ecliptic and th
   

   In this diagram the orbital plane (yellow) intersects a reference plane.For Earth-orbiting satellites, the reference plane is usually the earths equatorial plane and for satellites in the solar orbits it is the elliptic plane.The intersection is called the line of nodes as it connects the center of mass to ascending and descending nodes.This plane together with the vernal point establishes a reference frame. source:1

   e plane of the orbit.

2. Longitude of the Ascending Node (o)states the position in the orbit where the elliptical path of the planet passes through the plane of the ecliptic, from below the plane to above the plane.
3. Longitude of Perihelion (p)states the position in the orbit where the planet is
   closest to the Sun.
4. Mean distance (a)the value of the semi-major axis of the orbit – measured in Astronomical Units for the major planets.
5. Daily motion (n)states how far in degrees the planet moves in one (mean solar) day. This figure can be used to find the mean anomaly of the planet for a given number of days either side of the date of the elements. The figures quoted in the Astronomical Almanac do not tally with the period of the planet as calculated by applying Kepler’s 3rd Law to the semi-major axis.
6. Eccentricity (e)eccentricity of the ellipse which describes the orbit.
7. Mean Longitude (L)Position of the planet in the orbit on the date of the elements.source:2

All the above parameters are change constantly (but very very slowly) due to gravitational perturbations by other objects and the effects of relativity because of that we need to take the latest orbital elements into our calculation for least possible error or highest accuracy.I just googled and got orbital elements for the year 2013 (text) .You can try different methods or go to the following sites to gather data.

• http://asa.usno.navy.mil/SecE/Osculating_Elements.html
• http://ssd.jpl.nasa.gov/ 
•  http://www.imcce.fr/langues/en/
• http://naif.jpl.nasa.gov/naif/
• Research gate forum

Now that you have the latest osculating (orbital) elements for a date lets say 16 August 2013, we will start our calculations but before that we need to know few terms

1. JDN – An Integer(Julian Day No.) that would give the day number count from the beginning of Julian year which will help us locating the planets/celestial bodies in their orbits for a particular date and time relative to a particular date.
2. Equatorial Coordinate System – The coordinates Right Ascension (RA) and Declination (DEC) will be used frequently here.Look at the GIF  to know what is RA and DEC.

The positions of objects in the sky as viewed from Earth are referred to a coordinate system whose alignment is changing with time in a complex way. A few of the important motions and effects are summarised below;

• The Earth is rotating on its axis once every siderial day
• The rotation axis is moving in a circle with a period of roughly 26,000 years (precession)
• The axis is ‘nodding’ up and down with a period of roughly 19 years (nutation)
• The finite speed of light (sometimes referred to as ‘aberration’ in some books)

The ‘fixed’ stars provide a reference system which allows us to account for the daily rotation of the Earth on its axis. We use the Equatorial coordinate system to refer positions to a frame in which the stars appear still, and the right ascension (RA) and declination (DEC) are used to give the coordinates of the planet with respect to the fixed stars. The ‘zero’ of RA is refered to the ‘vernal equinox’, in the same way that the ‘zero’ of longitude is taken as the Greenwich Meridian.

The presession of the equinoxes means that the ‘zero’ of RA is changing slowly with time, which means that star coordinates must always be referred to an epoch, or date. By using orbital elements referred to the fundamental epoch J2000, the orbits of the planets are described in a coordinate system which is based on the position the vernal equinox will have at J2000. A further advantage of this dodge is that our positions for the planets will correspond exactly to the positions found in most recent star charts. You should be able to plot the path of Mars directly onto a star chart such asThe Cambridge Star Atlas.

Nutation (which is a small effect anyway) can also be spirited away by referring our positions to the ‘mean ecliptic of J2000’. The word ‘mean’ indicates that no allowance for nutation has been made. Our observation platform (the Earth) is nodding, so the stars and planets will appear to nod together. Our J2000 elements will give is positions which match the co-ordinates of the stars found in star maps.

There is a problem with this use of J2000 equinox and mean ecliptic. If I just dial the values for RA and DEC into a computerised telescope, then the planet will not appear in the centre of the field of view – as the RA and DEC will not be referred to the ‘equinox and true ecliptic of date’. The effect will be very small for 10 years either side of J2000. (source:2)

So we would increase the accuracy by using elements from 2013.

Orbital Elements for 16th August 2013

Calculation

The sections below deal with calculating the RA and DEC of a planet from the osculating elements. Lets find the position of Mars at UT on the 1st of March 2016. The main steps in the calculation are;

• Finding the position of the planet in its orbit
  • Find the number of days since the date of the elements
  • Find the mean anomaly from the Mean Longitude and the daily motion
  • Find the true anomaly using the Equation of Centre
  • Find the radius vector of the planet
• Refer that position to the Ecliptic – hence find the heliocentric ecliptic coordinates of the planet
• Repeat most of above to find the heliocentric coordinates of the Earth
• Transform the heliocentric coordinates to geocentric coordinates by a change of origin
• Transform the geocentric ecliptic coordinates to geocentric equatorial coordinates by a rotation about the X axis
• Calculate the RA and DEC and Earth – planet distance from the rectangular coordinates

The method used here was adapted from Paul Schlyter’s page ‘How to compute planetary positions’ at source:3

Notation

Position of the planet in its orbit

Number of days from date of elements (d)

What you need:

• ‘day number’ (dele) of the elements
• ‘day number’ you want the position for (dpos)

The ‘day number’ can be the

1. Julian day, or
2. The number of days since the fundamental epoch J2000. I use the second alternative as less precision is needed for the numbers!

The following tables show the days from the beginning of the year to the beginning of each month, and the days from J2000 to the beginning of each year.

we can find the day number corresponding to the date of the elements (16th August 2013) as follows;

and the day number of the date we want the position for (1st Aprih 2016) is;

so the number of days after the date of the elements is

i.e. 928 days after the elements. (You must take dates before an epoch as negative in the calculations below.)

For fast moving planets such as Mercury and Mars, you need to include the time of day which you want the position for. Just add the Universal Time in decimal hours divided by 24 to the day number of your position (dele above);

Finding the Mean Anomaly of the planet

The Mean Anomaly of the planet is given by the very simple formula;

For our case of Mars and 928 days since the date of the elements;

Finding the true anomaly of the planet

Mean Anomaly is calculated considering the orbit to be circular, but True Anomaly gives the actual position of the planet as it considers orbit to be elliptical. We convert Mean Anomaly to True Anomaly using the following formula:

This is manually  calculated using the Equation of Centre from “The Astronomical Almanac (page E4)”.

For our Mars position, we have

Finding the radius vector of the planet

The distance from the planet to the focus of the ellipse is given by a simple formula based on the geometry of the ellipse;

In our Mars calculation we have;

Heliocentric coordinates of the planet

Having found the true anomaly and the radius vector of the planet, we can go on to find the position of the planet with respect to the plane of the ecliptic. The formulas below are a combination of ‘resolving’ to find components and rotations around various axes to transform the coordinates to the Ecliptic frame. We might expect the formulas to involve the inclination of the planet’s orbit (i), and various angles within the plane of the orbit, as well as the longitude of the ascending node (o).

In the case of Mars we have;

Heliocentric coordinates of Earth

Similarly find M(mean anomaly),V(true anomaly) and r(radius vector) for earth

We simply the equations for Earth, as the inclination of the Earth’s orbit is very small.so…

For the problem in hand we have;

Geocentric ecliptic coordinates of the planet

Just subtract  Earth’s coordinates from those of the planet and we get geocentric (earth as a centre ) coordinates.

We then have the geocentric ecliptic coordinates of the planet. For the case of Mars (1st March 2016) we have,

Geocentric equatorial coordinates of the planet

To change the coordinate system from geocentric ecliptic to geocentric equatorial is just a matter of a rotation around the X axis by an angle equal to the ‘obliquity of the Ecliptic. The X axis points towards the ‘First point of Aries’, which is the direction in space associated with the equinox. As we are using elements referred to the equinox of J2000.0, we use the obliquity for that epoch, which is 23.439292 degrees. the formulas are given below;

For Mars, we have

rectangular coordinates are not much use with star charts, so we calculate the familiar right ascension and declination using the formulas;

For Mars on the 1stMarch 2016 we have

With this one of the hardest part of our project is done! 😀

Azimuth and Elevation (Altitude) conversion from RA and DEC [equitorial coordinates to horizontal coordinates]

Now that we have RA and DEC it gives us the value from the vernal equinox at 1st March 2016,therefore we neede to give our position in order to see with respect to us (the observer as center ).So we first take the latitude(lat) and longitude(long) using a GPS module  of our current location.

For those who are wondering what is Azimuth and Elevation have a look of this image

Note : We will be connecting a GPS module to our arduino setup to constantly refresh and update our coordinates …Now lets proceed with the calculation.

We need to know the

• Local Sidereal Time – First we need the time (hour and minute in UT) and the day no of the date 1st march 2016 [i.e. dpos= 5903.5 ] .Now we use this formula to get LST.  LST = (100.46 + 0.985647 * day + Long + 15 * (hour + minute / 60) + 360) – (((int)                ((100.46 + 0.985647 * day + Long + 15 * (hour + minute / 60) + 360)/360))*360);

• Hour Angle       HA = (LST – RA + 360)- ((int)((LST – RA + 360)/360))*360 ;

• HA, DEC, Lat to Alt, AZ
  x = cos(HA * (pi / 180)) * cos(Dec * (pi / 180));
  y = sin(HA * (pi / 180)) *cos(Dec * (pi / 180));
  z = sin(Dec * (pi / 180));

• Horizontal coordinates :                                                                                                                 xhor = x * cos((90 – Lat) * (pi / 180)) – z *sin((90 – Lat) * (pi / 180));
  yhor = y;
  zhor = x * sin((90 – Lat) * (pi / 180)) + z *cos((90 – Lat) * (pi / 180));

az = atan2(yhor, xhor) * (180 / pi) + 180;
alt = asin(zhor) * (180 / pi);

If you want to know more about these formulae  then do some research.

So finally we have the Result as

Here is a rough Code : Altaz_RaDec.cpp

This is my Repo link : github

We can feed these values to our pan-tilt servo mechanism with proper mapping Link.

References :

1. http://www.stargazing.net/kepler/ellipse.html
2. http://www.stjarnhimlen.se/comp/tutorial.html

RTPT (Real Time Planet Tracking System and Trajectory Prediction) Copyright © 2016 Shubham Paul , Samhita Ganguly ,Rohit Kumar  GNU GPL3+