Horizontal coordinate system

From Wikipedia, the free encyclopedia

 HORIZONTAL COORDINATES. Azimuth, from the North point (red) -also from the South point toward the West (blue). Altitude, green.
HORIZONTAL COORDINATES. Azimuth, from the North point (red) -also from the South point toward the West (blue). Altitude, green.

The horizontal coordinate system is a celestial coordinate system that uses the observer's local horizon as the fundamental plane. This conveniently divides the sky into the upper hemisphere that you can see, and the lower hemisphere that you cannot (because the Earth is in the way). The pole of the upper hemisphere is called the zenith. The pole of the lower hemisphere is called the nadir.

The horizontal coordinates are:

  • altitude (Alt), sometimes referred to as elevation, that is the angle between the object and the observer's local horizon.
  • azimuth (Az), that is the angle of the object around the horizon (measured from the south point towards the west in astronomy and from north point towards the east in geodesy).

The horizontal coordinate system is sometimes also called the az/el[1] or Alt/Az coordinate system.

Contents

[edit] General observations

The horizontal coordinate system is fixed to the Earth, not the stars. Therefore, the altitude and azimuth of an object changes with time, as the object appears to drift across the sky. In addition, because the horizontal system is defined by your local horizon, the same object viewed from different locations on Earth at the same time will have different values of altitude and azimuth.

Horizontal coordinates are very useful for determining the rise and set times of an object in the sky. When an object's altitude is 0°, it is on the horizon, if at that moment its altitude is increasing, it is rising, if its altitude is decreasing it is setting. However all objects on the celestial sphere are subject to the diurnal motion, which is always from east to west, so the inherent cumbersome determination whether altitude is increasing or decreasing can be easily found by considering the azimuth of the celestial object instead:

  • if the azimuth is between 0° and 180° (south—west—north), it is setting.
  • if the azimuth is between 180° and 360° (north—east—south), it is rising.

There are the following special cases:

  • At the northpole all directions are south, so the azimuth is undefined there; no problem, the altitude of all celestial objects remains the same all the time, none ever rises or sets.
  • At the southpole all directions are north; same consideration.
  • At the equator objects on the celestial poles stay at one point on the horizon

Note that the above considerations are strictly speaking true for the geometric horizon only: the horizon as it would appear for an observer on sealevel on a perfect smooth Earth without atmosphere. In practice the apparent horizon, which you see, has a negative altitude, which absolute value gets larger when you come higher, due to the curvature of the Earth. In addition the atmospheric refraction adds another 0.5° to that value.

[edit] Transformation of coordinates

It is possible to pass from the equatorial coordinate system to the horizontal coordinate system and back, once the observer's geographic latitude φ is known (+90° on the northpole, 0° on the equator, -90° on the southpole).

We will use A for the azimuth, a for the altitude.

We will use δ for the declination, H for the hour angle.

[edit] equatorial to horizontal

\sin a = \sin \phi \cdot \sin \delta + \cos \phi \cdot \cos \delta \cdot \cos H

\cos A \cdot \cos a = \cos \phi \cdot \sin \delta - \sin \phi \cdot \cos \delta \cdot \cos H

\sin A \cdot \cos a = - \cos \delta \cdot \sin H

One may be tempted to 'simplify' the last two equations by dividing out the cosa leaving one expression in tanA only. But the tangent cannot distinguish between, let say an azimuth of 45° and 225°, while these two values are very different, they are opposite directions, NE and SW respectively. Only do this when the situation is such that one knows in advance in which quadrant the wanted azimuth has to be.

If the calculation is done with an electronic pocket calculator, it is best not to use the functions arcsin and arccos when possible, because of their limited 180° only range, and also because of the low accuracy the former gets around ±90° and the latter around 0° and 180°. Most scientific calculators have a rectangular to polar (R->P) and polar to rectangular (P->R) function, which avoids that problem and gives us an extra sanity check as well.

The algorithm then becomes as follows.

  • Calculate the terms right of the = sign of the 3 equations given above
  • Apply a R->P conversion taking the \cos A \cdot \cos a as the X value and the \sin A \cdot \cos a as the Y value
  • The angle part of the answer is the azimuth, an angle over the full range of 0° to 360° (or -180° to +180° etc)
  • Apply a second R->P conversion taking the radius part of the last answer as the X and the sina of the first equation as the Y value
  • The angle part of the answer is the altitude, an angle between -90° and +90°
  • The radius part of the answer must be 1 exactly, if not you have made an error.

[edit] horizontal to equatorial

\sin \delta = \sin \phi \cdot \sin a + \cos \phi \cdot \cos a \cdot \cos A

\cos \delta \cdot \cos H = \cos \phi \cdot \sin a - \sin \phi \cdot \cos a \cdot \cos A

\cos \delta \cdot \sin H = - \sin A \cdot \cos a

Same considerations as for the first set also holds for this set of formulas.

[edit] The position of the Sun

There are several ways to compute the apparent position of the Sun in horizontal coordinates.

Complete and accurate algorithms to obtain precise values can be found in Jean Meeus's book Astronomical Algorithms.

Instead a simple approximate algorithm is the following:

Given:

You have to compute:

  • The Sun declination of the corresponding day of the year, which is given by the following formula:

\delta = -23.45^\circ \cdot \cos \left ( \frac{360^\circ}{365} \cdot \left ( N + 10 \right ) \right )

where N is the number of days spent since January 1.

  • The true hour angle that is the angle which the earth should rotate to take the observer's location directly under the sun.
    • Let hh:mm be the time the observer reads on the clock.
    • Merge the hours and the minutes in one variable T = hh + mm/60 measured in hours.
    • hh:mm is the official time of the time zone, but it is different from the true local time of the observer's location. T has to be corrected adding the quantity + (Longitude/15 - Time Zone), which is measured in hours and represents the difference of time between the true local time of the observer's location and the official time of the time zone.
    • If it is summer and Daylight Saving Time is used, you have to subtract one hour in order to get Standard Time.
    • The value of the Equation of Time in that day has to be added. Since T is measured in hours, the Equation of Time must be divided by 60 before being added.
    • The hour angle can be now computed. In fact the angle which the earth should rotate to take the observer's location directly under the sun is given by the following expression: H = (12 - T) * 15. Since T is measured in hours and the speed of rotation of the earth 15 degrees per hour, H is measured in degrees. If you need H measured in radians you just have to multiply by the factor 2π/360.

[edit] References

This article's initial version originated from 'Jason Harris' Astroinfo which comes along with KStars, a Desktop Planetarium for Linux/KDE. See http://edu.kde.org/kstars/index.phtml