river_to_shapefile: Export river network as shapefile

In rivnet: Extract and Analyze Rivers from Elevation Data

Export river network as shapefile

Description

Export river network as shapefile. Reach attributes can be added.

Usage

river_to_shapefile(river, filename, atts = NULL, 
                    EPSG = NULL, ...)

Arguments

river	river object. It must have been aggregated (via aggregate_river).
filename	Character. Output file name. It must contain the ".shp" extension.
atts	Attributes at AG level that can be exported. This should be a character vector, with entries corresponding to the names of fields within the AG sub-list of river. See example.
EPSG	EPSG code. Default is NULL, which does not produce a .prj file (i.e., the shapefile does not contain projection information).
...	Additional arguments to be passed to writeVector (e.g., overwrite = TRUE allows overwriting an existing shapefile).

Value

No output is produced. This function is used for its side effetcs.

See Also

terra::writeVector

Examples

library(terra) # to use "vect"
fp <- system.file("extdata/wigger.tif", package="rivnet")
river <- extract_river(outlet=c(637478,237413), DEM=fp)
river <- aggregate_river(river)

tmpname <- paste0(tempfile(), ".shp")
river_to_shapefile(river, tmpname, overwrite = TRUE)

# read output
vv <- vect(tmpname)
vv
plot(vv)


# add attributes to shapefile (drainage area, stream order)
river_to_shapefile(river, tmpname,
					atts = c("A", "streamOrder"), # same names as in river$AG
					overwrite = TRUE) 
vv <- vect(tmpname)
vv

# add projection 
river_to_shapefile(river, tmpname, 
					EPSG = 21781, 
					overwrite = TRUE)
vv <- vect(tmpname)
vv					

rivnet documentation built on April 4, 2025, 1:25 a.m.