Подмножество пространственных точек с многоугольником

Question

У меня есть SpatialPolygonsDataFrame (spolydf) и SpatialPointsDataFrame (spointdf). Слои имеют разную протяженность, но перекрываются.

Я могу выбрать точки, которые попадают в многоугольник, используя

fall.within.poly <- spointdf[spolydf,]

Как выбрать точки, которые выходят за пределы многоугольника? пытался

fall.outside.poly <- spointdf[-spolydf,]

но не работает. Я делаю что-то простое - помогите, пожалуйста.

Answer 1

Уже немного поздно, но у меня была та же проблема сегодня, поэтому я решил опубликовать свое решение, используя 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)

Answer 2

Самый простой способ, который я нашел:

      library(splancs)
set.seed(123)
my.shape <- matrix(runif(10), 5)
my.points <- data.frame(x=runif(500), y=runif(500))
my.points$in.shape <- 1:500 %in% inpip(my.points, my.shape)

plot(my.points[1:2], col=1 + my.points$in.shape)
polygon(my.shape)