我想从降雨时间系列中提取降雨事件,同时在同一事件中允许X干小时(作为参数).因此,通过降雨事件,我的意思是大约连续降雨(RF> 0),内部最大X连续干小时(RF = 0).
我实际上不想用迭代器和增量的方式来做它,我寻找可以放心的pandas或numpy/scipy工具.
这是我的数据帧的示例.RF
是原始降雨量,RFfill
RF.interpolate()填充nodata.evtId
是为了存储事件唯一ID而创建的字段.
TS RF RFfill evtId
0 1997-11-27 14:00:00 0.3 0.3 NaN
1 1997-11-27 15:00:00 1.1 1.1 NaN
2 1997-11-27 16:00:00 0.2 0.2 NaN
3 1997-11-27 17:00:00 0.0 0.0 NaN
4 1997-11-27 18:00:00 0.0 0.0 NaN
5 1997-11-27 19:00:00 1.1 1.1 NaN
6 1997-11-27 20:00:00 0.6 0.6 NaN
7 1997-11-27 21:00:00 0.0 0.0 NaN
8 1997-11-27 22:00:00 0.0 0.0 NaN
9 1997-11-27 23:00:00 0.0 0.0 NaN
10 1997-11-28 00:00:00 0.0 …
Run Code Online (Sandbox Code Playgroud)