compute_iFLS_weights {rdwplus} | R Documentation |
Compute iFLS weights
Description
Compute an iFLO weight raster outside of the compute_metrics()
function.
Usage
compute_iFLS_weights(
streams,
null_streams,
flow_dir,
out_flow_length,
out_iFLS,
out_iFLS_no_stream,
watershed,
idwp,
remove_streams,
...
)
Arguments
streams |
Pour point raster containing a single pour point (i.e., the outlet). |
null_streams |
A streams raster with NoData for the stream cells and 1s everywhere else |
flow_dir |
A flow direction raster. |
out_flow_length |
Name of the output flow length raster. |
out_iFLS |
Name of the output weights raster. |
out_iFLS_no_stream |
Name of the output weights raster excluding cells on the stream line (ignored inf |
watershed |
Watershed raster to use as a mask for the flow-path calculations. This is optional for the iFLS weight calculations. |
idwp |
An inverse distance weighting power. This should be negative. The value |
remove_streams |
A logical indicating whether cells corresponding to the stream line should be removed from the weights raster. Defaults to |
... |
Optional extra arguments to |
Value
Nothing.
Examples
if(check_running()){
# Retrieve paths to data sets
dem <- system.file("extdata", "dem.tif", package = "rdwplus")
lus <- system.file("extdata", "landuse.tif", package = "rdwplus")
sts <- system.file("extdata", "site.shp", package = "rdwplus")
stm <- system.file("extdata", "streams.shp", package = "rdwplus")
# Set environment
set_envir(dem)
# Get other data sets (stream layer, sites, land use, etc.)
raster_to_mapset(lus)
vector_to_mapset(c(stm, sts))
# Reclassify streams
out_stream <- paste0(tempdir(), "/streams.tif")
rasterise_stream("streams", out_stream, TRUE)
reclassify_streams("streams.tif", "streams01.tif", overwrite = TRUE)
# Burn in the streams to the DEM
burn_in("dem.tif", "streams01.tif", "burndem.tif", overwrite = TRUE)
# Fill dem
fill_sinks("burndem.tif", "filldem.tif", "fd1.tif", "sinks.tif", overwrite = TRUE)
# Derive flow direction and accumulation grids
derive_flow("dem.tif", "fd.tif", "fa.tif", overwrite = T)
# Derive a new stream raster from the FA grid
derive_streams("dem.tif", "fa.tif", "new_stm.tif", "new_stm", min_acc = 200, overwrite = T)
# Recode streams
reclassify_streams("new_stm.tif", "null_stm.tif", "none")
# Snap sites to streams and flow accumulation
snap_sites("site", "new_stm.tif", "fa.tif", 2, "snapsite", T)
# Get watersheds
get_watersheds("snapsite", "fd.tif", "wshed.tif", T)
# Get iFLS weights
compute_iFLS_weights(
"new_stm.tif",
"null_stm.tif",
"fd.tif",
"fl_streams.tif",
"iFLS_weights.tif",
idwp = -1,
watershed = "wshed.tif",
remove_streams = FALSE,
overwrite = T
)
plot_GRASS("iFLS_weights.tif", col = topo.colors(12))
}