Scripts

Instructions:
Please provide software codes including the datasets used for testing and the detailed explanation of the processing steps and algorithms used. You are welcome to use our reference dataset [1].

geomorphometry.org Scripts are intended for the free exchange. Please alert the web-administrator if this script is a registered product with commercial value. The author, who claims creation of this work, expressly publishes it to the public domain. This work is free for the taking and cannot be appropriated by a single author even though it may be included in a copyrighted work. It may be freely used and redistributed and is provided "AS-IS" without warranty of any kind. No technical support is provided.

Average:

Your rating: None Average: 4 (1 vote)

Characteristic scale

Short title:

characteristicScale.lsc

Inputs: A DEM (the Baranja Hill 25m DEM is used as an example). The name of the parameter to measure and the minimum and maximum window sizes over which to measure it.
Outputs: Two rasters, one containing the measured parameter, the other the window size at which the parameter is most extreme.

Purpose and use:

Finds the scale at which a geomorphometric parameter is most extreme for each cell in a DEM. Part of this script appears in Hengl and Reuter (2008) [the Geomorphometry book].
Script to measure surface parameter at characteristic scales. It is designed to incorporate scale-based analysis into surface parameterisation. It measures the given parameter (e.g. slope, profile curvature etc.) at a range of scales, and finds the scale at which that parameter is most extreme. Can be used to explore scale sensitivity of a surface.

Programming environment:

Landserf

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
characteristicScale.lsc_.zip [4]	1.02 KB

Average:

Your rating: None

Conver to Landserf format

Short title:

importData.lsc

Inputs: Baranja Hill data (DEM25m.asc, DEM25srtm.asc, orthophoto.asc).
Outputs: LandSerf files of the input files plus a difference map of the two DEMs.

Purpose and use:

LandScript to import Baranja Hill data and convert to LandSerf format. This script is also included in Chapter 14 of Hengl and Reuter (2008) [the Geomorphometry book]

Programming environment:

Landserf

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
importData.lsc_.zip [5]	1.07 KB

Average:

Your rating: None

Digital soil mapping for Baranja Hill

Short title:

BARSOIL

Inputs: baranja.txt - 59 field profile observations; 15 raster maps (predictors);
Outputs: Predictions and simulations of SOLUM and GLEY_P.

Purpose and use:

Script to interpolate soil thickness (SOLUM) and occurrence of gley horizon (GLEY_P) using land-surface parameters and mapping units. Prepared for the needs of book 'Geomorphometry: concepts, software, applications'.

Programming environment:

R / S language

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
BARSOIL.zip [6]	798.73 KB

Average:

Your rating: None

Flow magnitude and drainage basins

Short title:

flowMag.lsc

Inputs: A DEM (the Baranja Hill 25m DEM is used as an example).
Outputs: A raster containing flow magnitude values for the DEM and a raster containing drainage basins.

Purpose and use:

Script to calculate flow magnitude and drainage basins from a DEM. It demonstrates how recursive function calls in LandScript can be used to calculate zonal geomorphometric parameters using map algebra. To use this script, start LandSerf 2.3 or above and then open the LandScript editor (menu: Edit->LandSctipt Editor).
Part of this script appears in Chapter 14 of Hengl and Reuter (2008) [the Geomorphometry book]. This script runs slowly so is better as a demonstration of how recursive function calls can be made in LandScript rather than a production quality flow magnitude calculator.

Programming environment:

Landserf

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
flowMag.lsc_.zip [7]	1.29 KB

Average:

Your rating: None

GDEM assessment

Short title:

GDEMerr

A methodological framework for assessment of accuracy of a DEM product is described using four case studies (Booschord in the Netherlands, Calabria in Italy, Fishcamp in USA and Zlatibor in Serbia). Focus is put on evaluating the true accuracy of ASTER GDEM using LiDAR data aggregated to 30 m resolution. Three aspects of accuracy have been evaluated: (a) absolute accuracy of elevations (goodness of fit between true and GDEM elevations), (2) accuracy of stream networks (goodness of fit for buffer distance maps for stream networks), and (3) accuracy of surface roughness parameters (goodness of representation of nugget variation and residual errors).

Purpose and use:

Procedures are explained how to assess:

Accuracy of absolute elevations (absolute error);
Positional and attribute accuracy of hydrological features (streams, watersheds, landforms etc);
Accuracy of surface roughness (i.e. representation fo the short-range variation);

The script is attached to the paper indicated below.

Programming environment:

R / S language

Status of work:

Public Domain

Reference:

How accurate and usable is GDEM? A statistical assessment of GDEM using LiDAR data [8]

Data set name:

Fishcamp [9]

Attachment	Size
GDEM_assessment.R [10]	9.48 KB

Average:

Your rating: None Average: 4 (1 vote)

Geomorphological mapping

Short title:

geomorph

For a complete description of the processing steps, see the original publication.

Purpose and use:

Automated extraction of geomorphological features using digital elevation data: case study Drente; outputs: extracted classes and summary statistics; various plots and images

Programming environment:

R / S language

Status of work:

Public Domain

Reference:

{Semi-automated identification and extraction of geomorphological features using digital elevation data} [11]

Data set name:

Boschoord case study [12]

Attachment	Size
maps_KML.zip [13]	84.14 KB
Boschoord.R [14]	20.58 KB

Average:

Your rating: None

Geostatistical simulations of topography

Short title:

DEMsim

Inputs: control.txt - 1020 precise measurements (photogrametric + spot heights); elevations.txt - 2051 points (contours + spot heights); dem10m_tin.asc - 100x150 pixels 30m DEM.
Outputs: simulated DEMs, simulated error surfaces, error assessment statistics.

Purpose and use:

Script to generate and simulate DEMs and assess the error of the height measurements. Prepared for the needs of a research paper 'Geostatistical modelling of topography using auxiliary maps'. Please consider testing the script before you use it with large datasets.

Programming environment:

R / S language

Status of work:

Public Domain

Reference:

{Geostatistical modelling of topography using auxiliary maps} [15]

Data set name:

Zlatibor [16]

Attachment	Size
zlatibor.R [17]	8.26 KB

Average:

Your rating: None

Landform classes (Pennock and Corre, 2001)

Short title:

landform

Inputs: Digital Elevation Model, classification parameters.
Outputs: Map showing landform classes according to Pennock and Corre (2001).

Purpose and use:

Landform Classification based on method of Pennock and Corre (2001). Classify landform units based on relief parameters provided using the topo.aml. Original Source are papers by Pennock et al., Rewritten in aml as closely as possible. Requirements: topo.aml, killgrids.aml, logoff.aml

Programming environment:

Arc AML

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Average:

Your rating: None

Landslide probabilities

Short title:

landprob

Inputs: DEM in metric projection, Source Landslide areas, H/L value as stopping parameter.
Outputs: pq_limi and all other parameters as described in Gruber et al., 2008 article in the Geomorphometry book.

Purpose and use:

Computes Landslide probabilities --- generate landslide probabilities areas using a DEM according to Gruber et al 2008.
&r landslide1 <dem> <sourcearea_dem> {h/l stopping value}

Programming environment:

Arc AML

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
landslide_probabilities.zip [18]	841 bytes

Average:

Your rating: None

Openess

Short title:

openness.lsc

Inputs: A DEM (uses the Baranja Hill DEM25m.asc as an example).
Outputs: A raster of 'openness' that describes how exposed a surface is at each location. The measure is an angle where values close to 0 are very sheltered and 180 is an exposed plane.

Purpose and use:

Script to calculate 'openness' as defined by Yokoyama, Shirasawa and Pike (2002). It is designed to demonstrate how LandScript can be used to calculate new geomorphometric parameters using map algebra. To use this script, start LandSerf 2.3 or above and then open the LandScript editor (menu: Edit->LandSctipt Editor).

Programming environment:

Landserf

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
openness.lsc_.zip [19]	1.02 KB

Average:

Your rating: None

TOPO parameters

Short title:

TOPO.AML

Inputs: INPUTDEM - a grid representing a continuous surface.
Outputs: a set of geomorphometric parameters.
&r topo INPUTDEM {streamflow threshold} {streamcover} (i.e. &r topo ls10 150)

Purpose and use:

It's a compilation of grid commands to calculate some geomorphometric parameters in one go. To avoid mistyping, different approaches between different projects, time saving, etc.
More description in topomanual.pdf included in file.

Programming environment:

Arc AML

Status of work:

Public Domain

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Data set name:

Baranja hill [3]

Attachment	Size
TOPO.AML_.zip [20]	289.98 KB

Average:

Your rating: None

Uncertainty of stream networks derived from elevation data

Short title:

streams_error

A Sextante implementation of these algorithms called FlowTools can be obtained here [21] (contributed by Daniel Nüst [22]). The project report can be downloaded as a PDF [23]. You can access the Eclipse [24] project "flowTools" in the following public Subversion [25] repository:

http://svn.xp-dev.com/svn/FlowTools/

You can also download an initial version (might be outdated!) of the project folder as a zip file [26].

Purpose and use:

Extraction of stream networks from a DEM using error propagation technique.

Programming environment:

R / S language

Status of work:

Public Domain

Reference:

On the uncertainty of stream networks derived from elevation data: the error propagation approach [27]

Data set name:

Baranja hill [3]

Attachment	Size
streams_error.R [28]	12.53 KB
stream_sims.zip [29]	1.2 MB

Average:

Your rating: None Average: 3 (1 vote)

geomorphometry terrain toolbox

Short title:

amltools

ArcInfo Arc Marco Language Code (AML) Geomorphometry Toolbox contains a collection of terrain analysis scripts which can be used with a installation of ArcInfo's Workstation. Please note: They are not working with ARCVIEW and also not directly with ArcGIS, even if it is possible to run AML Code from inside ArcGIS. Still, you need ArcInfo Work Station installed. For a quick start please refer to readme.txt for some command calls or to the outdated topomanual2.pdf.

The ArcGIS Geomorphomerty Python toolbox (DEMO Version) works up to a cell limit of 5000 cells. It contains a significant larger amount of functionality. As of Version 1.0.6 you do not need ArcInfo WS to be installed anymore if AG10 is installed. I rewrote IFTHENELSE and DOCELL in python. Previous Versions(<AG9.4) will still need to install ArcInfo Workstation on your machine as ESRI programmers have not been able to port certain GRID functions as DOCELL. The gtb_demo zip contains now three different zip files for the respective Python Version (2.4-2.6) you might be using. If you get the magic number error, the installed tollbox does not match your installed python version. Current Version is 1.0.6.

For further information and FULL version please refer to www.ai-relief.org [30].

Programming environment:

Arc AML

Status of work:

Shareware

Reference:

Geomorphometry: Concepts, Software, Applications [2]

Attachment	Size
geomorphometry_aml.zip [31]	1.18 MB
gtb_demoV1.0.6.zip [32]	392.81 KB

Average:

Your rating: None Average: 4 (1 vote)