Netcdf 数据集从开始时间到UTC小时的秒转换

Question

我正在使用netcdf格式代码，我需要将时间从开始时间(2016-01-01 00：00：00.0(的秒转换为UTC时间。我对这一切都很陌生，所以我真的很挣扎！

我试过使用netCDF4的num2date。

from netCDF4 import date2num , num2date, Dataset
time=f.variables['time'][:]
dates=netCDF4.num2date(time[:],time.units)
print(dates.strftime('%Y%m%d%H') for date in dates)

属性

错误："屏蔽数组"对象没有属性"单位">

Answer 1

由于您从time=f.variables['time'][:]中的变量中提取time，它将丢失其关联的单位(正如错误所说，时间只是一个屏蔽数组(。 你必须喂给num2date()的是variables['time'].units，例如

from netCDF4 import date2num, num2date, Dataset
file = ... # your nc file
with Dataset(file) as root:
time = root.variables['time'][:]
dates = num2date(time, root.variables['time'].units)
## directly get UTC hours here:
# unit_utchours = root.variables['time'].units.replace('seconds', 'hours')
## would e.g. be 'hours since 2019-08-15 00:00:00'
# utc_hours = date2num(dates, unit_utchours)
# check:
print(dates[0].strftime('%Y%m%d%H'))
# e.g. prints 2019081516

。要将日期作为数字获取，您可以例如

num_dates = [int(d.strftime('%Y%m%d%H')) for d in dates]
# replace int with float if you need floating point numbers etc.

。若要获取以 UTC 小时为单位的日期，请参阅第一个代码块中的注释部分。由于日期数组包含 datetime.datetime 类型的对象，因此您也可以这样做

utc_hours = [d.hour+(d.minute/60)+(d.second/3600) for d in dates]