다각형을 자르고 데이터를 유지 하시겠습니까?

이 두 개의 다각형이 있습니다.

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를 다음과 같이 시도했는데, 이것으로 spdf1교차되는 부분을 ​​잘라 내고 SpatialPolygons spdf2로 변환 spdf1한다 (즉,에 포함 된 정보를 버림 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)

다각형 오버레이를 다루는 더 많은 함수는 '래스터 패키지'(섹션 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")