更改 python 中 netcdf 变量的单位

Question

我在Python中有一个netCDF文件，其中时间变量的单位是days_since_Jan11900包无法读取的xarray。它抛出错误

unable to decode time units 'days_since_Jan11900' with 'the default calendar'. Try opening your dataset with decode_times=False or installing cftime if it is not installed.

我不想使用decode_times，只是将其更改为days_since_0. 这是自 0000 年 1 月 1 日以来的天数。

如何将 netCDF 中的单位更改为可读的单位？

<class 'netCDF4._netCDF4.Variable'>

float64 days(time)

    units: days_since_Jan11900

    long_name: calendar_days

unlimited dimensions: 

current shape = (87600,)

filling on, default _FillValue of 9.969209968386869e+36 used

Answer 1

按照我的要求，我自己用 Python 解决了这个问题，而不是通过 cdo。这是我的做法，非常简单：ds['days'] = ds.days - 730966ds.days.attrs['units'] = 'days since 1860-01-01'ds.to_netcdf(path + 'name_of_file.nc')这730966是我在数组中的起始值days，因此我减去了它并将我的单位更改为适当的单位，在我的例子中days since 1860-01-01。

Answer 2

如果您使用的是 Unix，您可以尝试使用 CDO 更改日历。以下内容应该可以解决您的问题：cdo&nbsp;-L&nbsp;&nbsp;-setreftime,0001-01-01&nbsp;&nbsp;-settaxis,0001-01-01,0:00:00,1day&nbsp;&nbsp;-setcalendar,standard&nbsp;infile.nc&nbsp;outfile.nc我已将时间设置为第 1 年，因为我认为任何有效日历中都不存在第 0 年。