6
如何在SpatialPointsDataFrame上覆盖多边形并保留SPDF数据?
我有SpatialPointsDataFrame一些其他数据。我想提取多边形内的那些点,同时保留SPDF对象及其对应的数据。 到目前为止,我运气不佳,只能通过一个通用ID进行匹配和合并,但这只能工作,因为我已经使用单个IDS对数据进行了网格划分。 这是一个简单的示例,我正在寻找红场内的点。 library(sp) set.seed(357) pts <- data.frame(x = rnorm(100), y = rnorm(100), var1 = runif(100), var2 = sample(letters, 100, replace = TRUE)) coordinates(pts) <- ~ x + y class(pts) plot(pts) axis(1); axis(2) ply <- matrix(c(-1,-1, 1,-1, 1,1, -1,1, -1,-1), ncol = 2, byrow = TRUE) ply <- SpatialPolygons(list(Polygons(list(Polygon(ply)), ID …