Convertir un objet polygonal spatial en bloc de données à l'aide de R

Mon objectif est de modifier un fichier de formes existant en fusionnant certains polygones.

Après avoir importé le fichier de formes et utilisé la commande UnionSpatialPolygons, j'obtiens le contour de polygone souhaité.

Cependant, il s'agit maintenant d'un objet SpatialPolygons et non d'un SpatialPolygonsDataFrame, donc je ne peux pas l'exporter vers un fichier de formes à l'aide de writeOGR.

Comment contourner ce problème?

Comme son nom l'indique, un SpatialPolygonsDataFrame est fondamentalement juste un objet SpatialPolygons avec des données attachées (la table attributaire). Les données doivent avoir au moins autant de lignes qu'il y a de fonctionnalités

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

----ÉDITER----

Si vous souhaitez convertir votre SpatialPolygonsDataFramedos en SpatialPolygonsobjet, il vous suffit d'adresser la structure de l'objet dans R

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

Problèmes:

1: le résultat de UnionSpatialPolygons est un polygone spatial

2: reconvertir le résultat en une trame de données de polygone spatial est une vraie douleur

-une. vous avez besoin d'un bloc de données très précis à attacher à un polygone spatial

-b. les données que vous avez utilisées pour UnionSpatialPolygons ont plus de lignes que la sortie et ne sont pas formatées de la manière requise.

Ma (moche) solution:

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