在R中提取具有多边形的点 [英] Extracting points with polygon in R
本文介绍了在R中提取具有多边形的点的处理方法,对大家解决问题具有一定的参考价值,需要的朋友们下面随着小编来一起学习吧!
问题描述
我正在尝试使用'sp'包函数'over'通过多边形提取点
I'm trying to extract points by a polygon using the 'sp' package function '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)),]
不幸的是,我收到一个错误消息
Unfortunately, I'm getting an ERROR
Error in d[!is.na(over(d, plgn)), ] :
matrix argument not supported in SpatialPointsDataFrame selection
有什么主意吗?是因为我的多边形包含100个要素吗?
Any idea?? Is it because my polygon contains 100s of features?
推荐答案
您的plgn
是SpatialPolygonsDataFrame
,因此,is.na(over(d, plgn))
返回逻辑矩阵.这不能用于子集SpatialPoints*
.您可以执行以下操作将逻辑矩阵转换为子设置操作可以容纳的向量:
Your plgn
is a SpatialPolygonsDataFrame
, and as such, is.na(over(d, plgn))
returns a logical matrix. This cannot be used to subset your SpatialPoints*
. You can do the following to convert the logical matrix to a vector that the subsetting operation can accommodate:
d[complete.cases(over(d, plgn)), ]
这篇关于在R中提取具有多边形的点的文章就介绍到这了,希望我们推荐的答案对大家有所帮助,也希望大家多多支持IT屋!
查看全文