Creates a new SpatialLinesNetwork (for SpatialLines) or sfNetwork (for sf) object that can be used for routing analysis within R.
SpatialLinesNetwork(sl, uselonglat = FALSE, tolerance = 0)A SpatialLines or SpatialLinesDataFrame containing the lines to use to create the network.
A boolean value indicating if the data should be assumed
to be using WGS84 latitude/longitude coordinates. If FALSE or not
set, uses the coordinate system specified by the SpatialLines object.
A numeric value indicating the tolerance (in the units of the coordinate system) to use as a tolerance with which to match nodes.
This function is used to create a new SpatialLinesNetwork from an existing SpatialLines or SpatialLinesDataFrame object. A typical use case is to represent a transport network for routing and other network analysis functions. This function and the corresponding SpatialLinesNetwork class is an implementation of the SpatialLinesNetwork developed by Edzer Pebesma and presented on RPubs. The original implementation has been rewritten to better support large (i.e., detailed city-size) networks and to provide additional methods useful for conducting transport research following on from the initial examples provided by Janoska(2013).
Pebesma, E. (2013). Spatial Networks, URL:http://rpubs.com/edzer/6767.
Janoska, Z. (2013). Find shortest path in spatial network, URL:http://rpubs.com/janoskaz/10396.
Other rnet: calc_catchment_sum,
calc_catchment,
calc_moving_catchment,
calc_network_catchment,
find_network_nodes, gsection,
islines, lineLabels,
overline2, overline,
plot,SpatialLinesNetwork,ANY-method,
plot,sfNetwork,ANY-method,
sln2points,
sum_network_links,
sum_network_routes
# NOT RUN {
sln <- SpatialLinesNetwork(route_network)
class(sln)
weightfield(sln) # field used to determine shortest path
plot(sln)
points(sln2points(sln)[1, ], cex = 5)
points(sln2points(sln)[50, ], cex = 5)
shortpath <- sum_network_routes(sln, 1, 50, sumvars = "length")
plot(shortpath, col = "red", lwd = 4, add = TRUE)
points(sln2points(sln)[35, ], cex = 5)
shortpath <- sum_network_routes(sln, 1, 35, sumvars = "length")
plot(shortpath, col = "red", lwd = 4, add = TRUE)
library(sf)
sln_sf <- SpatialLinesNetwork(route_network_sf)
plot(sln_sf)
shortpath <- sum_network_routes(sln_sf, 1, 50, sumvars = "length")
plot(shortpath$geometry, col = "red", lwd = 4, add = TRUE)
# }
Run the code above in your browser using DataLab