ch_wbt_channels: Generate stream network

In CSHS-hydRology/CSHShydRology: Canadian Hydrological Analyses

Generate stream network

Description

Generate stream network

Usage

ch_wbt_channels(
  fn_flowacc,
  fn_flowdir,
  fn_channel_ras,
  fn_channel_vec,
  threshold = NULL,
  ...
)

Arguments

fn_flowacc	File name for flow accumulation grid.
fn_flowdir	File name for flow direction grid.
fn_channel_ras	File name for raster version of channel network.
fn_channel_vec	File name for vector version of channel networks.
threshold	Threshold for channel initiation.
...	Other parameters for whitebox function wbt_extract_streams

Value

Returns a sf vector object of the stream channels.

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)

# get flow directions
flow_dir_file <- tempfile("flow_dir", fileext = c(".tif"))
flow_dir <- ch_wbt_flow_direction(no_sink_raster_file, flow_dir_file)
channel_raster_file <- tempfile("channels", fileext = c(".tif"))
channel_vector_file <- tempfile("channels", fileext = c(".shp"))
channels <- ch_wbt_channels(flow_acc_file, flow_dir_file, channel_raster_file,
channel_vector_file, 1)
plot(channels)

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