if (FALSE) {
library(sf)
nc = st_read(system.file("shape/nc.shp", package="sf"))
my_grid = gpat_st_make_grid(nc)
my_grid$id = 1:100
grid_centroids = st_centroid(my_grid) %>%
st_coordinates(grid_centroids) %>%
as_data_frame() %>%
mutate(id = 1:100)
ggplot() +
geom_sf(data = my_grid) +
geom_text(data = grid_centroids, aes(x = X, y = Y, label = id)) +
theme_void()
}
Run the code above in your browser using DataLab