geoaxe
- split geospatial objects into pieces
Stable CRAN version
install.packages("geoaxe")
Development version
pak::pkg_install("ropensci/geoaxe")
Works for only SpatialPolygons
for now, but aim to include related classes soon.
library("rgeos")
wkt <- "POLYGON((-180 -20, -140 55, 10 0, -140 -60, -180 -20))"
poly <- rgeos::readWKT(wkt)
polys <- chop(x = poly)
Plot original polygon
plot(poly, lwd = 6)
Add chopped up polygon bits
plot(polys, add = TRUE)
plots go left to right, then down, and repeat
layout(matrix(c(1,2,3,4), 2, 2, byrow = TRUE))
par(mar = c(1, 0, 1, 0))
plot(chop(wkt, n = 10))
plot(chop(wkt, n = 15))
plot(chop(wkt, n = 20))
plot(chop(wkt, n = 50))
plots go left to right, then down, and repeat
layout(matrix(1:8, 4, 2, byrow = TRUE))
par(mar = c(1, 0, 1, 0))
plot(chop(wkt, size = 2))
plot(chop(wkt, size = 4))
plot(chop(wkt, size = 8))
plot(chop(wkt, size = 15))
plot(chop(wkt, size = 25))
plot(chop(wkt, size = 50))
plot(chop(wkt, size = 100))
plot(chop(wkt, size = 200))