Преобразовать объект пространственного многоугольника во фрейм данных, используя R

Моя цель - изменить существующий шейп-файл путем объединения определенных полигонов.

После импорта шейп-файла и использования команды UnionSpatialPolygons я получаю желаемый контур многоугольника.

Однако теперь это объект SpatialPolygons, а не SpatialPolygonsDataFrame, поэтому я не могу экспортировать его в шейп-файл с помощью writeOGR.

Как я могу обойти эту проблему?

Как следует из названия, SpatialPolygonsDataFrame - это просто объект SpatialPolygons с прикрепленными данными (таблица атрибутов). Данные должны содержать как минимум столько строк, сколько имеется объектов

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

----РЕДАКТИРОВАТЬ----

Если вы хотите преобразовать свою SpatialPolygonsDataFrameспину в SpatialPolygonsобъект, вам просто нужно обратиться к структуре объекта в R

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

Проблемы:

1: результат UnionSpatialPolygons - это пространственный многоугольник

2: преобразование результата обратно во фрейм данных пространственного многоугольника - реальная боль

-a. вам нужен очень точный фрейм данных для присоединения к пространственному многоугольнику

-b. Данные, которые вы использовали для UnionSpatialPolygons, имеют больше строк, чем выходные данные, и не отформатированы так, как нужно.

Мое (некрасивое) решение:

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