I'm looking to intersect 2 spatial layers, keeping all the non-intersecting features as well.
My first layer is the SA2 from NSW, Australia, which look like
# Join all back
SA2_end <- dplyr::bind_rows(
SA2_diff,
intersection
)
ggplot(SA2_end)
geom_sf(aes(fill = koala))