1. External links

2. References

Map algebra is a set-based algebra for manipulating geographic data, proposed by Dr. Dana Tomlin in the early 1980s. It is a set of primitive operations in a geographic information system (GIS) which allows two or more raster layers ("maps") of similar dimensions to produce a new raster layer (map) using algebraic operations such as addition, subtraction etc.

Depending on the spatial neighborhood, GIS transformations are categorized into four classes: local, focal, global, and zonal. Local operations works on individual raster cells, or pixels. Focal operations work on cells and their neighbors, whereas global operations work on the entire layer. Finally, zonal operations work on areas of cells that share the same value. The input and output for each operator being map, the operators can be combined into a procedure or script, to perform complex tasks.^[1]

When map algebra is performed in cells from local operations, different types of operations can be used:

• Arithmetic operations uses basic mathematical functions like addition, subtraction, multiplication and division.
• Statistical operations uses statistical operations such as minimum, maximum, average and median.
• Relational operations compares cells using functions such as greater than, smaller than or equal to.
• Trigonometric operations uses sine, cosine, tangent, arcsine between two or more raster layers.
• Exponential and logarithmic operations use exponent and logarithm functions.^[2]

Several major GIS systems use map algebra concepts, including ERDAS Imagine and ArcGIS. ArcGIS 10 implements Map Algebra in Python; functions are imported Python methods and Python's overloading capability^[3] is used for operators. For example, rasters can be multiplied using the "*" arithmetic operator.^[4]

Here are some examples, in MapBasic:

1. demo for Brown's Pond data set
2. Give layers
3. altitude
4. development – 0: vacant, 1: major, 2: minor, 3: houses, 4: buildings, 5 cement
5. water – 0: dry, 2: wet, 3: pond

1. calculate the slope at each location based on altitude

slope = IncrementalGradient of altitude

1. identify the areas that are too steep

toosteep = LocalRating of slope

  where 1 replaces 4 5 6  where VOID replaces ...

1. create layer unifying water and development

occupied = LocalRating of development

notbad = 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 ... 10

aspect = IncrementalAspect of altitude

southface = LocalRating of aspect

  where 1 replaces 135 ... 225  where VOID replaces ...

sites = LocalMinimum of nearroad and southface and notbad

sitenums = FocalInsularity of sites at 0 ... 1

sitesize = ZonalSum of 1 within sitenums

bestsites = LocalRating of sitesize

  where sitesize replaces 100 ... 300  where VOID replaces ...

External links

• osGeo-RFC-39 about Layer Algebra

References

• Half Brother
• Half Brother (album)
• Half Brother (band)
• Half (Canadian coin)
• Half-century (cricket)
• Half-chord
• Half-circle
• Halfcircle
• Half circle
• Half-circumferentor
• Half-cloverleaf interchange
• Half-collared gnatwren
• Half-colon
• Half colon
• Half cousin
• Half-cousin
• Halfdan
• Halfdan Haneborg Hansen
• Halfdan Jebe
• Halfdan Jonsson
• Half-demon
• Half-dimininished seventh chord
• Half Dome (Aye Nako song)
• Half Dome Granite
• Half Dome Granodiorite