Planar partition from disparate polygon inputs. Overlaps aggregate to n
.
Value
Set of sf-format polygons with additional column, n
, denoting
number of overlaps contributing to each of the resultant polygons.
Examples
g <- sf::st_sfc (list (
sf::st_point (cbind (0, 0)),
sf::st_point (cbind (0, 1)),
sf::st_point (cbind (1, 0))
))
pts <- sf::st_sf (a = 1:3, geometry = g)
overlapping_polys <- sf::st_buffer (pts, 0.75)
## decompose and count space-filling from overlapping polygons
x <- ms_aggregate_polys (overlapping_polys)
plot (x)
if (FALSE) {
library (ggplot2)
ggplot (x) +
geom_sf () +
facet_wrap (~n)
}
library (sf)
#> Linking to GEOS 3.10.2, GDAL 3.4.1, PROJ 8.2.1; sf_use_s2() is TRUE
set.seed (6)
pts <- expand.grid (x = 1:8, y = 1:10) %>% st_as_sf (coords = c ("x", "y"))
xsf <- sf::st_buffer (pts, runif (nrow (pts), 0.2, 1.5))
if (FALSE) {
out <- ms_aggregate_polys (xsf)
}