我正在尝试从多多边形 shapefile 生成大量二进制栅格。
我的快照栅格有一个大像素，0.5 x 0.5 度。我在栅格化大多边形时没有大问题，但是对于小的多边形，我正在写空栅格（全为 0）。我想知道 rasterize 函数中是否有一个容差参数，通过该参数我可以将 1 分配给多边形接触的每个像素（即使接触的多边形的百分比非常小）。
这是代码的一部分
for (i in 1:length(shape)) {
shape.r
谢谢！ 哈维尔。
【问题讨论】：
标签： r gis shapefile r-raster rasterize
这是一些示例数据
f <- system.file("ex/lux.shp", package="terra")
v <- vect(f)
x <- lapply(1:nrow(v), \(i)rescale(v[i,], 0.2))
vv <- vect(x)
r <- rast(v, ncols=10, nrows=10)
b <- as.lines(r)
标准的rasterize 方法会遗漏很多多边形，如果它们相对于像元大小来说很小。
x <- rasterize(vv, r, "ID_2")
plot(x)
lines(b, col="light gray")
lines(vv)
touches=TRUE 的参数在这种情况下会有所帮助
xx <- rasterize(vv, r, "ID_2", touches=TRUE)
plot(xx)
lines(b, col="light gray")
lines(vv)
但如果它们的多边形更小，它们仍然可能被遗漏。解决这个问题的一种方法是也栅格化多边形的质心（栅格化可能应该有一个参数来自动执行此操作）。
cnt <- centroids(vv)
# leaving out touches=T otherwise the example is not interesting
# with these data
x <- rasterize(vv, r, "ID_2")
y <- rasterize(cnt, r, "ID_2")
z <- cover(x, y)
plot(z)
lines(b, col="light gray")
lines(vv)
你也可以考虑使用
x <- rasterize(vv, r, cover=TRUE)
x <- x > 0
【讨论】：