Adjust the encountered distance using a vertical factor

Available with Image Server

Once the adjusted straight-line distance is calculated, you can use the vertical factor to control the rate at which the distance is encountered. You can use the cost surface, characteristics of the mover, and horizontal factor to control the rate as well.

The vertical factor accounts for the effort to move through the slopes of the landscape. This effort affects how the distance is encountered. Going uphill may take more effort, going downhill requires less effort, and crossing slopes falls somewhere in between. Modifying the adjusted straight-line distance for this effort helps capture the rate at which distance is encountered by the traveler.

A hiker going downhill
Less energy is expended going downhill so the mover can cover the distance at a faster rate.

Do not confuse the vertical factor, which represents the effort to overcome slopes, with the surface distance, which is an adjustment to the straight-line distance for the actual distance the traveler covers as they move over the ups and downs in the landscape.

Slope is often relevant to cost distance analysis. Intuitively, it is more costly to overcome steep slopes than flatter slopes. Typically, the Surface Parameters tool is used to generate a slope raster; however, sometimes that slope raster is incorrectly input into the cost surface.

The Slope option in the Surface Parameters tool calculates the rate of change of elevation for each digital elevation model (DEM) cell. It is the first derivative of a DEM. But, as discussed above, it matters how the slope is encountered when moving through a cell. The traveler may avoid a cell assigned a steep slope in a cost surface. This may be effective if the traveler is moving up the slope when traveling through the cell. However, if the traveler is moving down or across the slope when moving through the cell, the cell may be preferred for travel.

You can account for the effort to overcome the slope using the elevation raster in the vertical factor. Do not include the slope raster in the cost surface when the direction the slope is encountered matters.

The direction that slopes are calculated in the vertical factor can also be altered by the Travel direction source characteristic parameter. That is, moving toward or away from a source will change the direction the traveler will enter a cell and as a result, how the slope is encountered.

Incorporating the vertical factor (VF) is a multiplier modifier to the adjusted straight-line distance calculations. Details about how the vertical factor is calculated are provided in the Distance accumulation algorithm topic.

Vertical factor use examples

The vertical factor can be used in various scenarios, such as the following:

  • Locate a new hiking trail between two campgrounds that is longer but easier to traverse than hiking directly up the shortest path between them.
  • Examine the affects of salt spread on the roads in the winter on the health of the surrounding vegetation. Vegetation downslope from the road will be more impacted by the runoff.
  • Determine the movement of a sea creature that is dependent on the change in the concentration of salinity.

Incorporate a vertical factor

Distance analysis can be divided conceptually into the following related functional areas:

From the second functional area, determine the rate at which the distance will be encountered through a vertical factor as illustrated below. The scenario involves a collection of four forest ranger stations (purple dots), and some rivers (blue lines).

Map of accumulative cost distance from four ranger stations
The least-cost distance from each cell to the closest ranger stations incorporates a barrier, surface raster, and cost surface.

To incorporate the effort for the rangers to overcome the slopes, a vertical factor is specified. The elevation surface is used as the vertical raster.

Map of accumulative cost distance with a vertical factor added
The least-cost distance from each cell to the closest ranger station, incorporating the vertical factor that has a barrier, surface raster, and cost surface. The vertical factor influences the center of the study area in particular.

Create a distance raster using a vertical factor

To create a raster that incorporates a vertical factor, complete the following steps:

  1. Open the Distance Accumulation tool.
  2. Provide a source for the Input raster or feature source data parameter.
  3. Name the output distance accumulation raster.
  4. Expand the Costs relative to vertical movement category.
  5. Provide a vertical factor raster for the Input vertical raster parameter.

    This input is used to calculate the slope that is encountered when moving between cells. Typically, an elevation raster should be specified.

    The Vertical factor parameter appears.

  6. Specify the settings for the Vertical factor parameter.

    This parameter identifies the multiplier to be applied to the cost to adjust for the effort to move through the encountered slopes.

  7. Click Run.

Vertical factor affects the rate distance is encountered

To modify the rate distance that is encountered and account for the effort of the traveler to overcome slopes, the tool does the following:

  • Calculates how the slope is encountered when moving from one cell to the next. This is referred to as the vertical relative moving angle (VRMA).
  • Identifies how the VRMA will modify the rate the distance that is encountered.

Calculate the VRMA

The VRMA is the angle of slope from the processing cell (the From cell) to the cell that the traveler is moving into (the To cell). The distance is being calculated for the To cell. The heights from which the slopes are calculated are defined by the input vertical factor raster.

The slope is calculated using the Pythagorean theorem formula of rise/run. The base of the triangle necessary for determining the slope is derived from the adjusted straight-line distance. The height is established by subtracting the From cell value from the To cell value. The resultant angle is the VRMA.

VRMA calculation

The VRMA is specified in degrees. The range of values for the VRMA is -90 to +90 degrees, accounting for both positive and negative slopes.

Identify the VF multiplier

The VRMA value is then plotted on the specified vertical factor graph to obtain the vertical factor multiplier that will be used in the calculations that determine the cost to reach the To cell. The distance value to move through the cell is multiplied by the identified vertical factor. The larger the vertical factor, the more difficult the movement. A VF greater than 1 increases the cost distance encountered. A VF less than 1 but greater than 0 allows the mover to encounter the distances at a faster rate.

For example, the following graph shows the relationship of the VF and the VRMA for a linear VF function:

VF and VRMA in a linear type graph

The vertical factor functions that allow you to capture the interaction of the traveler to the slopes they are encountering are Binary, Linear, Inverse linear, Symmetric linear, Symmetric inverse linear, Cos, Sec, Cos-Sec, and Sec-Cos. See the Additional information section below for details of each function.

Note:

The vertical factor is a multiplier. Take care when specifying the units when combining the vertical factor with a cost surface, source characteristics, or a horizontal factor. Generally, when a cost surface is entered, the vertical factor should be a multiplier adjustment of the rate of cost surface units. If time is the unit for the rate of the cost surface, the vertical factor should be a modifier of time. Only one of these factors can define the units for the rate. The other factors are unitless, and their values are multiplier modifiers of the specified units.

Sample applications that use a vertical factor

Sample applications that use a vertical factor are described below.

Create a downslope buffer to understand the effect of winter salt on vegetation

You want to identify areas that are downslope and within 50 meters of a road because these areas may be impacted by winter saltwater runoff. You want to measure the distance along the surface of the terrain. You can use the Binary vertical factor setting to prevent the Distance Accumulation tool from identifying cells that are higher than road cells. Examples of the resulting downslope buffers are shown below.

Road with a 50 meter downslope buffer around it
The areas that can be impacted by winter salt runoff below roads are identified. Using the Distance Accumulation tool with the binary vertical factor setting, the downslope cells within a surface distance of 50 meters from roads are determined (orange shading).

For comparison, a different section of the road is used to show the difference in the straight-line and downslope-only buffers.

Road with a 50 meter downslope buffer overlaid on a 50 meter straight-line buffer
A straight-line distance 50 meter buffer (blue shading) overlaid on a downslope-only surface distance buffer (orange shading). The latter is narrower in locations where you cannot travel downhill for 50 meters.

To create a downslope buffer, complete the following steps:

  1. Open the Distance Accumulation tool.
  2. Provide the roads in the Input raster or feature source data parameter.
  3. Name the Output distance accumulation raster parameter value.
  4. Expand the Costs relative to vertical movement category.
  5. Provide an elevation raster in the Input vertical raster parameter.
  6. Set Binary for the Vertical factor parameter.
  7. Expand the Characteristics of the source category.
  8. Set the Maximum accumulation distance parameter to 50 meters.
  9. Click Run.

Tobler’s hiking function

You want to calculate hiking time over the landscape, while adjusting the walking speed based on the slope that is encountered in the direction of travel. Tobler’s hiking function is an empirical model that is used to perform that adjustment. In the model, it is assumed that a base walking speed is 6 km/h, which is achieved when traveling slightly downslope (about -3 degrees).

Tobler’s hiking function

Where d is the angle that slope is calculated from.

When graphed, the speed function W looks like the following:

Graph of Tobler's speed function
Tobler's speed function as a function of slope in degrees is expressed in km/h. The maximum travel speed of 6 km/h occurs when walking slightly downslope.

You want to know how long it takes to travel a given distance (one cell), not how far you can travel in a specified amount of time, so you need to work with the reciprocal of speed, called pace. The pace in hours per meter (because the horizontal distance analysis units are meters) is expressed instead of hours per kilometer:

Pace function

Where d is the angle that slope is calculated from.

The pace function, called hiking time, looks like the following:

Graph of Tobler's speed function converted to pace function
The reciprocal of Tobler's speed function is the pace function, expressed here in hours per meter, plotted as a function of slope in degrees.

To create a hiking time map, complete the following steps:

  1. Open the Distance Accumulation tool.
  2. Provide the hiking starting location as the Input raster or feature source data parameter value.
  3. Name the Output distance accumulation raster parameter value.
  4. Provide an elevation surface with the Input surface raster parameter.
  5. Expand the Costs relative to vertical movement parameter category, and do the following
    1. Provide an elevation raster in the Input vertical raster parameter.
    2. Set the Vertical factor parameter value to Hiking time.
    3. Set the Low cut angle and High cut angle parameters as the default values.
  6. Expand the Characteristics of the source category, and set the Travel Direction parameter value to Travel from source.
  7. Click Run.

On the resulting accumulation raster you can click anywhere to find the hiking time, in hours, from the source to the location you clicked, while taking the most optimal route. Including the direction of travel is important, since hiking up a gradient is done at a slower speed than hiking down the same gradient.

An additional option is available for the Vertical factor parameter that is also based on Tobler's hiking function, the Bidirectional Hiking Time option. This setting finds the average time to hike in one direction on an out-and-back hiking trail, using the most optimal route. If you want to find the most optimal out-and-back trail this is the best way to determine the optimal route. When you use this option, the Travel Direction parameter does not have any affect on the result, because it is an average of both directions.

Additional information

The following sections provide additional information regarding vertical factors.

Vertical factors

To define the vertical factor function, a list of graphs is provided or you can create a custom function using an ASCII file. The following vertical factor functions are available in the Distance Accumulation tool:

Vertical factor options, modifiers, and default values

FunctionZero factorLow cut angleHigh cut angleSlopePowerCos powerSec power
Binary

1

-30

30

N/A

N/A

N/A

N/A

Linear

1

-90

90

1.111E-02

N/A

N/A

N/A

Inverse linear

1

-45

45

-2.222E-02

N/A

N/A

N/A

Symmetric linear

1

-90

90

1.111E-02

N/A

N/A

N/A

Symmetric inverse linear

1

-45

45

-2.222E-02

N/A

N/A

N/A

Cos

N/A

-90

90

N/A

1

N/A

N/A

Sec

N/A

-90

90

N/A

1

N/A

N/A

Cos – Sec

N/A

-90

90

N/A

N/A

1

1

Sec – Cos

N/A

-90

90

N/A

N/A

1

1

Hiking Time

N/A

-70

70

N/A

N/A

N/A

N/A

Bidirectional Hiking Time

N/A

-70

70

N/A

N/A

N/A

N/A

Binary

When the VRMA is larger than the low-cutting angle and smaller than the high-cutting angle, the VF for moving between the two cells is set to the value associated with the zero factor. If the VRMA is greater than the cutting angle, the VF is set to infinity. The default cutting angle is 30 degrees if none is specified.

Default Binary vertical factor graph

Linear

The VFs are determined by a straight line in the VRMA-VF coordinate system. The line intercepts the y-axis, equitable to the VF factor at the zero factor value. The slope of the line can be specified using the Slope modifier. If no slope is identified, the default is 1/90 (specified as 0.01111). The default low-cutting angle is -90 degrees, and the default high-cutting angle is 90 degrees.

Default Linear vertical factor graph

Inverse linear

The VFs are determined by the inverse values from a straight line in the VRMA-VF coordinate system. The line intercepts the y-axis, equitable to the VF factor at the zero factor value. The slope of the line can be specified using the Slope modifier. If no slope is identified, the default is -1/45 (specified as 0.02222). The default low-cutting angle is -45 degrees, and the default high-cutting angle is 45.

Default Inverse linear vertical factor graph

Symmetric linear

This vertical factor is composed of two linear functions relative to the VRMAs that are symmetrical to the VF (y) axis. Both lines intercept the y-axis at the VF value associated with the zero factor. The slope of the lines is defined as a single slope relative to the positive VRMA using the Slope vertical factor modifier, which mirrors the negative VRMAs. The default slope is 1/90 (specified as 0.01111). The default low-cutting angle is -90, and the default high-cutting angle is 90.

Default Symmetric linear vertical factor graph

Symmetric inverse linear

This vertical factor is the inverse of the Symmetric linear vertical factor keyword. It is composed of two inverse linear functions relative to the VRMAs that are symmetrical to the VF (y) axis. Both lines intercept the y-axis at the VF value of 1. The slope of the lines is defined as a single slope relative to the positive VRMA using the Slope vertical factor modifier, which mirrors the negative VRMAs. The default slope is -1/45 (specified as 0.02222). The default low-cutting angle is -45, and the default high-cutting angle is 45.

Default Symmetric inverse linear vertical factor graph

Cos

The VF is determined by the cosine function of the VRMA. The default low-cutting angle is -90 degrees, and the default high-cutting angle is 90 degrees. The default Cos power value is 1.0.

Default Cos vertical factor graph

Sec

The VF is determined by the secant function of the VRMA. The default low-cutting angle is -90 degrees, and the default high-cutting angle is 90 degrees. The default Sec power value is 1.0.

Default Sec vertical factor graph

Cos-Sec

When the VRMA is a negative degree value, the VF is determined by the cosine function of the VRMA. If the VRMA is a positive degree value, the VF is determined by the secant function of the VRMA. The default low-cutting angle is -90 degrees, and the default high-cutting angle is 90 degrees. The default Cos power and Sec power values are both 1.0.

Default Cos-Sec vertical factor graph

Sec-Cos

When the VRMA is a negative degree value, the VF is determined by the secant function of the VRMA. If the VRMA is a positive degree value, the VF is determined by the cosine function of the VRMA. The default low-cutting angle is -90 degrees, and the default high-cutting angle is 90 degrees. The default Sec power and Cos power values are both 1.0.

Default Sec-Cos vertical factor graph

Hiking time

The VF is the reciprocal of Tobler's hiking function and the result is the pace in hours. The lowest cost is when the VRMA is about -3 degrees. The default low-cutting angle is -70 degrees, and the default high-cutting angle is 70 degrees. The formula will be adjusted for map units.

The formula for Tobler's hiking function is as follows:

Tobler’s hiking function

The hiking time function, which is the reciprocal of Tobler's hiking function, is as follows:

Pace or hiking time function

Where d is the angle that slope is calculated from.

Hiking time vertical factor graph

Bidirectional hiking time

This function determines the average time to hike one direction on an out-and-back hiking trail, using the most optimal route. The VF is based on the reciprocal of Tobler's hiking function, but it uses an average and the result is the pace in hours. The formula will be adjusted for map units.

The formula for Tobler's hiking function is as follows:

Tobler’s hiking function

The hiking time function, which is the reciprocal of Tobler's hiking function is as follows:

Pace or hiking time function

The Bidirectional hiking function time is the average of the hiking time using the positive and negative angles.

Bidirectional hiking function

Where d is angle that slope is calculated from.

Bidirectional hiking time vertical factor graph

Table

The table is an ASCII file with two columns on each line.

The first column identifies the VRMA in degrees, and the second column identifies the VF. Each line specifies a point. Two consecutive points produce a line segment in the VRMA-VF coordinate system. The angles must be input in ascending order and range from -90 to 90. The VF factor for any VRMA angle less than the first (lowest) input value or larger than the final (largest) input will be set to infinity. An infinite VF is represented by -1 in the ASCII table.

The following is a sample vertical factor ASCII table. The units of the first column are degrees and the units of the second column are hours per meter.

    -90  -1
    -80  -1
    -70   2.099409721
    -60   0.060064462
    -50   0.009064613
    -40   0.00263818
    -30   0.001055449
    -20   0.000500142
    -10   0.00025934
      0   0.000198541
     10   0.000368021
     20   0.000709735
     30   0.001497754
     40   0.003743755
     50   0.012863298
     60   0.085235529
     70   2.979204206
     80  -1
     90  -1

Vertical factor modifiers

You can further control the VRMA function using modifiers that allow refinement of the vertical factors. There may be a threshold angle such that if the VRMA exceeds it, the cost is so great that it becomes a barrier to travel. This threshold is referred to as the cut angle. The VF is assigned to infinity when the VRMA exceeds this value.

The vertical factor graph will have both lower and upper cut angles in contrast to the horizontal factor graph, which will only have a single cut angle.

Using these modifiers, cut angles can be specified for each of the functions, the trigonometric curves can be raised by a power, the zero factor can alter the y-axis intercept for the nontrigonometric functions, and the slope of the line in the linear functions can be defined.

Zero factor

This modifier establishes the vertical factor used when the VRMA is zero. This factor positions the y-intercept of the specified function.

Low cut angle

This modifier is the VRMA degree that defines the lower threshold, below which (less than) the VFs are set to infinity regardless of the specified vertical factor keywords.

High cut angle

This modifier is the VRMA degree that defines the upper threshold, beyond which (larger than) the VFs are set to infinity regardless of the specified vertical factor keywords.

Low and High cut angle vertical factor modifiers example

Slope

This modifier identifies the slope of the straight lines in the VRMA-VF coordinate system for the Linear, Inverse linear, Symmetric linear, and Symmetric inverse linear keywords. Slope is specified as the rise over the run (for example, a 30-degree slope is 1/30 specified as 0.03333). See the Linear VRMA diagram for an example of a linear function with a slope of 1/90.

Power

This modifier is the power to which the values will be raised.

Cos power

This modifier is the power to which the nonnegative values in the Sec-Cos VRMA function and the negative values in the Cos-Sec VRMA function will be raised. The VF is determined by the following:

VF = cos(VRMA)power

Sec power

This modifier is the power to which the nonnegative values in the Cos-Sec VRMA function and the negative values in the Sec-Cos VRMA function will be raised. The VF is determined by the following:

VF = sec(VRMA)power

Table name

This modifier identifies the name of the ASCII file that will be used with the Table vertical factor keyword.

Additional resources

Tobler, Waldo. 1993. Three Presentations on Geographical Analysis and Modeling: Non-Isotropic Geographic Modeling; Speculations on the Geometry of Geography; and Global Spatial Analysis (93-1). Retrieved from https://escholarship.org/uc/item/05r820mz