Konwertuj przestrzenny obiekt wielokąta na ramkę danych za pomocą R.

Moim celem jest zmodyfikowanie istniejącego pliku kształtu przez scalenie niektórych wielokątów.

Po zaimportowaniu pliku kształtu i użyciu polecenia UnionSpatialPolygons otrzymuję kontur wielokąta, który chcę.

Jednak jest to teraz obiekt SpatialPolygons, a nie SpatialPolygonsDataFrame, więc nie jestem w stanie wyeksportować go do pliku kształtu za pomocą writeOGR.

Jak mogę obejść ten problem?

Jak sama nazwa wskazuje, SpatialPolygonsDataFrame jest w zasadzie tylko obiektem SpatialPolygons z dołączonymi danymi (tabela atrybutów). Dane muszą mieć co najmniej tyle wierszy, ile jest funkcji

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

----EDYTOWAĆ----

Jeśli chcesz przekonwertować SpatialPolygonsDataFrameplecy na SpatialPolygonsobiekt, musisz po prostu zająć się strukturą obiektu w R.

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

Problemy:

1: wynikiem UnionSpatialPolygons jest wielokąt przestrzenny

2: przekształcenie wyniku z powrotem w ramkę danych wielokąta przestrzennego jest prawdziwym bólem

-za. potrzebujesz bardzo dokładnej ramki danych, aby dołączyć do wielokąta przestrzennego

-b. dane użyte w UnionSpatialPolygons zawierają więcej wierszy niż dane wyjściowe i nie są sformatowane w sposób, który jest potrzebny.

Moje (brzydkie) rozwiązanie:

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