(This is a test post)
library(ggplot2)
nc <- sf::st_read(system.file("shape/nc.shp", package = "sf"), quiet = TRUE)
ggplot(nc) +
geom_sf(aes(fill = AREA))
## although coordinates are longitude/latitude, st_intersection assumes that they are planar
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.