В R намерете многоъгълника, съдържащ точка - Shapefile

Имам шейпфайл и бих искал да намеря многоъгълника, съдържащ списък с точки. Например,

rio <- readShapeSpatial("setores_rio.shp")  
bairrorio.fort<- fortify(rio , region = "neighborhood")

head(bairrorio.fort)
   long       lat order  hole piece          group           id
1 -43.17769 -22.91814     1 FALSE     1 330455705001.1 330455705001
2 -43.17771 -22.91814     2 FALSE     1 330455705001.1 330455705001
3 -43.17771 -22.91808     3 FALSE     1 330455705001.1 330455705001
4 -43.17793 -22.91811     4 FALSE     1 330455705001.1 330455705001
5 -43.17811 -22.91768     5 FALSE     1 330455705001.1 330455705001
6 -43.17802 -22.91766     6 FALSE     1 330455705001.1 330455705001

Да приемем, че p = c(long, lat) е точка с lat long локализация. Бих искал да намеря id(neighborhood) (вижте bairrorio.fort), съдържащ точка p.


r shapefile
person Vasco    schedule 27.05.2015    source източник
comment
bairrorio.fort[bairrorio.fort$long==long & bairrorio.fort$lat==lat,"id"] ще даде това, което искате   -  person user227710    schedule 28.05.2015
comment
проверете over() в пакета sp.   -  person cory    schedule 28.05.2015
comment
Сигурни ли сте, че точката p= c(long, lat) може да бъде всяка точка? Не работи... Мисля, че вашето предложение може да работи само за точка в bairrorio.fort, въпреки това, когато стартирам с p = c(-43.17769 -22.91814) (вижте първия ред на bairrorio.fort) резултатът е › bairrorio.fort[bairrorio.fort$long == -43.17769 & bairrorio.fort$lat == -22.91814, id] символ (0)   -  person Vasco    schedule 28.05.2015
comment
Не използвайте readShapeSpatial за четене на шейпфайл. Използвайте readOGR от rgdal пакет.   -  person    schedule 28.05.2015

Отговори (1)


arrow_upward
1
arrow_downward

Можете да използвате point.in.polygon от sp:

library(sp)
library(magrittr)

bairrorio.fort %>% 
  split(.$id) %>% 
  sapply(function(x) point.in.polygon(p[1], p[2], x$long, x$lat) > 0) %>%  
  names(.)[.]
person bergant    schedule 27.05.2015