Map algebra is analgebra for manipulatinggeographic data, primarily fields. Developed by Dr.Dana Tomlin and others in the late 1970s, it is a set of primitive operations in ageographic information system (GIS) which allows one or morerasterlayers ("maps") of similar dimensions to produce a new raster layer (map) using mathematical or other operations such as addition, subtraction etc.

Prior to the advent of GIS, the overlay principle had developed as a method of literally superimposing different thematic maps (typically anisarithmic map or achorochromatic map) drawn on transparent film (e.g.,cellulose acetate) to see the interactions and find locations with specific combinations of characteristics.^[1] The technique was largely developed by landscape architects and city planners, starting withWarren Manning and further refined and popularized byJaqueline Tyrwhitt,Ian McHarg and others during the 1950s and 1960s.^[2][3][4]

In the mid-1970s, landscape architecture student C. Dana Tomlin developed some of the first tools for overlay analysis in raster as part of the IMGRID project at theHarvard Laboratory for Computer Graphics and Spatial Analysis, which he eventually transformed into the Map Analysis Package (MAP), a popular raster GIS during the 1980s. While a graduate student atYale University, Tomlin and Joseph K. Berry re-conceptualized these tools as a mathematical model, which by 1983 they were calling "map algebra."^[5][6] This effort was part of Tomlin's development ofcartographic modeling, a technique for using these raster operations to implement the manual overlay procedures of McHarg. Although the basic operations were defined in his 1983 PhD dissertation, Tomlin had refined the principles of map algebra and cartographic modeling into their current form by 1990.^[7][8] Although the term cartographic modeling has not gained as wide an acceptance as synonyms such assuitability analysis, suitability modeling and multi-criteria decision making, "map algebra" became a core part of GIS. Because Tomlin released the source code to MAP, its algorithms were implemented (with varying degrees of modification) as the analysis toolkit of almost every raster GIS software package starting in the 1980s, including GRASS, IDRISI (nowTerrSet), and the GRID module of ARC/INFO (later incorporated into the Spatial Analyst module of ArcGIS).

This widespread implementation further led to the development of many extensions to map algebra, following efforts to extend the raster data model, such as adding new functionality for analyzing spatiotemporal and three-dimensional grids.^[9][10]

Like otheralgebraic structures, map algebra consists of a set of objects (thedomain) and a set ofoperations that manipulate those objects withclosure (i.e., the result of an operation is itself in the domain, not something completely different). In this case, the domain is the set of all possible "maps," which are generally implemented as raster grids. A raster grid is a two-dimensional array ofcells (Tomlin called themlocations orpoints), each cell occupying a square area of geographic space and being coded with a value representing the measured property of a given geographic phenomenon (usually a field) at that location. Each operation 1) takes one or more raster grids as inputs, 2) creates an output grid with matching cell geometry, 3) scans through each cell of the input grid (or spatially matching cells of multiple inputs), 4) performs the operation on the cell value(s), and writes the result to the corresponding cell in the output grid.^[7] Originally, the inputs and the output grids were required to have the identical cell geometry (i.e., covering the same spatial extent with the same cell arrangement, so that each cell corresponds between inputs and outputs), but many modern GIS implementations do not require this, performing interpolation as needed to derive values at corresponding locations.^[11]

Tomlin classified the many possible map algebra operations into three types, to which some systems add a fourth:^[12]

Local Operators

Operations that operate on one cell location at a time during the scan phase. A simple example would be an arithmetic operator such as addition: to compute MAP3 = MAP1 + MAP2, the software scans through each matching cell of the input grids, adds the numeric values in each using normal arithmetic, and puts the result in the matching cell of the output grid. Due to this decomposition of operations on maps into operations on individual cell values, any operation that can be performed on numbers (e.g., arithmetic, statistics, trigonometry, logic) can be performed in map algebra. For example, a LocalMean operator would take in two or more grids and compute thearithmetic mean of each set of spatially corresponding cells. In addition, a range of GIS-specific operations has been defined, such as reclassifying a large range of values to a smaller range of values (e.g., 45 land cover categories to 3 levels of habitat suitability), which dates to the original IMGRID implementation of 1975.^[13] A common use of local functions is for implementingmathematical models, such as an index, that are designed to compute a resultant value at a location from a set of input variables.

Focal Operators

Functions that operate on a geometricneighborhood around each cell. A common example is calculating slope from a grid of elevation values. Looking at a single cell, with a single elevation, it is impossible to judge a trend such as slope. Thus, the slope of each cell is computed from the value of the corresponding cell in the input elevation grid and the values of its immediate neighbors. Other functions allow for the size and shape of the neighborhood (e.g. a circle or square of arbitrary size) to be specified. For example, a FocalMean operator could be used to compute the mean value of all the cells within 1000 meters (a circle) of each cell.

Zonal Operators

Functions that operate onregions of identical value. These are commonly used with discrete fields (also known as categorical coverages), where space is partitioned into regions of homogeneous nominal or categorical value of a property such asland cover,land use, soil type, or surface geologic formation. Unlike local and focal operators, zonal operators do not operate on each cell individually; instead, all of the cells of a given value are taken as input to a single computation, with identical output being written to all of the corresponding cells. For example, a ZonalMean operator would take in two layers, one with values representing the regions (e.g., dominant vegetation species) and another of a related quantitative property (e.g., percent canopy cover). For each unique value found in the former grid, the software collects all of the corresponding cells in the latter grid, computes thearithmetic mean, and writes this value to all of the corresponding cells in the output grid.

Global Operators

Functions that summarize the entire grid. These were not included in Tomlin's work, and are not technically part of map algebra, because the result of the operation is not a raster grid (i.e., it is notclosed), but a single value or summary table. However, they are useful to include in the general toolkit of operations. For example, a GlobalMean operator would compute thearithmetic mean of all of the cells in the input grid and return a single mean value. Some also consider operators that generate a new grid by evaluating patterns across the entire input grid as global, which could be considered part of the algebra. An example of these are the operators for evaluating cost distance.^[14]

Several GIS software packages implement map algebra concepts, includingPostGIS,ERDAS Imagine,QGIS,GRASS GIS,TerrSet,PCRaster, andArcGIS.

In Tomlin's original formulation of cartographic modeling in the Map Analysis Package, he designed a simple procedural language around the algebra operators to allow them to be combined into a complete procedure with additional structures such as conditional branching and looping.^[8] However, in most modern implementations, map algebra operations are typically one component of a general procedural processing system, such as a visual modeling tool or a scripting language. For example, ArcGIS implements Map Algebra in both its visual ModelBuilder tool and inPython. Here, Python's overloading capability^[15] allows simple operators and functions to be used for raster grids. For example, rasters can be multiplied using the same "*" arithmetic operator used for multiplying numbers.^[16]

Here a modern MapAlgebra implementation, embedding map algebra expressions into SQL (of PostGIS and others), seefunctionST_MapAlgebra() guide:

WITHfooAS(SELECTST_AddBand(ST_MakeEmptyRaster(10,10,0,0,1,1,0,0,0),'32BF'::text,1,-1)ASrast)SELECTST_MapAlgebra(rast,1,NULL,'ceil([rast]*[rast.x]/[rast.y]+[rast.val])')FROMfoo;

Here are some examples inMapBasic, the scripting language forMapInfo Professional:

# demo for Brown's Pond data set# Give layers#  altitude#  development – 0: vacant, 1: major, 2: minor, 3: houses, 4: buildings, 5 cement#  water – 0: dry, 2: wet, 3: pond# calculate the slope at each location based on altitudeslope = IncrementalGradient of altitude# identify the areas that are too steeptoosteep = LocalRating of slope  where 1 replaces 4 5 6  where VOID replaces ...# create layer unifying water and developmentoccupied = LocalRating of development  where water replaces VOIDnotbad = LocalRating of occupied and toosteep  where 1 replaces VOID and VOID  where VOID replaces ... and ...roads = LocalRating of development  where 1 replaces 1 2  where VOID replaces ...nearread = FocalNeighbor of roads at 0 ... 10aspect = IncrementalAspect of altitudesouthface = LocalRating of aspect  where 1 replaces 135 ... 225  where VOID replaces ...sites = LocalMinimum of nearroad and southface and notbadsitenums = FocalInsularity of sites at 0 ... 1sitesize = ZonalSum of 1 within sitenumsbestsites = LocalRating of sitesize  where sitesize replaces 100 ... 300  where VOID replaces ...

• Mathematical morphology
• Field (geography)

• osGeo-RFC-39 about Layer Algebra

• B. E. DavisGIS: A Visual Approach (2001 Cengage Learning) pp. 249ff.

• Geographic information systems
• Applied mathematics
• Algebra
• Geographic information science
• Spatial analysis

• CS1 maint: location missing publisher
• Articles with short description
• Short description matches Wikidata