gpt4 book ai didi

r - R中的相交多边形,在Spatstat中导入多边形进行空间聚类分析

转载 作者:行者123 更新时间:2023-12-04 05:47:40 27 4
gpt4 key购买 nike

我正在尝试使用 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 中的修复几何函数清理多边形,但没有帮助。我能做什么?

最佳答案

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

例如:

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

原答案如下:

结束更新

您可能需要更好的工具来修复几何体,或者自己手动调查问题。 rgeos包可用于描述现有问题,您可以直观地识别它们发生的位置 - 然后直接在 R 中编辑它们或返回到您的 GIS 并在那里进行。这是 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()

Self-intersection polygon

由于您正在从 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
...

关于r - R中的相交多边形,在Spatstat中导入多边形进行空间聚类分析,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/10474299/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com