Planetary-scale answers, unlocked.
A Hands-On Guide for Working with Large-Scale Spatial Data. Learn more.
Authors
A map projection is a set of formulas that turns positions on the curved Earth, given as latitude and longitude, into x and y positions on a flat map. A globe cannot be flattened without stretching or tearing, so every projection distorts something. This article tests one claim: no flat map keeps both area and shape. An equal-area projection keeps area and stretches shapes. A conformal projection keeps local shapes and inflates area. A compromise projection such as Robinson gets both wrong, by smaller amounts.
The meaning of projection in cartography is a rule that transfers the graticule from a globe onto a flat surface. The graticule is the grid of latitude and longitude lines. John P. Snyder's Map Projections: A Working Manual (US Geological Survey (USGS) Professional Paper 1395, 1987) is the standard reference for projection formulas. It notes that a few projections are geometric, like a light casting the globe's lines onto paper. Most are mathematical, with spacing chosen to keep a property such as area.
A projection is one part of a coordinate reference system, which also fixes the datum, units and axis order. That article covers datums, EPSG codes, the registry numbers for coordinate systems, and Universal Transverse Mercator (UTM) zones. This one covers the projection step and the distortion it causes.
Snyder defines a developable surface as one that can be turned into a plane without distortion. He names three: the cylinder, the cone and the plane. A globe is not developable. Snyder states the result: "no map projection shows scale correctly throughout the map". He also writes: "No map can be both equal-area and conformal."
Every geospatial dataset drawn or measured on a flat grid has passed through a projection. A GeoTIFF records the projection of its pixel grid in its tags. Landsat and Sentinel-2 scenes used in satellite imagery analysis arrive on UTM grids.
Snyder lists the properties a mapmaker weighs when choosing a projection. Four of them describe distortion.
A projection that keeps none of the four exactly is a compromise projection. Snyder's table of suggested world maps lists Miller Cylindrical and Robinson under "compromise distortion".
Nicolas Auguste Tissot published his analysis of projection distortion in 1859 and 1881, as Snyder recounts. He showed that an infinitely small circle on the Earth projects as a perfect ellipse on any map. On a conformal map, that ellipse is a circle, larger or smaller than the original. On an equal-area map, the ellipse keeps the circle's area and changes its shape. The ellipse's longest and shortest axes measure the scale distortion and the largest angle error at that point. Mulcahy and Clarke's Symbolization of Map Projection Distortion: A Review (Cartography and Geographic Information Science, 2001) surveys this and other ways to draw distortion on a map.
The Wherobots test below uses a finite version of the indicatrix: a circle 500 km in radius, drawn on the ellipsoid and then projected.
Snyder groups most projections by the developable surface they are built on. National Geographic's guide to selecting a map projection names the same three basic kinds: planar, conical and cylindrical.
A cylinder or cone can also cut through the globe along two parallels, which gives two standard parallels of true scale. Turning the surface sideways or tilting it gives transverse and oblique aspects.
An equal-area projection draws every region at its true size relative to every other region, as defined in the four properties above. Common examples are the Albers equal-area conic for mid-latitude regions, the Mollweide and the Equal Earth projection for world maps, and the Gall-Peters cylindrical projection. Each keeps area by stretching shapes, since no map can be both equal-area and conformal.
Area calculations on projected x and y need an equal-area projection. ST_Area on projected coordinates returns the area the map draws. In the Wherobots test below, Mercator returned 16.4 times Greenland's true area and Robinson 1.599 times. Equal Earth (EPSG:8857) and the equal-area cylindrical projection returned exactly 1 for all eight countries. ST_AreaSpheroid measures on the ellipsoid and skips the projection step.
ST_Area
ST_AreaSpheroid
The table lists what eight common projections keep and give up. The next section tests all of them except Goode homolosine.
Gerardus Mercator presented his projection in 1569 to aid navigation, according to Snyder. On a Mercator map, a rhumb line, a course of constant compass bearing, is a straight line. The poles lie at infinity. PROJ, the open-source projection library, calls Mercator "inappropriate for world maps because of the gross distortions in area" in its Mercator documentation. It gives Greenland as the example: drawn larger than South America, with one eighth of its area. The EPSG:3395 World Mercator record lists its scope as small-scale conformal mapping between 80°S and 84°N.
Many online maps use a variation of Mercator, as the National Geographic guide notes. Battersby, Finn, Usery and Yamamoto examined that choice in Implications of Web Mercator and Its Use in Online Mapping (Cartographica, 2014).
The PROJ documentation for Robinson classes it as pseudocylindrical and implements it as a spherical projection only. Snyder lists it as a compromise projection for world maps. National Geographic used the Van der Grinten projection from 1922 until 1988. In 1995, the Winkel tripel replaced the Robinson projection on the Society's world maps, according to its guide.
Bojan Šavrič, Tom Patterson and Bernhard Jenny introduced the Equal Earth map projection in the International Journal of Geographical Information Science in 2018. The PROJ documentation for Equal Earth says it was inspired by Robinson and keeps the relative size of areas. It has an ellipsoidal form. EPSG registered it as EPSG:8857 on WGS 84, the World Geodetic System 1984. Its stated scope is small-scale equal-area mapping.
Carl B. Mollweide presented his equal-area projection of the Earth inside an ellipse in 1805, Snyder records. National Geographic used an interrupted Mollweide for its September 2012 map of the ocean floors. The PROJ documentation for Mollweide lists only a spherical form.
Johann Heinrich Lambert described the cylindrical equal-area projection in 1772. James Gall's version, with true scale along the 45th parallels north and south, dates from 1855. He described it in Use of cylindrical projections for geographical, astronomical, and scientific purposes (Scottish Geographical Magazine, 1885). Snyder reports that Arno Peters of Germany presented a near-identical projection in 1967 as an original concept. The PROJ documentation for the equal-area cylindrical projection lists Gall-Peters as the setting with true scale at 45°.
J. Paul Goode of the University of Chicago interrupted the sinusoidal projection in 1916 to reduce its distortion, Snyder writes. He then joined the sinusoidal projection between about 40°N and 40°S with the Mollweide projection beyond. He presented the result in The homolosine projection: a new device for portraying the earth's surface entire (Annals of the Association of American Geographers, 1925). The PROJ documentation for the interrupted Goode homolosine describes six sinusoidal and six Mollweide regions. They meet at 40° 44' 11.8", the latitude where the two scales are equal. The lobes are placed to show the land.
An azimuthal projection keeps directions from its center. Snyder notes that of the common azimuthal projections, one is equal-area, one is conformal and one is equidistant. Lambert presented the azimuthal equal-area projection in 1772. On the azimuthal equidistant projection, distances measured from the center are true. Snyder adds that the center is its only point without distortion. Mercator's 1569 world map carried polar azimuthal equidistant insets.
The test measures area and shape in the same projections, on the same reference.
ST_Transform reprojects a geometry, with the target given as an EPSG code or a PROJ string. The reference is ST_AreaSpheroid and ST_DistanceSpheroid, which measure on the WGS 84 ellipsoid with no projection.
ST_Transform
ST_DistanceSpheroid
ST_Transform accepted every projection below. It rejected the interrupted Goode homolosine (+proj=igh) and the Winkel tripel (+proj=wintri) with "Unknown projection", so neither is in the test.
+proj=igh
+proj=wintri
WITH c AS ( SELECT country, subtype, geometry AS g FROM wherobots_open_data.overture_maps_foundation.divisions_division_area WHERE subtype IN ('country','dependency') AND class = 'land' AND country IN ('CD','BR','IN','AU','DZ','DE','CA','GL') ), m AS ( SELECT country, subtype, g, ST_AreaSpheroid(g) AS true_m2, ST_Area(ST_Transform(g, 'EPSG:4326', 'EPSG:3395')) AS merc, ST_Area(ST_Transform(g, 'EPSG:4326', '+proj=eqc +lon_0=0 +datum=WGS84 +units=m')) AS eqc, ST_Area(ST_Transform(g, 'EPSG:4326', '+proj=robin +lon_0=0 +datum=WGS84 +units=m')) AS robin, ST_Area(ST_Transform(g, 'EPSG:4326', 'EPSG:8857')) AS eqearth, ST_Area(ST_Transform(g, 'EPSG:4326', '+proj=moll +lon_0=0 +datum=WGS84 +units=m')) AS moll, ST_Area(ST_Transform(g, 'EPSG:4326', '+proj=cea +lat_ts=45 +lon_0=0 +datum=WGS84 +units=m')) AS gallpeters FROM c ) SELECT country, subtype, ROUND(ST_Y(ST_Centroid(g)), 1) AS centroid_lat, ROUND(true_m2 / 1e6, 0) AS km2_ellipsoid, ROUND(merc / true_m2, 3) AS mercator, ROUND(eqc / true_m2, 3) AS plate_carree, ROUND(robin / true_m2, 3) AS robinson, ROUND(eqearth / true_m2, 4) AS equal_earth, ROUND(moll / true_m2, 4) AS mollweide, ROUND(gallpeters / true_m2, 4) AS gall_peters FROM m ORDER BY ABS(ST_Y(ST_Centroid(g)))
Each value is the projected area divided by the area on the ellipsoid, so 1 means the projection kept area. The query ran in 30.2 seconds.
Mercator's area ratio climbs with latitude, from 1.008 for the Democratic Republic of the Congo to 16.4 for Greenland. Overture lists Greenland as a dependency, so the query selects that subtype too. Plate carrée, which plots longitude and latitude directly as x and y, reaches 3.8. Robinson shrinks the tropics to 0.823 and enlarges Greenland to 1.599. Equal Earth and Gall-Peters read 1 for all eight.
Mollweide reads 0.9946 to 1.0066. PROJ implements Mollweide, and Robinson, as spherical projections. Applied to WGS 84 coordinates, which sit on a flattened ellipsoid, the Mollweide formulas come within 0.66% of equal area without reaching it. Equal Earth and the equal-area cylindrical projection have ellipsoidal forms in PROJ and read exactly 1.
Area alone does not show shape. This query draws a circle of 500 km radius on the ellipsoid with ST_Buffer(..., true) at four latitudes and projects it. Roundness is the shape's area divided by the area of the smallest circle that encloses it. A circle scores 1. An ellipse twice as wide as it is tall scores 0.5.
ST_Buffer(..., true)
WITH pts AS ( SELECT * FROM VALUES (0.0), (30.0), (60.0), (75.0) AS t(lat) ), c AS ( SELECT lat, ST_Buffer(ST_Point(20.0, lat), 500000, true) AS g FROM pts ), p AS ( SELECT lat, g, ST_AreaSpheroid(g) AS true_m2, proj, ST_Transform(g, 'EPSG:4326', code) AS pg FROM c LATERAL VIEW explode(map( 'mercator', 'EPSG:3395', 'plate_carree', '+proj=eqc +lon_0=0 +datum=WGS84 +units=m', 'robinson', '+proj=robin +lon_0=0 +datum=WGS84 +units=m', 'equal_earth', 'EPSG:8857', 'gall_peters', '+proj=cea +lat_ts=45 +lon_0=0 +datum=WGS84 +units=m', 'azimuthal_eqd_20E_0N', '+proj=aeqd +lat_0=0 +lon_0=20 +datum=WGS84 +units=m' )) e AS proj, code ) SELECT proj, lat, ROUND(ST_Area(pg) / true_m2, 3) AS area_ratio, ROUND(ST_Area(pg) / ST_Area(ST_MinimumBoundingCircle(pg)), 3) AS roundness, ROUND((ST_XMax(pg) - ST_XMin(pg)) / (ST_YMax(pg) - ST_YMin(pg)), 3) AS width_to_height FROM p ORDER BY proj, lat
The query ran in 2.2 seconds.
Mercator keeps the circle round at every latitude, at 0.961 or higher, and inflates its area to 15.515 times at 75°N. Its roundness dips at 75°N because the circle spans about 70.5°N to 79.5°N. Mercator's scale changes across that span, and conformal maps keep shape only for small areas.
Equal Earth and Gall-Peters keep the area at 1 at all four latitudes. At 75°N, Equal Earth draws the circle 4.568 times wider than tall, and Gall-Peters 7.441 times. Gall-Peters also stretches the circle at the equator, to half as wide as tall.
Robinson's area runs from 0.821 to 1.648, and its roundness from 0.936 down to 0.337.
The third query checks distance. It compares the ellipsoidal distance from Chicago to seven cities with the straight-line distance in three projections. The azimuthal equidistant projection is centered on Chicago.
WITH cities AS ( SELECT * FROM VALUES ('London', -0.1276, 51.5072), ('New York', -74.0060, 40.7128), ('Cairo', 31.2357, 30.0444), ('Tokyo', 139.6917, 35.6895), ('Sao Paulo', -46.6333, -23.5505), ('Sydney', 151.2093, -33.8688), ('Anchorage', -149.9003, 61.2181) AS t(city, lon, lat) ), o AS (SELECT ST_Point(-87.6298, 41.8781) AS chicago) SELECT city, ROUND(ST_DistanceSpheroid(chicago, ST_Point(lon, lat)) / 1e3, 0) AS km_ellipsoid, ROUND(ST_Distance(ST_Transform(chicago, 'EPSG:4326', '+proj=aeqd +lat_0=41.8781 +lon_0=-87.6298 +datum=WGS84 +units=m'), ST_Transform(ST_Point(lon, lat), 'EPSG:4326', '+proj=aeqd +lat_0=41.8781 +lon_0=-87.6298 +datum=WGS84 +units=m')) / 1e3, 0) AS km_azimuthal_equidistant, ROUND(ST_Distance(ST_Transform(chicago, 'EPSG:4326', 'EPSG:3395'), ST_Transform(ST_Point(lon, lat), 'EPSG:4326', 'EPSG:3395')) / 1e3, 0) AS km_mercator, ROUND(ST_Distance(ST_Transform(chicago, 'EPSG:4326', 'EPSG:8857'), ST_Transform(ST_Point(lon, lat), 'EPSG:4326', 'EPSG:8857')) / 1e3, 0) AS km_equal_earth FROM cities CROSS JOIN o ORDER BY km_ellipsoid
The query ran in 1.9 seconds.
The azimuthal equidistant map matches all seven distances to the kilometer. Mercator reads 25,321 km to Tokyo, against 10,161 km on the ellipsoid. Both world maps are cut at 180° longitude, so their straight line to Tokyo or Sydney runs east the long way round. Query B shows the cost of the azimuthal map away from its center. A circle 75° from the center reads 1.354 times its area, with roundness 0.749.
The results support the claim. No projection in the test kept both area and shape at 75°N.
Snyder writes that "it cannot be said that there is one 'best' projection for mapping." The best map projection to use depends on what the map has to keep. His table of suggested projections gives a starting point:
The International Cartographic Association's book Choosing a Map Projection (Springer, 2017) treats the same question at book length. Snyder and Voxland's An Album of Map Projections (USGS Professional Paper 1453, 1989) shows the candidates side by side.
For geospatial analysis, measure distance and area on the ellipsoid, as the test above did, and use a projection to draw the result. A spatial join needs both layers in the same coordinate system before it compares them. Reprojecting a raster such as a digital elevation model resamples its grid, which gives cells new values. Reprojecting vector data moves each vertex and resamples nothing, one of the differences between raster and vector data. Global grids face the same trade-off: H3 maps its hexagons from an icosahedron with a gnomonic projection, and their area varies by about a factor of two worldwide.
According to Snyder, Mercator probably set the spacing of his parallels by drawing. Edward Wright published the tables that fix the spacing in 1599. Johann Heinrich Lambert's 1772 Beiträge added seven projections, including the conformal conic, the cylindrical equal-area and the azimuthal equal-area.
In 1962, Waldo Tobler classified projections by whether each map coordinate depends on latitude, longitude or both. In his scheme, a pseudocylindrical projection has y depending on latitude alone and x on both.
Goldberg and Gott's Flexion and Skewness in Map Projections of the Earth (Cartographica, 2007) proposed two more measures for scoring world maps. Flexion measures how much a map bends straight lines, and skewness how lopsided it makes them.
Reproject and measure geometries in any of these projections with a Wherobots free trial at cloud.wherobots.com.
A map projection is a set of formulas that converts latitude and longitude on the curved Earth into x and y positions on a flat map. Because a globe cannot be flattened without stretching, every projection distorts area, shape, distance or direction somewhere. Most projections are designed to keep one or two of these properties. Compromise projections such as Robinson keep none exactly and limit the error in all of them.
No single projection is best for every map. Use an equal-area projection such as Equal Earth or Mollweide to compare sizes or map densities. Use a conformal projection such as Mercator for navigation and local shapes. Use an azimuthal equidistant projection centered on one place to show distances from it. Use a compromise projection such as Robinson for a general world reference map.
By construction, the main types are cylindrical, conic, azimuthal (also called planar) and pseudocylindrical. By the property they keep, projections are equal-area, conformal, equidistant, azimuthal, or compromise projections that keep none of these exactly and limit the distortion of all of them.
In geography, projection means transferring the grid of latitude and longitude lines from a globe onto a flat surface by a rule. A few projections are geometric, like a light casting the globe’s lines onto a sheet of paper. Most are mathematical formulas chosen to keep a property such as area or shape.
Mercator keeps local angles, so it stretches the map north-south by the same factor it stretches it east-west, and that factor grows toward the poles. In Wherobots, Greenland’s land area measured 16.4 times its true area in World Mercator (EPSG:3395), while the Democratic Republic of the Congo, on the equator, measured 1.008 times.
The Robinson projection is a pseudocylindrical compromise projection for world maps. It is neither equal-area nor conformal. In Wherobots, it drew eight countries at 0.823 to 1.599 times their true area, and it flattened a 500 km circle at 75 degrees north to a roundness of 0.337. National Geographic replaced it with the Winkel tripel on its world maps in 1995.
The Peters projection, also called Gall-Peters, is a cylindrical equal-area projection with true scale along the 45th parallels. James Gall described a near-identical projection in 1855, and Arno Peters presented his version in 1967. It shows every country at its true relative area and stretches shapes: a 500 km circle at 75 degrees north becomes an ellipse 7.441 times wider than it is tall.
The Goode homolosine projection is an interrupted equal-area world map that J. Paul Goode published in 1925. It joins the sinusoidal projection between about 40 degrees north and south with the Mollweide projection toward the poles. It is interrupted at chosen meridians, with lobes placed to show the land.
An azimuthal projection maps the globe onto a plane that touches it at one point. Directions from that center point to every other point are shown correctly. The azimuthal equidistant version also shows true distances from the center, the Lambert azimuthal equal-area version keeps area, and the stereographic version is conformal.
Tissot’s indicatrix is a way to show distortion that Nicolas Auguste Tissot published in 1859 and 1881. A tiny circle on the Earth becomes an ellipse on the map. On a conformal map the ellipse stays a circle and changes size. On an equal-area map the ellipse keeps the same area and changes shape.
No. John P. Snyder’s USGS manual on map projections states that no map can be both equal-area and conformal. In a Wherobots test, every projection that kept area exactly flattened a circle near 75 degrees north, and Mercator, which kept the circle round, inflated its area 15.515 times.
What is an Isochrone? Isochrone Maps and Travel Time
An isochrone is the area reachable from a place within a set travel time. How isochrones are computed, isochrone vs radius, and a King County test.
What is a Geospatial Database? Definition, Examples
A geospatial database stores geometry, indexes it by location and answers spatial SQL. The definition, examples, and how to choose one by workload.
What is Geospatial Intelligence? GEOINT Defined, Examples
Geospatial intelligence (GEOINT) analyzes imagery and map data to describe activity on earth. Definition, history, examples, and a New Orleans walkthrough.
share this article
Awesome that you’d like to share our articles. Where would you like to share it to: