2017-03-01 75 views
0

我绘制了R中的地图,对应于我拥有的纬度长值,一些纬度长值落在特定区域(比如印度) ,有没有办法知道哪些纬度长的值是那些?在R中绘制地图(轨迹值落在特定区域之外)

+0

欢迎来到SO。请将鼠标悬停在R标签上 - 它要求提供一个最小可重现的例子,就像我在答案中提供的那样。 [这是一个指南](http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example#answer-5963610);也请看R帮助文件中的示例部分。一个很好的问题通常提供最少的输入数据,所需的输出数据,代码尝试包括所需的软件包 - 全新的/干净的R会话中的所有复制 - 粘贴 - 运行。 *为什么?*它让所有人都可以更轻松地跟随和参与,而无需猜测。 – lukeA

回答

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) 

enter image description here

(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%")所有的变化。

+0

谢谢@lukeA ...它会做需要的。 –