ch_wbt_flow_accumulation {CSHShydRology}

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

# Only proceed if Whitebox executable is installed
library(whitebox)
if (check_whitebox_binary()){
  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)
} else {
  message("Examples not run as Whitebox executable not found")
}