Sight reduction

Last updated

In astronavigation, sight reduction is the process of deriving from a sight (in celestial navigation usually obtained using a sextant) the information needed for establishing a line of position, generally by intercept method.

Contents

Sight is defined as the observation of the altitude, and sometimes also the azimuth, of a celestial body for a line of position; or the data obtained by such observation. [1]

The mathematical basis of sight reduction is the circle of equal altitude. The calculation can be done by computer, or by hand via tabular methods and longhand methods.

Algorithm

Steps for measuring and correcting Ho using a sextant. Corrections for Sextant Altitude.en.jpg
Steps for measuring and correcting Ho using a sextant.
Using Ho, Z, Hc in intercept method. MarcqSaintHilaire.en.jpg
Using Ho, Z, Hc in intercept method.

Given:

First calculate the altitude of the celestial body using the equation of circle of equal altitude:

The azimuth or (Zn=0 at North, measured eastward) is then calculated by:

These values are contrasted with the observed altitude . , , and are the three inputs to the intercept method (Marcq St Hilaire method), which uses the difference in observed and calculated altitudes to ascertain one's relative location to the assumed point.

Tabular sight reduction

The methods included are:

Longhand haversine sight reduction

This method is a practical procedure to reduce celestial sights with the needed accuracy, without using electronic tools such as calculator or a computer. And it could serve as a backup in case of malfunction of the positioning system aboard.

Doniol

The first approach of a compact and concise method was published by R. Doniol in 1955 [4] and involved haversines. The altitude is derived from , in which , , .

The calculation is:

n = cos(LatDec) m = cos(Lat + Dec) a = hav(LHA) Hc = arcsin(na ⋅ (m + n))

Ultra compact sight reduction

Haversine Sight Reduction algorithm Haversine Sight Reduction.jpg
Haversine Sight Reduction algorithm

A practical and friendly method using only haversines was developed between 2014 and 2015, [5] and published in NavList.

A compact expression for the altitude was derived [6] using haversines, , for all the terms of the equation:

where is the zenith distance,

is the calculated altitude.

The algorithm if absolute values are used is:

if same name for latitude and declination (both are North or South)   ''n'' = hav(|<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Lat</span>| − |<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Dec</span>|)   ''m'' = hav(|<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Lat</span>| + |<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Dec</span>|) if contrary name (one is North the other is South)   ''n'' = hav(|<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Lat</span>| + |<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Dec</span>|)   ''m'' = hav(|<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Lat</span>| − |<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Dec</span>|) ''q'' = ''n'' + ''m'' ''a'' = hav(''LHA'') hav(''ZD'') = ''n'' + ''a'' · (1 − ''q'') ''ZD'' = archav() -> inverse look-up at the haversine tables ''Hc'' = 90° − ''ZD''

For the azimuth a diagram [7] was developed for a faster solution without calculation, and with an accuracy of 1°.

Azimuth diagram by Hanno Ix Azimuth diagram by Hanno Ix.jpg
Azimuth diagram by Hanno Ix

This diagram could be used also for star identification. [8]

An ambiguity in the value of azimuth may arise since in the diagram . is E↔W as the name of the meridian angle, but the N↕S name is not determined. In most situations azimuth ambiguities are resolved simply by observation.

When there are reasons for doubt or for the purpose of checking the following formula [9] should be used:

The algorithm if absolute values are used is:

if same name for latitude and declination (both are North or South)   ''a'' = hav(90° − |<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Dec</span>|) if contrary name (one is North the other is South)   ''a'' = hav(90° + |<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Dec</span>|) ''m'' = hav(|<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Lat</span>| + ''Hc'') ''n'' = hav(|<span class="nowrap" style="padding-left:0.1em; padding-right:0.1em;">Lat</span>| − ''Hc'') ''q'' = ''n'' + ''m'' hav(''Z'') = (''a'' − ''n'') / (1 − ''q'') ''Z'' = archav() -> inverse look-up at the haversine tables if Latitude ''N'':   if ''LHA'' > 180°, ''Zn'' = ''Z''   if ''LHA'' < 180°, ''Zn'' = 360° − ''Z'' if Latitude ''S'':   if ''LHA'' > 180°, ''Zn'' = 180° − ''Z''   if ''LHA'' < 180°, ''Zn'' = 180° + ''Z''

This computation of the altitude and the azimuth needs a haversine table. For a precision of 1 minute of arc, a four figure table is enough. [10] [11]

An example

Data:   ''Lat'' = 34° 10.0′ N (+)   ''Dec'' = 21° 11.0′ S (−)   ''LHA'' = 57° 17.0′ Altitude ''Hc'':   ''a'' = 0.2298   ''m'' = 0.0128   ''n'' = 0.2157   hav(''ZD'') = 0.3930   ''ZD'' = archav(0.3930) = 77° 39′   ''Hc'' = 90° - 77° 39′ = 12° 21′ Azimuth ''Zn'':   ''a'' = 0.6807   ''m'' = 0.1560   ''n'' = 0.0358   hav(''Z'') = 0.7979   ''Z'' = archav(0.7979) = 126.6°   Because ''LHA'' < 180° and Latitude is ''North'': ''Zn'' = 360° - ''Z'' = 233.4°

See also

Related Research Articles

<span class="mw-page-title-main">Pauli matrices</span> Matrices important in quantum mechanics and the study of spin

In mathematical physics and mathematics, the Pauli matrices are a set of three 2 × 2 complex matrices which are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

<i>n</i>-sphere Generalized sphere of dimension n (mathematics)

In mathematics, an n-sphere or hypersphere is an n-dimensional generalization of the 1-dimensional circle and 2-dimensional sphere to any non-negative integer n. The n-sphere is the setting for n-dimensional spherical geometry.

<span class="mw-page-title-main">Residue (complex analysis)</span> Attribute of a mathematical function

In mathematics, more specifically complex analysis, the residue is a complex number proportional to the contour integral of a meromorphic function along a path enclosing one of its singularities. Residues can be computed quite easily and, once known, allow the determination of general contour integrals via the residue theorem.

In calculus, and more generally in mathematical analysis, integration by parts or partial integration is a process that finds the integral of a product of functions in terms of the integral of the product of their derivative and antiderivative. It is frequently used to transform the antiderivative of a product of functions into an antiderivative for which a solution can be more easily found. The rule can be thought of as an integral version of the product rule of differentiation.

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols , (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δf (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p).

Integration is the basic operation in integral calculus. While differentiation has straightforward rules by which the derivative of a complicated function can be found by differentiating its simpler component functions, integration does not, so tables of known integrals are often useful. This page lists some of the most common antiderivatives.

<span class="mw-page-title-main">Rhumb line</span> Arc crossing all meridians of longitude at the same angle

In navigation, a rhumb line, rhumb, or loxodrome is an arc crossing all meridians of longitude at the same angle, that is, a path with constant bearing as measured relative to true north.

<span class="mw-page-title-main">Great-circle distance</span> Shortest distance between two points on the surface of a sphere

The great-circle distance, orthodromic distance, or spherical distance is the distance along a great circle.

In mathematics, the Jacobi elliptic functions are a set of basic elliptic functions. They are found in the description of the motion of a pendulum, as well as in the design of electronic elliptic filters. While trigonometric functions are defined with reference to a circle, the Jacobi elliptic functions are a generalization which refer to other conic sections, the ellipse in particular. The relation to trigonometric functions is contained in the notation, for example, by the matching notation for . The Jacobi elliptic functions are used more often in practical problems than the Weierstrass elliptic functions as they do not require notions of complex analysis to be defined and/or understood. They were introduced by Carl Gustav Jakob Jacobi (1829). Carl Friedrich Gauss had already studied special Jacobi elliptic functions in 1797, the lemniscate elliptic functions in particular, but his work was published much later.

<span class="mw-page-title-main">Versine</span> 1 minus the cosine of an angle

The versine or versed sine is a trigonometric function found in some of the earliest trigonometric tables. The versine of an angle is 1 minus its cosine.

<span class="mw-page-title-main">Haversine formula</span> Formula for the great-circle distance between two points on a sphere

The haversine formula determines the great-circle distance between two points on a sphere given their longitudes and latitudes. Important in navigation, it is a special case of a more general formula in spherical trigonometry, the law of haversines, that relates the sides and angles of spherical triangles.

In astronomical navigation, the intercept method, also known as Marcq St. Hilaire method, is a method of calculating an observer's position on Earth (geopositioning). It was originally called the azimuth intercept method because the process involves drawing a line which intercepts the azimuth line. This name was shortened to intercept method and the intercept distance was shortened to 'intercept'.

<span class="mw-page-title-main">Longitude by chronometer</span>

Longitude by chronometer is a method, in navigation, of determining longitude using a marine chronometer, which was developed by John Harrison during the first half of the eighteenth century. It is an astronomical method of calculating the longitude at which a position line, drawn from a sight by sextant of any celestial body, crosses the observer's assumed latitude. In order to calculate the position line, the time of the sight must be known so that the celestial position i.e. the Greenwich Hour Angle and Declination, of the observed celestial body is known. All that can be derived from a single sight is a single position line, which can be achieved at any time during daylight when both the sea horizon and the sun are visible. To achieve a fix, more than one celestial body and the sea horizon must be visible. This is usually only possible at dawn and dusk.

Ex-meridian is a celestial navigation method of calculating an observer's position on Earth. The method gives the observer a position line on which the observer is situated. It is usually used when the Sun is obscured at noon, and as a result, a meridian altitude is not possible. The navigator measures the altitude of the Sun as close to noon as possible and then calculates where the position line lies.

<span class="mw-page-title-main">Sine and cosine</span> Fundamental trigonometric functions

In mathematics, sine and cosine are trigonometric functions of an angle. The sine and cosine of an acute angle are defined in the context of a right triangle: for the specified angle, its sine is the ratio of the length of the side that is opposite that angle to the length of the longest side of the triangle, and the cosine is the ratio of the length of the adjacent leg to that of the hypotenuse. For an angle , the sine and cosine functions are denoted simply as and .

<span class="mw-page-title-main">Denavit–Hartenberg parameters</span> Convention for attaching reference frames to links of a kinematic chain

In mechanical engineering, the Denavit–Hartenberg parameters are the four parameters associated with a particular convention for attaching reference frames to the links of a spatial kinematic chain, or robot manipulator.

<span class="mw-page-title-main">Circle of equal altitude</span> Real line of position in celestial navigation

The circle of equal altitude, also called circle of position (CoP), is defined as the locus of points on Earth on which an observer sees a celestial object such as the sun or a star, at a given time, with the same observed altitude. It was discovered by the American sea-captain Thomas Hubbard Sumner in 1837, published in 1843 and is the basis of an important method in celestial navigation.

The isoazimuth is the locus of the points on the Earth's surface whose initial orthodromic course with respect to a fixed point is constant.

References

  1. The American Practical Navigator (2002)
  2. Pub. 249 Volume 1. Stars; Pub. 249 Volume 2. Latitudes 0° to 39°; Pub. 249 Volume 3. Latitudes 40° to 89°
  3. Pub. 229 Volume 1. Latitudes 0° to 15°; Pub. 229 Volume 2. Latitudes 15° to 30°; Pub. 229 Volume 3. Latitudes 30° to 45°; Pub. 229 Volume 4. Latitudes 45° to 60°; Pub. 229 Volume 5. Latitudes 60° to 75°; Pub. 229 Volume 6. Latitudes 75° to 90°.
  4. Table de point miniature (Hauteur et azimut), by R. Doniol, Navigation IFN Vol. III Nº 10, Avril 1955 Paper
  5. Rudzinski, Greg (July 2015). Ix, Hanno. "Ultra compact sight reduction". Ocean Navigator. Portland, ME, USA: Navigator Publishing LLC (227): 42–43. ISSN   0886-0149 . Retrieved 2015-11-07.
  6. Altitude haversine formula by Hanno Ix http://fer3.com/arc/m2.aspx/Longhand-Sight-Reduction-HannoIx-nov-2014-g29121
  7. Azimuth diagram by Hanno Ix. http://fer3.com/arc/m2.aspx/Gregs-article-havDoniol-Ocean-Navigator-HannoIx-jun-2015-g31689
  8. Hc by Azimuth Diagram http://fer3.com/arc/m2.aspx/Hc-Azimuth-Diagram-finally-HannoIx-aug-2013-g24772
  9. Azimuth haversine formula by Lars Bergman http://fer3.com/arc/m2.aspx/Longhand-Sight-Reduction-Bergman-nov-2014-g29441
  10. "NavList: Re: Longhand Sight Reduction (129172)".
  11. Natural-Haversine 4-place Table; PDF; 51kB