# NOT RUN {
library(ggplot2)
library(sf)
# Convert stems to sf object
census_1_ex_sf <- census_1_ex %>%
st_as_sf(coords = c("gx", "gy"))
# Plot stems with plot boundary
ggplot() +
geom_sf(data = study_region_ex) +
geom_sf(data = study_region_bw, color = "red", fill = "transparent")
# }
Run the code above in your browser using DataLab