我正在尝试将霜冻天数添加到 xarray 教程数据集中。
airtemps = xr.tutorial.load_dataset('air_temperature')
# The set spans more than a year, let's take only one
airtemps = airtemps.sel(time=slice('2013-01-01', '2013-12-31'))
airtemps['air'] = airtemps.air - 273.15
数据在时间上是相当高分辨率的,所以为了更容易处理,我对它们重新采样 air_day = airtemps.resample('1D', 'time', how='mean') air_month = air_day.resample('1M', 'time' , 如何='意思')
结果,我想在air_month
所有三个维度上添加一个额外的变量,其中包含平均值低于零的天数。
我非常天真的尝试将是……。像这样
air_month['frost'] = sum(air_day.air < 0)
但是通过sum()
ing,我在这里失去了时间维度。我被困在这里,仍然没有将 xarray 概念包裹在我的脑海中。
感谢帮助!