Slope Algorithms

Znw

Zn

Zne

Zw

Z

Ze

Zsw

Zs

Zse

The slope is taken as a dZ value divided by a horizontal distance.  The horizontal distance is the data spacing (east-west, north-south, or diagonal) for the one point and nine point methods, and twice that distance for the four and eight point methods.  For DEMs like the USGS 10 m and 30, the NS and EW spacings are the same; they differ for data like the USGS NED, USGS 1:250K DEM, SRTM, and DTED with geographic spacing.

Aspect, the downhill direction, can be computed as a byproduct of the slope computation.

Slope and aspects, the magnitude and direction of the vector tangent to the topographic surface pointing downhill at a point, have been computed with a multitude of methods; Eyton (1991) and Carter (1992) list a number of references. The following discussion centers on the methods implemented in MICRODEM; the references point to clear expositions of the method, and not necessarily the first use of the method, many of which are in hard to location publications in the gray literature. The algorithms fall into several categories, depending on the number of neighboring points considered to compute the slope and aspect. 

In general the choice of slope algorithm does not make much difference, and you should just stick with the default.


Two neighbors:

Three neighbors:

Four neighbors to N, S, E, and W (excluding point itself) (rook's case):

Four neighbors to NW, SE,N E, and SW (excluding point itself) (bishop's case):

Eight neighbors (excluding point itself) (Queen's case):

Nine points:


Table. Selected Slope and Aspect Algorithms expanded and adapted from Guth (1995).

Name Short Name Neighbors Used Neighbors Comments References
Steepest Down Hill SDN 9   Also called D8; only 8 aspects O'Callahan and Mark, 1987
Steepest Adjacent Neighbor SAN 9   Only 8 aspects Sharpnack & Akin, 1969
Average Neighbor AVN 9   Only 8 aspects  
Guth Hybrid (Steepest + 8 even for aspect) HYB 9   Preferred in Guth, 1995 Guth, 1995
Four Neighbors N4 4 zn,zs,ze,zw "Rook's case": Preferred in Jones, 1998

Second order finite difference

Fleming and Hoffer, 1979; Ritter, 1987; Zevenbergen and Thorne, 1987; O'Neill & Mark (1987)
3 Neighbors N3 3 z,zn,zne In a grid, two choices for the triangles O'Neill and Mark, 1987
8 neighbors, even weighting N8E 8   "Queen's case": Equivalent to fitting a second order trend surface. Preferred by Evans, 1998; Florinsky, 1998.

Identical weights. Third order finite difference

Sharpnack and Akin, 1969; Evans, 1980; Horn, 1981; Heerdegen and Beran, 1982;  Wood, 1996
8 neighbors weighted N82 8   Third order finite difference weighted by reciprocal of twice distance Differential weights (Horn, 1981); Sobel operator (Richards, 1986)
8 neighbors weighted N8R 8   Third order finite difference weighted by reciprocal of sqrt2 * distance Unwin, 1981
Frame Finite Difference FFD 4   Difference above, below, right, and left  
Simple Difference SD 3     The other choice of points from O'Neill and Mark, 1987
Two neighbors   2 z,ze No aspect Batson & others (1975)

Maps showing where slope algorithms differ.

Assess the effect of region size on computed slope.

Slope references

Slope units, in percent.

Slope algorithms for arc second DEMs


Three slope methods appear to have the broadest support:

The Steepest Adjacent Neighbor (SAN): requires adjustments to get reliable aspect distribution, and does not give a smooth slope distribution. It was the preferred method in MICRODEM until the summer of 2003, and there are still reasons to prefer it.  The Eight Neighbors Unweighted method became the default in MICRODEM the summer of 2003, primarily because it produces the smoothest, most realistic slope histograms, and reasonable aspect distributions  It still has problems with smoothing of slopes in valleys, peaks, and ridges.


Guth (1995) reviewed all the independent slope and aspect algorithms he could find. He suggested that the steepest adjacent neighbor provided the best estimate of slope, and the 8 neighbors with even weights provided the best estimate of aspect.  More recent comparisons of slope algorithms include Hodgson (1998) and  Jones (1998), but neither cited Guth (1995).

Justification to use largest slope of 8 to adjacent elevations in grid:

  1. Each point is surrounded by 8 neighbors. Four are to the N,S,E,& W, and 4 approximately NE,SE,SW,& NW. How close those are depend upon the type of DEM and latitude; it is exact at the equator where x and y spacing is equal or using a UTM based DEM. Withe DTED, the worst cases will be in the vicinity of 50, where just to the south the spacing is 60 and 90 m, and just to the north the spacing is 90 and 120 m. There the diagonal values are at 53 and 34 instead of 45 (arc tan (120/90) and arc tan (60/90)).
  2. The direction of maximum slope can thus occur no more than 28 away from the largest of the eight directional slopes to nearest neighbors in the DTED grid (half of 56, the complement of 34).
  3. If we assume a plane for the surface in the vicinity of the greatest slope, we can use vector arithmetic to project the slope in various directions. The projection will be a function of the cosine. For the worst case, the estimated slope to the DTED value in the closest direction will be 87.9% of the true slope (cos(28) = 0.879), and that occurs only when the true maximum lies directly between two grid values. At the equator, with maximum error in estimating, the estimated slope is 92.4% of its true value.
  4. Since maximum error is only about 10%, and given the overall accuracy of the data to start with, we do not feel it necessary to obtain the added rigor of fitting a surface to the data points in the vicinity and then computing the slope of the tangent surface.
  5. The other methods effectively compute the slope are a region twice the size of the data spacing.
  6. See Guth (1995) for additional details and a comparison on the effect of algorithm on calculated slope.

 

Discussions of slope on arc-second DEMs


Last revised 2/9/2021