Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It only takes a minute to sign up.
Sign up to join this community
I'm trying to make a buffer around each polygon in my data. At some point I get:
Error in RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td, unaryUnion_if_byid_false, : TopologyException: Input geom 1 is invalid: Ring Self-intersection at or near point
I tried to repair the shapefile with
land1 <- gBuffer(land1, byid=TRUE, width=0)
but didn't work.
How can I find the polygons that are problematic?
I tried which(!gIsValid(land1, byid=TRUE)) but didn't get the right answer.
I've got a draft solution for this problem but it relies on a dev package only available on Github (rangl).
A recent use of it is here: https://github.com/mtennekes/tmap/issues/96#issuecomment-271091194
The key is to triangulate the polygon (with constraints imposed by the polygon boundary on the mesh) and then find ambiguous regions by point-in-poly tests on the triangle centroids. The ambiguous regions are inside overlaps where the rings wind in an invalid way. The approach is described here, but the implementation is all work in progress etc.
https://github.com/r-gris/toposhop/issues/2#issuecomment-264831912
I don't know any other way in R, though rmapshaper and/or polyclip might work.
