在R中用多边形提取点答案

【问题标题】：Extracting points with polygon in R在R中用多边形提取点
【发布时间】：2014-03-01 06:44:11
【问题描述】：

我正在尝试使用“sp”包函数“over”通过多边形提取点

library(sp)
library(rgeos)
#my polygon plgn (many polygon features in one)
plot(plgn)
proj4string(plgn) = CRS("+proj=utm +zone=46 +datum=WGS84 +units=m +no_defs")
#giving spatial reference to point data d
coordinates(d) <- ~X+Y
proj4string(d) = CRS("+proj=utm +zone=46 +datum=WGS84 +units=m +no_defs")
#USE overlay (there are many NAs)
overlay=d[!is.na(over(d, plgn)),]

很遗憾，我遇到了一个错误

Error in d[!is.na(over(d, plgn)), ] : 
matrix argument not supported in SpatialPointsDataFrame selection

有什么想法吗？？是因为我的多边形包含 100 多个特征吗？

【问题讨论】：

您是否尝试过从内向外运行该生产线？例如，检查 over(d, plgn) 返回的内容，然后查看 is.na(over(d, plgn))。他们的结果有意义吗？

标签： r extract polygon geospatial

【解决方案1】：

您的plgn 是SpatialPolygonsDataFrame，因此，is.na(over(d, plgn)) 返回一个逻辑矩阵。这不能用于子集您的SpatialPoints*。您可以执行以下操作将逻辑矩阵转换为子集操作可以容纳的向量：

d[complete.cases(over(d, plgn)), ]

【讨论】：