ポリゴンをクリップしてデータを保持しますか？

次の2つのポリゴンがあります。

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)

とspdf1交差する部分を切り取りたいspdf2。ただし、spdf1SpatialPolygonsDataFrameのままにして、に含まれる情報を保持したいですspdf1@data。

次のようにgDifferenceを試しました。gDifferenceは、とspdf1交差する部分を切り取りますspdf2が、その後spdf1SpatialPolygonsに変換します（つまり、に含まれる情報を破棄しますspdf1@data）。

gDifference(spdf1, spdf2, byid=T)

どのように私はにカットできるspdf1とのspdf2データが含まれているに保持しているがspdf1@data？SpatialPointsDataFrameにポリゴンをオーバーレイしてSPDFデータを保存する方法なしに、この同様の質問を確認してみました。

最も簡単なアプローチは

  library(raster)
  x <- spdf1 - spdf2

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

ポリゴンオーバーレイを処理するその他の関数については、「 'raster-package'（セクションXIV）」を参照してください。これらの関数は、「ユーザーレベル」（「開発者レベル」とは対照的に）関数で、フードの下でrgeosの基本関数を使用します。

回避策は、クリップの実行後にからSpatialPolygonsに変換しながら属性を再追加すること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")