2017-03-01 77 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 ...它會做需要的。 –