如何使用 R 获取距离点最近的关联多边形的信息？

Question

我正在研究如何在形状文件中的点和多边形之间进行交集（空间连接）。我的想法是获得最近的点和多边形内完全匹配的点。在 ARGIS 中，有一个名为 CLOSEST 的匹配选项函数，它们的定义如下：“连接要素中最接近目标要素的要素被匹配。两个或多个联接要素与目标的距离可能相同当这种情况发生时，会随机选择一个连接特征作为匹配特征。”

我有一个将点相交成多边形的函数，它是由 Lyndon Estes 在 r-sig-geo 列表中慷慨贡献的，当所有多边形都填充了所有区域时，代码运行得很好。第二种情况称为空间连接距离，在 ArcGIS 中，当 match_option 为 CLOSEST 时，称为 INTERSECT，正如 ArcGIS 所做的那样。因此，当区域没有被所有多边形填充时，您可以修改点与多边形之间的最小距离。

这是第一个 INTERSECT 的数据和函数：

library(rgeos)
library(sp) 
library(maptools)
library(rgdal)
library(sp)
xy.map <- readShapeSpatial("http://www.udec.cl/~jbustosm/points.shp")
manzana.map <- readShapeSpatial("http://www.udec.cl/~jbustosm/manzanas_from.shp" )

IntersectPtWithPoly <- function(x, y) { 
# Extracts values from a SpatialPolygonDataFrame with SpatialPointsDataFrame, and appends table (similar to 
# ArcGIS intersect)
# Args: 
#   x: SpatialPoints*Frame
#   y: SpatialPolygonsDataFrame
# Returns:
# SpatialPointsDataFrame with appended table of polygon attributes

  # Set up overlay with new column of join IDs in x
  z <- overlay(y, x)

  # Bind captured data to points dataframe
  x2 <- cbind(x, z)

  # Make it back into a SpatialPointsDataFrame 
  # Account for different coordinate variable names 
  if(("coords.x1" %in% colnames(x2)) & ("coords.x2" %in% colnames(x2))) {
    coordinates(x2) <- ~coords.x1 + coords.x2  
  } else if(("x" %in% colnames(x2)) & ("x" %in% colnames(x2))) {
    coordinates(x2) <- ~x + y 
  }

  # Reassign its projection if it has one
  if(is.na(CRSargs(x@proj4string)) == "FALSE") {
    x2@proj4string <- x@proj4string  
  }
  return(x2)
}


test<-IntersectPtWithPoly (xy.map,manzana.map)

与 Lyndon 分享一些想法，他告诉我这一点：

---

我认为最简单的方法是在每个点周围放置一个缓冲区（如果在，则可以指定 50 m）投影坐标），将它们转换为多边形，然后您的任务就变成两个不同多边形对象的交集。

我还没有在 R 中完成过这种类型的操作，但我怀疑您可以使用以下函数找到答案：

library(sp)
?over

library(rgeos)
?gBuffer
?gIntersects

我建议放置数据的子集来说明问题，然后也许其他人对多边形有更好的想法多边形相交/叠加可以建议该方法。

应该在形状文件中的点半径中进行设置，以便使它们进入最近的多边形。

我知道这个功能可以帮助实现它。

library(sp)
?over

library(rgeos)
?gBuffer
?gIntersects

---

我正在努力，因此任何评论或帮助，将不胜感激！

Answer 1

我知道可以使用 sp 和 rgeos 进行多边形到多边形的叠加。加载 sp 后，您需要加载 rgeos。

library(rgeos)
over(polygon1, polygon2)