Jump to content

Contour line

From Wikipedia, the free encyclopedia
(Redirected from Aclinic line)
The bottom part of the diagram shows some contour lines with a straight line running through the location of the maximum value. The curve at the top represents the values along that straight line.
A three-dimensional surface, whose contour graph is below.
A two-dimensional contour graph of the three-dimensional surface in the above picture.

A contour line (also isoline, isopleth, isoquant or isarithm) of a function of two variables is a curve along which the function has a constant value, so that the curve joins points of equal value.[1][2] It is a plane section of the three-dimensional graph of the function parallel to the -plane. More generally, a contour line for a function of two variables is a curve connecting points where the function has the same particular value.[2]

In cartography, a contour line (often just called a "contour") joins points of equal elevation (height) above a given level, such as mean sea level.[3] A contour map is a map illustrated with contour lines, for example a topographic map, which thus shows valleys and hills, and the steepness or gentleness of slopes.[4] The contour interval of a contour map is the difference in elevation between successive contour lines.[5]

The gradient of the function is always perpendicular to the contour lines. When the lines are close together the magnitude of the gradient is large: the variation is steep. A level set is a generalization of a contour line for functions of any number of variables.

Contour lines are curved, straight or a mixture of both lines on a map describing the intersection of a real or hypothetical surface with one or more horizontal planes. The configuration of these contours allows map readers to infer the relative gradient of a parameter and estimate that parameter at specific places. Contour lines may be either traced on a visible three-dimensional model of the surface, as when a photogrammetrist viewing a stereo-model plots elevation contours, or interpolated from the estimated surface elevations, as when a computer program threads contours through a network of observation points of area centroids. In the latter case, the method of interpolation affects the reliability of individual isolines and their portrayal of slope, pits and peaks.[6]

History

[edit]
Edmond Halley's New and Correct Chart Shewing the Variations of the Compass (1701)

The idea of lines that join points of equal value was rediscovered several times. The oldest known isobath (contour line of constant depth) is found on a map dated 1584 of the river Spaarne, near Haarlem, by Dutchman Pieter Bruinsz.[7] In 1701, Edmond Halley used such lines (isogons) on a chart of magnetic variation.[8] The Dutch engineer Nicholas Cruquius drew the bed of the river Merwede with lines of equal depth (isobaths) at intervals of 1 fathom in 1727, and Philippe Buache used them at 10-fathom intervals on a chart of the English Channel that was prepared in 1737 and published in 1752. Such lines were used to describe a land surface (contour lines) in a map of the Duchy of Modena and Reggio by Domenico Vandelli in 1746, and they were studied theoretically by Ducarla in 1771, and Charles Hutton used them in the Schiehallion experiment. In 1791, a map of France by J. L. Dupain-Triel used contour lines at 20-metre intervals, hachures, spot-heights and a vertical section. In 1801, the chief of the French Corps of Engineers, Haxo, used contour lines at the larger scale of 1:500 on a plan of his projects for Rocca d'Anfo, now in northern Italy, under Napoleon.[9][10][11]

By around 1843, when the Ordnance Survey started to regularly record contour lines in Great Britain and Ireland, they were already in general use in European countries. Isobaths were not routinely used on nautical charts until those of Russia from 1834, and those of Britain from 1838.[9][12][13]

As different uses of the technique were invented independently, cartographers began to recognize a common theme, and debated what to call these "lines of equal value" generally. The word isogram (from Ancient Greek ἴσος (isos) 'equal' and γράμμα (gramma) 'writing, drawing') was proposed by Francis Galton in 1889 for lines indicating equality of some physical condition or quantity,[14] though isogram can also refer to a word without a repeated letter. As late as 1944, John K. Wright still preferred isogram, but it never attained wide usage. During the early 20th century, isopleth (πλῆθος, plethos, 'amount') was being used by 1911 in the United States, while isarithm (ἀριθμός, arithmos, 'number') had become common in Europe. Additional alternatives, including the Greek-English hybrid isoline and isometric line (μέτρον, metron, 'measure'), also emerged. Despite attempts to select a single standard, all of these alternatives have survived to the present.[15][16]

When maps with contour lines became common, the idea spread to other applications. Perhaps the latest to develop are air quality and noise pollution contour maps, which first appeared in the United States in approximately 1970, largely as a result of national legislation requiring spatial delineation of these parameters.

Types

[edit]

Contour lines are often given specific names beginning with "iso-" according to the nature of the variable being mapped, although in many usages the phrase "contour line" is most commonly used. Specific names are most common in meteorology, where multiple maps with different variables may be viewed simultaneously. The prefix "'iso-" can be replaced with "isallo-" to specify a contour line connecting points where a variable changes at the same rate during a given time period.

An isogon (from Ancient Greek γωνία (gonia) 'angle') is a contour line for a variable which measures direction. In meteorology and in geomagnetics, the term isogon has specific meanings which are described below. An isocline (κλίνειν, klinein, 'to lean or slope') is a line joining points with equal slope. In population dynamics and in geomagnetics, the terms isocline and isoclinic line have specific meanings which are described below.

Equidistant points

[edit]

A curve of equidistant points is a set of points all at the same distance from a given point, line, or polyline. In this case the function whose value is being held constant along a contour line is a distance function.

Isopleths

[edit]

In 1944, John K. Wright proposed that the term isopleth be used for contour lines that depict a variable which cannot be measured at a point, but which instead must be calculated from data collected over an area, as opposed to isometric lines for variables that could be measured at a point; this distinction has since been followed generally.[16][17] An example of an isopleth is population density, which can be calculated by dividing the population of a census district by the surface area of that district. Each calculated value is presumed to be the value of the variable at the centre of the area, and isopleths can then be drawn by a process of interpolation. The idea of an isopleth map can be compared with that of a choropleth map.[18][19]

In meteorology, the word isopleth is used for any type of contour line.[20]

Meteorology

[edit]
Isohyetal map of precipitation

Meteorological contour lines are based on interpolation of the point data received from weather stations and weather satellites. Weather stations are seldom exactly positioned at a contour line (when they are, this indicates a measurement precisely equal to the value of the contour). Instead, lines are drawn to best approximate the locations of exact values, based on the scattered information points available.

Meteorological contour maps may present collected data such as actual air pressure at a given time, or generalized data such as average pressure over a period of time, or forecast data such as predicted air pressure at some point in the future.

Thermodynamic diagrams use multiple overlapping contour sets (including isobars and isotherms) to present a picture of the major thermodynamic factors in a weather system.

Barometric pressure

[edit]
Video loop of isallobars showing the motion of a cold front

An isobar (from Ancient Greek βάρος (baros) 'weight') is a line of equal or constant pressure on a graph, plot, or map; an isopleth or contour line of pressure. More accurately, isobars are lines drawn on a map joining places of equal average atmospheric pressure reduced to sea level for a specified period of time. In meteorology, the barometric pressures shown are reduced to sea level, not the surface pressures at the map locations.[21] The distribution of isobars is closely related to the magnitude and direction of the wind field, and can be used to predict future weather patterns. Isobars are commonly used in television weather reporting.

Isallobars are lines joining points of equal pressure change during a specific time interval.[22] These can be divided into anallobars, lines joining points of equal pressure increase during a specific time interval,[23] and katallobars, lines joining points of equal pressure decrease.[24] In general, weather systems move along an axis joining high and low isallobaric centers.[25] Isallobaric gradients are important components of the wind as they increase or decrease the geostrophic wind.

An isopycnal is a line of constant density. An isoheight or isohypse is a line of constant geopotential height on a constant pressure surface chart. Isohypse and isoheight are simply known as lines showing equal pressure on a map.

[edit]
The 10 °C (50 °F) mean isotherm in July, marked by the red line, is commonly used to define the border of the Arctic region

An isotherm (from Ancient Greek θέρμη (thermē) 'heat') is a line that connects points on a map that have the same temperature. Therefore, all points through which an isotherm passes have the same or equal temperatures at the time indicated.[26][2] An isotherm at 0 °C is called the freezing level. The term lignes isothermes (or lignes d'égale chaleur) was coined by the Prussian geographer and naturalist Alexander von Humboldt, who as part of his research into the geographical distribution of plants published the first map of isotherms in Paris, in 1817.[27][28]

An isocheim is a line of equal mean winter temperature, and an isothere is a line of equal mean summer temperature.

An isohel (ἥλιος, helios, 'Sun') is a line of equal or constant solar radiation.

An isogeotherm is a line of equal temperature beneath the Earth's surface.

Rainfall and air moisture

[edit]

An isohyet or isohyetal line (from Ancient Greek ὑετός (huetos) 'rain') is a line on a map joining points of equal rainfall in a given period. A map with isohyets is called an isohyetal map.

An isohume is a line of constant relative humidity, while an isodrosotherm (from Ancient Greek δρόσος (drosos) 'dew' and θέρμη (therme) 'heat') is a line of equal or constant dew point.

An isoneph is a line indicating equal cloud cover.

An isochalaz is a line of constant frequency of hail storms, and an isobront is a line drawn through geographical points at which a given phase of thunderstorm activity occurred simultaneously.

Snow cover is frequently shown as a contour-line map.

Wind

[edit]

An isotach (from Ancient Greek ταχύς (tachus) 'fast') is a line joining points with constant wind speed. In meteorology, the term isogon refers to a line of constant wind direction.

Freeze and thaw

[edit]

An isopectic line denotes equal dates of ice formation each winter, and an isotac denotes equal dates of thawing.

Physical geography and oceanography

[edit]

Elevation and depth

[edit]
Topographic map of Stowe, Vermont. The brown contour lines represent the elevation. The contour interval is 20 feet.

Contours are one of several common methods used to denote elevation or altitude and depth on maps. From these contours, a sense of the general terrain can be determined. They are used at a variety of scales, from large-scale engineering drawings and architectural plans, through topographic maps and bathymetric charts, up to continental-scale maps.

"Contour line" is the most common usage in cartography, but isobath for underwater depths on bathymetric maps and isohypse for elevations are also used.

In cartography, the contour interval is the elevation difference between adjacent contour lines. The contour interval should be the same over a single map. When calculated as a ratio against the map scale, a sense of the hilliness of the terrain can be derived.

Interpretation
[edit]

There are several rules to note when interpreting terrain contour lines:

  • The rule of Vs: sharp-pointed vees usually are in stream valleys, with the drainage channel passing through the point of the vee, with the vee pointing upstream. This is a consequence of erosion.
  • The rule of Os: closed loops are normally uphill on the inside and downhill on the outside, and the innermost loop is the highest area. If a loop instead represents a depression, some maps note this by short lines called hachures which are perpendicular to the contour and point in the direction of the low.[29] (The concept is similar to but distinct from hachures used in hachure maps.)
  • Spacing of contours: close contours indicate a steep slope; distant contours a shallow slope. Two or more contour lines merging indicates a cliff. By counting the number of contours that cross a segment of a stream, the stream gradient can be approximated.

Of course, to determine differences in elevation between two points, the contour interval, or distance in altitude between two adjacent contour lines, must be known, and this is normally stated in the map key. Usually contour intervals are consistent throughout a map, but there are exceptions. Sometimes intermediate contours are present in flatter areas; these can be dashed or dotted lines at half the noted contour interval. When contours are used with hypsometric tints on a small-scale map that includes mountains and flatter low-lying areas, it is common to have smaller intervals at lower elevations so that detail is shown in all areas. Conversely, for an island which consists of a plateau surrounded by steep cliffs, it is possible to use smaller intervals as the height increases.[30]

Electrostatics

[edit]

An isopotential map is a measure of electrostatic potential in space, often depicted in two dimensions with the electrostatic charges inducing that electric potential. The term equipotential line or isopotential line refers to a curve of constant electric potential. Whether crossing an equipotential line represents ascending or descending the potential is inferred from the labels on the charges. In three dimensions, equipotential surfaces may be depicted with a two dimensional cross-section, showing equipotential lines at the intersection of the surfaces and the cross-section.

The general mathematical term level set is often used to describe the full collection of points having a particular potential, especially in higher dimensional space.

Magnetism

[edit]
Isogonic lines for the year 2000. The agonic lines are thicker and labeled with "0".

In the study of the Earth's magnetic field, the term isogon or isogonic line refers to a line of constant magnetic declination, the variation of magnetic north from geographic north. An agonic line is drawn through points of zero magnetic declination. An isoporic line refers to a line of constant annual variation of magnetic declination .[31]

An isoclinic line connects points of equal magnetic dip, and an aclinic line is the isoclinic line of magnetic dip zero.

An isodynamic line (from δύναμις or dynamis meaning 'power') connects points with the same intensity of magnetic force.

Oceanography

[edit]

Besides ocean depth, oceanographers use contour to describe diffuse variable phenomena much as meteorologists do with atmospheric phenomena. In particular, isobathytherms are lines showing depths of water with equal temperature, isohalines show lines of equal ocean salinity, and isopycnals are surfaces of equal water density.

Geology

[edit]

Various geological data are rendered as contour maps in structural geology, sedimentology, stratigraphy and economic geology. Contour maps are used to show the below ground surface of geologic strata, fault surfaces (especially low angle thrust faults) and unconformities. Isopach maps use isopachs (lines of equal thickness) to illustrate variations in thickness of geologic units.

Environmental science

[edit]

In discussing pollution, density maps can be very useful in indicating sources and areas of greatest contamination. Contour maps are especially useful for diffuse forms or scales of pollution. Acid precipitation is indicated on maps with isoplats. Some of the most widespread applications of environmental science contour maps involve mapping of environmental noise (where lines of equal sound pressure level are denoted isobels[32]), air pollution, soil contamination, thermal pollution and groundwater contamination. By contour planting and contour ploughing, the rate of water runoff and thus soil erosion can be substantially reduced; this is especially important in riparian zones.

Ecology

[edit]

An isoflor is an isopleth contour connecting areas of comparable biological diversity. Usually, the variable is the number of species of a given genus or family that occurs in a region. Isoflor maps are thus used to show distribution patterns and trends such as centres of diversity.[33]

Social sciences

[edit]
From economics, an indifference map with three indifference curves shown. All points on a particular indifference curve have the same value of the utility function, whose values implicitly come out of the page in the unshown third dimension.

In economics, contour lines can be used to describe features which vary quantitatively over space. An isochrone shows lines of equivalent drive time or travel time to a given location and is used in the generation of isochrone maps. An isotim shows equivalent transport costs from the source of a raw material, and an isodapane shows equivalent cost of travel time.

A single production isoquant (convex) and a single isocost curve (linear). Labor usage is plotted horizontally and physical capital usage is plotted vertically.

Contour lines are also used to display non-geographic information in economics. Indifference curves (as shown at left) are used to show bundles of goods to which a person would assign equal utility. An isoquant (in the image at right) is a curve of equal production quantity for alternative combinations of input usages, and an isocost curve (also in the image at right) shows alternative usages having equal production costs.

In political science an analogous method is used in understanding coalitions (for example the diagram in Laver and Shepsle's work[34]).

In population dynamics, an isocline shows the set of population sizes at which the rate of change, or partial derivative, for one population in a pair of interacting populations is zero.

Statistics

[edit]

In statistics, isodensity lines [35] or isodensanes are lines that join points with the same value of a probability density. Isodensanes are used to display bivariate distributions. For example, for a bivariate elliptical distribution the isodensity lines are ellipses.

Thermodynamics, engineering, and other sciences

[edit]

Various types of graphs in thermodynamics, engineering, and other sciences use isobars (constant pressure), isotherms (constant temperature), isochors (constant specific volume), or other types of isolines, even though these graphs are usually not related to maps. Such isolines are useful for representing more than two dimensions (or quantities) on two-dimensional graphs. Common examples in thermodynamics are some types of phase diagrams.

Isoclines are used to solve ordinary differential equations.

In interpreting radar images, an isodop is a line of equal Doppler velocity, and an isoecho is a line of equal radar reflectivity.

In the case of hybrid contours, energies of hybrid orbitals and the energies of pure atomic orbitals are plotted. The graph obtained is called hybrid contour.

Other phenomena

[edit]

Algorithms

[edit]

Graphical design

[edit]

To maximize readability of contour maps, there are several design choices available to the map creator, principally line weight, line color, line type and method of numerical marking.

Line weight is simply the darkness or thickness of the line used. This choice is made based upon the least intrusive form of contours that enable the reader to decipher the background information in the map itself. If there is little or no content on the base map, the contour lines may be drawn with relatively heavy thickness. Also, for many forms of contours such as topographic maps, it is common to vary the line weight and/or color, so that a different line characteristic occurs for certain numerical values. For example, in the topographic map above, the even hundred foot elevations are shown in a different weight from the twenty foot intervals.

Line color is the choice of any number of pigments that suit the display. Sometimes a sheen or gloss is used as well as color to set the contour lines apart from the base map. Line colour can be varied to show other information.

Line type refers to whether the basic contour line is solid, dashed, dotted or broken in some other pattern to create the desired effect. Dotted or dashed lines are often used when the underlying base map conveys very important (or difficult to read) information. Broken line types are used when the location of the contour line is inferred.

Numerical marking is the manner of denoting the arithmetical values of contour lines. This can be done by placing numbers along some of the contour lines, typically using interpolation for intervening lines. Alternatively a map key can be produced associating the contours with their values.

If the contour lines are not numerically labeled and adjacent lines have the same style (with the same weight, color and type), then the direction of the gradient cannot be determined from the contour lines alone. However, if the contour lines cycle through three or more styles, then the direction of the gradient can be determined from the lines. The orientation of the numerical text labels is often used to indicate the direction of the slope.

Plan view versus profile view

[edit]

Most commonly contour lines are drawn in plan view, or as an observer in space would view the Earth's surface: ordinary map form. However, some parameters can often be displayed in profile view showing a vertical profile of the parameter mapped. Some of the most common parameters mapped in profile are air pollutant concentrations and sound levels. In each of those cases it may be important to analyze (air pollutant concentrations or sound levels) at varying heights so as to determine the air quality or noise health effects on people at different elevations, for example, living on different floor levels of an urban apartment. In actuality, both plan and profile view contour maps are used in air pollution and noise pollution studies.

Contour map labeled aesthetically in an "elevation up" manner.

Labeling contour maps

[edit]

Labels are a critical component of elevation maps. A properly labeled contour map helps the reader to quickly interpret the shape of the terrain. If numbers are placed close to each other, it means that the terrain is steep. Labels should be placed along a slightly curved line "pointing" to the summit or nadir, from several directions if possible, making the visual identification of the summit or nadir easy.[36][37] Contour labels can be oriented so a reader is facing uphill when reading the label.

Manual labeling of contour maps is a time-consuming process, however, there are a few software systems that can do the job automatically and in accordance with cartographic conventions, called automatic label placement.

See also

[edit]

References

[edit]
  1. ^ Courant, Richard, Herbert Robbins, and Ian Stewart. What Is Mathematics?: An Elementary Approach to Ideas and Methods. New York: Oxford University Press, 1996. p. 344.
  2. ^ a b c Hughes-Hallett, Deborah; McCallum, William G.; Gleason, Andrew M. (2013). Calculus : Single and Multivariable (6 ed.). John wiley. ISBN 978-0470-88861-2.
  3. ^ "Definition of contour line". Dictionary.com. Retrieved 2022-04-04.
  4. ^ "Definition of CONTOUR MAP". Merriam-Webster. Retrieved 2022-04-04.
  5. ^ Tracy, John C. Plane Surveying; A Text-Book and Pocket Manual. New York: J. Wiley & Sons, 1907. p. 337.
  6. ^ Davis, John C., 1986, Statistics and data analysis in geology, Wiley ISBN 0-471-08079-9
  7. ^ Morato-Moreno, Manuel (2017). "Orígenes de la representación topográfica del terreno en algunos mapas hispanoamericanos del s. XVI". Boletín de la Asociación de Geógrafos Españoles. doi:10.21138/bage.2414.
  8. ^ Thrower, N. J. W. Maps and Civilization: Cartography in Culture and Society, University of Chicago Press, 1972, revised 1996, page 97; and Jardine, Lisa Ingenious Pursuits: Building the Scientific Revolution, Little, Brown, and Company, 1999, page 31.
  9. ^ a b R. A. Skelton, "Cartography", History of Technology, Oxford, vol. 6, pp. 612–614, 1958.
  10. ^ Colonel Berthaut, La Carte de France, vol. 1, p. 139, quoted by Close.
  11. ^ C. Hutton, "An account of the calculations made from the survey and measures taken at Schehallien, in order to ascertain the mean density of the Earth", Philosophical Transactions of the Royal Society of London, vol. 68, pp. 756757
  12. ^ C. Close, The Early Years of the Ordnance Survey, 1926, republished by David and Charles, 1969, ISBN 0-7153-4477-3, pp. 141–144.
  13. ^ T. Owen and E. Pilbeam, Ordnance Survey: Map Makers to Britain since 1791, HMSO, 1992, ISBN 0-11-701507-5.
  14. ^ Galton, Francis (1889). "On the Principle and Methods of Assigning Marks for Bodily Efficiency". Nature. 40 (1044): 651. Bibcode:1889Natur..40..649.. doi:10.1038/040649a0. S2CID 3996216.
  15. ^ Wright, John K. (Apr 1930). "Isopleth as a Generic Term". Geographical Review. 20 (2): 341. JSTOR 208890.
  16. ^ a b Wright, John K. (Oct 1944). "The Terminology of Certain Map Symbols". Geographical Review. 34 (4): 653–654. Bibcode:1944GeoRv..34..653W. doi:10.2307/210035. JSTOR 210035.
  17. ^ Robinson AH (1971). "The genealogy of the isopleth". Cartographic Journal. 8 (1): 49–53. Bibcode:1971CartJ...8...49R. doi:10.1179/caj.1971.8.1.49.
  18. ^ T. Slocum, R. McMaster, F. Kessler, and H. Howard, Thematic Cartography and Geographic Visualization, 2nd edition, Pearson, 2005, ISBN 0-13-035123-7, p. 272.
  19. ^ ArcGIS, Isopleth: Contours, 2013.
  20. ^ NOAA's National Weather Service, Glossary.
  21. ^ Edward J. Hopkins, Ph.D. (1996-06-10). "Surface Weather Analysis Chart". University of Wisconsin. Retrieved 2007-05-10.
  22. ^ World Meteorological Organisation. "Isallobar". Eumetcal. Archived from the original on 16 April 2014. Retrieved 12 April 2014.
  23. ^ World Meteorological Organisation. "Anallobar". Eumetcal. Archived from the original on 24 September 2015. Retrieved 12 April 2014.
  24. ^ World Meteorological Organisation. "Katallobar". Eumetcal. Archived from the original on 5 February 2008. Retrieved 12 April 2014.
  25. ^ "Forecasting weather system movement with pressure tendency". Chapter 13 – Weather Forecasting. Lyndon State College Atmospheric Sciences. Retrieved 12 April 2014.
  26. ^ DataStreme Atmosphere (2008-04-28). "Air Temperature Patterns". American Meteorological Society. Archived from the original on 2008-05-11. Retrieved 2010-02-07.
  27. ^ Daum, Andreas W. (2024). Alexander von Humboldt: A Concise Biography. Trans. Robert Savage. Princeton, N.J.: Princeton University Press. pp. 106–107. ISBN 978-0-691-24736-6.
  28. ^ Munzar, Jan (1967-09-01). "Alexander Von Humboldt and His Isotherms". Weather. 22 (9): 360–363. Bibcode:1967Wthr...22..360M. doi:10.1002/j.1477-8696.1967.tb02989.x. ISSN 1477-8696.
  29. ^ Leveson, David J. (2002). "Depression Contours – Getting Into and Out of a Hole". City University of New York.
  30. ^ Sark (Sercq), D Survey, Ministry of Defence, Series M 824, Sheet Sark, Edition 4 GSGS, 1965, OCLC OCLC 27636277. Scale 1:10,560. Contour intervals: 50 feet up to 200, 20 feet from 200 to 300, and 10 feet above 300.
  31. ^ "isoporic line". historicalcharts.noaa.gov. 1946. Archived from the original on 2015-07-21. Retrieved 2015-07-20.
  32. ^ "Isobel". www.sfu.ca. 2005-01-05. Retrieved 2010-04-25.
  33. ^ Specht, Raymond. Heathlands and related shrublands: Analytical studies. Elsevier. pp. 219–220.
  34. ^ Laver, Michael and Kenneth A. Shepsle (1996) Making and breaking governments pictures.
  35. ^ Fernández, Antonio (2011). "A Generalized Regression Methodology for Bivariate Heteroscedastic Data" (PDF). Communications in Statistics – Theory and Methods. 40 (4): 598–621. doi:10.1080/03610920903444011. S2CID 55887263.
  36. ^ Imhof, E., "Die Anordnung der Namen in der Karte," Annuaire International de Cartographie II, Orell-Füssli Verlag, Zürich, 93–129, 1962.
  37. ^ Freeman, H., "Computer Name Placement," ch. 29, in Geographical Information Systems, 1, D.J. Maguire, M.F. Goodchild, and D.W. Rhind, John Wiley, New York, 1991, 449–460.
[edit]