在R中绘制地图（轨道值落在特定区域之外）

Question

我在R中绘制了一张与我所拥有的纬度值相对应的地图，一些纬度较长的值落在特定区域之外（让我们说印度），有没有办法知道哪个lat长值是那些？

Answer 1

这是一个例子

library(sp)
spoly = SpatialPolygons(list(Polygons(list( Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2))) ), "s1")), 1L)
spoints <- SpatialPoints(matrix(c(1:4, 2:5), ncol=2))
plot(spoly)
plot(spoints, add=T)
axis(1);axis(2)

(ov <- spoints %over% spoly)
#  1  2  3  4 
# NA  1  1  1 
coordinates(spoints)[is.na(ov),]
# coords.x1 coords.x2 
#         1         2 

有关所有变体，请参阅help("%over%")。