我用这样的熊猫创建了一个时间序列:
date_series_min = pd.date_range(yy+'-'+mm, periods=44640, freq='T')
yy 是文件名中的年份(例如:2001),mm 是月份(例如:01)。 然后我使用上面指出的方法每分钟创建一个时间序列的日期。 我看到这个系列的变量类型是:
<class 'pandas.core.indexes.datetimes.DatetimeIndex'>
如何在 netCDF 文件中写入这个时间序列?
提前致谢
我做了一个这样的脚本,我使用 netCDF4 库来保存 netCDf 文件:
#!/usr/bin/env ipython
# ---------------------
import pandas as pd
from netCDF4 import Dataset,num2date,date2num
import datetime
# ----------------------
yy = '2001';mm='01';
date_series_min = pd.date_range(yy+'-'+mm, periods=44640, freq='T') # This is DatetimeIndex, perhaps ordinary datetime objects are preferred...
datesout = date_series_min.to_pydatetime()
# ----------------------
unout = 'days since 2001-01-01 00:00:00'
fileout = 'test.nc'
# ----------------------
with Dataset(fileout,'w','NETCDF3') as ncout:
ncout.createDimension('time',None);
timevar = ncout.createVariable('time','float32',('time'));timevar.setncattr('units',unout);
timevar[:]= date2num(datesout,unout);
控制台的输出(
ncdump -h $fileout
)是这样的:
netcdf test {
dimensions: time = UNLIMITED ; // (44640 currently)
variables:
float time(time) ;
time:units = "days since 2001-01-01 00:00:00" ;
}