Skip to content
Prev 23210 / 29559 Next

densify geometries

This is a nice idea with gBuffer, but testing on a simple case shows you
need explicit curvature to get any extra verts.

library(rgeos)
library(raster) ## just for the "cheat's" polygon creator
 x <- as(extent(0, 1, 0, 1), "SpatialPolygons")
plot(as(as(gBuffer(x, width = 0.1, quadsegs = 55), "SpatialLines"),
"SpatialPoints"))


I want the extra vertices for curvature that's in the space traversed, not
defined by the segment.

Ultimately we need it to be driven by the right curve (great circle,
loxodrome, etc.) so geosphere is a good start - but still the simplest case
of inserting explicitly seems like an obvious gap here (pardon the pun),
even for GEOS?

Cheers, Mike.


On Fri, 7 Aug 2015 at 15:23 chris english <englishchristophera at gmail.com>
wrote: