在R中相交多边形，在Spatstat中导入多边形以进行空间聚类分析

Question

我正在尝试使用readShapePoly命令将一个多边形导入到R中，并使用as(filename,"owin")将该多边形转换为边界。我一直收到这样的响应:多边形包含重复的顶点和自交的多边形：

Polygon 1 contains duplicated vertices
Polygon 1 is self-intersecting
2, 3, 4, 5, 6, 7, 8, 9, 10.
done.
Polygon 1 contains duplicated vertices
Polygon 1 is self-intersecting
Checking for cross-intersection between 10 polygons...1, 2, 3, 4, 5, 6, 7, 8, 9.
done.
Error in owin(poly = opls) : 
  Polygon data contain duplicated vertices and self-intersection

我试着用ArcGIS中的修复几何体功能清理多边形，但没有帮助。我能做什么?

Answer 1

更新:您可以使用零宽度缓冲区以自动方式修复大多数此类问题。使用风险自负，如果不对所涉及的数据进行调查，就无法知道这是否是明智的方式。

例如：

 gIsValid(gBuffer(p4, width = 0, byid = TRUE))
 [1] TRUE

最初的答案如下：

结束更新

您可能需要更好的工具来修复几何图形，或者自己手动调查问题。rgeos包可以用来描述现有的问题，您可以直观地识别它们出现的位置-然后直接在R中编辑它们，或者返回到您的地理信息系统并在那里进行编辑。下面是gIsValid函数帮助中的一个简单示例：

library(rgeos)
p4 = readWKT("POLYGON ((0 40, 0 0, 40 40, 40 0, 0 40))")
gIsValid(p4, reason = TRUE)
[1] "Self-intersection[20 20]"

更改reason返回值以获得更具描述性的警告：

gIsValid(p4)
[1] FALSE
Warning message:
In RGEOSUnaryPredFunc(spgeom, byid, "rgeos_isvalid") :
  Self-intersection at or near point 20 20

因此，现在绘制多边形并突出显示攻击点：

plot(p4, col = "grey")
points(20, 20, pch = 3, cex = 3, col = "red")
axis(1);axis(2);box()

​

​

由于您是从shapefile读取数据，因此可以使用byid参数来测试哪些多边形需要调查(尽管在本例中它们都很好)：

library(maptools)
xx <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],
  IDvar="FIPSNO", proj4string=CRS("+proj=longlat +ellps=clrk66"))
gIsValid(xx, byid = TRUE)
37001 37003 37005 37007 37009 37011 
 TRUE  TRUE  TRUE  TRUE  TRUE  TRUE
... 

等。