Berechnen von globalen Statistiken in Raster-Bricks in R ignorieren Pixel Wert 0

Question

Ich versuche globale Stats (SD) im Stapel von 40 Raster Brick zu berechnen, ohne 0 bis NA zu ersetzen. Gibt es dafür einen schnelleren Weg als unten? Ich kann eine Schleife so erstellen, dass die globalSD für jeden Baustein separat berechnet wird. Aber es löst immer noch nicht das Problem, 0 in meiner Berechnung zu ignorieren.

rc <- function(x) { ifelse(x >0, cellStats(fSDstack, stat="sd"), ifelse(x <= 0,x,NA)) } 

test <- calc(fSDstack,fun=rc) 

Dank TiB

Answer 1

Wie wäre das?

library(raster) 
b <- brick(system.file("external/rlogo.grd", package="raster")) 

# here using 255 instead of 0, as that goes better with the example data 
rc <- function(x) { sd(x[x != 255], na.rm=TRUE) } 
test <- calc(b,fun=rc)