我有一个由3个多边形组成的spatialPolygonsDataFrame。第三多边形具有与第一多边形相同的形状，但是在第二多边形所在的位置具有孔。
我使用另一个问题（How to add a hole to a polygon within a SpatialPolygonsDataFrame?）的答案来弥补这个缺陷。
library(raster)
library(sp)
# create rasters and store them in a list
r1 <- raster(xmn=1, xmx=5, ymn=1, ymx=5, nrows=4, ncols=4)
r1[] <- 1:length(r1)
# create SpatialPolygonsDataFrame
Sr1 = Polygon(cbind(c(1,5,4,1,1),c(1,2,5,4,1)))
Sr2 = Polygon(cbind(c(2,4,3,2),c(3,2,4,3)))
SpP = SpatialPolygons(list(Polygons(list(Sr1), "s1"), Polygons(list(Sr2), "s2")),
1:2)
dat = data.frame(ID = c("s1", "s2"), value = c("a", "b"))
row.names(dat) <- c("s1", "s2")
p <- SpatialPolygonsDataFrame(SpP, data = dat,
match.ID = TRUE)
AddHoleToPolygon <-function(poly,hole){
# invert the coordinates for Polygons to flag it as a hole
coordsHole <- hole@polygons[[1]]@Polygons[[1]]@coords
newHole <- Polygon(coordsHole,hole=TRUE)
# punch the hole in the main poly
listPol <- poly@polygons[[1]]@Polygons
listPol[[length(listPol)+1]] <- newHole
punch <- Polygons(listPol,poly@polygons[[1]]@ID)
# make the polygon a SpatialPolygonsDataFrame as the entry
new <- SpatialPolygons(list(punch),proj4string=poly@proj4string)
new <- SpatialPolygonsDataFrame(new,data=as(poly,"data.frame"))
return(new)
}
punchedPoly <-AddHoleToPolygon(p[1,],p[2,])
p1 <- rbind(p, punchedPoly, makeUniqueIDs = TRUE)
p1 <- p1[2:3,]
当我使用mask()“修剪”栅格r1时，便会创建孔，尽管三角形多边形具有值并且实际上不是真正的孔。但是它被带有孔的第三个多边形“覆盖”：
masked_hole <- mask(r1, p1)
plot(masked_hole)
当我更改多边形的顺序时，则不会创建孔：
m3 <- mask(r1, p1[c(2,1),])
plot(m3)
光栅化功能以相同的方式受到影响：
r2 <- rasterize(p1, r1, field = "value")
plot(r2)
r3 <- rasterize(p1[c(2,1),], r1, field = "value")
plot(r3)
在我的真实数据中，我有一些孔，其中没有“填充”多边形，而我想保留这些多边形。
如何为在没有孔的地方创建多边形的多边形修复spatialPolygonsDataFrame？
如何解决此问题而无需重新排序，而是“转换”创建孔的多边形？