我正在尝试使用'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个特征吗?
答案 0 :(得分:4)
您的plgn
是SpatialPolygonsDataFrame
,因此,is.na(over(d, plgn))
会返回逻辑矩阵。这不能用于SpatialPoints*
的子集。您可以执行以下操作将逻辑矩阵转换为子集化操作可以容纳的向量:
d[complete.cases(over(d, plgn)), ]