Wczytanie i zapisz Shapefile z wykorzystanie biblioteki SF

Tym razem pokażemy jak wczytać i zapisać plik warstwy wektorowej z użyciem nowej biblioteki sf.

Kod:

library(sf)

shp = read_sf("d:/GIS_in_R/powiaty.shp")
plot(shp$geometry)

shp
shp[5,]
select = which(shp$jpt_nazwa_ == "powiat pruszkowski")
select
single_polygon = shp[select,]
plot(single_polygon$geometry)
write_sf(single_polygon,"d:/GIS_in_R/single_pol.shp")