How Curvature works

The Curvature tool calculates the second derivative value of the input surface on a cell-by-cell basis.

For each cell, a fourth-order polynomial of the form:

 Z = Ax²y² + Bx²y + Cxy² + Dx² + Ey² + Fxy + Gx + Hy + I
is fit to a surface composed of a 3x3 window. The coefficients a, b, c, and so on, are calculated from this surface.

The relationships between the coefficients and the nine values of elevation for every cell numbered as shown on the diagram are as follows:

Curvature values diagram
Curvature values diagram
 A = [(Z1 + Z3 + Z7 + Z9) / 4  - (Z2 + Z4 + Z6 + Z8) / 2 + Z5] / L4
 B = [(Z1 + Z3 - Z7 - Z9) /4 - (Z2 - Z8) /2] / L3
 C = [(-Z1 + Z3 - Z7 + Z9) /4 + (Z4 - Z6)] /2] / L3
 D = [(Z4 + Z6) /2 - Z5] / L2
 E = [(Z2 + Z8) /2 - Z5] / L2
 F = (-Z1 + Z3 + Z7 - Z9) / 4L2
 G = (-Z4 + Z6) / 2L
 H = (Z2 - Z8) / 2L
 I = Z5

The output of the Curvature tool is the second derivative of the surface—for example, the slope of the slope—such that:

 Curvature = -2(D + E) * 100

From an applied viewpoint, the output of the tool can be used to describe the physical characteristics of a drainage basin in an effort to understand erosion and runoff processes. The slope affects the overall rate of movement downslope. Aspect defines the direction of flow. The profile curvature affects the acceleration and deceleration of flow and, therefore, influences erosion and deposition. The planform curvature influences convergence and divergence of flow.

Displaying contours over a raster may help with understanding and interpreting the data resulting from the execution of the Curvature tool. An example of the process follows:

Interpreting results from Curvature

Displaying contours over a raster may help with understanding and interpreting the data resulting from the execution of the tool. An example of the process follows.

  1. Create a curvature raster:

    Input raster : elev_ras

    Output curvature raster : curv_ras

    Z factor : 1

    Output profile curve raster : profile_ras

    Output plan curve raster : plan_ras

  2. Create contours of the surface raster:

    Input raster : elev_ras

    Output polyline features : cont_lines

    Contour interval : 100

    Base contour : ""

    Z factor : 1

  3. Create a slope raster:

    Input raster : elev_ras

    Output raster : slope_ras

    Output measurement : DEGREE

    Z factor : 1

  4. Then create contours of the slope:

    Input raster : elev_ras

    Output polyline features : cont_lines

    Contour interval : 100

    Base contour : ""

    Z factor : 1

  5. Add the curvature raster as a layer in ArcMap. Overlay the two contour feature datasets just created and apply different color symbology for each.

References

Moore, I. D., R. B. Grayson, and A. R. Landson. 1991. Digital Terrain Modelling: A Review of Hydrological, Geomorphological, and Biological Applications. Hydrological Processes 5: 3–30.

Zeverbergen, L. W., and C. R. Thorne. 1987. Quantitative Analysis of Land Surface Topography. Earth Surface Processes and Landforms 12: 47–56.

Related Topics


6/29/2011