Module Morphometric Features
Uses a multi-scale approach by fitting quadratic parameters to any size window (via least squares) to derive slope, aspect and curvatures (optional output) for subsequent classification of morphometric features (peaks, ridges, passes, channels, pits and planes). This is the method as proposed and implemented by Jo Wood (1996) in LandSerf and GRASS GIS (r.param.scale).
Optional output is described in the following. Generalised elevation is the smoothed input DEM. Slope is the magnitude of maximum gradient. It is given for steepest slope angle and measured in degrees. Aspect is the direction of maximum gradient. Profile curvature is the curvature intersecting with the plane defined by the Z axis and maximum gradient direction. Positive values describe convex profile curvature, negative values concave profile. Plan curvature is the horizontal curvature, intersecting with the XY plane. Longitudinal curvature is the profile curvature intersecting with the plane defined by the surface normal and maximum gradient direction. Cross-sectional curvature is the tangential curvature intersecting with the plane defined by the surface normal and a tangent to the contour - perpendicular to maximum gradient direction. Minimum curvature is measured in direction perpendicular to the direction of of maximum curvature. The maximum curvature is measured in any direction.
References:
Wood, J. (1996): The Geomorphological characterisation of Digital Elevation Models. Diss., Department of Geography, University of Leicester, U.K. online.
Wood, J. (2009): Geomorphometry in LandSerf. In: Hengl, T. and Reuter, H.I. [Eds.]: Geomorphometry: Concepts, Software, Applications. Developments in Soil Science, Elsevier, Vol.33, 333-349.
LandSerf Homepage.
- Author: O.Conrad (c) 2013
- Menu: Terrain Analysis|Terrain Classification
Parameters
| Name | Type | Identifier | Description | Constraints |
Input | Elevation | Grid (input) | DEM | - | - |
Output | Morphometric Features | Grid (output) | FEATURES | - | - |
Generalized Surface (*) | Grid (optional output) | ELEVATION | - | - |
Slope (*) | Grid (optional output) | SLOPE | - | - |
Aspect (*) | Grid (optional output) | ASPECT | - | - |
Profile Curvature (*) | Grid (optional output) | PROFC | - | - |
Plan Curvature (*) | Grid (optional output) | PLANC | - | - |
Longitudinal Curvature (*) | Grid (optional output) | LONGC | - | - |
Cross-Sectional Curvature (*) | Grid (optional output) | CROSC | - | - |
Maximum Curvature (*) | Grid (optional output) | MAXIC | - | - |
Minimum Curvature (*) | Grid (optional output) | MINIC | - | - |
Options | Grid system | Grid system | PARAMETERS_GRID_SYSTEM | Grid system | - |
Scale Radius (Cells) | Integer | SIZE | Size of processing window (= 1 + 2 * radius) given as number of cells | Minimum: 1 Default: 5 |
Slope Tolerance | Floating point | TOL_SLOPE | Slope tolerance that defines a 'flat' surface (degrees) | Default: 1.000000 |
Curvature Tolerance | Floating point | TOL_CURVE | Curvature tolerance that defines 'planar' surface | Default: 0.000100 |
Distance Weighting Exponent | Floating point | EXPONENT | Exponent for distance weighting (0.0-4.0) | Minimum: 0.000000 Maximum: 4.000000 Default: 0.000000 |
Vertical Scaling | Floating point | ZSCALE | Vertical scaling factor | Default: 1.000000 |
Constrain | Boolean | CONSTRAIN | Constrain model through central window cell | Default: 0 |
(*) optional |
Command-line
Usage: saga_cmd ta_morphometry 23 [-DEM <str>] [-FEATURES <str>] [-ELEVATION <str>] [-SLOPE <str>] [-ASPECT <str>] [-PROFC <str>] [-PLANC <str>] [-LONGC <str>] [-CROSC <str>] [-MAXIC <str>] [-MINIC <str>] [-SIZE <num>] [-TOL_SLOPE <double>] [-TOL_CURVE <double>] [-EXPONENT <double>] [-ZSCALE <double>] [-CONSTRAIN <str>]
-DEM:<str> Elevation
Grid (input)
-FEATURES:<str> Morphometric Features
Grid (output)
-ELEVATION:<str> Generalized Surface
Grid (optional output)
-SLOPE:<str> Slope
Grid (optional output)
-ASPECT:<str> Aspect
Grid (optional output)
-PROFC:<str> Profile Curvature
Grid (optional output)
-PLANC:<str> Plan Curvature
Grid (optional output)
-LONGC:<str> Longitudinal Curvature
Grid (optional output)
-CROSC:<str> Cross-Sectional Curvature
Grid (optional output)
-MAXIC:<str> Maximum Curvature
Grid (optional output)
-MINIC:<str> Minimum Curvature
Grid (optional output)
-SIZE:<num> Scale Radius (Cells)
Integer
Minimum: 1
Default: 5
-TOL_SLOPE:<double> Slope Tolerance
Floating point
Default: 1.000000
-TOL_CURVE:<double> Curvature Tolerance
Floating point
Default: 0.000100
-EXPONENT:<double> Distance Weighting Exponent
Floating point
Minimum: 0.000000
Maximum: 4.000000
Default: 0.000000
-ZSCALE:<double> Vertical Scaling
Floating point
Default: 1.000000
-CONSTRAIN:<str> Constrain
Boolean
Default: 0