Map Projections

The Earth is a spheroid - it is round like a ball or sphere but flattened by slightly by the centrifugal force of rotation. The Earth has a circumference of around 24,900 miles, but is around 27 miles wider than it is tall.

The Earth exists in three-dimensions but, other than globes, most representations of the earth, such as printed maps or web maps, are two dimensional. A projection is a set of mathematical transformations used to represent the three-dimensional world in two dimensions.

This tutorial introduces some fundamental concepts about projections that will be useful when choosing projections for your GIS projects.

Coordinate Systems

The process of representing the three-dimensional world as two-dimensional visualizations requires a sequence of steps.

  1. The lumpy surface of the earth is approximated with a simplified representation called a geoid.
  2. A smooth mathematical ellipsoid or spheroid representation of that geoid that is used as the reference for latitude and longitude degrees. The geoid and ellipsoid together are referred as a geodetic datum, often just called a datum for convenience.
  3. A datum defines a geographic coordinate system of latitudes and longitudes that indicates where locations are on the surface of the planet. Because different datums place their ellipsoids in slightly different places, a single latitude and longitude can refer to different places depending on which datum is used.
  4. Finally, a set of mathematical transformations are performed on the three-dimensional latitudes and longitudes to create a projected coordinate system that specifies the two-dimensional X and Y locations where features are drawn on the paper map or computer display.
Figure
Figure
Relationship between terrain, an datum, and an ellipsoid (NASA 2007, European Environment Agency 2003)
Figure
The relationship between geographic and projected coordinates

The science of creating geoids is called geodesy or geodetics. Geodetic engineers measure three fundamental properties of the earth. Because these properties are constantly changing (albeit in often subtle ways), geodetic engineers always have something to do in support of the earth sciences.

There are dozens of different datums that have been developed over the past century to best fit a datum to specific parts of the earth using the technology of the time. In the United States you will commonly encounter these datums:

GIS software can transform data created with one geographic coordinate system into another, and when you are just visualizing data on a simple map, geographic coordinate systems are generally not a concern. However, there are times you need to pay attention to geographic coordinate systems:

Coordinate System ID Numbers

Projections are often named after the cartographers that developed them. While some common projections have names that are generally understood (like spherical Mercator), a name is usually not specific enough to indicate exactly which projection you want to use.

For example, although the Mercator is a well known general type of projection, different Mercator projections have different parameters depending on how they wrap the world in different ways (transverse, oblique) at different places (central meridian) and use different coordinate units.

Most coordinate systems can be specifically identified with a WKID (well-known identifier), which consists of the name of the authority followed by a four- to six-digit ID number. Some examples:

The authority for a WKID is the organization that decides which projection goes with which ID number. There are two primary authorities:

You can look up WKID numbers at:

Although projections are usually specified with a single WKID for the projected coordinate systems, because projected coordinate systems are built on top of datums with geographic coordinate systems, a projection will have a separate WKID for the geographic coordinate system.

For example, the details for the world Robinson projection listed in ArcGIS Pro (ESRI 54030) show that projection is based on a WGS 1984 geographic coordinate system (EPSG 4326).

Figure
Details on the Robinson Projection in ArcGIS Pro

Common World Projected Coordinate Systems

There are dozens of different types of projections. When you consider that many of those projections have parameters that can be varied depending on the part of the earth you wish to focus on, the number of possible projections is infinite.

However, there are a much more limited number of projections that are commonly used. The following are some North-America-centric suggestions that can give you a good starting point in seeking an appropriate projection for your map. A much longer list of world projections is available here.

World Mercator

Mercator is one of the oldest and most commonly seen projections. It was devised by the Flemish cartographer Gerardus Mercator (1512-1594) during the great age of European exploration, and as such it is conformal (preserves direction) in order to facilitate naval navigation.

A primary disadvantage of the Mercator is that it makes areas look larger the further they are from the Equator. This is best exemplified by the highly exaggerated size of Greenland, something called the Greenland Problem (Block 2004).

Because we visually associate size with power, this characteristic makes countries at northern latitudes (which have also historically been imperial powers) look more powerful than countries in the Global South. In the late 1960s and early 1970s, German historian Arno Peters began actively crusading against use of the Mercator projection (and for the use of his own Gall-Peters projection) as an antidote to what he asserted was the European colonial domination of other peoples of the world reflected in the distortions of the Mercator projection. Peters' critique of the Mercator projection remains active in the cartographic community despite the continued ubiquity of maps using that projection (Peters 1983, Robinson 1985, Monmonier 2010).

Figure
World Mercator (EPSG 3395)

The Mercator projection is a cylindrical projection, with a developable surface equivalent to wrapping the map as a cylinder around the earth.

Figure
Cylindrical Projection

Spherical Mercator (Web Mercator)

Spherical Mercator is a variation of the Mercator projection that mathematically represents the world as a sphere rather than an ellipsoid to make calculations easier. This is the projection used by almost all web maps.

Because the world is slightly flattened, this spherical representation causes the map to be inaccurate when mapping large areas closer to the poles. However, that is not usually a major problem with a web mapping app like Google Maps that is used in lower latitudes for mundane tasks like finding sushi restaurants.

Google Maps Spherical Mercator (EPSG 3857, ESRI 102100)

Robinson Projection

The Robinson projection is a compromise projection developed by Arthur Robinson in 1963 (ESRI 2023, Robinson 1974). The Robinson projection solved the "Greenland Problem" by representing the country closer to its actual size than the Mercator projection. The Robinson projection is considered by many to be visually appealing and in 1988 was adopted by the National Geographic Society for some of its world maps in place of the older Van der Grinten projection.

Figure
The World Robinson Projection (ESRI 54030)

Winkel Tripel

Winkel Tripel is another compromise projection devised by Oswald Winkel in 1921 that strives to represent the whole world while minimizing distortion of area, direction and distance. This was the third of a group of world projections devised by Winkel, and tripel means three in German. This is the projection used for maps by the National Geographic Society and many educational publishers.

Figure
Winkel Tripel Projection (Wikipedia - ESRI 54042)

McArthur's Universal Corrective Map

McArthur's Universal Corrective Map is one of a number of alternative projections intended to provoke viewers to critically examine assumptions about maps. It is a variant on a Mercator projection by Australian cartographer Stuart McArthur that debuted in 1979 with a reversal of north and south.

There is, of course, no geological or astrophysical reason why north should be up, and the reversal of the poles asks us to call into question the biases inherent to a traditional north-up orientation.

Figure
McArthur's Universal Corrective Map

Unprojected Longitude / Latitude

Since longitude and latitude represent x/y coordinates, they can be directly graphed on a grid, and this is the default when you open an unprojected feature class in ArcGIS Pro. However, because longitudes are closer together the further you move away from the equator, the shapes of areas are flattened and expanded the further you move away from the equator.

Equirectangular (unprojected) longitude and latitude is useful for diagnosing problems with data, but should not be used in any professional presentation, as that implies carelessness or ignorance of cartographic principles and standards.

Figure
Unprojected Latitude and Longitude (EPSG 4326)

Common Large Country and Continent Projected Coordinate Systems

Albers Equal Area Conic

Albers Equal Area Conic is a common projection for maps of the continental United States that has historically been used by the United States Geological Society (USGS). It was developed by Heinrich C. Albers in 1805.

Figure
North America Albers Equal Area Conic Projection (ESRI 102008)

This is an example of a conic projection that is equivalent to wrapping the developable surface as a cone around the planet in a way that intersects the surface of the earth in one or two places, which are called standard parallels. This projection is an equal-area projection preserves area while slightly distorting shape. The selection of standard parallels can be based on which part of the earth is the point of focus.

Figure
Conic Projection

Lambert Conformal Conic

The Lambert Conformal Conic is another conic projection commonly used for the continental United States. While it distorts area, it preserves shapes better, and since it preserves angles (conformal) it is commonly used for aeronautical charts since a straight line on the map approximates the shortest-path great-circle distance between two points. It was first developed an published by Swiss cartographer Johann Heinrich Lambert in 1772.

Figure
North America Lambert Conformal Conic Projection (ESRI 102009)

Common Local Projected Coordinate Systems

While projections for world or country maps are primarily focused on presentation or navigation, projections at the local level are more important for the ability to accurately find specific locations and measure distances, such as when surveying property lines or engaging in military activity.

State Plane Coordinate System

The State Plane Coordinate System (SPCS) is a set of planar projections commonly used by local and state governments in the USA.

There are 124 different SPCS zones in different parts of the US.

State Plane Coordinate Systems express their coordinates in different length units as determined by the policies of the different state governments:

There are different SPCS variants for different datums. Unless your organization or project uses an older datum or length unit, the safest option is to use the most recent NAD 83 (2011) datum with lengths in meters.

State Plane Coordinate System Zones (click for zone and WKID)

Universal Transverse Mercator

Universal Transverse Mercator (UTM) is a coordinate system based on the Mercator projection and used as an alternative to latitude and longitude.

UTM divides the world is divided into 120 zones with zone coordinates specified in meters. Zones are numbered 1 - 60, with separate zones (N or S) for the northern and southern hemisphere.

Universal Transverse Mercator Zones (click for zone and WKID)

Transverse Mercator projections are based on a cylindrical wrapping of the developable surface around the planet like the Mercator. However, the wrapping is around the poles (transverse) rather than around the equator (equatorial), so the zones of minimal distortion cover broad swaths of area that can be used as battle spaces.

Figure
Transverse Cylindrical Projection

Map Projections in ArcGIS Pro

Projections can be changed on an ArcGIS Pro map by right-clicking on the layer in the Contents pane and selecting Properties, Coordinate Systems.

The dialog box will show you the current projection and allow you to search or browse for a new projection.

The projections dialog gives you two sets of coordinate systems to choose from:

The video below demonstrates how to change a world map from its existing projection (Web Mercator - EPSG 3857) to the common Robinson projection (EPSG 54030). Because different projections often have similar names, if you want a specific projection and you have a WKID number to search on, that will help you find that specific projection. In this case, we use 54030 for the Robinson projection.

Changing a projection in ArcGIS Pro

Getting Projection Information in ArcGIS Pro

While the software gives you the option to display a map in different geographic coordinate systems, you will normally only create maps using projected coordinate systems. The geographic coordinate system used to create the projected coordinate system is implied when you choose a projected coordinate system.

You can find the geographic coordinate system used in a projection by selecting the Details option in the ArcGIS Pro coordinate systems dialog.

For example, this video shows that the Robinson projected coordinate system (ESRI 54030) uses the WGS 1984 datum (EPSG 4326) as its geographic coordinate system.

Details on the Robinson projection

Grids and Graticules in Map Layout

A common feature of maps is a grid of lines that show map coordinates.

A graticule is a grid of lines that shows intervals of latitude and longitude on a map.

The video below shows how to add a graticule to a map frame.

Displaying a graticule in a map layout

By contrast, in ArcGIS Pro terms, a measured grid is a grid of lines that displays distance intervals in the projected coordinate system used in the map.

Graticules and measured grids add clutter to a map, especially if there is a base map or complex labeling, so they should generally be avoided unless absolutely needed to convey useful information to the map reader.

The video below shows how to add a measured grid to a map layout for a map displayed in the Illinois East State Plane Coordinate System.

Displaying a grid in a map layout

Euclidian Distance with Planar Coordinates

The State Plane Coordinate Systems treat the earth as a flat plane (hence, the name) and can therefore be conveniently used for calculating distance.

The example below uses two points in Chicago: The Sears Tower (41.8789, -87.6358) and the Cloud Gate sculpture in Millennium Park (41.882686, -87.623325). Chicago is in the Illinois East State Plane Coordinate System zone, and the WKID for the coordinate system in US feet is EPSG 3435.

The easting in a State Plane projection is the position east-west (similar to longitude), and the northing is the location north-south (similar to latitude). Because this is a US Feet State Plane variant, those numbers are in feet.

Using ArcGIS Pro as shown in the video, you can change the map projection to State Plane and then get the northing and easting for each point. By subtracting to get the east-west distance and the north-south distance between the points, you get the two sides of a triangle. You can then use the Pythagorean theorem to calculate the Euclidean distance (hypotenuse) between the two points in feet.

Note that you cannot do this with latitudes and longitudes because those are angular measurements, and the ground distance represented by a degree is different on different parts of the planet.

Sears Tower: 1174208 East, 1899233 North
Cloud Gate: 1177593 East, 1900640 North

East/west distance = 1177593 - 1174208 = 3385 feet apart east/west
North/south distance = 1900640 - 1899233 = 1407 feet apart north/south

square_root(3385^2 + 1407^2) = square_root(11458225 + 1979649) = 3666 feet
Figure
Distance calculation using northings and eastings
Getting SPCS eastings and northings in ArcGIS Pro

Projections in Python

By default, the GeoPandas GeoPandas plot() method will display a GeoDataFrame in its current coordinate system. In this case, data read from the GeoJSON file is in unprojected WGS 1984 lat/long (EPSG 4326).

import geopandas
import matplotlib.pyplot as plt

countries = geopandas.read_file("https://michaelminn.net/tutorials/data/2022-world-bank-indicators.geojson")
graticule = geopandas.read_file("https://michaelminn.net/tutorials/data/2023-graticule.geojson")

axis = countries.plot('Energy MM BTU per Capita', legend=True, scheme='quantiles', cmap='coolwarm')
graticule.plot(facecolor='none', edgecolor='#00000020', ax=axis)
axis.set_axis_off()
plt.show()
Figure
GeoDataFrame plot() unprojected WGS 1984

You can transform a GeoDataFrame to a new projection by passing the WKID to the to_crs() method. This example uses the world Robinson projection (ESRI 54030)

countries = countries.to_crs("ESRI:54030")
graticule = graticule.to_crs("ESRI:54030")

axis = countries.plot('Energy MM BTU per Capita', legend=True, scheme='quantiles', cmap='coolwarm')
graticule.plot(facecolor='none', edgecolor='#00000020', ax=axis)
axis.set_axis_off()
plt.show()
Figure
to_crs() to the world Robinson projection

Projections can also be passed as PROJ strings (proj.4). This example demonstrates use of a PROJ string to create a MacArthur's Universal Corrective map.

macarthur = '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +axis=esu +units=m +no_defs'
countries = countries.to_crs(macarthur)
graticule = graticule.to_crs(macarthur)

axis = countries.plot('Energy MM BTU per Capita', legend=True, scheme='quantiles', cmap='coolwarm')
graticule.plot(facecolor='none', edgecolor='#00000020', ax=axis)
axis.set_axis_off()
plt.show()
Figure
to_crs() to the MacArthur's Universal Corrective projection

Projections in R

By default, the sf plot() function will display an sf object in its current coordinate system. In this case, data read from the GeoJSON file is in unprojected WGS 1984 lat/long (EPSG 4326).

library(sf)

countries = st_read("https://michaelminn.net/tutorials/data/2022-world-bank-indicators.geojson")
graticule = st_graticule(lat=seq(-80, 80, 20), lon=seq(-180, 180, 20))

palette = colorRampPalette(c("red", "gray", "navy"))

plot(countries["Energy.MM.BTU.per.Capita"], breaks="quantile", pal=palette, reset=F)
plot(graticule$geometry, col="#00000040", add=T)
Figure
sf plot() unprojected WGS 1984

You can transform an sf object to to a new projection by passing the WKID to the st_transform() function. This example uses the world Robinson projection (ESRI 54030)

countries = st_transform(countries, 'ESRI:54030')
graticule = st_transform(graticule, 'ESRI:54030')

plot(countries["Energy.MM.BTU.per.Capita"], breaks="quantile", pal=palette, reset=F)
plot(graticule$geometry, col="#00000040", add=T)
Figure
st_transform() to the world Robinson projection

Projections can also be passed as PROJ strings (proj.4). This example demonstrates use of a PROJ string to create a MacArthur's Universal Corrective map.

macarthur = '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +axis=esu +units=m +no_defs'
countries = st_transform(countries, macarthur)
graticule = st_transform(graticule, macarthur)

plot(countries["Energy.MM.BTU.per.Capita"], breaks="quantile", pal=palette, reset=F)
plot(graticule$geometry, col="#00000040", add=T)
Figure
st_transform() to the MacArthur's Universal Corrective projection