R中的多边形交点

我有一个包含很多多边形，点和线的列表。现在我想检查每个多边形之间的交集区域。我知道我可以使用gIntersect，但我的数据是在data.frame中，而不是在spatialpolygon文件中。有没有一个聪明的方法来做到这一点？R中的多边形交点

lon  lat ID 
127.0167 29.14449 1 
127.0148 29.14507 1 
123.8638 17.63341 2 
123.8593 17.62754 2 
123.8568 17.63601 2 
123.8528 17.65023 2 
127.6790 49.01934 3 
127.2719 49.12513 6 
127.0249 49.14633 7 
127.6763 49.02139 8 
127.6710 49.02426 8 
127.6684 49.02668 8 
127.4648 49.07208 13 
127.3757 49.08205 13 
127.4198 30.04310 14 
127.4259 40.04974 14 
127.3136 39.09050 15 
127.3197 29.09516 15 
127.2360 39.16492 16 
127.2099 49.16787 16 
127.2266 29.01800 17 
128.2771 48.36411 18 
128.1930 28.44411 18 
128.1925 18.44530 18 
128.1928 28.44553 18 
128.1932 48.44598 18 
128.1953 38.44774 18 
128.1978 28.44532 18 
125.7947 28.71272 19 
125.7982 28.72078 19 
125.8402 18.74029 19 
125.8572 18.74141 19

这里是列表。它包含原始文件中的更多数据点。具有相同ID的坐标属于同一个多边形，点或线。我需要检查每个多边形与列表中所有其他多边形的交集。

来源

2015-11-21 quallenjäger

有你在这里考虑这个例子：https://stat.ethz.ch/R-manual/R-devel/library/mgcv/html/in.out.html – BigDataScientist

是的，我的问题是这种转换列表到空间数据，所以我可以使用gIntersection。我可以使用多边形（f，ID = ID）来转换它吗？与上面的列表。我不知道如何识别每个矩阵。 –

您可以将data.frame转换为SpatialPolygons。

library(sp) 
polys <- lapply(unique(coords$ID), function(i) { 
    Polygons(list(Polygon(coords[coords$ID==i, 1:2])), ID=i) 
}) 
spa_polys <- SpatialPolygons(polys)

另外，我还发现，gIntersects不会为一个点（一个点的空间多边形）或线（二点空间多边形）的工作。

library(rgeos) 
> gIntersects(spa_polys[13], spa_polys[1]) 
Error in RGEOSBinPredFunc(spgeom1, spgeom2, byid, func) : 
    IllegalArgumentException: Invalid number of points in LinearRing found 3 - must be 0 or >= 4 
> gIntersects(spa_polys[13], spa_polys[2]) 
[1] FALSE 
> gIntersects(spa_polys[13], spa_polys[3]) 
Error in RGEOSBinPredFunc(spgeom1, spgeom2, byid, func) : 
    IllegalArgumentException: point array must contain 0 or >1 elements 
> gIntersects(spa_polys[13], spa_polys[13]) 
[1] TRUE

来源

2015-11-21 16:23:11 kitman0804

谢谢！多数民众赞成我的第二个问题，我如何检查数据是线或点或多边形，我怎么能过滤出我的设置？ –

@quallenjäger我不知道我们如何检查'SpatialPolygons'对象是否是一条线或一个点或多边形，但我们可以在将'coords'从'data.frame'转换为'SpatialPolygons'之前过滤掉它们。 – kitman0804

R中的多边形交点

回答

相关问题