InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Earth and Planetary Sciences » Geology and Geophysics » "Hydro-Geomorphology - Models and Trends", book edited by Dericks P. Shukla, ISBN 978-953-51-3574-6, Print ISBN 978-953-51-3573-9, Published: October 18, 2017 under CC BY 3.0 license. © The Author(s).

Chapter 5

Digital Elevation Models in Geomorphology

By Bartłomiej Szypuła
DOI: 10.5772/intechopen.68447

Article top

Overview

Example of the GRID model (A), and TIN model (B).
Figure 1. Example of the GRID model (A), and TIN model (B).
Elevation map (A), and simplified Hammond’s landform classification (B).
Figure 2. Elevation map (A), and simplified Hammond’s landform classification (B).
Classification of form elements by plan and profile curvature (after Dikau [124]).
Figure 3. Classification of form elements by plan and profile curvature (after Dikau [124]).
Elevation map (A), and TPI 3-category slope classification (B).
Figure 4. Elevation map (A), and TPI 3-category slope classification (B).

Digital Elevation Models in Geomorphology

Bartłomiej Szypuła
Show details

Abstract

This chapter presents place of geomorphometry in contemporary geomorphology. The focus is on discussing digital elevation models (DEMs) that are the primary data source for the analysis. One has described the genesis and definition, main types, data sources and available free global DEMs. Then we focus on landform parameters, starting with primary morphometric parameters, then morphometric indices and at last examples of morphometric tools available in geographic information system (GIS) packages. The last section briefly discusses the landform classification systems which have arisen in recent years.

Keywords: geomorphometry, DEM, DTM, LiDAR, morphometric variables and parameters, landform classification, ArcGIS, SAGA

1. Introduction

Geomorphology, the study of the Earth’s physical land-surface features, such as landforms and landscapes and on-going creation and transformation of the Earth’s surface, is one of the most important research disciplines in Earth science. The term geomorphology was first used to describe the morphology of the Earth’s surface in the end of nineteenth century [1]. Geomorphological studies have focused on the description and classification of landforms (geometric shape, topologic attributes, and internal structure), on the dynamical processes characterizing their evolution and existence and on their relationship to and association with other forms and processes [2].

Geomorphology dates back to sixth to fifth century BC, when Xenophanes of Colophon (580–480 BC) speculated that, as seashells are found on the top of mountains, the surface of the Earth must have risen and fallen, or Herodotus (484–420 BC) thought that the lower part of Egypt was a former marine bay referring to the year-by-year accumulation of river-borne silt in the Nile delta region [3]. Geomorphology as an independent scientific discipline developed in the late nineteenth century [4]. The first modern theory of landscape evolution was the ‘geographical cycle’ expounded by Davis [57]. Another important theory, or rather a variation on Davis’ scheme was offered by Penck [8, 9] who stated that according to the Davis model, uplift and planation take place alternately.

Geomorphology in the twentieth century experienced rapid evolution and growth, and many overlapping phases of development can be identified. Contemporary geomorphology combines and consists of many individual fields of science, that is, geology, hydrology, meteorology, cartography, geographic information system (GIS), engineering, biology, archaeology, etc. This complexity and the long tradition associated with its development helped geomorphologist to use many research techniques and measurement tools according to the research requirements. It is often said about the four main research directions, which are historical, dynamic, structural, and climatic, in geomorphology. However, as rightly observed by Migoń [10], each area has its own history, dynamic nature of the processes, and is present in specific geological structures and with the participation of specific climatic conditions.

Herein, we will focus on research approach to landform as the main study subject in geomorphology. In general, over the last 50 years four primary approaches were distinguished in geomorphology.

  1. Morphography: It describes physical appearance of landforms and is qualitative approach to the form. It is linked with the direct observation of forms in situ which allows specifying the appearance of the form and its morphographic classification (plain, hill, valley, ridge, etc.). These terms do not indicate the way of creation of forms rather only determine their external expressions.

  2. Morphogenesis: It focuses on explaining the origin of the forms and determine the mechanisms of their contemporary development. Geomorphologists use different methods to determine the nature of the process in the past and the present form.

  3. Morphochronology: It aims to specify age of the forms and the age-relationships between adjacent landforms. Geomorphologists examine both absolute as well as relative age between the forms.

  4. Morphometry: It deals with establishing geometric features of the landforms on the basis of measurements. This chapter is just dedicated to morphometry that at the beginning was an element of quantitative geomorphology and later became an independent discipline in the Earth sciences—geomorphometry.

The quantitative phase in geomorphology was developed in 1940–1970 and reflected a broader trend within many of the Earth sciences disciplines towards enhanced use of sophisticated technologies (often derived from military purposes) to measure, describe and analyse the Earth’s surface features in number of categories. Horton’s publications [11, 12] on stream networks and drainage basin processes are classically identified as the precursor to this quantitative movement. Strahler [13] has limited a method for the quantitative analysis of the forms modelled by flowing water and gravitational movements over a longer period of time. The data were obtained mainly from measurements on detailed topographic maps and then applying statistical and morphometric analysis for various calculations and indicators of the fluvial relief within the drainage basin [1418], karst [19], or glacial relief [20, 21] were made.

Further development of the geomorphometry resulted in working out of new theoretical and methodological basis. Lustig [22] distinguished two approaches for quantitative analysis and Earth’s surface characteristics which are (1) the characteristics of individual forms based largely on field surveys and (2) analysis of the surface area as a whole based on analysis of the map. Then Evans [23] proposed the division of geomorphometry into two forms: general and specific. In general, geomorphometry applies to and describes the continuous land surface. It provides a basis for the quantitative comparison of qualitatively different landscapes and it can adapt to the methods of surface analysis used outside geomorphology. Whereas specific geomorphometry applies to discrete landforms, describes selected relief or landform types as well as their geometry and laws of formation and development. Although such divisions may seem somewhat artificial, since the geometry of individual forms consists of the geometry of the surface as a whole, but they are useful to define more closely the principles of various types of morphometric analysis.

At the end of the 1980s and the beginning of the 1990s, with the personal computer revolution, algorithms have been implemented in many raster-based GIS packages (ArcInfo, MicroStation, MicroDEM, etc.) and it was possible to process digital elevation models (DEMs) over fairly large areas. Then began the new era of geomorphometry with new opportunities to visualize and compute land-surface parameters. Now-a-days, geomorphometry is an important component of terrain analysis and surface modelling including measurements of morphometry of continental ice surfaces, characterizing glacial troughs, mapping sea-floor terrain types, guiding missiles, assessing soil erosion, analysing wildfire propagation, and mapping ecoregions [24, 25]; it is a broad field that is important not only in various aspects of Earth sciences but also in engineering, biology, and medicine [26].

2. Digital elevation models (DEMs)

2.1. History and definition

As Pike [27] noted, a numerical description of the ground surface is helpful in addressing many geomorphological problems. In the definition of the subject of this chapter appears the word ‘model’, that is, a representation, generally in miniature, to show the construction or appearance of something. Meyer [28] neatly expressed it—reality scaled down and converted to a form which we can comprehend. In this case it is used to represent the original situation—approximation of topography. The term and concept of ‘terrain model’ was first described by Miller and Laflamme [29] but it did not come into general use until the 1960s and even later because of the technological limitations (computers, processors, etc). They defined digital terrain model (DTM) as just a statistical representation of continuous Earth’s surface which consists of many points of known co-ordinates x, y, z in the arbitrary co-ordinate system. In the following years, a number of terms related to the presentation of the Earth’s surface by numerical methods, that is, Burrough [30] described DEM as a regular gridded matrix representation of the continuous variation of relief over space. Among many terms, the most common and accepted in geomorphometric and GIS terminologies is digital elevation model (DEM) or digital terrain model (DTM). Some authors make distinctions between them, that is, DEM is ‘an ordered array of numbers that represent the spatial distribution of elevations above some arbitrary datum in a landscape’. DTM is ‘an ordered array of numbers that represent the spatial distribution of terrain attributes’; therefore DEM is a subset of DTM. Moore et al. [31, 32], Weibel and Heller [33], and Li et al. [34] defined DTM as a digital (numerical) representation of the terrain. Herein both terms will be treated as synonyms but, in fact, the concept of DEM is broader and more universal. In summary, we can say that digital elevation model is the set of digital data describing elevation values of Earth’s ground surface (or any other surface) which contains additional information about the character of this surface (i.e. structural lines, break lines, water bodies, etc.) and interpolation algorithm, which is the best for approximation (modelling) of the real topography. A DEM is a complete representation of a land surface which means that heights are available at each point in the area of interest [35].

2.2. Types

Due to geometry and data organization there are several basic types of models: raster (GRID), vector triangulated irregular network (TIN), point and hybrid [34]. Herein we will skip point and hybrid models because they are rarely used and absent in the most popular GIS software packages (i.e. ArcGIS, QGIS, MapInfo, Surfer).

Raster model (GRID) is the most widely used digital height data structure during the last years because of its simplicity [36]. The simplicity consists, in fact, that elevations are stored using a regular square grid that is consistent in each part of the study area. All square grids form regular matrix of heights for which plain co-ordinates (x, y) can be easily calculated due to the regular spacing of the grid points (Figure 1). This kind of DEM has many advantages, such as simple elevation matrices that record topological relations between data points implicitly and ease of computer implementation [31, 37]. Moreover, DEM is considerably easier to design land-surface parameters and objects using grids because simpler algorithms can be used; grids have a uniform spatial structure and almost all properties of gridded DEMs are defined by a single characteristic which is cell size; a grid model is more suited to the computer models used in image processing and for printing [35]. For this reason, some DEM software packages accept only grid data.

media/F1.png

Figure 1.

Example of the GRID model (A), and TIN model (B).

Raster DEMs also have disadvantages, such as grids present under-sampled topography in areas where the topography is complex, and they over-sample smooth topography; re-projection of a grid is slow and sometimes leads to a loss of accuracy (because the initial grid loses its regular structure in a new projection and so it has to be re-calculated); the different distances between grid centres in cardinal and diagonal directions have a negative impact on the precision of many hydrological modelling [35]; the computed upslope flow paths will tend to have zigzag pattern across the landscape and increase the difficulty of calculating specific catchment areas accurately [38, 31]; square grids cannot handle abrupt changes in elevation easily and they will often skip important details of the land surface in flat areas [39].

The second most popular and widely used is a vector model named triangulated irregular network (TIN). TIN is proposed by Hormann [40] which devised a TIN idea, linking selected points on divides, drainage lines and breaks in slope to interrelate height, slope gradient and aspect. TIN is based on triangular elements (facets that vary in shape and size) with vertices at the sample height points [31]. These facets consist of planes joining the three adjacent points in the network and are usually constructed using Delaunay triangulation [33]. TIN can be interpolated directly from surveyed points or discrete features that are extracted manually from maps or by computer from a grid or contour DEM [27]. The triangle may be regarded as the most basic and universal unit in all geometrical patterns (because of their great flexibility in terms of shape and size), since a regular grid of square or rectangular cells or any polygon with any shape can be decomposed into a series of triangles [34].

The advantages of such a model are: relevance to gravitational movements, and especially to hydrological applications [41]; TINs are widely used in perspective representations of surfaces, especially in dynamic fly-through displays where the foreground is represented in full detail but the background can be simplified as larger triangles; for areas with high relief or rougher surfaces, irregular DEMs can use smaller spacing between points, and larger spacing where relief is lower or where the surface is smoother and in this way they can more accurately describe geological faults and other sharp elevation changes using the same number of points as grids [35]; the TIN structure gives best reflecting processes of erosion and deposition mimics paths of steepest gradient [41]. Of course we can use the algorithms to convert the TIN into grid and vice versa, but usually at the loss on the quality and accuracy of the model.

2.3. Data sources

DEM source means the data acquired by different techniques and from different sources. In general, such data can be acquired from topo maps and from field surveys.

Manual or semi-automated digitizing of contour lines and every single elevation point on topographic maps on the computer screen or digitizer can be done relatively cheaper without the need for resurveying. Derived DEMs represent the underlying terrain without surface vegetation and buildings. Despite many critical voices regarding the use of this method [4243], it is a quick and effective method to get quite good DEM and till the end of 1990s this was the most common source of elevation data for DEMs.

The second group of data sources are surveys, such as ground survey techniques, digital photogrammetry, and remote sensing. Field surveys are carried out using total stations, theodolite, levelling instruments, global positioning system (GPS), and differential global positioning system (DGPS). These kinds of surveys have high accuracy (sometimes even less than 1 cm), flexibility (the measurement density can be varied, depending on the terrain), and very little processing is required after the measurements have been taken, but especially suited for measuring small areas. Digital photogrammetry relies on the stereoscopic interpretation of aerial photographs or satellite imagery using manual or automatic stereoplotters [39, 33]. Remote sensing surveys consist of Airborne laser scanning (ALS) or satellite platforms and include laser-ranging altimetry (LiDAR—light detection and ranging), synthetic aperture radar interferometry (InSAR or IfSAR). The final results will often include tree-top canopies and buildings. This gives higher elevation values, rough surfaces, and high slope values [43]. The advantages of the LiDAR method are production time that is typically shorter than that for photogrammetrically generated DEMs [44] and great spatial and vertical accuracy (<0.5 m). The main disadvantage of LiDAR data is point cloud which produces a very dense and detailed land-surface model that could be difficult to handle during the production process and sometimes the accuracy of the readings vary according to the characteristics of the terrain (i.e. very steep slopes) [45]. However, these days LiDAR is definitely the best method of the DEM production. Several countries have already produced national LiDAR DEM/DSM (e.g. Belgium, the Netherlands at resolutions of 2–5 m, Poland at 1 m).

2.4. Free global DEMs

Currently, many elevation data with better accuracy and parameters are freely available in the world including LiDAR data. Individual countries or institutions offer different kinds of models. In Table 1, there are models (and their properties) available for the whole world, completely free of charge.

ModelSpatial resolution (grid size)Accuracy vertical/horizontalCo-ordinate SystemData formatData extentSourceInstitution, publication
GTOPO3030 × 30″ (~1000 m)10–300/150 mWGS-84DTED, USGS DEM, DCW90°N to 90°Shttps://lta.cr.usgs.gov/GTOPO30USGS (USA) and NASA, UNEP/GRID, USAID, INEGI (Mexico), GSI (Japan), MWLR (New Zealand), SCAR 1996
GLOBE DEM30 × 30″ (~1000 m)250/160 mWGS-84DTED90°N to 90°Shttps://www.ngdc.noaa.gov/mgg /topo/gltiles.htmlNOAA, 1999
SRTM v31 × 1″ (~30 m)
3″ × 3″ (~90 m)
10/13 mWGS-84/EGM96DTED, BIL, GeoTIFF, hgt60°N to 54°Shttps://earthexplorer.usgs.gov/NASA, NGA, 2013
ETOPO11 × 1″ (~30 m)?WGS-84netCDF, GRD98, xyz, GeoTIFF90°N to 90°Shttps://www.ngdc.noaa.gov/mgg/global/global.htmlNOAA, 2009
ASTGTM v02
AST14DEM v03
1 × 1″ (~30 m)
1 × 1″ (~30 m)
20/30 mWGS-84/EGM96
WGS-84/UTM
GeoTIFF83°N to 83°Shttps://search.earthdata.nasa.gov/search?q=ASTGTM+V002NASA, METI, 2009, 2011
GMTED20107.5′ × 7.5″ (~250 m)
15 × 15′ (~500 m)
30″ × 30′ (~1000 m)
26–30 m/?
29–32 m/?
25–42 m/?
WGS-84ESRI ArcGrid,
GeoTIFF
84°N to 56°Shttps://topotools.cr.usgs.gov/gmted_viewer/viewer.htmUSGS, NGA (USA), 2011
AW3D301 × 1″ (~30 m)5/5 mITRF97/GRS80GeoTIFF82°N to 82°Shttp://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htmJAXA, 2016

Table 1.

Basic information about the free global DEMs.

GTOPO30 (Global 30 arc-second Elevation): This is a global DEM with a horizontal grid spacing of 30 arc seconds (approx. 1 km) resulting in a product having dimension of 21,600 rows and 43,200 columns. GTOPO30 was derived from eight several raster and vector sources of topographic information which are digital terrain elevation data (DTED), Digital Chart of the World, USGS 1-degree DEM’s, army map service (AMS) 1:1,000,000, International Map of the World 1:1,000,000, Peru Map 1:1,000,000, New Zealand DEM, and Antarctic Digital Database [46]. The horizontal accuracy is ±150 m linear error at 90% confidence and vertical accuracy is 10–300 m at the 90% confidence level [4749].

GLOBE DEM (Global Land 1-km base elevation DEM): It is one of the first global DEM, which was initially spearheaded in 1990. This data set covers 180°W to 180°E longitude and 90°N to 90°S latitude. The horizontal resolution is 0.5 arc-minute in latitude and longitude, resulting in dimensions of 21,600 rows and 43,200 columns [50]. GLOBE version 1.0 has 11 broad sources of information: 6 gridded DEMs and 5 cartographic sources, which were adapted for use in GLOBE. These were digital terrain elevation data (DTED), Digital Chart of the World, Australian DEM, Antarctic Digital Database, Brazil Maps 1:1,000,000, DEM for Greenland, army map service (AMS) Maps 1:1,000,000, DEM for Japan, DEM for Italy, DEM for New Zealand, and Peru Map 1:1,000,000 [51]. Vertical accuracy expressed as ±30 m linear error at 90% confidence can also be described as a root mean square error (RMSE) of 18 m. Linear error distribution is 97 m for RMSE and it can be expressed as ±160 m linear error at 90% confidence.

Shuttle radar topography mission v3 (SRTM): It was flown to aboard the space shuttle Endeavour February 11–22, 2000. The National Aeronautics and Space Administration (NASA) and the National geospatial-intelligence agency (NGA) participated in an international project to acquire radar data that was used to create near-global set of land elevations. SRTM successfully collected radar data over 80% of the Earth’s land surface between 60°N and 56°S latitude with data points posted every 1 arc-second (approximately 30 m). Absolute height error of SRTM data sets is 5–10 m and absolute geolocation error is 7–13 m [52]. The level of processing and the resolution of the data vary by SRTM data set:

1. SRTM Non-Void Filled. This version was edited or finalized by the NGA to delineate and flatten water bodies, better define coastlines, remove spikes and wells, and fill small voids. Data were sampled at 1 arc-seconds (USA) and 3 arc-seconds (rest of world);

2. SRTM Void Filled. This elevation data are the result of additional processing to address areas of missing data or voids in the SRTM Non-Void Filled collection. The resolution for SRTM is 1 arc-seconds (USA) and 3 arc-seconds (rest of world);

3. SRTM 1 Arc-Second Global. This elevation data offer worldwide coverage of void filled data at 1 arc-second (30 m) resolution and provides open distribution of this high-resolution global data set. Some tiles may still contain voids. Please note that tiles above 50°N and below 50°S latitude are sampled at a resolution of 2 by 1 arc-second [53].

ETOPO1 (Global Relief Model at 1 arc-minute resolution): It was made in 2008 by the National Geophysical Data Centre (NGDC), an office of the National Oceanic and Atmospheric Administration (NOAA). This model was developed as an improvement to the ETOPO2v2 Global Relief Model. ETOPO1 is available in two versions: ‘Ice Surface’ (top of Antarctic and Greenland ice sheets) and ‘Bedrock’ (base of the ice sheets). These versions of ETOPO1 were generated from diverse (regional and global) digital data sets, which were shifted to common horizontal and vertical datum. Next steps were evaluated and edited as needed [54]. Shoreline, bathymetric, topographic, integrated bathymetric–topographic, and bedrock digital data sets (13 sources) were obtained from several U.S. government agencies, international agencies, and academic institutions.

ASTGTM (ASTER global digital elevation model v002 and AST14DEM—ASTER digital elevation model v003): These DEMs were developed jointly by the US NASA and Japan’s Ministry of Economy, Trade and Industry (METI). ASTER is capable of collecting in-track stereo using nadir and aft-looking near infrared cameras. Since 2001, these stereo pairs have been used to produce single-scene (60 × 60 km) DEMs having vertical RMSE accuracies generally between 10 and 25 m [55]. The ASTER GDEM covers land surfaces between 83°N and 83°S and is comprised of 22,702 tiles. This model is distributed as georeferenced tagged image file format (GeoTIFF) files. The data have resolution at 1 arc-second (approximately 30 m) grid and referenced to the 1984 World Geodetic System (WGS84)/1996 Earth Gravitational Model (EGM96) geoid. Although the ASTER GDEM v. 002 is better model than ASTER GDEM v. 001, users have to know that the data still may contains anomalies and artefacts. One should know that these mistakes can introduce large elevation errors on local scales [56].

GMTED2010 (Global Multi-resolution Terrain Elevation Data 2010): This model was made by collaborating USGS and the NGA, which replaced GTOPO30 model. The new elevation data set has been generated at three separate horizontal resolutions of 30 (about 1 km), 15 (about 500 m), and 7.5 arc-seconds (about 250 m). The global aggregated vertical accuracy of GMTED2010 can be summarized in root mean square error (RMSE). At 30 arc-seconds, the RMSE range is 25–42 m; at 15 arc-seconds, range is 29–32 m, and at 7.5 arc-seconds, range is in between 26 and 30 m. This new product suite provides global coverage of all land areas from latitude 84°N to 56°S for most products, and coverage from 84°N to 90°S for several products [57]. An additional advantage of the new multi-resolution model over GTOPO30 is that at each resolution seven new raster elevation data sets are available. The new models have been produced using the various aggregation methods: minimum elevation, maximum elevation, mean elevation, median elevation, standard deviation of elevation, systematic subsample, and breakline emphasis. GMTED2010 is based on data derived from different raster-based elevation sources, such as SRTM DTED, non-SRTM DTED, Canadian digital elevation data (CDED), Satellite Pour l’Observation de la Terre (SPOT 5), Reference3D, National elevation dataset (NED), GEODATA, an Antarctica satellite radar and laser altimeter DEM, and Greenland satellite radar altimeter DEM [57].

AW3D30 (ALOS Global Digital Surface Model): The global digital surface model (DSM) at 1 arc-second (approximately 30 m) resolution that was released by Japan aerospace exploration agency (JAXA). This model has been compiled with images acquired by the advanced land observing satellite (ALOS). The elevation data are published based on the DSM data set (5-m mesh version) of the ‘World 3D Topographic Data’ [58]. The source data were a huge amount of stereo-pairs images derived from satellite mission in the years 2006–2011. Next, they were processed semi-automatically to provide digital surface model (DSM). The height accuracy of the data set is approximately <5 m from the evaluation with ground control points (GCPs) or reference DSMs derived from the LiDAR [59].

In addition to the above global DEMs, there are many elevation data with much better accuracy, but for smaller areas, for example, for Europe EU-DEM with grid size 25 m, for Spain MDT05/MDT05-LIDAR (5 m), for Netherlands AHN3 DTM (0.5 m), for all Slovenia (LiDAR data), or LiDAR data for Haiti (1 m).

3. Morphometric landform properties

The main object of the studies in geomorphology is the relief (of the Earth’s surface). The term relief of the Earth’s surface was used to describe the vertical dimension or amplitude of topography [60] or as the elevation difference over a pre-determined area or collective elevations and their inequalities of a land surface [61]. One can also say that the relief is the complexity of the Earth’s surface shapes, and these shapes form the landforms. If we would present a continuous Earth’s surface as a matrix of discrete points with a defined distance interval from one another, the shape dimensions would express in the changes of height, distance, and direction between these points. These three vectors and the relationships between them are the basis for most morphometric landform properties.

3.1. Morphometric parameters

Elementary data used in geomorphometry are DEM. Strength of DEM is a suggestive (plastic and three-dimensional) visual message and the ability of quantifying the topography. Quantification of the Earth’s surface is expressed by topographic attributes which can be determined from the derivatives of the topographic surface. Generally, one can say, these derivatives measure rate at which elevation changes in response to changes in location. Deng et al. [62] noted that local terrain shape (as the continuous variation of elevation values over the terrain surface from point to point) has an enormous impact on local terrain attributes, but this role is influenced by data and computational factors.

Theoretical assumptions concerning the morphometric properties of the surface come from the 1950s [6365], and even earlier. Stone and Dugundji [66] and Hobson [67] stated that measures of landforms can be considered as a kind of the roughness of the surface. In general, roughness refers to the irregularity of a topographic surface and cannot be completely defined by any single measure but must be represented by a roughness vector or set of parameters. With roughness, concept of wavelength and amplitude ideas was related. The significant wavelengths of topography are termed grain or texture, while amplitudes associated with these wavelengths correspond to the concept of relief [60]. Texture and grain are terms that have been used to indicate in some way the scale of horizontal variations in the topography. Texture is used to refer the shortest significant wavelength in the topography and grain is used for the longest significant wavelength. Texture is related to the smallest landform elements, while the grain is related to the size of area over which one measures other parameters. Wood and Snell [68] defined grain as the size of area over which the other factors are to be measured. It is dependent on the spacing of major ridges and valleys and thus indicates texture of topography. Texture refers to the shortest significant topographic wavelength.

The systems proposed by Evans [23] and Krcho [69] for field variables are local-based, relating to a vanishingly small area around each point. They divided topographic properties into primary geomorphometric parameters (height, slope, aspect, profile, and planar curvature) and statistical measures derived in square-matrices out of the gridded DEMs (such as relief, standard deviation, and skewness of heights).

Mark [60] stated, that probably the most important single class of processes which has shaped the Earth’s surface could be divided into geometrical properties that involve the relationships among dimensional properties, such as elevation, lengths, areas, volumes, and topological properties which relate numbers of objects in the drainage network (e.g. the bifurcation ratio). Then Moore and Thornes [70] developed LEAP-land erosion analysis programs to examine the spatial distribution of slope length, slope steepness, and the plan curvature for assessing topographic erosion potential.

Pike [71] noted that terrain can be abstracted using geometric signature—a set of measurements that describes topographic form sufficiently well to distinguish geomorphically disparate landscapes. Then he developed the concept of a geometric signature, a multi-variate description of topography using a suite of measures, and later expanded the concept with a listing of 49 variables that could be grouped into 22 attributes [72, 73]. He considered roughness and height, the two most important attributes, with two measures of texture at seventh and eleventh position. Fifteen different variables contribute to roughness. Pike et al. [74] also referred to grain concept. Relief at topographic grain is an estimate of local relief optimized by varying unit-cell size. In homogeneous terrain, local relief within nested circles increases with circle size and then levels off at a diameter termed grain, a measure of characteristic local ridgeline-to-channel spacing. Topographic grain is the characteristic horizontal spacing of major ridges and valleys and is an important descriptor of meso-scale topographic texture. The grain concept arose from the need for a variable and non-arbitrary unit-cell size and also enables to calculate local relief parameter.

Moore and Grayson [41] and Moore et al. [31, 32] described terrain attributes adapted from Speight [75, 76], and then often mentioned and cited [3, 36, 7780]. They divided terrain attributes into primary and secondary (compound) attributes. Primary attributes, hydrologically related, are directly calculated from elevation data and include variables, such as altitude, upslope height, aspect, slope, upslope slope, dispersal slope, catchment slope, upslope area, dispersal area, catchment area, specific catchment area, flow path length, upslope length, dispersal length, catchment length, profile curvature, and plan curvature. It is interesting that Bork and Rohdenburg [81] have developed a digital relief model for estimating and displaying the distribution of these morphographic parameters.

Schmidt and Dikau [78] subdivided primary geomorphometric parameters into three types: simple, complex, and combined. Simple primary geomorphometric parameters (usually derived through a filter operation within a 3 × 3 moving window) are height, slope, aspect, profile curvature, contour curvature, drainage direction, and real area of pixel. Complex geomorphometric parameters are derived through the analysis of the whole matrix of a DEM. They contain structural information about the surrounding morphometry and consist of contributing area, mean slope of contributing area, average, and variance of primary parameters in contributing area, length of flowpath to outlet, average and variance of primary parameters in flowpath to outlet, length of flowpath to stream, average and variance of primary parameters in flowpath to stream, x- and y-co-ordinates of corresponding stream point, height of corresponding stream point, height distance to corresponding stream point, length of minimum flowpath to watershed, relative slope position after minimum, slope length, relative slope position after maximum, and minimum and maximum slope length.

Interesting approach to the morphometric parameters was presented by Shary et al. [82], who stated that morphometric variables describe not the land surface itself, but rather the system, land surface + vector field, where vector fields of common interest are gravitational field and solar irradiation. He divided morphometric variables and concepts into field-specific (may refer to this system description) and field-invariant (invariant with respect to any vector field, that is, describing the land surface geometrical form). On the other hand, morphometric variables may divide into local, regional, or global (when height data of all the Earth is needed for their determination).

Basso [79] divided topographic attributes into: (i) local (calculated from a small neighbouring area surrounding the DEM cell, usually 3 × 3), (ii) regional (calculated using considerably larger geometric area than the local attributes, less sensitive to the DEM resolution), (iii) catchment oriented (related to the whole catchment area, and are the measurement of certain catchment characteristics), and (iv) process oriented (describe or characterize the spatial variability of a simple representation of specific processes that occur on the landscape).

Goodwin and Tarboton [84] described five categories of the morphometric parameters of drainage basin which are size properties (provide measures of scale that can be used to compare the magnitudes of two or more drainage basins), surface properties (quantities depicted by fields comprising a value at each point within a domain, drainage basin), shape properties (i.e. length, width, perimeter, and more complex function), relief properties (total basin relief, relief ratio [85], and hypsometric curve), and texture properties (amount of landscape dissection by a channel network).

Olaya [86] presented division of the morphometric parameters into local and regional. Local parameters consist of geometric (slope, aspect, curvatures, visibility, visual exposure, and visibility index) and statistical (i.e. average value, standard deviation, skewness coefficient, kurtosis coefficient, range of values, etc.). Regional parameters are connected with hydrological properties (catchment area, height, slope, proximity, etc. [86]).

The latest approach to the classification system of the fundamental geomorphometric variables is presented by Evans and Minár [87]. They proposed field and object variables. Field variables include specific to gravity field (local point-based, local area-based, and regional), specific to other fields and field-invariant variables (local point- and regional-based). Object variables differ between areal, linear, and point features.

Thus the optimal number of variables depends on spatial scale, resolution of the source data (DEM), and the requirements of the research problem; and many measures describe the same attribute of surface form and thus are redundant; for this reason new geomorphometric parameters are very rare.

3.2. Morphometric indices

Morphometric parameters, which were discussed in the previous section, showed that there are many different classifications of these parameters. In this section, we want to look at some popular morphometric indices called secondary or composite attributes [3, 32, 36, 62, 80, 89], combined or compound geomorphometric parameters [78, 83], statistical parameters [23, 86] or process oriented [79]. These indices are combinations of the primary morphometric attributes and describe or characterize the spatial variability of specific processes occurring on the landscape. Sometimes these morphometric indices can be derived empirically but it is preferable to develop them through the application and simplification of the underlying physics of the processes. With the index approach we simplify some physical sophistication to allow improved estimates of spatial patterns in the landscape [31].

In Table 2, some selected geomorphometric indices commonly used and their definitions were listed.

ParameterFormulaDescriptionSource
Drainage densityDD= L/A
where: L is sum of the channel lengths and A is basin area
The sum of the channel lengths divided by basin area. It is important indicator of the linear scale of landform elements in a drainage basin, indicates the closeness of spacing of channels, thus providing a quantitative measure of the average length of stream channel for the whole basin.Horton [12]
Form factorRf= A(lb)2
where: A is area of basin, Lb is the basin length
The ratio of the basin area to the square of the basin length. Indicates the flow intensity of a basin of a defined area. The form factor value should be always less than 0.7854 (the value corresponding to a perfectly circular basin).Horton [12]
Terrain ruggedness index (TRI)TRI = Y[Σ(xijx00)2]1/2
where: xij is elevation of neighbour cell to cell (0,0)
Measure of topographic heterogeneity (amount of elevation difference between adjacent cells of a DEM).Riley et al. [90]
Stream power indexΩ = pgq.tanβ)
where: pg is the unit weight of water, q is the discharge per unit width and β is the representative slope angle
This is the time rate of energy expenditure and has been used extensively in studies of erosion, sediment transport and geomorphology as a measure of the erosive power of flowing water.Moore et al. [31]
Elevation-relief ratioE = (Hmean Hmin)/(Hmax Hmin)
where: H is elevation
Expresses relative proportion of upland to lowland within a sample region. Usually ranges from 0.15 to 0.85. Low values occur in terrains characterized by isolated relief features standing above extensive level surfaces, and high values describe broad and level surfaces, broken by occasional depressions.Pike and Wilson [91] after Wood and Snell [68]
Topographic openness algorithmPositive openness
φL = (0φL+45φL+…+315φL)/8
negative openness
L =(0ⵖL+45ⵖL+…+315ⵖL)/8
where: L is specified distance
Describes the degree of enclosure of a location on an irregular Earth’s surface.
Positive values expressing openness over surface topography is high for convex forms, whereas negative values expressing and describing attribute below the surface topography and are high concave forms.
Yokoyama et al. [92]
Terrain shape index_
TSI = Z/R
where: Z is mean elevation of the sample plot boundary, and R is plot radius measured in the units used for elevation
This index is equivalent to the mean slope gradient of the plot boundary as viewed from the plot centre, with units of meters change in elevation per meter of plot radius.
Typical TSI values for mountain landforms ranged from −0.24 to −0.12, on convex upper slopes near ridge tops, and from +0.09 to +0.17 on concave lower slopes.
McNab [93]
Landform indexLI = H°/100
where: H° is vertical gradient to the horizon
The landform index is the average vertical gradient to the topographic horizon, divided by 100 to convert percent to a decimal value.
The index is dimensionless and the effects of height and distance to the landform are compensating factors.
McNab [94]
Compound topographic index (CTI)CTI = ln (Af/tanβ)
where: Af is the specific catchment area draining through the point, and β is the representative slope angle
Ratio between slope and catchment area; quantification of catenary topographic convergence represented by slope angle and catchment. For the same contributing area CTI values are higher for pixels with lower slopes —this means that CTI primarily reflects accumulation processes.Moore et al. [32]
Heat load indexHLI = [1 − cos(θ−45)]/2
where: θ is aspect in degrees east of north
Quantitative measure of aspect and steepness of slope. The equations applied to 0–60° north latitude, slopes from 0 to 90° and all aspects.McCune and Keon [95]
Anisotropy indexANI = Ȓmin/Ȓmax
where: Ȓmin is smallest estimated range parameter, and Ȓmax is highest estimated range parameter in various directions
Ratio between the minimum and maximum range parameter of spatial dependence, fitted for various directionsBishop and Minasny [96]
Shape complexity index​SCI = P / (2rπ ) , r = ​√ __ ​ A __ π ​ ​​
where: P is the perimeter of the polygon, A its area, r is the radius of the circle with the same area
Index which is used to describe polygons on DEM slices. Indicates how compact (or oval) a feature is.Hengl et al. [97]
Basin relief ratioRh = H/L where: H is total basin relief, and L is basin lengthRatio between total basin relief (difference in elevation of basin mouth and summit) and basin length, measured as the longest dimension of the drainage basin. Indicates overall slope of the watershed surface. It is a dimensionless number, readily correlated with other measures that do not depend on total drainage basin dimensions.Schumm [98]
Relative reliefRp = H/P where H is total basin relief, and P is basin perimeterRatio between total basin relief and drainage basin perimeter.Melton [99, 100]
Drainage basin compactnessBc = P/A
where: P is drainage basin perimeter, and A is drainage basin area
Ratio between perimeter and area of drainage basin. Higher values correspond to the basins of developing the long-term share erosion running in conditions of relative peace tectonic or are typical for catchment formed in low resistance rocks.Engstrom [101]
Drainage basin shape ratioBs = Bl/Bw
where: Bl is max length of the drainage basin, and Bw is max width of the drainage basin
Ratio between maximum length and maximum width of drainage basin. Higher values correspond to more elongated basins and also indicate a relatively higher tectonic activity of the area.Cannon [102], Ramírez-Herrera [103]
Surface area ratioSAR = A/Ap
where: A is cell’s surface area, and Ap is planimetric area of that cell
Ratio between surface area and planimetric area; surface ratios will always be greater than or equal to 1.Jennes [104]
Valley height-width ratioVf = 2Vfw/[(Eld Esc) + (Erd Esc)]
where: Vfw is width of the valley floor, Eld/Erd height of the left/right watershed and Esc height of the valley floor
Ratio of the width of the valley bottom and the average height of the slopes; allows comparison of erosional patterns between watersheds. Low index value Vf characteristics are deeply cut river valleys developing under high uplift area, while high values correspond to the valleys in the wide areas of thin tectonic activity.Bull and McFadden [105]
Dissection indexDI = RR/AR
where: RR is relative relief, and AR is absolute relief
Ratio between relative and absolute relief which always varies between zero (complete absence of dissection) and one (extreme case, vertical cliff). Assess the degree of incision of a landscape.Sharma [106]
Relief index (RI)RI = CL/AP
where: CL is total length of contour lines, and AP is planar surface area
Ratio of the summary length of the contour lines and the planar surface area at which they occur; shows relief variability. RI is based on a combination of local relief (number of contour lines and elevational changes) and degree of surface cut (length and shape of the contour lines) with reference to the planar surface area.Szypuła [107]

Table 2.

Selected geomorphometric indices derived from DEMs.

3.3. Morphometric tools available in GIS packages

The current availability of high speed computing platforms and high-resolution (less than 10 m spatial resolution) DEMs now provides the opportunity to perform quantitative analyses and calculating morphometric indices on new level. Many GIS packages offer tools to work with DEMs. On one hand we have the comWmercial software (ArcGIS, MapInfo, Surfer, Global Mapper, Terra Solid), and on other one, great free programs (SAGA, QGIS, MicroDEM), which offer a lot of interesting tools. In Table 3, we have presented useful tools for geomorphometric (or geomorphological) analysis. As examples we chose one commercial package (ArcGIS) and one free (SAGA). One should remember that rapid change of the GIS applications cause, next new versions a large number of greater tools.

SAGA 3.0.0ArcGIS 10.5
Terrain analysis—Channels
Channel network
Channel network and drainage basins
Overland flow distance to channel network
Strahler order
Valley depth
Vertical distance to channel network
Watershed basins
Watershed basins (Extended)
3D analyst tools—Raster Surface
Aspect
Contour
Contour list
Contour with barriers
Curvature
Cut fill
Hillshade
Slope
Terrain analysis—compound analyses
Elevation
Analytical hillshading
Slope
Aspect
Plan curvature
Profile curvature
Convergence index
Closed depressions
Total catchment area
Topographic wetness index
LS-factor
Channel network
Drainage basins
Channel network base level
Channel network distance
Valley depth
Relative slope position
3D analyst tools—visibility
An overview of the visibility toolset
Construct sight lines
Intervisibility
Line of sight
Observer points
Skyline
Skyline barrier
Skyline graph
Sun shadow volume
Viewshed
Viewshed 2
Visibility
Spatial analyst tools—solar radiation
Area solar radiation
Points solar radiation
Solar radiation graphics
* Exemplary additional installed tools :
Terrain analysis—hydrology
Cell balance
Downslope area (interactive)
Edge contamination
Flow accumulation (Flow Tracing)
Flow accumulation (Mass-Flux Method)
Flow accumulation (Recursive)
Flow accumulation (Top-Down)
Flow depth (interactive)
Flow path length
Flow sinuosity (interactive)
Flow width and specific catchment area
Isochrones constant speed (interactive)
Isochrones variable speed (interactive)
LS factor
LS-factor, field-based
Lake flood
Lake flood (interactive)
Melton ruggedness number
SAGA wetness index
Slope length
Slope limited flow accumulation
Stream power index
TCI low
Topographic wetness index (TWI)
Upslope area
Upslope area (interactive)
Hypsometry tools
Hypsometric curve
Area gradient tools
Max watershed
Surface area analysis
Surface area
Surface ratio
Flat area
Geomorphometric and gradient metrics
Classify aspect
Linear aspect
Mean slope
Slope impedance
Dissection
Hierarchical slope position
Landform curvature
Roughness
Slope position
Surface/area ratio
Surface relief ratio
Compound topographic index
Heat load index
Integrated moisture index
Site exposure index
Topographic radiation aspect index
Terrain analysis—lighting, visibility
Analytical hillshading
Potential incoming solar radiation
Sky view factor
Topographic correction
Topographic openness
Visibility (points)
Visibility (single point) (interactive)
Topography tools
Beer’s aspect
McCune and Keon heat load index
Landform classification
PRISM data helper
Slope position classification
Solar illumination index
Topographic convergence/wetness index
Topographic position index
Terrain analysis—morphometry
Convergence index
Convergence index (Search radius)
Curvature classification
Diurnal anisotropic heating
Downslope distance gradient
Effective air flow heights
Fuzzy landform element classification
Hypsometry
Land-surface temperature
Mass balance index
Morphometric features
Morphometric protection index
Multi-scale topographic position index (TPI)
Multi-resolution index of valley bottom flatness (MRVBF)
Real surface area
Relative heights and slope positions
Slope, aspect, curvature
Surface specific points
TPI-based landform classification
Terrain ruggedness index (TRI)
Terrain surface classification (Iwahashi and Pike)
Terrain surface convexity
Terrain surface texture
Topographic position index (TPI)
Upslope and downslope curvature
Valley and ridge detection (Top hat approach)
Vector ruggedness measure (VRM)
Wind effect (Windward/Leeward Index)
Wind exposition index
Local relief model
Terrain tools
Terrain ruggedness (VRM)
Riparian topography tools
Calculate flooding height
Calculate inundation area
Height above river
Prepare HAR for flooding
Flow accumulation for both positive and negative values
Land facet analysis
Calculate density raster
Shannon’s diversity index
Identify termini polygons
Topographic position index tools
Mahalanobis distance tools

Table 3.

The list of exemplary morphometric tools available in SAGA and ArcGIS software (after [108110]).

4. Landform classifications

Geomorphology studies the relief. If we want to understand the relief of the Earth’s surface (which is highly complex) we need to simplify and subdivide it into landforms. We have to focus on describing landforms, their spatial arrangement and the processes which led to their formation. Of course, many landforms can be delineated manually using photo-interpretation to assess their form, size, scale, adjacency, surface roughness, hydrological and contextual position but there is always problem with the boundary of the landform. And herein DEMs are helpful. Availability of global DEMs and high accuracy national LiDAR data made new possibilities of analysing these data to extract and classify geomorphic entities. The landform elements can be extracted automatically from DEMs by using land-surface parameters, such as slope, curvatures, catchment area, distance to streams, peaks and depression depth, etc. [111]. The goal of automated extraction of landforms and landform elements using semi-automated or fully-automated algorithms is to find their geometric signature, which Pike [72] defined as a set of measures that describe topographic form well enough to distinguish among geomorphically disparate landscapes.

There are many landform classification systems, and herein we take this issue very briefly and in general. The pioneer works in quantitative systems for landform classification were conducted by Hammond [112], Wood and Snell [68] and Anstey [113, 114] by using topographic contour maps. These studies were aimed primarily towards the systematization and logical interpretation of terrain data to assist in the determination of design criteria for materiel, and secondarily towards the development of a universal system for the quantification of landform data.

In the mid-1960s, Hammond [115] devised the three-level system of regional landform classification, based purely on geomorphometric parameters calculated in the chosen window size (Figure 2). For each window position the following parameters were calculated: (1) percentage of area where the ground is flat or gentle (less than 8% slope); (2) local relief (maximum minus minimum elevation); and (3) profile type (relative proportion of flat or gently sloping terrain that occurs in lowlands or uplands).

media/F2.png

Figure 2.

Elevation map (A), and simplified Hammond’s landform classification (B).

Peucker and Douglas [116] showed a few methods designed to detect pits, peaks, passes, ridges, ravines, and breaks, given an array of sampled, quantized terrain heights. Described methods used local analysis which means the results at each point do not depend on the results already achieved at other points.

As the hillslopes constitute a basic element of all landscapes and a fundamental component of geomorphologic systems, many subjective classifications of slope profiles that were intended to create conceptual classifications of hillslopes have been proposed. Ruhe [117] divided hillslopes into summits, shoulders, backslopes, footslopes, and toeslopes. Dalrymple et al. [118] and Conacher and Dalrymple [119] proposed nine unit classification of hillslopes: (i) interfluve, (ii) seepage slope, (iii) convex creep slope, (iv) fall face, (v) transportational midslope, (vi) colluvial footslope, (vii) alluvial toeslope, (viii) channel wall, and (ix) channel bed (after [120]).

Next, there were several approaches to standardize hillslope units using qualitative terms [121]. Most commonly, a hillslope was described by a series of basic units describing changes in slope, curvature, and processes along the hillslope profile [122].

As Young [123] noted the curvature can be classified into convex, concave, and rectilinear surfaces, Dikau [124] defined basic form elements of the landscape as the combination of three slope profile curvature characteristics and three plan curvature characteristics which lead to nine possible hillslope units (Figure 3). They deliver a disjunctive description of the hillslope surface into units of homogeneous curvature characteristics. Then Dikau [124] proposed digital geomorphological relief model (DGRM) which will generate form facets and elements (as basic relief units) for geomorphological mapping and simulations for the derivation of more complex relief units. Dikau et al. [125] developed and applied an automated method for classifying macro landform types from DEM that was based on analysis of variation in topographic measures within areas defined by moving windows.

media/F3.png

Figure 3.

Classification of form elements by plan and profile curvature (after Dikau [124]).

The original classification of Pennock et al. [126] explicitly assumed that surface form, as described by curvature, could be directly related to surface processes and to relative landform position (divergent/convergent shoulder, backslope footslope and level). Thus, strong profile convexity was assumed to be indicative of upper, water-shedding slope positions; whereas strong profile concavity was associated with lower, water-receiving landform positions, and planar surfaces were associated with backslopes or flat areas. Of course, this pattern is not always adhered to and there are many instances where convex-concave patterns repeat over short distances along a longer hill slope [121].

Speight [121] described 10 morphological types of topographic landform positions: (i) crest, (ii) depression (open, closed), (iii) flat, (iv) slope, (v) simple slope, (vi) upper slope, (vii) mid-slope, (viii) lower slope, (ix) hillock, and (x) ridge (after [120]).

Very simple and interesting method for classifying relief on the base DEM is topographic position index (TPI). TPI is the difference between the elevation at a cell and the average elevation in a neighbourhood surrounding that cell (Figure 4). Positive values means that the cell is higher than its neighbours (indicate ridges, hills, etc.) while negative values means the cell is lower (indicate canyons, valleys, etc). TPI is a simplification of the landscape position index (LPI) described by Fels and Zobel [127] and was developed in detail by Weiss [128]. TPI values provide a simple and powerful means to classify the landscape into morphological classes [129]. TPI is naturally very scale-dependent, it means neighbourhood size (and shape) and DEM resolutions are critical to the final analysis, so the work with this index should be based on experiments with different threshold values.

media/F4.png

Figure 4.

Elevation map (A), and TPI 3-category slope classification (B).

Drăguţ and Blaschke [130] presented an automated classification system of landform elements based on object-oriented image analysis. Firstly, one need to derive elevation, profile curvature, plan curvature, and slope gradient from DEM. Next, relatively homogenous objects are determined at several levels through segmentation of the image. These object primitives are classified as landform elements using a relative classification model, built both on the surface shape and on the altitudinal position of objects. The classification has nine classes which are peaks and toe slopes, steep and flat/gentle slopes, shoulders and negative contacts, head slopes, and side and nose slopes. Classes are defined using flexible fuzzy membership functions.

Iwahashi and Pike [88] developed an iterative procedure that classifies topography automatically into terrain types, grid cell by grid cell, on the basis of three morphometric variables, such as slope gradient, local convexity, and surface texture. They applied an unsupervised nested-means algorithm for 16 topographic types of the world. They noted the procedure is unsupervised and reflects frequency distributions of the input variables but not defined criteria. It causes that resulting classes are undefined and have to be calibrated empirically by subsequent analysis.

Within the context of defining landform units that maximize internal homogeneity and external differences, Minár and Evans [131] presented the concept of elementary forms (segments and units) defined by constant values of fundamental morphometric properties and limited by discontinuities of the properties. The basic system of form-defining properties represents altitude and its derivatives, constant values of which provide elementary forms with various types of homogeneity.

Drăguţ et al. [132] presented an algorithm to derive elementary forms from DEMs. Elementary forms were defined by constant values of fundamental morphometric properties and limited by discontinuities of these properties. A multi-resolution segmentation technique was customized to partition the layers of altitude derivatives into homogeneous divisions through a self-scalable procedure, which reveals the pattern encoded within the data. Layers were segmented successively, following the order of elevation derivatives (i.e. gradient, aspect, profile curvature, and plan curvature).

Jasiewicz and Stepinski [133] introduced a novel method for unsupervised classification and mapping of landforms from a DEM. This method involves the pattern recognition, not the differential geometry. The foundation of this idea is the concept of geomorphon (geomorphologic phonotypes), which is a simple ternary pattern that serves as an archetype of a particular terrain morphology. A finite number of 498 geomorphons constitute a comprehensive and exhaustive set of all possible morphological terrain types including standard elements of landscape as well as unfamiliar forms rarely found in natural terrestrial surfaces.

And one should remember that horizontal and vertical resolution of the elevation data used to present a terrain surface has a significant influence on the level of detail and the accuracy of portrayal of surface features and on the values of land-surface parameters that are computed from a DEM. One should first test out predictive efficiency for various DEM resolutions and neighbourhood sizes, and then objectively derive the most suitable resolution and search size [120].

Acknowledgements

The publication has been (partially) financed from the funds of the Leading National Research Centre (KNOW) received by the Centre for Polar Studies of the University of Silesia, Poland.

References

1 - de Margerie E. Geologie. Polybiblion Revue Bibliographique Universelle. Partie litteraire. 1886;2(24):310-330
2 - Bauer BO. Geomorphology. In: Goudie AS, editor. Encyclopedia of Geomorphology. London & New York: Taylor & Francis e-Library; 2006. pp. 428-435
3 - Huggett RJ. Fundamentals of Geomorphology. Routledge Fundamentals of Physical Geography. New York: Taylor & Francis e-Library; 2005. p. 386
4 - Klimaszewski M. Geomorfologia. 6th ed. Warszawa: PWN; 1980. p. 1063
5 - Davis WM. The rivers and valleys of Pennsylvania. National Geographical Magazine. 1889;1:183-253. (Also in Geographical Essays)
6 - Davis WM. The geographical cycle. Geographical Journal. 1899;14:481-504
7 - Davis WM. Geographical Essays. Boston, Mass.: Ginn; 1909
8 - Penck W. Die morphologische Analyse, ein Kapitel der physikalischen Geologie. Stuttgart: Engelhorn; 1924
9 - Penck W. Morphological Analysis of Landforms. London: Macmillan; 1953
10 - Migoń P. Geomorfologia. 1st ed. Warszawa: PWN; 2006. p. 461
11 - Horton RE. Drainage basin characteristics. Transactions, American Geophysical Union. 1932;14:350-361
12 - Horton RE. Erosional development of streams and their drainage basins; hydrophysical approach to quantitative morphology. Bulletin of the Geological Society of America. 1945;56:275-370
13 - Strahler AN. Quantitative analysis of watershed geomorphology. Transactions, American Geophysical Union. 1957;38:913920
14 - Strahler AN. Hypsometric (area-altitude) analysis of erosional topography. Geological Society of America Bulletin. 1952;63:1117-1142. DOI: 10.1130/0016-7606(1952)63[1117:HAAOET]2.0.CO;2
15 - Strahler AN. Quantitative slope analysis. Geological Society of America Bulletin. 1956;67:571-596
16 - Strahler AN. Quantitative geomorphology of drainage basins and channel networks. In: Chow VT, editor. Handbook of Applied Hydrology. New York: McGraw-Hill; 1964. pp. 439-476
17 - Schumm S. The relation of drainage basin relief to sediment loss. International Association of Hydrological Sciences Publications. 1954;36:216-219
18 - Chorley RJ. The drainage basin as the fundamental geomorphic unit. In: Chorley RJ, editor. Water, Earth and Man. London: Methuen; 1969. pp. 77-100
19 - Williams PW. Morphometric analysis of polygonal karst in New Guinea. Geological Society of America Bulletin. 1972;83:761-796
20 - Evans IS. The geomorphology and morphometry of glacial and nival areas. In: Chorley RJ, editor. Water, Earth and Man. London: Methuen; 1969. pp. 369-380
21 - Evans IS. The geomorphometry and asymmetry of glaciated mountains [Unpubl. Ph.D. Thesis]. University of Cambridge; 1974
22 - Lustig LK. Quantitative analysis of desert topography. Arid Lands in Perspective. 1969;45-48
23 - Evans IS. General geomorphometry, derivatives of altitude, and descriptive statistics. In: Chorley RJ, editor. Spatial Analysis in Geomorphology. Harper & Row, New York. 1972. pp. 17-90
24 - Pike RJ. Geomorphometry – progress, practice, and prospect. Zeitschrift fur Geomorphologie, Supplementband. 1995;101:221-238
25 - Pike RJ. A Bibliography of Terrain Modeling (Geomorphometry), the Quantitative Representation of Topography - Supplement 4.0. (USGS, Open File Report 02-465). Menlo Park, Calif.: US Department of the Interior, United States Geological Survey; 2002
26 - Evans IS. Geomorphometry. In: Goudie AS, editor. Encyclopedia of Geomorphology. London & New York: Taylor & Francis e-Library; 2006. pp. 435-439
27 - Pike RJ. Digital elevation model. In: Goudie AS, editor. Encyclopedia of Geomorphology. London & New York: Taylor & Francis e-Library; 2006. pp. 260-263
28 - Meyer W. Concepts of Mathematical Modelling. New York: McGraw-Hill Book Company; 1995
29 - Miller CL, Laflamme RA. The digital terrain model — Theory and application. Photogrammetric Engineering. 1958;24(3):433-442
30 - Burrough PA. Methods of interpolation. In: Principles of Geographical Information Systems for Land Resources Assessment. Oxford: OUP; 1986. pp. 147-166
31 - Moore ID, Grayson RB, Ladson AR. Digital terrain modelling: A review of hydrological, geomorphological, and biological applications. Hydrological Processes. 1991;5(1):3-30. DOI: 10.1002/hyp.3360050103
32 - Moore ID, Turner AK, Wilson JP, Jenson S, Band L. GIS and land-surface–subsurface process modelling. In: Goodchild MF, Parkes BO, Steyaert LT, editors. Environmental Modeling with GIS. New York: Oxford University Press; 1993. pp. 196-230
33 - Weibel R., Heller M. Digital terrain modeling. In: Maguire D, Goodchild M, Rhind D, editors. Geographical Information Systems. London: Longman Scientific and Technical; 1991. pp. 269-297
34 - Li Z, Zhu Q, Gold C. Digital terrain modeling: Principles and methodology. Boca Raton: CRC Press; 2005. p. 319
35 - Hengl T, Evans IS. Mathematical and digital models of the land surface. In: Hengl T, Reuter HI, editors. Geomorphometry. Concepts, Software, Applications. Elsevier; 2009. pp. 31-64. DOI: 10.1016/S0166-2481(08)00002-0
36 - Wilson JP, Gallant JC. Digital terrain analysis. In: Wilson JP, Gallant JC, editors. Terrain Analysis: Principles and Applications. John Wiley & Sons, Inc.; 2000. pp. 1-28
37 - Wise SM. The effect of GIS interpolation errors on the use of digital elevation models in geomorphology. In: Lane SN, Richards KS, Chandler JH, editors. Landform Monitoring, Modelling and Analysis. Wiley; 1998. pp. 139-164
38 - Zevenbergen LW, Thorne CR. Quantitative analysis of land surface topography. Earth Surface Processes and Landforms. 1987;12:47-56
39 - Carter JR. Digital representations of topographic surfaces. Photogrammetric Engineering and Remote Sensing. 1988;54(11):1577-1580
40 - Hormann K. Geomorphologische Kartenanalyse mit Hilfe elektronischer Rechenanlagen. Zeitschrift für Geomorphologie. 1969;133(1):75-98
41 - Moore ID, Grayson RB. Terrain-based catchment partitioning and runoff prediction using vector elevation data. Water Resources Research. 1991;27:1177-1191
42 - Wood J. Digital elevation model (DEM). In: Kemp KK, editor. Encyclopedia of Geographic Information Science. SAGE Publications, Inc; 2008. Los Angeles. pp. 107-109. DOI: 10.4135/9781412953962
43 - Nelson A, Reuter HI, Gessler P. DEM production methods and sources. In: Hengl T, Reuter HI, editors. Geomorphometry. Concepts, Software, Applications. Elsevier, Amsterdam; 2009. pp. 65-86. DOI: 10.1016/S0166-2481(08)00003-2
44 - Baltsavias EP. A comparison between photogrammetry and laser scanning. ISPRS Journal of Photogrammetry and Remote Sensing, 1999;54:83-94
45 - Smith SE. Topographic mapping. In: Grunwald S, editor. Environmental Soil–Landscape Modeling: Geographic Information Technologies and Pedometrics. Vol. 1. New York: CRC Press; 2005. pp. 155-182
46 - Gesch DB, Greenlee S. GTOPO30 Documentation. U.S. Geological Survey; 1997. p. 19
47 - DMA (Defense Mapping Agency). Defense Mapping Agency Product Specifications for Digital Terrain Elevation Data (DTED). 2nd ed. St. Louis, Missouri: Defense Mapping Agency Aerospace Center; 1986. p. 26
48 - DMA (Defense Mapping Agency). Digitizing the Future. 3d ed. Washington, D.C.: Defense Mapping Agency; 1990. p. 105
49 - USGS (United States Geological Survey). Digital Elevation Models, data User Guide 5. Reston, Virginia, USGS; 1993. p. 50
50 - Hastings DA, Dunbar PK. Global Land One-kilometer Base Elevation (GLOBE) Digital Elevation Model, Documentation, Volume 1.0. Key to Geophysical Records Documentation (KGRD) 34. 325 Broadway, Boulder, Colorado 80303, U.S.A: National Oceanic and Atmospheric Administration, National Geophysical Data Center; 1999
51 - Hastings DA. Global Land One-km Base Elevation (GLOBE) DEM. Southern Africa Subset; 2000. p. 6
52 - Rodriguez E, et al. An assessment of the SRTM topographic products [Technical Report JPL D-31639]. Pasadena, California: Jet Propulsion Laboratory; 2005. p. 143
53 - Farr TG, et al. The shuttle radar topography mission. Journal of Geophysical Research. 2007;45:43. RG2004, DOI: 10.1029/2005RG000183
54 - Amante C, Eakins BW. ETOPO1 1 arc-minute global relief model: Procedures, data sources and analysis. Boulder, Colorado: National Geophysical Data Center Marine Geology and Geophysics Division; 2009. p. 19
55 - NASA JPL. ASTER Global Digital Elevation Model [Data set]. NASA JPL; 2009. DOI: 10.5067/ASTER/ASTGTM.002
56 - Tachikawa T, Hato M, Kaku M, Iwasaki A. The characteristics of ASTER GDEM version 2; IGARSS, July 2011
57 - Danielson JJ, Gesch DB. Global multi-resolution terrain elevation data 2010 (GMTED2010): U.S. Geological Survey Open-File Report 2011-1073. 2011. p. 26
58 - Tadono T, Ishida H, Oda F, Naito S, Minakawa K, Iwamoto H. Precise global DEM generation by ALOS PRISM. ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences. 2014;II-4:71-76. DOI: 10.5194/isprsannals-II-4-71-2014
59 - Takaku J, Tadono T, Tsutsui K. Generation of high resolution Global DSM from ALOS PRISM. The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Vol. XL-4, ISPRS TC IV Symposium, Suzhou, China; 2014. p. 243-248. DOI: 10.5194/isprsarchives-XL-4-243-2014.
60 - Mark DM. Geomorphometric parameters: a review and evaluation. Geografiska Annaler. 1975;57A(3-4):165-177
61 - Montgomery DR. Relief. In: Goudie AS, editor. Encyclopedia of Geomorphology. London & New York: Taylor & Francis e-Library; 2006. pp. 840-842
62 - Deng Y, Wilson JP, Gallant JC. Terrain analysis. In: Wilson JP, Fotheringham AS, editors. The Handbook of Geographic Information Science. 2008. pp. 417-435. DOI: 10.1002/9780470690819.ch23
63 - Chorley RJ. Climate and morphometry. Journal of Geology. 1957;65:628-638
64 - Chorley RJ. Group operator variance in morphometric work with maps. American Journal of Science. 1958;256:208-218
65 - Chorley RJ, Morgan MA. Comparison of morphometric features, Unaka Mountains, Tennessee and North Carolina, and Dartmoor, England. Bulletin, Geological Society of America. 1962;73:17-34
66 - Stone RO, Dugundji J. A Study of Microrelief: Its Mapping, Classification, and Quantification by Means of a Fourier Analysis. Amsterdam: Elsevier; 1965. p. 97
67 - Hobson RD. Fortran IV programs to determine surface roughness in topography for the CDC 3400 computer. Kansas Geological Survey Computer Contribution. 1967;14:28
68 - Wood W.F, Snell JB. A Quantitative System for Classifying Landforms [Technical Report EP-124]. Natick, MA: U.S. Army Quartermaster Research and Engineering Center; 1960. p. 20
69 - Krcho J. Morphometric analysis of relief on the basis of geometric aspect of field theory. Acta geographica Universitatis Comenianae, Geographico-physica. 1973;1:11-233
70 - Moore RF. Thornes JB. Leap - a suite of Fortran IV programs for generating erosional potentials of land surfaces from topographic information. Computers & Geosciences. 1976;2:493-499
71 - Pike RJ. Information content of planetary terrain: varied effectiveness of parameters for the Earth. In: Lunar and Planetary Science Conference. vol. 18; 1987
72 - Pike RJ. The geometric signature: quantifying landslide-terrain types from digital elevation models: Mathematical Geology. 1988;20:491-512. DOI: 10.1007/BF00890333
73 - Pike RJ. Geometric signatures-experimental design, first results. In: Ohmori H, editor. DEMs and Geomorphometry, Special Publications of the Geographic Information Systems Association, Proceedings of the Symposia on New Concepts and Modeling in Geomorphology and Geomorphometry, DEMs and GIS; In: Fifth International Conference on Geomorphology; 24-26 august 2001; Tokyo. pp. 50-51
74 - Pike RJ, Acevedo W, Card DH. Topographic grain automated from digital elevation models. In: Proceedings of the Ninth International Symposium on Computer Assisted Cartogtraphy. Baltimore, MD: ASPRS/ASCM; 1989. pp. 128-137
75 - Speight JG. A parametric approach to landform regions. Special Publication Institute of British Geographers. 1974;7:213-230
76 - Speight JG. The role of topography in controlling throughflow generation: A discussion. Earth Surface Processes and Landforms. 1980;5:187-191
77 - Schmidt J, Merz B, Dikau R. Morphological structure and hydrological process modelling. Zeitschrift für Geomorphologie. Suppl.-Bd. 1998;112:55-66
78 - Schmidt J, Dikau R. Extracting geomorphometric attributes and objects from digital elevation models — semantics, methods, future needs. In: Dikau R, Saurer H, editors. GIS for Earth Surface Systems — Analysis and Modelling of the Natural Environment. Schweizbart’sche Verlagsbuchhandlung; 1999. pp. 153-173
79 - Basso B. Digital terrain analysis: Data source, resolution and applications for modeling physical processes in agroecosystems. Rivista Italiana di Agrometeorologia. 2005;2:5-14
80 - Wilson JP, Bishop MP. Geomorphometry. In: Shroder JF, Bishop MP, editors. Treatise on Geomorphology, vol 3, Remote Sensing and GIScience in Geomorphology. San Diego: Academic Press; 2013. pp. 162-186
81 - Bork HR, Rohdenburg H. Transferable parameterization methods for distributed hydrological and agroecological catchment models. Catena. 1986;13:99-117
82 - Shary PA, Sharaya SL, Mitusov AV. Fundamental quantitative methods of land surface analysis. Geoderma. 2002;107:1-32. DOI: 10.1016/S0016-7061(01)00136-7
83 - Pennock DJ. Terrain attributes, landform segmentation, and soil redistribution. Soil and Tillage Research, 2003;69(1-2):15-26.
84 - Goodwin CN, Tarboton DG. Morphometric properties. In: Goudie AS, editor. Encyclopedia of Geomorphology. London & New York: Taylor & Francis e-Library; 2006. pp. 696-698
85 - Hadley RF, Schumm SA. Sediment Sources and Drainage-Basin Characteristics in the Upper Cheyenne River Basin, Washington. US Geological Survey Water Supply Paper; 1961. p. 1,531
86 - Olaya V. Basic Land-Surface parameters. In: Hengl T, Reuter HI, editors. Geomorphometry. Concepts, Software, Applications. Elsevier, Amsterdam; 2009. pp. 141-169. DOI: 10.1016/S0166-2481(08)00006-8
87 - Evans IS, Minár J. A classification of geomorphometric variables. In: International Geomorphometry 2011; May 30 2011; Redlands, CA: Geomorphometry.org;. 2011. p. 105-108
88 - Iwahashi J, Pike RJ. Automated classifications of topography from DEMs by an unsupervised nested-means algorithm and a three-part geometric signature. Geomorphology. 2007;86(3-4):409-440. DOI: 10.1016/j.geomorph.2006.09.012
89 - Florinsky I.V.: Accuracy of local topographic variables derived from digital elevation models. International Journal of Geographical Information Science. 1998;12(1):47-62.
90 - Riley SJ, DeGloria SD, Elliot R. A terrain ruggedness index that quantifies topographic heterogeneity. Intermountain Journal of Sciences. 1999;5:23-27
91 - Pike RI, Wilson SE. Elevation-relief ratio, hypsometric integral and geomorphic area altitude analysis. Geological Society of America Bulletin. 1971;82:1079-1084
92 - Yokoyama R, Sirasawa M, Pike RJ. Visualizing topography by openness: A new application of image processing to digital elevation models. Photogrammetric Engineering & Remote Sensing. 2002;68(3):257-265
93 - McNab WH. Terrain shape index: quantifying effect of minor landforms on tree height. Forest Science. 1989;35(1):91-104
94 - McNab WH. A topographic index to quantify the effect of mesascale landform on site productivity. Canadian Journal of Forest Research. 1992;23:1100-1107
95 - McCune B, Keon D. Equations for potential annual direct incident radiation and heat load index. Journal of Vegetation Science. 2002;13:603-606
96 - Bishop TFA, Minasny B. Digital soil-terrain modelling: The predictive potential and uncertainty. In: Grunwald S, editor. Environmental Soil–Landscape Modeling: Geographic Information Technologies and Pedometrics. Boca Raton, FL: CRC Press; 2005. pp. 185-213
97 - Hengl T, Gruber S, Shrestha DP. Digital Terrain Analysis in ILWIS. Lecture Notes. Enschede: International Institute for Geo-Information Science & Earth Observation (ITC); 2003. p. 56
98 - Schumm SA. Evolution of drainage systems and slopes in badlands at Perth Amboy, New Jersey. Geological Society of America Bulletin. 1956;67(5):597-646. DOI: 10.1130/0016-7606(1956)67[597:EODSAS]2.0.CO;2
99 - Melton MA. An analysis of the relation among elements of climate, surface properties and Geomorphology [ONR Technical Report 11]. New York: Columbia University; 1957
100 - Melton MA. Correlation structure of morphometric properties of drainage systems and their controlling agents. Journal of Geology. 1958;66:442-460
101 - Engstrom WN. Morphometric analysis of mountain drainage basins in the Basin and Range Province, USA. Zeitschrift für Geomorphologie N. F. 1989;33:443-453
102 - Cannon PJ. Generation of explicite parameter for a quantitative geomorphic study of the Mill Creek drainage basin. Oklahoma Geology Notes. 1976;36(1):3-16
103 - Ramírez-Herrera MT. Geomorphic assessment of active tectonics in the Acambay graben, Mexican Volcanic Belt. Earth Surface Processess and Landforms. 1998;23:317-332
104 - Jenness J. Surface Areas and Ratios from Elevation Grid (surfgrids.avx) extension for ArcView 3.x, v. 1.2 Jenness Enterprises [Internet]. 2002. Available from: http://www.jennessent.com/arcview/surface_areas.htm (Accessed on: 1 November 2011)
105 - Bull WB, McFadden LD. Tectonic geomorphology north and south of the Garlock Fault, California. In: Doehring DO, editor. Geomorphology in Arid Regions: Annual Binghamton Conference. State University of New York at Binghamton; 1977. pp. 115-136
106 - Sharma HS, editor. Perspectives in Geomorphology. Essays on Indian Geomorphology. vol. 4. New Delhi: Concept Publishing; 1982. p. 358
107 - Szypuła B. Relief Index (RI) as a simple tool for geomorphometry. In: Jasiewicz J, Zwoliński Z, Mitasova H, Hengl T, editors. Geomorphometry for Geosciences. Poznań: Adam Mickiewicz University in Poznań - Institute of Geoecology and Geoinformation, International Society for Geomorphometry; 2015. pp. 127-128
108 - SAGA-GIS Tool Library Documentation v3.0.0 [Internet]. 2016. Available from: http://www.saga-gis.org/en/index.htm (Accessed on: 21 December 2016)
109 - Conrad O, Bechtel B, Bock M, Dietrich H, Fischer E, Gerlitz L, Wehberg J, Wichmann V, Böhner J. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. Geoscientific Model Development. 2015;8:1991-2007. DOI: 10.5194/gmd-8-1991-2015
110 - ESRI (Environmental Systems Resource Institute). ArcGIS Desktop 10.5. Redlands, California; 2016
111 - Evans IS, Hengl T, Gorsevski P. Applications in geomorphology. In: Hengl T, Reuter HI, editors. Geomorphometry. Elsevier, Amsterdam; 2009. pp. 497-526. DOI: 10.1016/S0166-2481(08)00022-6
112 - Hammond EH. Small-scale continental landform maps. Annals of the Association of American Geographers. 1954;44(1):33-42
113 - Anstey RL. Physical characteristics of alluvial fans. Technical Report ES-20, US Army Natick Labs; 1965
114 - Anstey RL. Application of Quantified Terrain Data to Vehicle Design Criteria. In: Proceedings of the 5th US Army Science Conference. West Point; 1966
115 - Hammond EH. Analysis of properties in land form geography: An application to broad-scale land form mapping. Annals of the Association of American Geographers. 1964;54:11-19
116 - Peucker TK, Douglas DH. Detection of Surface-Specific points by local parallel processing of discrete terrain elevation data. Computer Graphics and Image Processing. 1975;4(4):375-387. DOI: 10.1016/0146-664X(75)90005-2
117 - Ruhe RV. Elements of the soil landscape. In: Transactions of the 9th Congress of the International Society of Soil Science. Vol. 4. Madison, WI: International Soil Science Society; 1960. pp. 165-170
118 - Dalrymple JB, Blong RJ, Conacher AJ. A hypothetical nine unit landsurface model. Zeitschrift für Geomorphologie. 1968;12:60-76
119 - Conacher AJ, Dalrymple JB. The nine-unit landsurface model: an approach to pedogeomorphic research. Geoderma. 1977;18(1-2):1-154
120 - MacMillan RA, Shary PA. Landforms and landform elements in geomorphometry. In: Hengl T, Reuter HI, editors. Geomorphometry. Concepts, Software, Applications. Elsevier, Amsterdam; 2009. pp. 227-254. DOI: 10.1016/S0166-2481(08)00009-3
121 - Speight JG. Landform. In: McDonald RC, Isbell RF, Speight JG, Walker J, editors. Australian Soil and Land Survey: Field Handbook. Melbourne: Inkata Press; 1990. pp. 9-57
122 - Dikau R, Rasemann S, Schmidt J. Hillslope, Form. In: Goudie AS, editor. Encyclopedia of Geomorphology. London & New York: Taylor & Francis e-Library; 2006. pp. 516-521
123 - Young A. Slopes. London: Longman; 1972. p. 268
124 - Dikau R. The application of a digital relief model to landform analysis in geomorphology. In: Raper J, editor. Three-dimensional Applications in Geographical Information Systems. London: Taylor and Francis; 1989. pp. 51-77
125 - Dikau R, Brabb EE, Mark RM. Landform classification of New Mexico by computer. Open File Report 91-634. U.S. Geological Survey; 1991. p. 15
126 - Pennock DJ, Zebarth BJ, De Jong E. Landform classification and soil distribution in hummocky terrain, Saskatchewan, Canada. Geoderma. 1987;40(3-4):297-315. DOI: 10.1016/0016-7061(87)90040-1
127 - Fels JE, Zobel R. Landscape position and classified landtype mapping for statewide DRASTIC mapping project. North Carolina State University Technical Report VEL. 95.1; 1995
128 - Weiss A. Topographic positions and landforms analysis (Conference Poster). ESRI International User Conference; 2001; San Diego, CA. p. 9-13
129 - Jenness J. Topographic position index (tpi_jen.avx) extension for ArcView 3.x. Jenness Enterprises [Internet]. 2005. Available from: http://www.jennessent.com (Accessed: 11 December 2017)
130 - Drăguţ L, Blaschke T. Automated classification of landform elements using object-based image analysis. Geomorphology. 2006;81:330-344. DOI: 10.1016/j.geomorph.2006.04.013
131 - Minár J, Evans IS. Elementary forms for land surface segmentation: The theoretical basis of terrain analysis and geomorphological mapping. Geomorphology. 2008;95(3-4):236-259. DOI: 10.1016/j.geomorph.2007.06.003
132 - Drăguţ L, Csillik O, Minár J, Evans IS. Land-surface segmentation to delineate elementary forms from Digital Elevation Models. In: Geomorphometry 2013; 16-20 october 2013; Nanjing, China; 2013
133 - Jasiewicz J, Stepinski TF. Geomorphons — A pattern recognition approach to classification and mapping of landforms. Geomorphology. 2013;182(2013):147-156. DOI: 10.1016/j.geomorph.2012.11.005