st_network_blend  R Documentation 
Blending a point into a network is the combined process of first snapping the given point to its nearest point on its nearest edge in the network, subsequently splitting that edge at the location of the snapped point, and finally adding the snapped point as node to the network. If the location of the snapped point is already a node in the network, the attributes of the point (if any) will be joined to that node.
st_network_blend(x, y, tolerance = Inf)
x 
An object of class 
y 
The spatial features to be blended, either as object of class

tolerance 
The tolerance distance to be used. Only features that are
at least as close to the network as the tolerance distance will be blended.
Should be a nonnegative number preferably given as an object of class

There are two important details to be aware of. Firstly: when the
snap locations of multiple points are equal, only the first of these points
is blended into the network. By arranging y
before blending you can
influence which (type of) point is given priority in such cases.
Secondly: when the snap location of a point intersects with multiple edges,
it is only blended into the first of these edges. You might want to run the
to_spatial_subdivision
morpher after blending, such that
intersecting but unconnected edges get connected.
The blended network as an object of class sfnetwork
.
Due to internal rounding of rational numbers, it may occur that the intersection point between a line and a point is not evaluated as actually intersecting that line by the designated algorithm. Instead, the intersection point lies a tinybit away from the edge. Therefore, it is recommended to set the tolerance to a very small number (for example 1e5) even if you only want to blend points that intersect the line.
library(sf, quietly = TRUE) # Create a network and a set of points to blend. n11 = st_point(c(0,0)) n12 = st_point(c(1,1)) e1 = st_sfc(st_linestring(c(n11, n12)), crs = 3857) n21 = n12 n22 = st_point(c(0,2)) e2 = st_sfc(st_linestring(c(n21, n22)), crs = 3857) n31 = n22 n32 = st_point(c(1,1)) e3 = st_sfc(st_linestring(c(n31, n32)), crs = 3857) net = as_sfnetwork(c(e1,e2,e3)) pts = net %>% st_bbox() %>% st_as_sfc() %>% st_sample(10, type = "random") %>% st_set_crs(3857) %>% st_cast('POINT') # Blend points into the network. # > By default tolerance is set to Inf # > Meaning that all points get blended b1 = st_network_blend(net, pts) b1 # Blend points with a tolerance. tol = units::set_units(0.2, "m") b2 = st_network_blend(net, pts, tolerance = tol) b2 ## Plot results. # Initial network and points. oldpar = par(no.readonly = TRUE) par(mar = c(1,1,1,1), mfrow = c(1,3)) plot(net, cex = 2, main = "Network + set of points") plot(pts, cex = 2, col = "red", pch = 20, add = TRUE) # Blend with no tolerance plot(b1, cex = 2, main = "Blend with tolerance = Inf") plot(pts, cex = 2, col = "red", pch = 20, add = TRUE) # Blend with tolerance. within = st_is_within_distance(pts, st_geometry(net, "edges"), tol) pts_within = pts[lengths(within) > 0] plot(b2, cex = 2, main = "Blend with tolerance = 0.2 m") plot(pts, cex = 2, col = "grey", pch = 20, add = TRUE) plot(pts_within, cex = 2, col = "red", pch = 20, add = TRUE) par(oldpar)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.