重新启动 Rasterbrick 中的运行计数

问题描述

我正在使用 ERA5 再分析,并且我在某个地理范围内获得了一个价值光栅。每个时间层代表从 1.1.1950 到 28.2.2019 的一天,但非常重要的是,只有冬天的日子。我正在寻找至少 5 天的运行值大于 1(值范围从 0.5 到 1.5)。我正在使用 followin 函数这样做,效果很好:

 ff<-function(x){
 y<-x > 1
 n<- ave(y,cumsum(y == 0),FUN = cumsum)
 sum(n >= 5)
 }
 x1<-calc(x,ff)

然而,在冬季结束时有一个问题。在时间层 28.2.,1.12. 之后。接下来,可能会有一些运行,在季节之间的范围内,例如从 26.2 开始运行。结束 5.12。由我的函数 ff 检测到,但我不想检测到这一点,因为那是两个不同的季节。为方便起见,我将 2 月 29 日的所有日期都放在一边,这样每个季节都有 90 天(31、31、28)。总而言之,我想以某种方式每 90 天重新启动一次运行计数。我希望它是可以理解的,感谢每一个帮助!

解决方法

暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!

如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@)