（时区知道）netcdf中使用Python的datetime

Question

我试图在netcdf文件中保存一个timeserie。根据我发现的文档，这可以使用netCDF4模块中的date2num方法完成。我不能然而，它的工作（见下面的例子）：

from datetime import datetime as dt 
from netCDF4 import Dataset 
from netCDF4 import num2date, date2num 
import pytz 

filename = 'test.nc' 
root = Dataset(filename, 'w', format='NETCDF4_CLASSIC') 
root.name = 'test' 

# create dimension 
root.createDimension('datetime', None) 

# create variable 
timeserie = root.createVariable('timeserie', 'f4', ('datetime',)) 
timeserie.units = 'days since 1970-01-01 00:00:00 UTC' # reference: epoch 
timeserie.calendar = 'gregorian' 

# create testList 
listDT = [dt.now(pytz.utc), dt(1970, 1, 2, 12, 0, 0, 0, pytz.utc)] 
print date2num(listDT, units=timeserie.units, calendar=timeserie.calendar) 

root.close() 

错误它给了我：

Traceback (most recent call last): 
    File "test.py", line 20, in <module> 
    print date2num(listDT, units=timeserie.units, calendar=timeserie.calendar) 
    File "netCDF4\_netCDF4.pyx", line 5128, in netCDF4._netCDF4.date2num (netCDF4\_netCDF4.c:60367) 
TypeError: can't subtract offset-naive and offset-aware datetimes 

我将如何保存日期时间数据（最好时区感知）到netcdf使用python？

Answer 1

date2num() doc表示datetime对象必须是UTC。一个数字不能是时区感知的，除非它是UTC。标准/日光时间转换会导致严重问题。

尝试

listDT = [dt.now(), dt(1970, 1, 2, 12, 0, 0, 0)] 
print date2num(listDT, units=timeserie.units, calendar=timeserie.calendar) 

表明它是时区意识到，尝试东部为单位的偏移。

print date2num(listDT, units='days since 1970-01-01 00:00:00-04:00', calendar=timeserie.calendar)