Press 'o' to toggle the slide overview and 'f' for full-screen mode.

Choose the theme in which to view this presentation:

Black - White - League - Sky - Beige - Simple
Serif - Blood - Night - Moon - Solarized

GEOG*3480

GIS and Spatial Analysis


Least-cost Pathway Analysis



John Lindsay
Fall 2015

Background readings

  • Your textbook doesn't discuss LCPA.

  • For background materials, please see ESRI's online help topic "Understanding cost distance analysis" (http://pro.arcgis.com/en/pro-app/tool-reference/spatial-analyst/understanding-cost-distance-analysis.htm)

Distance and Cost

  • We're not always interested in the straight-line distance between spatial objects

  • There are times (and types of spatial problems) that require the estimation of cost-weighted distances, i.e. calculating the 'cost' of moving across the landscape

Distance and Cost

  • Cost = money, time, energy, fuel, suffering, resistance or other parameters

  • Very often we want to minimize costs

  • Euclidean distance has the limitation of assuming that the cost of crossing distances in any cell on the way to a source cell is the same as the cost of crossing distances in any other cell

Least-cost pathway analysis (LCPA)

  • LCPA is a way of finding an optimum, or set of optimum, pathways across a landscape with a non-uniform spatial distribution of resistance to travel
  • Many GIS offer tools for performing LCPA because it is so widely used...it's a common spatial problem
  • Measures the effective distance rather than the Euclidean distance
  • Type of spatial analysis based on the raster model
ArcGIS LCPA
ArcGIS LCPA

LCPA Steps

  1. Calculate the cost of traversing through each grid cell...the Cost Surface
  2. Specify a Source Cell(s) and/or Target Cell(s)
  3. Perform Cost Accumulation analysis
  4. Identify the optimum Least-Cost Pathway(s)

The Cost Surface

  • We can create a cost surface to account for the factors affecting cost

  • Each grid cell contains a number that indicates the cost of going through a unit of distance in that cell

  • If it is three times harder to go through the forest as it is to go through the pasture, we might assign values of 1 to the pasture cells and 3 to the forest cells.

The Cost Surface

Physical barriers (e.g. streams) can be represented in the cost surface by exceedingly high costs. Consider the following maze example:

maze cost surface

The Cost Surface

Care must be taken when adding linear features (roads, streams, etc.) to the cost surface. Diagonal raster connectivity may allow the least-cost path to incorrectly pass through an area with high resistance.

line barriers

Cost Accumulation

The accumulated cost found in the output grid is the sum of the products formed by:

  1. the distance across each cell, and;
  2. the value of each cell, i.e. the price per unit distance for going through, which will depend on the grid resolution.

Cost Accumulation

  • Cost accumulation calculates the SPREAD distance, with a friction (cost) determined by the cost surface.
  • The accumulation procedure works outward from specified source cells.
  • Each cell in the grid is then assigned the minimum cost-distance connecting it to the ‘nearest’ source cell.
  • Effectively creates a continuous cost accumulation surface.

Cost Accumulation

maze cost accumulation

The Least-Cost Pathway

  • Least cost pathways can be estimated by using a D8 flow routing algorithm on the cost accumulation surface in much the same way that we've seen flow algorithms used to trace flowpaths over the surface defined by a DEM
  • The source cells will be the lowest points in the cost accumulation surface
  • You can trace the pathway of steepest descent from any target cell(s) to the nearest source cell, working backwards

The Least-Cost Pathway

maze least-cost pathway

Let's consider an example...

  • Two major highways travel parallel through a mountainous environment. When there are traffic accidents, large sections of the highway must be closed. Unfortunately because no secondary road connects either highway during the dangerous section of the highway, detours are impossible during times of closure. This has become a serious problem of the area and as a result the authorities are considering creating a new secondary road to join the main highways.
  • What would be the best route for this new highway?
example map

Things to consider that affect the cost

  • The raster grid resolution is 90 m

  • [note all costs are made up for the example]

  • Privately owned land is prohibitively expensive to acquire (set cost to $100,000 / 90 m, i.e. an effective barrier)

  • Bridges over lakes are not feasible (within our budget)

Things to consider that affect the cost

  • There is a cost associated with crossing rivers and rivers.
    • Small rivers can be culverted ($500 / 90 m)
    • Small bridges cost between $5000 / 90 m and $10,000 / 90 m
    • Large bridges are prohibitively expensive (set to $1,000,000 / 90 m)
  • We can use the Stream Order tool to distinguish small, medium and large streams.
  • Streams, being single cell wide lines, need to be represented in the cost surface as 'fat' networks and we can use a 3 x 3 Maximum filter to do this.

Things to consider that affect the cost

  • We cannot build a road through wetlands due to environmental reasons (set to $1,000,000 / 90 m)
  • It costs more to build roads on steeper slopes because of the need to cut bedrock to grade the slope
    • 0°-5° $1,000 / 90 m
    • 5°-10° $2,000 / 90 m
    • 10°-15° $3,000 / 90 m
    • >15° $100,000 / 90 m (major grading required)
  • Slope can be measured from the DEM using the Slope tool.

Once we've combined each cost factor into a single cost surface, we can use Whitebox GAT's Cost Accumulation tool.

cost accumulation tool

Here is the cost surface and cost accumulation surface for our example:

cost surfaces

Now we use the back-link grid output from the cost accumulation as an input for Whitebox's Least-Cost Pathway tool. For our example, we would like the secondary road to connect to Road 2 (i.e. the eastern road) at Old Town. A raster of Old Town's location is therefore used as the destination input.

LCP tool

Here is the back-link raster and final least-cost pathway connecting Old Town on Road 2 to Road 1:

backlink and LCP