r - 从R中的栅格中提取缓冲区

Question

我是 R 和空间分析的初学者。我有一个栅格，其中栖息地多边形的单元格等于 1。

 img <- readPNG("Paysage.png")
 map <- raster(img[,,2],xmn=0, xmx=999, ymn=0, ymx=999)
 projection(map) <- "+proj=utm +zone=18 +ellps=GRS80 +datum=NAD83 +units=m +no_defs"
 map[] <- polygons[,c("id_polygon")] 
 map[ map > 0] <- 1

从这个栅格中，我如何为每个多边形构建一个栅格，其中位于给定多边形周围 1 公里缓冲区中的所有多边形单元都等于 1，所有其他矩阵单元都等于 0。

非常感谢您的帮助。

Answer 1

使用 raster 包中的缓冲区功能：只需确保将不希望缓冲区增长的所有单元格设置为 NA：

map[map==0]=NA
library(raster)
b <- buffer(map, width=1000)