## Load districts of Germany
load(system.file("shapes", "districtsD.RData", package = "surveillance"))
plot(districtsD, border = "gray", asp = 1)
## Union these districts using either "sf" or "polyclip"
if (requireNamespace("sf")) {
stateD <- unionSpatialPolygons(districtsD, method = "sf")
plot(stateD, add = TRUE, border = 2, lwd = 2)
}
if (requireNamespace("polyclip")) {
stateD_pc <- unionSpatialPolygons(districtsD, method = "polyclip")
plot(stateD_pc, add = TRUE, border = 1, lwd = 2, lty = 2)
}
Run the code above in your browser using DataLab