问题描述
我想连续 3 天计算高于某个阈值（例如 0.95）。我在 bwee 中编写了代码来计算。我创建的这段代码并不是每年都会产生。它只给出全年的一个结果。
djf <- stack(list.files(pattern='meanDJF.nc'))
sd <- brick(djf)
fun <- function(x) {seq <- rle(as.numeric(x)>as.numeric(x)[[nlayers(sd)]])
n = length(seq$lengths == 3 & seq$values == TRUE) + quantile(x,probs = c(0.95),na.rm = T)
return(n)
}
sd.dn <- calc(sd,fun)
如何计算月度或年度连续 3 天高于 95%？
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）