Объединение нескольких SpatialPolygonDataFrames в 1 SPDF в R?


22

Я создал 2 полигона в QGIS. Используя их в R, полигоны автоматически становятся SpatialPolygonsDataFrame (SPDF). Я хотел бы объединить их в один SPDF (как это легко сделать в ArcGis, используя инструмент Merge ). Я уверен, что должен быть простой способ, как это сделать в R, но я не могу найти как. слияния функция , кажется, сливаются только data.frames, совокупная функция растворения нескольких полигонов в один SHP, gIntersect (путем ввода функции соединение) возвращает логическое значение, а не на всех SPDF.

введите описание изображения здесь

данные доступны здесь: http://ulozto.cz/xpoo5jfL/ab-zip

library(sp)
library(raster)
library(rgeos)
library(spatstat)
library(rgdal)     
library(maptools)

setwd("C:/...")
a<-readOGR(dsn=getwd(), layer="pol.a")
b<- readOGR(dsn=getwd(), layer="pol.b")

ab<-merge(a, b)  # what tool if not "merge" to use??

2
Смотрите? Rgeos :: gUnion и / или?
Raster

Ответы:


21

Если вам не нужно объединять топологию, а просто добавлять новые полигоны, вы можете просто использовать:

ab <- rbind(a,b)

Если вы получаете сообщение об ошибке «Неуникальные значения слотов идентификатора полигона», это означает, что имена строк объектов совпадают. Чтобы это исправить, вы можете использовать spChFID для изменения имен строк и связанных с ними отношений слотов. Поскольку слоты в объекте используют имена строк для связывания объекта, вы не можете просто изменить row.names в слоте @data.

b <- spChFIDs(b, paste("b", row.names(b), sep="."))

Функция union (union_sp) в растровом пакете делает это и вызывает gIntersects из rgeos, негласно и является очень удобной вспомогательной функцией.

**** Edit 08-06-2018 Существует недокументированный аргумент, который можно использовать, чтобы пропустить проблему с дублирующимся идентификатором.

ab <- rbind(a, b, makeUniqueIDs = TRUE) 

Привет, спасибо, я попробовал это, но получил ошибку: Ошибка в validObject (res): недопустимый объект класса «SpatialPolygons»: неуникальные значения слота ID полигона. Как я могу справиться с этой ошибкой?
Maycca

3
Вы можете сделать: ab <- bind(a, b) чтобы избежать этой ошибки
Роберт Хейманс

растр :: union в настоящее время не работает с пространственнымиPOINTSdataframes
Mox

19

Супер простое решение, предоставляемое @mdsumner:

library(sp)
library(raster)
library(rgeos)
library(spatstat)
library(rgdal)     
library(maptools)

setwd("C:/...")
a<-readOGR(dsn=getwd(), layer="pol.a")
b<- readOGR(dsn=getwd(), layer="pol.b")

# use union in {raster} package ?raster::union
ab<-union(a, b)

привело к:

класс (аb)

[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"

введите описание изображения здесь


6
Супер простое решение, предоставленное Робертом Хеймансом, автором растра :)
mdsumner

«Union» (в настоящее время) не работает для пространственных точекdataframes, хотя мне сказали, что это будет в следующем выпуске. @RobertH предложил использовать rbind, хотя я не совсем понимаю, как это работает.
Мокс

Подробности здесь: gis.stackexchange.com/questions/274609/…
Mox

Похоже, raster::unionработает и для класса SpatialLinesDataFrame!
Филиппордо

1
library(sp)
data(meuse)
plot(meuse)
slotNames(meuse) #".Data"     "names"     "row.names" ".S3Class" 
coordinates(meuse) <- ~x+y #Add "ID" column to "meuse"
slotNames(meuse) #[1] "data"        "coords.nrs"  "coords"      "bbox"        "proj4string"
class(meuse) #[1] "SpatialPointsDataFrame"
names(meuse@data)
#[1] "cadmium" "copper"  "lead"    "zinc"    "elev"    "dist"    "om"      "ffreq"   "soil"    "lime"   
#[11] "landuse" "dist.m"
meuse@data <- data.frame(ID=1:nrow(meuse), meuse@data) #adds an ID field
names(meuse@data)
#[1] "ID"      "cadmium" "copper"  "lead"    "zinc"    "elev"    "dist"    "om"      "ffreq"   "soil"   
#[11] "lime"    "landuse" "dist.m" 
#Create a data.frame "df.new" with "IDS" (note different name) and "y" columns.
meuse_table.df <- data.frame(IDS=1:nrow(meuse), y=runif(nrow(meuse)))
class(meuse_table.df) #"data.frame"
#Now we can merge "df.new" to "meuse" (@data slot)
meuse <- merge(meuse, meuse_table.df, by.x = "ID", by.y = "IDS")
#create a new file named meuse, consisting of a merge of:
#   the meuse spatial points (from the original)
#   the dataframe created from the original, using the data.frame command
#   BY the field "ID" in the spatialpointsdataframe
#   By the field "IDS" in the tabular dataframe (df.new) 
head(meuse@data)
# I think the source of unease is that adding an ID field to both files 
#is based on them having the same number of rows in the same order. 
#in ArcGIS, this would be an unreasonable and dangerous assumption.
#R seems to have some sort of 'innate' key field, based on the order read it. 
#This is all great when splitting one file, and merging it back together.
#but what about two files? 
#I think it can be done, but it's a three-step process. 
#First, merge the polygons. Add an ID field, as above.
#Second, merge the tables (as dataframes), and add ID's. as above. 
#Third, attach the merged tables to the merged polygons. 
#For it to work, the order of things in the merge (polgyons, dataframe) needs be identfical. 
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.