Xarray 在时间变量中改变单位

Question

我有一个 nc 文件，其中的时间变量无法被 xarray.open_dataset() 解码，因为单位是“days_since_Jan11900”，这是来自 days 变量的完整信息：

<xarray.DataArray 'days' (time: 87600)>
array([679352., 679353., 679354., ..., 766949., 766950., 766951.])
Dimensions without coordinates: time
Attributes:
    units:      days_since_Jan11900
    long_name:  calendar_days

我已经用ds = xr.open_dataset('path/tmax_bcc-csm1-1.nc', decode_times = False)打开了文件，但现在我想把nc文件的单位改成Python可以轻松解码的东西。就我而言，时间单位是自 0000 年 1 月 1 日以来的现在天数。它最终将使我的分析变得更容易。

我该怎么做？

Answer 1

要解码时间，xarray 会搜索包含units 属性的变量，其形式为"{time_unit} since {reference_date}"。在您的情况下，我将覆盖您的 DataArray 上的单位属性，然后手动解码时间：

ds.days.attrs["units"] = "days since 0000-01-01"
result = xr.decode_cf(ds)

正如我所评论的，一个警告是我怀疑您拥有的参考日期无效，因此解码时间可能会引发错误。