Aggregate disparate polygons

Planar partition from disparate polygon inputs. Overlaps aggregate to n.

ms_aggregate_polys(px, ...)

Arguments

px	input polygons (assumed overlapping poly/mpolys)
...	unused

Details

Input is a single simple features polygon data frame. No attribute data is considered.

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)
#> Loading required namespace: gibble
#> Loading required namespace: polyclip
plot(x)
# library(ggplot2)
# ggplot(x) + geom_sf() + facet_wrap(~n)

library(sf)
#> Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1
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))
# out <- ms_aggregate_polys (xsf)