ch_wbt_flow_accumulation: Creates flow accumulation grid file

In CSHS-hydRology/CSHShydRology: Canadian Hydrological Analyses

Creates flow accumulation grid file

Description

Creates flow accumulation grid file

Usage

ch_wbt_flow_accumulation(fn_dem_ns, fn_flowacc, return_raster = TRUE)

Arguments

fn_dem_ns	File name of dem with sinks removed.
fn_flowacc	File name for flow accumulation grid to be created.
return_raster	If TRUE (the default), the flow accumulation grid will be returned as a raster object, in addition to being written to fn_flowacc. If FALSE, the output file will still be created but a NULL value is returned.

Value

If return_raster = TRUE, the flow accumulation grid will be returned as a raster object, otherwise NULL is returned.

Author(s)

Dan Moore

Examples

# Not tested automatically as requires installation of Whitebox
library(raster)
test_raster <- ch_volcano_raster()
dem_raster_file <- tempfile(fileext = c(".tif"))
no_sink_raster_file <- tempfile("no_sinks", fileext = c(".tif"))

# write test raster to file
writeRaster(test_raster, dem_raster_file, format = "GTiff")

# remove sinks
removed_sinks <- ch_wbt_removesinks(dem_raster_file, no_sink_raster_file, method = "fill")

# get flow accumulations
flow_acc_file <- tempfile("flow_acc", fileext = c(".tif"))
flow_acc <- ch_wbt_flow_accumulation(no_sink_raster_file, flow_acc_file)
plot(flow_acc)

CSHS-hydRology/CSHShydRology documentation built on Aug. 18, 2022, 4:44 p.m.