Konversi objek poligon spasial ke bingkai data menggunakan R

Tujuan saya adalah memodifikasi bentuk yang ada dengan menggabungkan poligon tertentu.

Setelah mengimpor shapefile dan menggunakan perintah UnionSpatialPolygons, saya mendapatkan garis besar poligon yang saya inginkan.

Namun, ini sekarang menjadi objek SpatialPolygons dan bukan SpatialPolygonsDataFrame, jadi saya tidak dapat mengekspornya ke shapefile menggunakan writeOGR.

Bagaimana saya bisa mengatasi masalah ini?

Seperti namanya, SpatialPolygonsDataFrame pada dasarnya hanyalah objek SpatialPolygons dengan data yang dilampirkan (tabel atribut). Data harus memiliki paling tidak baris sebanyak fitur

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

---- EDIT ----

Jika Anda ingin mengonversi SpatialPolygonsDataFramepunggung Anda ke SpatialPolygonsobjek, Anda hanya perlu membahas struktur objek dalam R

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

Masalah:

1: hasil dari UnionSpatialPolygons adalah poligon spasial

2: mengubah hasilnya kembali menjadi bingkai data poligon spasial sangat menyebalkan

-Sebuah. Anda membutuhkan bingkai data yang sangat tepat untuk dilampirkan ke poligon spasial

-b. data yang Anda gunakan untuk UnionSpatialPolygons memiliki lebih banyak baris daripada output dan tidak diformat dengan cara yang diperlukan.

Solusi (jelek) saya:

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