我有一个 SpatialPolygonsDataFrame (spolydf) 和一个 SpatialPointsDataFrame (spointdf)。这些层有不同的范围,但重叠。
我可以使用选择落在多边形内的点
fall.within.poly <- spointdf[spolydf,]
如何选择落在多边形之外的点?试过
fall.outside.poly <- spointdf[-spolydf,]
但不起作用。我想念一些简单的东西 - 请帮忙。
最佳答案
有点晚了,但我今天遇到了同样的问题,所以我想我会使用 gDifference()
发布我的解决方案来自 rgeos
包裹:
require(rgeos)
require(sp)
##create spdf
coords=expand.grid(seq(150,151,0.1),seq(-31,-30,0.1))
spdf=data.frame("lng"=coords[,1],"lat"=coords[,2])
coordinates(spdf) = ~lng+lat
proj4string(spdf)<- CRS("+init=epsg:4326")
plot(spdf)
##create poly
poly1 = SpatialPolygons(list(Polygons(list(Polygon(cbind(c(150.45,150.45,150.75,150.75,150.45),c(-30.75,-30.45,-30.45,-30.75,-30.75)))),ID=1)))
proj4string(poly1)<- CRS("+init=epsg:4326")
lines(poly1)
##get difference
out = gDifference(spdf,poly1)
points(out,col="red",pch=16)
关于r - 用多边形子集空间点,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/49929603/