我试图使用over函数(rgeos)在两个多边形("grd"和"range")之间进行空间叠加.
此代码适用于所有文件但只有一个:
range.grd.overlay <- grd[!is.na(over(grd, range)), ]@data$grdID # get gridcells that overlap with range
并收到此错误消息:
Error in createPolygonsComment(p) :
rgeos_PolyCreateComment: orphaned hole, cannot find containing polygon for hole at index 9
我知道shapefile中有一个没有周围多边形的洞,但我无法弄清楚如何解决它.我可以简单地删除索引9处的洞吗？怎么样 ？如果您有一个解决方案来修复ArcGIS中的shapefile,那对我来说也没问题.
非常感谢您的帮助!
此致,卡斯滕