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