Объединение двух SpatialPolygonsDataFrame объектов

Question

Использование шейп доступной here я пытаюсь объединить два многоугольники Судана и Южного Судана, так что я получаю границы Судана в 2010 году мой код, чтобы сделать шейп доступный в R является

library(rgdal) 
africa <- readOGR(dsn = "Data/Shapes", layer = "AfricanCountries") 
class(africa) 

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

Я пробовал различные пакеты и решения, как raster::intersect, rgeos::gIntersect или maptools::unionSpatialPolygons. Я всегда получаю пространственный объект, который потерял данные, принадлежащие двум полигонам, или вообще не объединялся.

До сих пор у меня есть следующее решение, которое, кажется, не очень удобно:

# Seperate Polygons and Data 
sp_africa <- SpatialPolygons(africa@polygons, proj4string = CRS(proj4string(africa))) 
dat_africa <- africa@data 

# Delete South Sudan from dataframe 
dat_africa <- dat_africa[-which(dat_africa$COUNTRY == "South Sudan"),] 

# Add row.names to data according to polygon ID'S 
rownames(dat_africa) <- dat_africa$OBJECTID 

# Get all ID's 
allIDs <- africa$OBJECTID 

# Get the ID for Sudan, Primary Land (so we only merge the 
# Sudan main land (no islands) with the South Sudan main land 
sudanID <- africa[which(africa$COUNTRY == "Sudan" & africa$Land_Type == "Primary land"),]$OBJECTID 

# Change ID of South Sudan to that of Sudan 
allIDs[which(africa$COUNTRY == "South Sudan")] <- sudanID 

# Now unite polygons and afterwards merge polygons with data 
tmp  <- unionSpatialPolygons(sp_africa, IDs = allIDs) 
africa2 <- SpatialPolygonsDataFrame(tmp, data = dat_africa) 

Если проще, более удобный способ, которым я был бы рад узнать об этом.

Answer 1

Вы можете использовать aggregate из пакета raster.

nsudan <- africa[africa$COUNTRY == "Sudan",] 
ssudan <- africa[africa$COUNTRY == "South Sudan",] 

plot(nsudan, axes = T, ylim = c(0, 25)) 
plot(ssudan, add = T) 

sudan <- aggregate(rbind(ssudan, nsudan)) 
plot(sudan) 

, так как вы создаете новый объект SpatialPolygons, будет трудно сохранить данные из всех функций. Вы, вероятно, следует добавить старую информацию из nsudan

# remove Sudan and South Sudan 
africa <- africa[!africa$COUNTRY %in% c("Sudan", "South Sudan"),] 

# Adjust the Polygon IDs 
africa <- spChFIDs(africa, as.character(africa$OBJECTID)) 
sudan <- spChFIDs(sudan, as.character(max(africa$OBJECTID) + 1)) 

library(maptools) 
africaNew <- spRbind(sudan, africa) 
plot(africaNew)