raster Pufferung mit Werten in einem gegebenen Zustand in R

Question

I eine Rasterdatei haben Temperaturwerte enthält, und eine Liste von Punkten als Zentren meiner gewünschten Puffer zu koordinieren, in rot gefärbt, wie in der Figur:

Wie kann ich Puffer um die Zentroide extrahieren, die alle Nachbarzellen enthalten, die der Bedingung Tmax= Tcentroid+6ºC entsprechen? Auch ich würde gerne Puffer überlappen und für den Fall, dass sie sich überschneiden, sie zu einem einzigen zusammenführen.

Answer 1

Hier ist ein Ansatz, um das zu bekommen. Es ist wahrscheinlich ein effizienter Weg, dies über Funktionen im gdistance Paket zu tun

Beispieldaten:

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

Algorithmus (beachten Sie, dass ich die default „Rook“ adjacency Regel bin mit):

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 
} 

inspizieren Ergebnisse:

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