0
我绘制了R中的地图,对应于我拥有的纬度长值,一些纬度长值落在特定区域(比如印度) ,有没有办法知道哪些纬度长的值是那些?在R中绘制地图(轨迹值落在特定区域之外)
我绘制了R中的地图,对应于我拥有的纬度长值,一些纬度长值落在特定区域(比如印度) ,有没有办法知道哪些纬度长的值是那些?在R中绘制地图(轨迹值落在特定区域之外)
下面是一个例子
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%")
所有的变化。
谢谢@lukeA ...它会做需要的。 –
欢迎来到SO。请将鼠标悬停在R标签上 - 它要求提供一个最小可重现的例子,就像我在答案中提供的那样。 [这是一个指南](http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example#answer-5963610);也请看R帮助文件中的示例部分。一个很好的问题通常提供最少的输入数据,所需的输出数据,代码尝试包括所需的软件包 - 全新的/干净的R会话中的所有复制 - 粘贴 - 运行。 *为什么?*它让所有人都可以更轻松地跟随和参与,而无需猜测。 – lukeA