python百分位数以15天窗口为中心



我对Python函数和一般Python(用于大数据集的xarray)有点陌生。是否存在一个函数,可以计算90个百分位集中在一个长时间序列的一个域的15天窗口?

即我有62年的北美夏季日最高温度(我在JJA内与cdo合并全天)。我想基于居中的15天窗口(为了计算热浪)计算该时期(夏季)tmax的第90个日百分位数

我想知道它是否存在一些可以直接计算这个长时间序列的东西?百分位数是多少?我在想也许是np。百分位数是np。roll之类的?但不知道如何正确使用

提前感谢!

如果你想使用python,你可以使用:

ds = xr.open('file.nc')
ds_90p = ds.rolling(time=10, min_periods=1).construct('tmp').quantile(0.9, dim='tmp')

然而,这将需要很长时间,取决于你的电脑,显然。我强烈建议使用cdo,它可以更快地产生相同的结果。

cdo runpctl,90,10 input_file.nc output_file.nc

最新更新