library(terra)
library(ggplot2)
v <- vect(system.file("extdata/cyl.gpkg", package = "tidyterra"))
# Grouped
gr_v <- v |>
mutate(start_with_s = startsWith(name, "S")) |>
group_by(start_with_s)
# Dissolve geometries.
diss <- gr_v |>
summarise(n = dplyr::n(), mean = mean(as.double(cpro)))
diss
autoplot(diss, aes(fill = start_with_s)) +
ggplot2::labs(title = "Dissolved")
# Keep geometries separate.
no_diss <- gr_v |>
summarise(n = dplyr::n(), mean = mean(as.double(cpro)), .dissolve = FALSE)
# Same statistic.
no_diss
autoplot(no_diss, aes(fill = start_with_s)) +
ggplot2::labs(title = "Not Dissolved")
Run the code above in your browser using DataLab