contributiontopeak: Calculates the contribution of each cell to peak flow

In ilyamaclean/ecohydrotools: hydrogical modelling for ecological applications with R

Description

contributiontopeak calculates the contribution of each cell of a digital elevation dataset to peak flow

Usage

contributiontopeak(ft, basins, bysize = TRUE)

Arguments

ft	a raster object, two-dimensional array or matrix of flow times as returned by flowtimes()
basins	a raster object, two-dimensional array or matrix with basins numbered as integers as returned by basindelin().
bysize	an optional logical indicating whether or not to weight values by basin size

Details

contribution to peak flow is calculated by fitting a log-normal distribution to flow times, estimating the probability density function (pdf) of this distribution, weighted by the maximum value of of the pdf. If bysize is TRUE, output values are multiplied by the size of the basin thereby making the assumption that larger basins have higher peak flow. Alternatively real peak flows can be estimated from meterological data.

Value

a raster object of the contribution of each cell to peak flow

Examples

basins <- basindelin(dtm100m) # Takes ~20 seconds to run
ft <- flowtimes(dtm100m)
plot(contributiontopeak(ft, basins))
plot(contributiontopeak(ft, basins, bysize = F))

ilyamaclean/ecohydrotools documentation built on June 10, 2019, 5:45 a.m.