Mercator projection

Last updated
Mercator projection of the world between 85degS and 85degN. Note the size comparison of Greenland and Africa Mercator projection Square.JPG
Mercator projection of the world between 85°S and 85°N. Note the size comparison of Greenland and Africa
The Mercator projection with Tissot's indicatrix of deformation. Mercator with Tissot's Indicatrices of Distortion.svg
The Mercator projection with Tissot's indicatrix of deformation.
Mercator 1569 world map (Nova et Aucta Orbis Terrae Descriptio ad Usum Navigantium Emendate Accommodata
) showing latitudes 66degS to 80degN. Mercator 1569.png
Mercator 1569 world map (Nova et Aucta Orbis Terrae Descriptio ad Usum Navigantium Emendate Accommodata) showing latitudes 66°S to 80°N.

The Mercator projection ( /mərˈktər/ ) is a cylindrical map projection presented by Flemish geographer and cartographer Gerardus Mercator in 1569. It became the standard map projection for navigation because of its unique property of representing any course of constant bearing as a straight segment. Such a course, known as a rhumb or, mathematically, a loxodrome, is preferred by navigators because the ship can sail in a constant compass direction to reach its destination, eliminating difficult and error-prone course corrections. Linear scale is constant on the Mercator in every direction around any point, thus preserving the angles and the shapes of small objects and fulfilling the conditions of a conformal map projection. As a side effect, the Mercator projection inflates the size of objects away from the equator. This inflation is very small near the equator, but accelerates with latitude to become infinite at the poles. So, for example, landmasses such as Greenland and Antarctica appear far larger than they actually are relative to landmasses near the equator, such as Central Africa.

Contents

History

There is some controversy over the origins of the Mercator. German polymath Erhard Etzlaub engraved miniature "compass maps" (about 10×8 cm) of Europe and parts of Africa that spanned latitudes 0°–67° to allow adjustment of his portable pocket-size sundials. The projection found on these maps, dating to 1511, was stated by Snyder [1] in 1987 to be the same projection as Mercator's. However, given the geometry of a sundial, these maps may well have been based on the similar central cylindrical projection, a limiting case of the gnomonic projection, which is the basis for a sundial. Snyder amends his assessment to "a similar projection" in 1994. [2]

Joseph Needham, a historian of China, wrote that the Chinese developed the Mercator projection hundreds of years before Mercator did, using it in star charts during the Song Dynasty. [3] However, this was a simple, and common, case of misidentification. The projection in use was the equirectangular projection.

Portuguese mathematician and cosmographer Pedro Nunes first described the mathematical principle of the loxodrome and its use in marine navigation. In 1537, he proposed constructing a nautical atlas composed of several large-scale sheets in the cylindrical equidistant projection as a way to minimize distortion of directions. If these sheets were brought to the same scale and assembled, they would approximate the Mercator projection.

In 1569, Gerhard Kremer, known by his trade name Gerardus Mercator, announced a new projection by publishing a large planispheric map measuring 202 by 124 cm (80 by 49 in) and printed in eighteen separate sheets. Mercator titled the map Nova et Aucta Orbis Terrae Descriptio ad Usum Navigantium Emendata: "A new and augmented description of Earth corrected for the use of sailors". This title, along with an elaborate explanation for using the projection that appears as a section of text on the map, shows that Mercator understood exactly what he had achieved and that he intended the projection to aid navigation. Mercator never explained the method of construction or how he arrived at it. Various hypotheses have been tendered over the years, but in any case Mercator's friendship with Pedro Nunes and his access to the loxodromic tables Nunes created likely aided his efforts.

English mathematician Edward Wright published the first accurate tables for constructing the projection in 1599 and, in more detail, in 1610, calling his treatise "Certaine Errors in Navigation". The first mathematical formulation was publicized around 1645 by a mathematician named Henry Bond (c. 1600–1678). However, the mathematics involved were developed but never published by mathematician Thomas Harriot starting around 1589. [4]

The development of the Mercator projection represented a major breakthrough in the nautical cartography of the 16th century. However, it was much ahead of its time, since the old navigational and surveying techniques were not compatible with its use in navigation. Two main problems prevented its immediate application: the impossibility of determining the longitude at sea with adequate accuracy and the fact that magnetic directions, instead of geographical directions, were used in navigation. Only in the middle of the 18th century, after the marine chronometer was invented and the spatial distribution of magnetic declination was known, could the Mercator projection be fully adopted by navigators.

Despite those position-finding limitations, the Mercator projection can be found in many world maps in the centuries following Mercator's first publication. However, it did not begin to dominate world maps until the 19th century, when the problem of position determination had been largely solved. Once the Mercator became the usual projection for commercial and educational maps, it came under persistent criticism from cartographers for its unbalanced representation of landmasses and its inability to usefully show the polar regions.

The criticisms leveled against inappropriate use of the Mercator projection resulted in a flurry of new inventions in the late 19th and early 20th century, often directly touted as alternatives to the Mercator. Due to these pressures, publishers gradually reduced their use of the projection over the course of the 20th century. However, the advent of Web mapping gave the projection an abrupt resurgence in the form of the Web Mercator projection.

Today, the Mercator can still be found in marine charts, occasional world maps, and Web mapping service, but commercial atlases have largely abandoned it, and wall maps of the world can be found in many alternative projections. Google Maps, which relied on it since 2005, still uses it for local-area maps but dropped the projection from desktop platforms in 2017 for maps that are zoomed out of local areas. Many other online mapping services still exclusively use the Web Mercator.

Properties

As in all cylindrical projections, parallels and meridians on the Mercator are straight and perpendicular to each other. In accomplishing this, the unavoidable east-west stretching of the map, which increases as distance away from the equator increases, is accompanied in the Mercator projection by a corresponding north-south stretching, so that at every point location the east-west scale is the same as the north-south scale, making it a conformal map projection. Conformal projections preserve angles around all locations.

Because the linear scale of a Mercator map increases with latitude, it distorts the size of geographical objects far from the equator and conveys a distorted perception of the overall geometry of the planet. At latitudes greater than 70° north or south the Mercator projection is practically unusable, because the linear scale becomes infinitely large at the poles. A Mercator map can therefore never fully show the polar areas (as long as the projection is based on a cylinder centered on the Earth's rotation axis; see the transverse Mercator projection for another application).

All lines of constant bearing (rhumbs or loxodromes—those making constant angles with the meridians) are represented by straight segments on a Mercator map. The two properties, conformality and straight rhumb lines, make this projection uniquely suited to marine navigation: courses and bearings are measured using wind roses or protractors, and the corresponding directions are easily transferred from point to point, on the map, with the help of a parallel ruler (for example).

Distortion of sizes

As on all map projections, shapes or sizes are distortions of the true layout of the Earth's surface. The Mercator projection exaggerates areas far from the equator.

For example:

Uses

Practically every marine chart in print is based on the Mercator projection due to its uniquely favorable properties for navigation. It is also commonly used by street map services hosted on the Internet, due to its uniquely favorable properties for local-area maps computed on demand. [5]

On the other hand, because of great land area distortions, it is not well suited for general world maps. Therefore, Mercator himself used the equal-area sinusoidal projection to show relative areas. However, despite such distortions, the Mercator projection was, especially in the late 19th and early 20th centuries, perhaps the most common projection used in world maps, despite being much criticized for this use. [6] [7] [8] [9]

Because of its very common usage, the Mercator projection has been supposed to have influenced people's view of the world, [10] and because it shows countries near the Equator as too small when compared to those of Europe and North America, it has been supposed to cause people to consider those countries as less important. [11] As a result of these criticisms, modern atlases no longer use the Mercator projection for world maps or for areas distant from the equator, preferring other cylindrical projections, or forms of equal-area projection. The Mercator projection is still commonly used for areas near the equator, however, where distortion is minimal. It is also frequently found in maps of time zones.

Arno Peters stirred controversy beginning in 1972 when he proposed what is now usually called the Gall–Peters projection to remedy the problems of the Mercator. The projection he promoted is a specific parameterization of the cylindrical equal-area projection. In response, a 1989 resolution by seven North American geographical groups disparaged using cylindrical projections for general-purpose world maps, which would include both the Mercator and the Gall–Peters. [12]

Web Mercator

Many major online street mapping services (Bing Maps, Google Maps, MapQuest, OpenStreetMap, Yahoo! Maps, and others) use a variant of the Mercator projection for their map images [13] called Web Mercator or Google Web Mercator. Despite its obvious scale variation at small scales, the projection is well-suited as an interactive world map that can be zoomed seamlessly to large-scale (local) maps, where there is relatively little distortion due to the variant projection's near-conformality.

The major online street mapping services' tiling systems display most of the world at the lowest zoom level as a single square image, excluding the polar regions by truncation at latitudes of φmax = ±85.05113°. (See below.) Latitude values outside this range are mapped using a different relationship that does not diverge at φ = ±90°.[ citation needed ]

Mathematics of the Mercator projection

The spherical model

Although the surface of Earth is best modelled by an oblate ellipsoid of revolution, for small scale maps the ellipsoid is approximated by a sphere of radius a. Many different methods exist for calculating a. The simplest include (a) the equatorial radius of the ellipsoid, (b) the arithmetic or geometric mean of the semi-axes of the ellipsoid, and (c) the radius of the sphere having the same volume as the ellipsoid. [14] The range for a amongst the possible choices is about 35 km, but for small scale (large region) applications this variation may be ignored, and mean values of 6,371 km and 40,030 km may be taken for the radius and circumference respectively. These are the values used for numerical examples in later sections. Only high-accuracy cartography on large scale maps requires an ellipsoidal model.

Cylindrical projections

The spherical approximation of Earth with radius a can be modelled by a smaller sphere of radius R, called the globe in this section. The globe determines the scale of the map. The various cylindrical projections specify how the geographic detail is transferred from the globe to a cylinder tangential to it at the equator. The cylinder is then unrolled to give the planar map. [15] [16] The fraction R/a is called the representative fraction (RF) or the principal scale of the projection. For example, a Mercator map printed in a book might have an equatorial width of 13.4 cm corresponding to a globe radius of 2.13 cm and an RF of approximately 1/300M (M is used as an abbreviation for 1,000,000 in writing an RF) whereas Mercator's original 1569 map has a width of 198 cm corresponding to a globe radius of 31.5 cm and an RF of about 1/20M.

Cylindrical Projection basics2.svg

A cylindrical map projection is specified by formulae linking the geographic coordinates of latitude φ and longitude λ to Cartesian coordinates on the map with origin on the equator and x-axis along the equator. By construction, all points on the same meridian lie on the same generator [17] of the cylinder at a constant value of x, but the distance y along the generator (measured from the equator) is an arbitrary [18] function of latitude, y(φ). In general this function does not describe the geometrical projection (as of light rays onto a screen) from the centre of the globe to the cylinder, which is only one of an unlimited number of ways to conceptually project a cylindrical map.

Since the cylinder is tangential to the globe at the equator, the scale factor between globe and cylinder is unity on the equator but nowhere else. In particular since the radius of a parallel, or circle of latitude, is R cos φ, the corresponding parallel on the map must have been stretched by a factor of 1/cos φ = sec φ. This scale factor on the parallel is conventionally denoted by k and the corresponding scale factor on the meridian is denoted by h. [19]

Small element geometry

The relations between y(φ) and properties of the projection, such as the transformation of angles and the variation in scale, follow from the geometry of corresponding small elements on the globe and map. The figure below shows a point P at latitude φ and longitude λ on the globe and a nearby point Q at latitude φ + δφ and longitude λ + δλ. The vertical lines PK and MQ are arcs of meridians of length Rδφ. [20] The horizontal lines PM and KQ are arcs of parallels of length R(cos φ)δλ. [21] The corresponding points on the projection define a rectangle of width δx and height δy.

CylProj infinitesimals2.svg

For small elements, the angle PKQ is approximately a right angle and therefore

The previously mentioned scaling factors from globe to cylinder are given by

parallel scale factor    
meridian scale factor  

Since the meridians are mapped to lines of constant x, we must have x = R(λλ0) and δx = Rδλ, (λ in radians). Therefore, in the limit of infinitesimally small elements

Derivation of the Mercator projection

The choice of the function y(φ) for the Mercator projection is determined by the demand that the projection be conformal, a condition which can be defined in two equivalent ways:

  • Equality of angles. The condition that a sailing course of constant azimuth α on the globe is mapped into a constant grid bearing β on the map. Setting α = β in the above equations gives y(φ) = R sec φ.
  • Isotropy of scale factors. This is the statement that the point scale factor is independent of direction so that small shapes are preserved by the projection. Setting h = k in the above equations again gives y(φ) = R sec φ.
Mercator y plot.svg

Integrating the equation

with y(0) = 0, by using integral tables [22] or elementary methods, [23] gives y(φ). Therefore,

In the first equation λ0 is the longitude of an arbitrary central meridian usually, but not always, that of Greenwich (i.e., zero). The difference (λ  λ0) is in radians.

The function y(φ) is plotted alongside φ for the case R = 1: it tends to infinity at the poles. The linear y-axis values are not usually shown on printed maps; instead some maps show the non-linear scale of latitude values on the right. More often than not the maps show only a graticule of selected meridians and parallels

Inverse transformations

The expression on the right of the second equation defines the Gudermannian function; i.e., φ = gd(y/R): the direct equation may therefore be written as y = R·gd−1(φ). [22]

Alternative expressions

There are many alternative expressions for y(φ), all derived by elementary manipulations. [23]

Corresponding inverses are:

For angles expressed in degrees:

The above formulae are written in terms of the globe radius R. It is often convenient to work directly with the map width W = 2πR. For example, the basic transformation equations become

Truncation and aspect ratio

The ordinate y of the Mercator projection becomes infinite at the poles and the map must be truncated at some latitude less than ninety degrees. This need not be done symmetrically. Mercator's original map is truncated at 80°N and 66°S with the result that European countries were moved towards the centre of the map. The aspect ratio of his map is 198/120 = 1.65. Even more extreme truncations have been used: a Finnish school atlas was truncated at approximately 76°N and 56°S, an aspect ratio of 1.97.

Much web based mapping uses a zoomable version of the Mercator projection with an aspect ratio of unity. In this case the maximum latitude attained must correspond to y = ±W/2, or equivalently y/R = π. Any of the inverse transformation formulae may be used to calculate the corresponding latitudes:

Scale factor

Mercator scale plot.svg

The figure comparing the infinitesimal elements on globe and projection shows that when α=β the triangles PQM and P′Q′M′ are similar so that the scale factor in an arbitrary direction is the same as the parallel and meridian scale factors:

This result holds for an arbitrary direction: the definition of isotropy of the point scale factor. The graph shows the variation of the scale factor with latitude. Some numerical values are listed below.

at latitude 30° the scale factor is  k = sec 30° = 1.15,
at latitude 45° the scale factor is  k = sec 45° = 1.41,
at latitude 60° the scale factor is  k = sec 60° = 2,
at latitude 80° the scale factor is  k = sec 80° = 5.76,
at latitude 85° the scale factor is  k = sec 85° = 11.5

Working from the projected map requires the scale factor in terms of the Mercator ordinate y (unless the map is provided with an explicit latitude scale). Since ruler measurements can furnish the map ordinate y and also the width W of the map then y/R = 2πy/W and the scale factor is determined using one of the alternative forms for the forms of the inverse transformation:

The variation with latitude is sometimes indicated by multiple bar scales as shown below and, for example, on a Finnish school atlas. The interpretation of such bar scales is non-trivial. See the discussion on distance formulae below.

World Scale from DMA Series 1150 map.png

Area scale

The area scale factor is the product of the parallel and meridian scales hk = sec2φ. For Greenland, taking 73° as a median latitude, hk = 11.7. For Australia, taking 25° as a median latitude, hk = 1.2. For Great Britain, taking 55° as a median latitude, hk = 3.04.

Distortion

Tissot's indicatrices on the Mercator projection Tissot mercator.png
Tissot's indicatrices on the Mercator projection

The classic way of showing the distortion inherent in a projection is to use Tissot's indicatrix. Nicolas Tissot noted that the scale factors at a point on a map projection, specified by the numbers h and k, define an ellipse at that point. For cylindrical projections, the axes of the ellipse are aligned to the meridians and parallels. [19] [24] [25] For the Mercator projection, h = k, so the ellipses degenerate into circles with radius proportional to the value of the scale factor for that latitude. These circles are rendered on the projected map with extreme variation in size, indicative of Mercator's scale variations.

Accuracy

One measure of a map's accuracy is a comparison of the length of corresponding line elements on the map and globe. Therefore, by construction, the Mercator projection is perfectly accurate, k = 1, along the equator and nowhere else. At a latitude of ±25° the value of sec φ is about 1.1 and therefore the projection may be deemed accurate to within 10% in a strip of width 50° centred on the equator. Narrower strips are better: sec  = 1.01, so a strip of width 16° (centred on the equator) is accurate to within 1% or 1 part in 100. Similarly sec 2.56° = 1.001, so a strip of width 5.12° (centred on the equator) is accurate to within 0.1% or 1 part in 1,000. Therefore, the Mercator projection is adequate for mapping countries close to the equator.

Secant projection

In a secant (in the sense of cutting) Mercator projection the globe is projected to a cylinder which cuts the sphere at two parallels with latitudes ±φ1. The scale is now true at these latitudes whereas parallels between these latitudes are contracted by the projection and their scale factor must be less than one. The result is that deviation of the scale from unity is reduced over a wider range of latitudes.

Cylindrical Projection secant.svg

An example of such a projection is

The scale on the equator is 0.99; the scale is k = 1 at a latitude of approximately ±8° (the value of φ1); the scale is k = 1.01 at a latitude of approximately ±11.4°. Therefore, the projection has an accuracy of 1%, over a wider strip of 22° compared with the 16° of the normal (tangent) projection. This is a standard technique of extending the region over which a map projection has a given accuracy.

Generalization to the ellipsoid

When the Earth is modelled by a spheroid (ellipsoid of revolution) the Mercator projection must be modified if it is to remain conformal. The transformation equations and scale factor for the non-secant version are [26]

The scale factor is unity on the equator, as it must be since the cylinder is tangential to the ellipsoid at the equator. The ellipsoidal correction of the scale factor increases with latitude but it is never greater than e2, a correction of less than 1%. (The value of e2 is about 0.006 for all reference ellipsoids.) This is much smaller than the scale inaccuracy, except very close to the equator. Only accurate Mercator projections of regions near the equator will necessitate the ellipsoidal corrections.

Formulae for distance

Converting ruler distance on the Mercator map into true (great circle) distance on the sphere is straightforward along the equator but nowhere else. One problem is the variation of scale with latitude, and another is that straight lines on the map (rhumb lines), other than the meridians or the equator, do not correspond to great circles.

The distinction between rhumb (sailing) distance and great circle (true) distance was clearly understood by Mercator. (See Legend 12 on the 1569 map.) He stressed that the rhumb line distance is an acceptable approximation for true great circle distance for courses of short or moderate distance, particularly at lower latitudes. He even quantifies his statement: "When the great circle distances which are to be measured in the vicinity of the equator do not exceed 20 degrees of a great circle, or 15 degrees near Spain and France, or 8 and even 10 degrees in northern parts it is convenient to use rhumb line distances".

For a ruler measurement of a short line, with midpoint at latitude φ, where the scale factor is k = sec φ = 1/cos φ:

True distance = rhumb distance ≅ ruler distance × cos φ / RF.   (short lines)

With radius and great circle circumference equal to 6,371 km and 40,030 km respectively an RF of 1/300M, for which R = 2.12 cm and W = 13.34 cm, implies that a ruler measurement of 3 mm. in any direction from a point on the equator corresponds to approximately 900 km. The corresponding distances for latitudes 20°, 40°, 60° and 80° are 846 km, 689 km, 450 km and 156 km respectively.

Longer distances require various approaches.

On the equator

Scale is unity on the equator (for a non-secant projection). Therefore, interpreting ruler measurements on the equator is simple:

True distance = ruler distance / RF     (equator)

For the above model, with RF = 1/300M, 1 cm corresponds to 3,000 km.

On other parallels

On any other parallel the scale factor is sec φ so that

Parallel distance = ruler distance × cos φ / RF     (parallel).

For the above model 1 cm corresponds to 1,500 km at a latitude of 60°.

This is not the shortest distance between the chosen endpoints on the parallel because a parallel is not a great circle. The difference is small for short distances but increases as λ, the longitudinal separation, increases. For two points, A and B, separated by 10° of longitude on the parallel at 60° the distance along the parallel is approximately 0.5 km greater than the great circle distance. (The distance AB along the parallel is (a cos φ) λ. The length of the chord AB is 2(a cos φ) sin λ/2. This chord subtends an angle at the centre equal to 2arcsin(cos φ sin λ/2) and the great circle distance between A and B is 2a arcsin(cos φ sin λ/2).) In the extreme case where the longitudinal separation is 180°, the distance along the parallel is one half of the circumference of that parallel; i.e., 10,007.5 km. On the other hand, the geodesic between these points is a great circle arc through the pole subtending an angle of 60° at the center: the length of this arc is one sixth of the great circle circumference, about 6,672 km. The difference is 3,338 km so the ruler distance measured from the map is quite misleading even after correcting for the latitude variation of the scale factor.

On a meridian

A meridian of the map is a great circle on the globe but the continuous scale variation means ruler measurement alone cannot yield the true distance between distant points on the meridian. However, if the map is marked with an accurate and finely spaced latitude scale from which the latitude may be read directly—as is the case for the Mercator 1569 world map (sheets 3, 9, 15) and all subsequent nautical charts—the meridian distance between two latitudes φ1 and φ2 is simply

If the latitudes of the end points cannot be determined with confidence then they can be found instead by calculation on the ruler distance. Calling the ruler distances of the end points on the map meridian as measured from the equator y1 and y2, the true distance between these points on the sphere is given by using any one of the inverse Mercator formulæ:

where R may be calculated from the width W of the map by R = W/2π. For example, on a map with R = 1 the values of y = 0, 1, 2, 3 correspond to latitudes of φ = 0°, 50°, 75°, 84° and therefore the successive intervals of 1 cm on the map correspond to latitude intervals on the globe of 50°, 25°, 9° and distances of 5,560 km, 2,780 km, and 1,000 km on the Earth.

On a rhumb

A straight line on the Mercator map at angle α to the meridians is a rhumb line. When α = π/2 or 3π/2 the rhumb corresponds to one of the parallels; only one, the equator, is a great circle. When α = 0 or π it corresponds to a meridian great circle (if continued around the Earth). For all other values it is a spiral from pole to pole on the globe intersecting all meridians at the same angle, and is thus not a great circle. [23] This section discusses only the last of these cases.

If α is neither 0 nor π then the above figure of the infinitesimal elements shows that the length of an infinitesimal rhumb line on the sphere between latitudes φ; and φ + δφ is a sec α δφ. Since α is constant on the rhumb this expression can be integrated to give, for finite rhumb lines on the Earth:

Once again, if Δφ may be read directly from an accurate latitude scale on the map, then the rhumb distance between map points with latitudes φ1 and φ2 is given by the above. If there is no such scale then the ruler distances between the end points and the equator, y1 and y2, give the result via an inverse formula:

These formulæ give rhumb distances on the sphere which may differ greatly from true distances whose determination requires more sophisticated calculations. [27]

See also

Notes

  1. Snyder, John P. (1987). Map Projections – A Working Manual. U.S. Geological Survey Professional Paper 1395. United States Government Printing Office, Washington, D.C. p. 38.
  2. Snyder, John P (1993). Flattening the Earth: Two Thousand Years of Map Projections. University of Chicago Press. p. 48. ISBN   0-226-76747-7.
  3. Needham, Joseph (1971). Science and Civilization in China. 4. Cambridge University Press. p. 359.
  4. Monmonier, Mark (2004). Rhumb Lines and Map Wars . University of Chicago Press. p.  72. ISBN   0-226-53431-6.
  5. https://www.mapthematics.com/forums/viewtopic.php?f=8&t=251&p=344&hilit=evil#p344
  6. Kellaway, G.P. (1946). Map Projections p. 37–38. London: Methuen & Co. LTD. (According to this source, it had been claimed that the Mercator projection was used for "imperialistic motives"
  7. Abelson, C.E. (1954). Common Map Projections s. 4. Sevenoaks: W.H. Smith & Sons.
  8. Chamberlin, Wellman (1947). The Round Earth on Flat Paper s. 99. Washington, D.C.: The National Geographic Society.
  9. Fisher, Irving (1943). "A World Map on a Regular Icosahedron by Gnomonic Projection." Geographical Review 33 (4): 605.
  10. "Mercator Projection vs. Peters Projection, part 2". Matt T. Rosenberg, about.com.
  11. "Mercator Projection vs. Peters Projection, part 1". Matt T. Rosenberg, about.com.
  12. American Cartographer. 1989. 16(3): 222–223.
  13. "Google Groups". groups.google.com. Retrieved 30 April 2018.
  14. Maling, pages 77–79.
  15. Snyder, Working manual pp 37—95.
  16. Snyder, Flattening the Earth.
  17. A generator of a cylinder is a straight line on the surface parallel to the axis of the cylinder.
  18. The function y(φ) is not completely arbitrary: it must be monotonic increasing and antisymmetric (y(−φ) = y(φ), so that y(0)=0): it is normally continuous with a continuous first derivative.
  19. 1 2 Snyder. Working Manual, page 20.
  20. R is the radius of the globe and φ is measured in radians.
  21. λ is measured in radians.
  22. 1 2 NIST. See Sections 4.26#ii and 4.23#viii
  23. 1 2 3 Osborne Chapter 2.
  24. Snyder, Flattening the Earth, pp 147—149
  25. More general example of Tissot's indicatrix: the Winkel tripel projection.
  26. Osborne, Chapters 5, 6
  27. See great-circle distance, the Vincenty's formulae, or Mathworld.

Related Research Articles

Latitude geographic coordinate specifying north–south position

In geography, latitude is a geographic coordinate that specifies the north–south position of a point on the Earth's surface. Latitude is an angle which ranges from 0° at the Equator to 90° at the poles. Lines of constant latitude, or parallels, run east–west as circles parallel to the equator. Latitude is used together with longitude to specify the precise location of features on the surface of the Earth. On its own, the term latitude should be taken to be the geodetic latitude as defined below. Briefly, geodetic latitude at a point is the angle formed by the vector perpendicular to the ellipsoidal surface from that point, and the equatorial plane. Also defined are six auxiliary latitudes which are used in special applications.

Geographic coordinate system Coordinate system

A geographic coordinate system is a coordinate system that enables every location on Earth to be specified by a set of numbers, letters or symbols. The coordinates are often chosen such that one of the numbers represents a vertical position and two or three of the numbers represent a horizontal position; alternatively, a geographic position may be expressed in a combined three-dimensional Cartesian vector. A common choice of coordinates is latitude, longitude and elevation. To specify a location on a plane requires a map projection.

Rhumb line 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 or magnetic north.

Transverse Mercator projection The transverse Mercator projection is the transverse aspect of the standard (or Normal) Mercator projection

The transverse Mercator map projection is an adaptation of the standard Mercator projection. The transverse version is widely used in national and international mapping systems around the world, including the UTM. When paired with a suitable geodetic datum, the transverse Mercator delivers high accuracy in zones less than a few degrees in east-west extent.

Scale (map) Ratio of distance on a map to the corresponding distance on the ground

The scale of a map is the ratio of a distance on the map to the corresponding distance on the ground. This simple concept is complicated by the curvature of the Earth's surface, which forces scale to vary across a map. Because of this variation, the concept of scale becomes meaningful in two distinct ways.

Equirectangular projection map projection that maps meridians and parallels to vertical and horizontal straight lines, respectively, producing a rectangular grid

The equirectangular projection is a simple map projection attributed to Marinus of Tyre, who Ptolemy claims invented the projection about AD 100. The projection maps meridians to vertical straight lines of constant spacing, and circles of latitude to horizontal straight lines of constant spacing. The projection is neither equal area nor conformal. Because of the distortions introduced by this projection, it has little use in navigation or cadastral mapping and finds its main use in thematic mapping. In particular, the plate carrée has become a standard for global raster datasets, such as Celestia and NASA World Wind, because of the particularly simple relationship between the position of an image pixel on the map and its corresponding geographic location on Earth.

Sinusoidal projection pseudocylindrical equal-area map projection

The sinusoidal projection is a pseudocylindrical equal-area map projection, sometimes called the Sanson–Flamsteed or the Mercator equal-area projection. Jean Cossin of Dieppe was one of the first mapmakers to use the sinusoidal, appearing in a world map of 1570.

Universal Transverse Mercator coordinate system coordinate system

The Universal Transverse Mercator (UTM) is a system for assigning coordinates to locations on the surface of the Earth. Like the traditional method of latitude and longitude, it is a horizontal position representation, which means it ignores altitude and treats the earth as a perfect ellipsoid. However, it differs from global latitude/longitude in that it divides earth into 60 zones and projects each to the plane as a basis for its coordinates. Specifying a location means specifying the zone and the x, y coordinate in that plane. The projection from spheroid to a UTM zone is some parameterization of the transverse Mercator projection. The parameters vary by nation or region or mapping system.

Tissots indicatrix

In cartography, a Tissot's indicatrix is a mathematical contrivance presented by French mathematician Nicolas Auguste Tissot in 1859 and 1871 in order to characterize local distortions due to map projection. It is the geometry that results from projecting a circle of infinitesimal radius from a curved geometric model, such as a globe, onto a map. Tissot proved that the resulting diagram is an ellipse whose axes indicate the two principal directions along which scale is maximal and minimal at that point on the map.

Lambert conformal conic projection map projection

A Lambert conformal conic projection (LCC) is a conic map projection used for aeronautical charts, portions of the State Plane Coordinate System, and many national and regional mapping systems. It is one of seven projections introduced by Johann Heinrich Lambert in his 1772 publication Anmerkungen und Zusätze zur Entwerfung der Land- und Himmelscharten.

Miller cylindrical projection Map projection

The Miller cylindrical projection is a modified Mercator projection, proposed by Osborn Maitland Miller in 1942. The latitude is scaled by a factor of ​45, projected according to Mercator, and then the result is multiplied by ​54 to retain scale along the equator. Hence:

Space-oblique Mercator projection map projection

Space-oblique Mercator projection is a map projection devised in the 1970s for preparing maps from Earth-survey satellite data. It is a generalization of the oblique Mercator projection that incorporates the time evolution of a given satellite ground track to optimize its representation on the map. The oblique Mercator projection, on the other hand, optimizes for a given geodesic.

Cassini projection map projection

The Cassini projection is a map projection described by César-François Cassini de Thury in 1745. It is the transverse aspect of the equirectangular projection, in that the globe is first rotated so the central meridian becomes the "equator", and then the normal equirectangular projection is applied. Considering the earth as a sphere, the projection is composed of the operations:

Cylindrical equal-area projection

In cartography, the cylindrical equal-area projection is a family of cylindrical, equal-area map projections.

In 1989 Bernard Russel Bowring gave formulas for the Transverse Mercator that are simpler to program but retain millimeter accuracy. Bowring rewrote the fourth order Redfearn series in a more compact notation by replacing the spherical terms, i.e. those independent of ellipticity, by the exact expressions used in the spherical transverse Mercator projection. There was no gain in accuracy since the elliptic terms were still truncated at the 1mm level. Such modifications were of possible use when computing resources were minimal.

The integral of the secant function of trigonometry was the subject of one of the "outstanding open problems of the mid-seventeenth century", solved in 1668 by James Gregory. In 1599, Edward Wright evaluated the integral by numerical methods – what today we would call Riemann sums. He wanted the solution for the purposes of cartography – specifically for constructing an accurate Mercator projection. In the 1640s, Henry Bond, a teacher of navigation, surveying, and other mathematical topics, compared Wright's numerically computed table of values of the integral of the secant with a table of logarithms of the tangent function, and consequently conjectured that

Central cylindrical projection

The central cylindrical projection is a perspective cylindrical map projection. It corresponds to projecting the Earth's surface onto a cylinder tangent to the equator as if from a light source at Earth's center. The cylinder is then cut along one of the projected meridians and unrolled into a flat map.

Gall stereographic projection Cylindrical map projection

The Gall stereographic projection, presented by James Gall in 1855, is a cylindrical projection. It is neither equal-area nor conformal but instead tries to balance the distortion inherent in any projection.

Web Mercator projection Mercator projection variant

Web Mercator, Google Web Mercator, Spherical Mercator, WGS 84 Web Mercator or WGS 84/Pseudo-Mercator is a variant of the Mercator projection and is the de facto standard for Web mapping applications. It rose to prominence when Google Maps adopted it in 2005. It is used by virtually all major online map providers, including Google Maps, Mapbox, Bing Maps, OpenStreetMap, Mapquest, Esri, and many others. Its official EPSG identifier is EPSG:3857, although others have been used historically.

Rectangular polyconic projection

The rectangular polyconic projection is a map projection was first mentioned in 1853 by the U.S. Coast Survey, where it was developed and used for portions of the U.S. exceeding about one square degree. It belongs to the polyconic projection class, which consists of map projections whose parallels are non-concentric circular arcs except for the equator, which is straight. Sometimes the rectangular polyconic is called the War Office projection due to its use by the British War Office for topographic maps. It is not used much these days, with practically all military grid systems having moved onto conformal projection systems, typically modeled on the transverse Mercator projection.

References