光柵緩衝帶中的R

Question

我有一個包含溫度值和座標點作爲我所需緩衝器的中心的一個列表中的光柵文件給定的條件內的值，在圖中以紅色着色像：

如何從包含匹配條件Tmax= Tcentroid+6ºC的所有鄰居單元的質心周圍提取緩衝區？ 另外我想重疊緩衝區，如果它們相交，合併成一個單一的。

Answer 1

這是一個獲得該方法的方法。有可能是一個更有效的方式通過函數在gdistance包來完成

示例數據：

library(raster) 
r <- raster(ncol=36, nrow=18) 
r[] <- 1:ncell(r) 
xy <- cbind(-55, seq(-25, 20, by=20)) 

算法（請注意，我使用的是默認的「車」鄰接規則）：

start <- cellFromXY(r, xy) 
svalues <- extract(r, xy) 

result <- list() 
for (i in 1:nrow(xy)) { 
    value <- svalues[i] 
    cells <- start[i] 
    allcells <- cells 
    while(TRUE) { 
     adj <- adjacent(r, cells, pairs=FALSE) 
     asel <- which(abs(r[adj] - value) < 5) 
     if (length(asel) == 0) break 
     cells <- adj[asel] 
     cells <- cells[!cells %in% allcells] 
     allcells <- c(allcells, cells) 
    } 
    result[[i]] <- allcells 
} 

檢查結果：

p <- xyFromCell(r, unlist(result)) 
plot(r) 
points(xy, pch=20) 
points(p, pch='+')