He creado 2 polígonos en QGIS. Usándolos en R, los polígonos se convierten automáticamente en SpatialPolygonsDataFrame (SPDF). Me gustaría combinarlos en un solo SPDF (como es muy fácil en ArcGis usando la herramienta Fusionar ). Estoy seguro de que debería haber una forma simple de completar eso en R, pero no puedo encontrar cómo. La función merge parece fusionar solo data.frames, la función aggregate disuelve múltiples polígonos en un shp, gIntersect (al escribir la función de unión) devuelve lógica valor, no en absoluto el SPDF.
datosestándisponiblesaquí:
library(sp)
library(raster)
library(rgeos)
library(spatstat)
library(rgdal)
library(maptools)
setwd("C:/...")
a<-readOGR(dsn=getwd(), layer="pol.a")
b<- readOGR(dsn=getwd(), layer="pol.b")
ab<-merge(a, b) # what tool if not "merge" to use??