st_check_crossings: Check for crossing roads

In Jean-Romain/MFFProads: Road alignment and measure from airborne LiDAR data for MFFP

Check for crossing roads

Description

Check if some roads are crossing elsewhere than at their ends. These occurences should be fixed as roads crossing each other should be split at their intersection point in a valid topological network.

Usage

st_check_crossings(roads)

Arguments

roads

multiples lines (sf format)

Value

Point object (sf format) corresponding to intersection points of crossing roads.

Examples

library(sf)

network <- system.file("extdata", "j53e_network.gpkg", package="ALSroads")
network <- st_read(network, "invalid_topology")

crossings <- st_check_crossings(network)

par(mfrow=c(1,2))
plot(st_geometry(network), axes = TRUE) # Inspected roads
plot(crossings, col = "blue", add = TRUE) # Crossings to fix
title("Full extent")
plot(st_geometry(network), axes = TRUE, xlim = c(260891, 260892), ylim=c(5250599, 5250599.5))
plot(crossings, col = "blue", add = TRUE)
title("Zoom at crossings")

Jean-Romain/MFFProads documentation built on Nov. 19, 2024, 5:12 p.m.