Lab: Weighted Overlay

A set of tools helping to solve multi-criteria problems like site selection or suitability analysis

A generic problem typically addressed with a Weighted Overlay method is the identification of agricultural land suitable for a particular crop. Below we will be working through a simplified exercise looking e.g. for open, agricultural land with a southern aspect, not overly steep and of course not built-up or covered by forest.

First, though, we want to familiarize ourselves with a sample study area in the Austrian alpine foreland:

Powered by Esri

An overview of the region including the study area - feel free to navigate and 'fly' around to get your bearings!


Data layers for the study area

Our study area is placed the northern edge of the Alps, straddling the Salzburg - Upper Austria border and characterized by a hilly landscape and several lakes.

Powered by Esri

Elevation across the study area has a moderate range from approx 500m a.s.l. to 1100m in the middle of our study area.

Slope is a factor in many physical processes and co-determines the suitability of land for most kinds of uses.

This sample is based on  Terrain: Slope in Degrees  (Living Atlas) and classified into 10° steps.

Land Cover is indispensable for not only understanding land surface dynamics, but also as a starting point for land use decision making.

This sample has been extracted from  World Land Cover 30m BaseVue 2013 , a data product derived from Landsat imagery.

Aspect determines the exposure of land to sunlight and solar radiation, therefore being an essential factor in many land use decisions.

Curvature is a geomorphometric descriptor quantifying local concavity vs convexity.

This and the above morphometric parameters have been resampled to the spatial resolution of the Land Cover layer - keeping in mind that any change in resolution would affect outcomes.


Case study

Equipped with these geospatial factor maps, we can start with finding suitable farming locations for the crop of our choice.

Before starting with an overlay analysis based on the above introduced factor layers, we can explore the study area in some more detail and make mental notes where we could be looking for suitable cropping areas:

Swipe between map and orthophoto, and zoom in to explore areas in detail


Setting up Weighted Overlay in ArcGIS Pro

For reproducing the visuals above in ArcGIS Pro, launch a New Map and Add the layers introduced above: L_WLCBaseVue2013, L_Elev250, L_Slope10, L_Aspect8 and L_Curv3 - and make sure you understand the classification of these layers. This is done easily when connected to https://zgis.maps.arcgis.com as an active portal and selecting these from the group AGI_SpatialAnalysis, or by using the links to the samples quoted above.

For completing the below outlined exercises, though, you would have to access or download and open a > layer package , which also can be found in your connection to https://zgis.maps.arcgis.com as an active portal and selecting this package from the group AGI_SpatialAnalysis.

(Again, you of course are more than welcome to use your own choice of problem and supporting data sets - just make sure these are integer rasters suitably reclassified and with a manageable resolution and extent)

If you have not worked with this toolset before, please read > How Weighted Overlay works , then open the Weighted Overlay tool from Geoprocessing (Spatial Analyst).

Rasters can be added to the analysis with the dropdown arrow, allowing the list of rasters to be populated directly from the table of contents.

As a next step, each (raster) layer is given a relative weight in %, with the total of weights enforced to be 100.

After choosing a rating scale (think 'points', in this example 1-3) each category per layer is rated by the analyst, NODATA categories are entirely excluded and 'Restricted' is given zero points (or more correctly, one point less than the minimum on the chosen rating scale).

Running the weighted overlay generates a raster layer with weighted points averages per cell:

The resulting scenario looks like this, with added labels and legend colors adjusted ('unsuitable' set to transparent in the Symbology tab).

Explore these results and compare with input layer values and rating of categories. Does this correspond with your expectations? Typically, several scenarios would be developed, either to reach a particular target (x hectares of farmland for a particular crop needed), or to test the sensitivity of parametrisation, or to develop several what-if scenarios.


Explore and assess the location of selected areas


Zonal analysis

Using Zonal Overlay, the areas in the 'best' category are then totaled per administrative area, demonstrating the cropping potential per municipality.


Matrix overlay

This kind of results from a Weighted Overlay (sometimes referred to as 'index overlay') will frequently satisfy the requirements and answer the questions and demands for information products put forward by decision makers. Still, we need to consider some of the constraints of this method:

  • Weighted averaging implies a metric scale of rating, a value of 4 means twice the value of 2.
  • High grades in one criterion (layer) can compensate lower grades in another criterion.

Sometimes, these points do not correspond with reality and multi-thematic (multi-layer) combinations of values only can be stated on nominal levels - this is the domain of Matrix Overlay. This name refers to a combinatorial matrix of grades given to tuples of values.

Matrix Overlay can be achieved in different ways, e.g. involving masking and NODATA encoding. Stopping short of scripting some conditional rules in Python, the  Raster Calculator  serves as a handy tool for evaluating non-metric ratings like this (fictitious) example:

The priority matrix for some action to be taken on different forest types in the elevation zones is rated in the scale 1-3

This translates into a clause like: Con(("L_Elev250" == 750) & ("L_WLCBaseVue2013" == 1), 3, Con(("L_Elev250" == 1000) & ("L_WLCBaseVue2013" == 1), 2, Con(("L_Elev250" == 750) & ("L_WLCBaseVue2013" == 2), 2, Con(("L_Elev250" == 1000) & ("L_WLCBaseVue2013" == 2), 1))))

And yes, a matrix of course can have more than 2 dimensions, this would be just harder to visualize and to code ...


Assignment

Think about and define a realistic multi-factor location or suitability problem and translate it into a Weighted Overlay application. Avoid getting started right away with 'doing it', but spend some time on deciding on the how-and-why:

  • which rating scale are you going to employ?
  • what is a suitable minimum mapping unit for the decision you're trying to support?
  • how do you choose and justify the weighting of layers?
  • how do you come with the rating ('points') for categories?
  • is your approach actually reasonably implemented in a 'metric' environment, or should you switch to a nominal domain, i.e. matrix overlay?
  • if you need a certain amount of farmland - e.g. 500ha aligned with the capacity of a processing factory - how could you come up with the best 500ha if more than that are identified as suitable?
  • there might be a minimum operational size of a plot, i.e. you do not really want to farm individual, widely distributed 'pixels' - how do you identify suitable plots with a minimum size? Hint - check out the > Region Group  function.
  • ...

Now it is your turn to set up your own case study, define a problem or question and independently establish a spatial priority or suitability selection in the domain of your choice, thereby creating a valuable information product for decision making.

The priority matrix for some action to be taken on different forest types in the elevation zones is rated in the scale 1-3