line2network: Create a River Network Object from a Shapefile

In mbtyers/riverdist: River Network Distance Computation and Applications

Description

Uses readOGR in package 'rgdal' to read a river shapefile, and establishes connectivity of segment endpoints based on spatial proximity.

Usage

line2network(sp = NA, path = ".", layer = NA, tolerance = 100,
  reproject = NULL, supplyprojection = NULL)

Arguments

sp	SpatialLinesDataFrame object. optional.
path	File path, default is the current working directory.
layer	Name of the shapefile, without the .shp extension.
tolerance	Snapping tolerance of segment endpoints to determine connectivity. Default is 100, therefore care should be exercised when working with larger units of distance, such as km.
reproject	A valid Proj.4 projection string, if the shapefile is to be re-projected. Re-projection is done using spTransform in package 'sp'.
supplyprojection	A valid Proj.4 projection string, if the input shapefile does not have the projection information attached.

Value

Returns an object of class "rivernetwork" containing all spatial and topological information. See rivernetwork-class.

Note

Since distance can only be calculated using projected coordinates, line2network() will generate an error if a non-projected input shapefile is detected. To resolve this, the shapefile can be re-projected in a GIS environment, or using reproject=, shown in the second example below.

Author(s)

Matt Tyers, Joseph Stachelek

Examples

filepath <- system.file("extdata", package="riverdist")

Gulk_UTM5 <- line2network(path=filepath, layer="Gulk_UTM5")
plot(Gulk_UTM5)

## Reading directly from a SpatialLinesDataFrame object

sp <- rgdal::readOGR(dsn = filepath, layer = "Gulk_UTM5", verbose = FALSE)
Gulk_UTM5 <- line2network(sp)
plot(Gulk_UTM5)

## Re-projecting in Alaska Albers Equal Area projection:

AKalbers <- "+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 
    +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0"
    
Gulk_AKalbers <- line2network(path=filepath, layer="Gulk_UTM5", reproject=AKalbers)
plot(Gulk_AKalbers)

mbtyers/riverdist documentation built on May 22, 2019, 12:58 p.m.