# NOT RUN {
# obtain data
library(sf)
x = st_read(system.file("gpkg/nc.gpkg", package="sf"), quiet = TRUE)[1, ] %>%
st_geometry %>%
st_transform(3395)
# snap to a grid of 5000 m
y = st_snap_to_grid(x, 5000)
# plot data for visual comparison
par(mfrow = c(1, 2))
plot(x, main = "orginal data")
plot(y, main = "snapped to 5000 m")
# }
Run the code above in your browser using DataLab