Converti un oggetto poligono spaziale in un frame di dati usando R

Il mio obiettivo è quello di modificare uno shapefile esistente unendo determinati poligoni.

Dopo aver importato lo shapefile e usando il comando UnionSpatialPolygons, ottengo il contorno poligonale che desidero.

Tuttavia, questo è ora un oggetto SpatialPolygons e non un oggetto SpatialPolygonsDataFrame, quindi non sono in grado di esportarlo in un file di forma utilizzando writeOGR.

Come posso aggirare questo problema?

Come dice il nome, uno SpatialPolygonsDataFrame è fondamentalmente solo un oggetto SpatialPolygons con i dati allegati (la tabella degli attributi). I dati devono avere almeno tante righe quante sono le funzionalità

library(rgdal)
ob <- SpatialPolygons(..)# Your SpatialPolygons Object
spp <-     SpatialPolygonsDataFrame(ob,data=as.data.frame("yourData"),proj4string=CRS("+proj=    aea > +ellps=GRS80 +datum=WGS84"))
writeOGR(spp,"shapes","testShape",driver="ESRI Shapefile",)

----MODIFICARE----

Se vuoi convertire la SpatialPolygonsDataFrameschiena in un SpatialPolygonsoggetto, devi solo indirizzare la struttura dell'oggetto all'interno di R

ob <- SpatialPolygons(spp@polygons,proj4string=spp@proj4string)

I problemi:

1: il risultato di UnionSpatialPolygons è un poligono spaziale

2: riconvertire il risultato in un frame di dati poligonale spaziale è un vero dolore

-un. è necessario un frame di dati molto preciso da collegare a un poligono spaziale

-b. i dati utilizzati per UnionSpatialPolygons hanno più righe dell'output e non sono formattati nel modo necessario.

La mia (brutta) soluzione:

### Coerce into spatial polygon data frame with id and row name of spatial polygon

# Make a data frame that meets the requirements above:

df<- data.frame(id = getSpPPolygonsIDSlots(your.spatialpolygon))
row.names(df) <- getSpPPolygonsIDSlots(your.spatialpolygon)

# Make spatial polygon data frame
spdf <- SpatialPolygonsDataFrame(your.spatialpolygon, data =df)

# Then don't forget to make sure the projection is correct
# XXXX is your SRID

proj4string(spdf) <- CRS("+init=epsg:XXXX");
spdf <- spTransform(spdf , CRS("+init=epsg:XXXX"));