build_layers: Construct basic landscape layer data for Dynamic TOPMODEL run

In dynatopmodel: Implementation of the Dynamic TOPMODEL Hydrological Model

Description

Given an elevation raster this function will create a basic multi-band raster that can be used to run Dynamic TOPMODEL after applying a suitable discretisation. It comprises the supplied elevations with the addition of upslope contributing area and topographic wetness index (TWI).

Usage

build_layers(dem, fill.sinks = TRUE, deg = 0.1)

Arguments

dem	Elevation raster using a projected coordinate system (e.g UTM) and regular grid spacing. Should have a resolution of a least 30m for the TWI to be meaningful.
fill.sinks	If TRUE (default) then run a sinkfill before calculating the upslope area and TWI.
deg	Threshold intercell slope to determine sinks (degrees).

Value

A multi-band raster (stack) comprising, in order, the elevations, upslope area and topographic wetness index values.

Author(s)

Peter Metcalfe

Examples

## Not run: 
require(dynatopmodel)
data("brompton")

# Upslope area and wetness index for Brompton catchment
layers <- build_layers(brompton$dem)

sp::plot(layers, main=c("Elevation AMSL (m)", "Upslope area (log(m^2/m))", "TWI ((log(m^2/m))"))

## End(Not run)

dynatopmodel documentation built on May 1, 2019, 7:32 p.m.