Konvertieren Sie ein räumliches Polygonobjekt mit R in einen Datenrahmen

Mein Ziel ist es, ein vorhandenes Shapefile durch Zusammenführen bestimmter Polygone zu ändern.

Nachdem ich das Shapefile importiert und den Befehl UnionSpatialPolygons verwendet habe, erhalte ich den gewünschten Polygonumriss.

Dies ist jedoch jetzt ein SpatialPolygons-Objekt und kein SpatialPolygonsDataFrame, sodass ich es nicht mit writeOGR in ein Shapefile exportieren kann.

Wie kann ich dieses Problem umgehen?

Wie der Name schon sagt, ist ein SpatialPolygonsDataFrame im Grunde nur ein SpatialPolygons-Objekt mit angehängten Daten (die Attributtabelle). Die Daten müssen mindestens so viele Zeilen enthalten, wie Features vorhanden sind

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

----BEARBEITEN----

Wenn Sie Ihren SpatialPolygonsDataFrameRücken in ein SpatialPolygonsObjekt umwandeln möchten, müssen Sie nur die Objektstruktur in R adressieren

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

Probleme:

1: Das Ergebnis von UnionSpatialPolygons ist ein räumliches Polygon

2: Das Ergebnis wieder in einen räumlichen Polygondatenrahmen umzuwandeln, ist ein echtes Problem

-ein. Zum Anhängen an ein räumliches Polygon benötigen Sie einen sehr genauen Datenrahmen

-b. Daten, die Sie für UnionSpatialPolygons verwendet haben, enthalten mehr Zeilen als die Ausgabe und sind nicht in der erforderlichen Weise formatiert.

Meine (hässliche) Lösung:

### 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"));