ラスターレンガのグローバル統計をピクセル無視の値0

Question

0に置き換えることなく、40ラスターブリックのスタックでグローバル統計（SD）を計算しようとしています。これを行うより速い方法が以下にありますか？ globalSDが別々に計算されるようなループを作成することができます。しかし、私の計算では0を無視するという問題は解決していません。これについて

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

test <- calc(fSDstack,fun=rc) 

おかげ のTiB

Answer 1

どのように？

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)