Gabungkan data titik spasial dengan data poligon banyak menggunakan R

Saya ingin bergabung spasial dengan SpatialPolygonsDataFrame dan SpatialPointsDataFrame.

dapatkan hanya satu data titik dalam banyak poligon, di bawah kode.

sp::over(meuse, srdf)

Beberapa titik data menjadi beberapa Poligon. Saya ingin mendapatkan data poin, menyimpan semua data poligon.

misalnya: titik A dalam poligon a dan b.

sekarang bisa mendapatkan:

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

Saya mendapatkan seperti ini:

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

Gunakan data:

    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)))

The st_joinfungsi dalam sfpaket ( versi pengembangan ) menyediakan sederhana dan solusi intuitif. Ini mengembalikan spasial sebagai data.frame, dan Anda dapat dengan mudah mengkonversi kembali ke 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")

EDIT:

st_joinsekarang dalam versi CRANsf , seperti yang sebelumnya ditunjukkan oleh SymbolixAU.

Menggunakan 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’.

Titik pertama melapisi satu poligon, titik ke 135 melapisi dua poligon:

> 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

Catatan jika Anda hanya ingin indeks poligon yang dilapis dan tidak semua atribut, cukup konversikan poligon ke SpatialPolygons dan Anda mendapatkan kembali daftar nomor indeks daripada bingkai data dengan semua atribut di:

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

Saya bisa mendapatkan data.frame di bawah kode:

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

Ini .idkolom yang ditambahkan secara otomatis .
Jika mendapatkan data meuse yang difilter atau diedit, lakukan di bawah kode.

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