Unir dados de pontos espaciais a vários dados poligonais usando R

Quero ingressar espacialmente com SpatialPolygonsDataFrame e SpatialPointsDataFrame.

obtenha apenas dados de um ponto em vários polígonos, sob o código

sp::over(meuse, srdf)

Alguns dados de pontos em vários polígonos. Eu quero obter dados de pontos manter todos os dados de polígono.

por exemplo: ponto A dentro do polígono a e b.

agora pode obter:

| point | polygon |
-------------------
| A     | a       |

Eu fico assim:

| point | polygon |
-------------------
| A     | a       |
| A     | b       |

Use dados:

    coordinates(meuse) = ~x+y
    sr1=Polygons(list(Polygon(cbind(c(180114, 180553, 181127, 181477, 181294, 181007, 180409,
                                      180162, 180114), c(332349, 332057, 332342, 333250, 333558, 333676,
                                                         332618, 332413, 332349)))),'1')
    sr2=Polygons(list(Polygon(cbind(c(180042, 180545, 180553, 180314, 179955, 179142, 179437,
                                      179524, 179979, 180042), c(332373, 332026, 331426, 330889, 330683,
                                                                 331133, 331623, 332152, 332357, 332373)))),'2')
    sr3=Polygons(list(Polygon(cbind(c(179110, 179907, 180433, 180712, 180752, 180329, 179875,
                                      179668, 179572, 179269, 178879, 178600, 178544, 179046, 179110),
                                    c(331086, 330620, 330494, 330265, 330075, 330233, 330336, 330004,
                                      329783, 329665, 329720, 329933, 330478, 331062, 331086)))),'3')
    sr4=Polygons(list(Polygon(cbind(c(180304, 180403,179632,179420,180304),
                                    c(332791, 333204, 333635, 333058, 332791)))),'4')
    sr5=Polygons(list(Polygon(cbind(c(179500, 180000, 180000, 179500),
                                    c(331000, 331000, 331500, 331500)))), '5')
    sr=SpatialPolygons(list(sr1,sr2,sr3,sr4, sr5))
    srdf=SpatialPolygonsDataFrame(sr, data.frame(row.names=c('1','2','3','4', '5'), PIDS=1:5, y=runif(5)))

A st_joinfunção no sfpacote ( versão de desenvolvimento ) fornece uma solução simples e intuitiva. Ele retorna a junção espacial como uma data.framee você pode facilmente converter novamente em Spatial*.

# Convert to sf-objects
srdf.sf <- st_as_sf(srdf)
meuse.sf <- st_as_sf(meuse)

# Keep all "meuse.sf", sort by row.names(meuse.sf). Default overlay is "intersects".
meuse_srdf <- st_join(meuse.sf, srdf.sf) 

# Keeps all "srdf.sf", sort by row.names(srdf.sf)
srdf_meuse <- st_join(srdf.sf, meuse.sf)

# Convert back to Spatial*
meuse_srdf <- as(meuse_srdf, "Spatial")
srdf_meuse <- as(srdf_meuse, "Spatial")

EDITAR:

st_joinestá agora na versão CRAN do sf, como apontado anteriormente pela SymbolixAU.

Usar over(meuse, srdf, returnList=TRUE)

 if ‘returnList’ is TRUE, a list of length
 ‘length(x)’, with list element ‘i’ the vector of all indices of
 the geometries in ‘y’ that correspond to the $i$-th geometry in
 ‘x’.

O primeiro ponto sobrepõe um polígono, o 135º ponto sobrepõe dois:

> rr = over(meuse,srdf,returnList=TRUE)
> rr[[1]]
  PIDS         y
1    1 0.2069365
> rr[[135]]
  PIDS          y
2    2 0.34809708
5    5 0.08942346

Observe que se você deseja apenas os índices de polígonos sobrepostos e não todos os atributos, basta converter os polígonos em SpatialPolygons e obter uma lista dos números de índice em vez dos quadros de dados com todos os atributos em:

> rr = over(meuse,as(srdf,"SpatialPolygons"),returnList=TRUE)
> rr[[1]]
[1] 1
> rr[[135]]
[1] 2 5

Eu posso obter data.frame sob o código:

library(tidyverse)
sp::over(meuse, srdf, returnList=TRUE) %>%
plyr::ldply(., data.frame) -> tmp

É .idcoluna adicionada automaticamente .
Se for filtrado ou editado os dados de uso, siga o código.

use_id <- 
  tmp %>% select(`.id`) %>% unlist(., use.names = F)
rownames(meuse) %in% use_id