rdwplus (version 0.1.0)

reclassify_streams: Reclassify streams into the format required for the land use metric calculations

Description

Given a streams raster, this function will either create a binary streams raster (0 for non-stream cells and 1 for stream cells) or a unary streams raster (1 for stream cells and NoData for all other cells). Another option is to reclassify the streams raster such that stream cells are given the value NoData and non-stream cells are given the value 1.

Usage

reclassify_streams(stream, out, out_type = "binary", overwrite = FALSE)

Value

Nothing. A file with the name out will be written to the current GRASS mapset. This raster will be in unsigned integer format.

Arguments

stream

Name of a streams raster in the current GRASS mapset. Typically this is the result of rasterise_stream. The raster should have NoData values for all non-stream cells. Stream cells can have any other value.

out

The output file.

out_type

Either 'binary', 'unary', or 'none'. See the Description above.

overwrite

A logical indicating whether the output should be allowed to overwrite any existing files. Defaults to FALSE.

Examples

# Will only run if GRASS is running
if(check_running()){

# Load data set
dem <- system.file("extdata", "dem.tif", package = "rdwplus")
stream_shp <- system.file("extdata", "streams.shp", package = "rdwplus")

# Set environment parameters and import data
set_envir(dem)
vector_to_mapset(vectors =  stream_shp)

# Create binary stream
out_name <- paste0(tempdir(), "/streams_rast.tif")
rasterise_stream("streams", out_name, overwrite = TRUE)
reclassify_streams("streams_rast.tif", "streams_binary.tif", out_type = "binary", overwrite = TRUE)
reclassify_streams("streams_rast.tif", "streams_unary.tif", out_type = "unary", overwrite = TRUE)
reclassify_streams("streams_rast.tif", "streams_none.tif", out_type = "none", overwrite = TRUE)

# Plot
plot_GRASS("streams_rast.tif", col = topo.colors(2), main = "Rasterized Streams")
plot_GRASS("streams_binary.tif", col = topo.colors(2), main = "Binary Streams")
plot_GRASS("streams_unary.tif", col = topo.colors(2), main = "Unary Streams")
plot_GRASS("streams_none.tif", col = topo.colors(2), main = "Null (none) Streams")
}