Selekcja przestrzenna z SF

Dzisiaj pokażemy Wam jak używać biblioteki sf do selekcji przestrzennej. Do naszego przykładu użyjemy dobrze znanej już Wam warstwy dróg i linii kolejowych. Przykład znowu w wersji wideo. Miłego oglądania.

Kod:

library(sf)

rail <- read_sf("d:/GIS_in_R/rail.shp")
road <- read_sf("d:/GIS_in_R/roads.shp")

plot(rail$geometry)
plot(road$geometry,col="red",add=T)

int <- st_intersection(rail,road)

plot(int$geometry,col="green",add=T,pch=20)

ind <- st_intersects(rail,road)
ind <- apply(ind, 1, any)
ind

plot(rail$geometry[ind])