Polygon abschneiden und Daten behalten?

Ich habe diese zwei Polygone:

library(sp); library(rgeos); library(maptools)

coords1 <- matrix(c(-1.841960, -1.823464, -1.838623, -1.841960, 55.663696,
                    55.659178, 55.650841, 55.663696), ncol=2)
coords2 <- matrix(c(-1.822606, -1.816790, -1.832712, -1.822606, 55.657887,
                    55.646806, 55.650679, 55.657887), ncol=2)
p1 <- Polygon(coords1)
p2 <- Polygon(coords2)
p1 <- Polygons(list(p1), ID = "p1")
p2 <- Polygons(list(p2), ID = "p2")
myPolys <- SpatialPolygons(list(p1, p2))
spdf1 = SpatialPolygonsDataFrame(myPolys, data.frame(variable1 = c(232,
                                                                   242), variable2 = c(235, 464), row.names = c("p1", "p2")))
proj4string(spdf1) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84
+towgs84=0,0,0")
plot(spdf1, col="red")

coords1a <- matrix(c(-1.830219, -1.833753, -1.821154, -1.830219, 55.647353,
                     55.656629, 55.652122, 55.647353), ncol=2)
p1a <- Polygon(coords1a)
p1a <- Polygons(list(p1a), ID = "p1a")
myPolys1 <- SpatialPolygons(list(p1a))
spdf2 = SpatialPolygonsDataFrame(myPolys1, data.frame(variable1 = c(2),
                                                      variable2 = c(3), row.names = c("p1a")))
proj4string(spdf2) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84
+towgs84=0,0,0")
plot(spdf2, col="yellow", add=T)

Ich möchte Teile davon ausschneiden, die von spdf1geschnitten werden spdf2. Ich möchte spdf1jedoch als SpatialPolygonsDataFrame bleiben und alle darin enthaltenen Informationen behalten spdf1@data.

Ich habe gDifference wie folgt ausprobiert, wobei Teile spdf1, die von geschnitten werden spdf2, ausgeschnitten werden , aber dann spdf1in SpatialPolygons konvertiert werden (dh die darin enthaltenen Informationen verworfen werden spdf1@data).

gDifference(spdf1, spdf2, byid=T)

Wie kann ich schneiden in spdf1mitspdf2 den darin enthaltenen Daten arbeiten, diese aber behalten spdf1@data? Ich habe diese ähnliche Frage geprüft und ausprobiert, ohne ein Polygon über SpatialPointsDataFrame zu legen und die SPDF-Daten beizubehalten.

Der einfachste Ansatz wäre

  library(raster)
  x <- spdf1 - spdf2

  # or, more formally
  y <- erase(spdf1,  spdf2)

Weitere Funktionen für die Überlagerung von Polygonen finden Sie unter? 'Raster-package' (Abschnitt XIV). Diese Funktionen verwenden die Basisfunktionen von RGEOS unter der Haube in Funktionen auf Benutzerebene (im Gegensatz zu Funktionen auf Entwicklerebene).

Eine Problemumgehung wäre, die Attribute nach dem Ausführen des Clips erneut hinzuzufügen, während Sie von SpatialPolygonsnach konvertieren SpatialPolygonsDataFrame.

sp3 <- gDifference(spdf1, spdf2, byid = TRUE)
row.names(sp3) <- row.names(spdf1)

spdf3 <- SpatialPolygonsDataFrame(sp3, data = spdf1@data)

spdf3@data

   variable1 variable2
p1       232       235
p2       242       464

plot(spdf3, col="red")