使用 NetCDF python 编写时间日期系列

问题描述 投票:0回答:1

我用这样的熊猫创建了一个时间序列:

date_series_min = pd.date_range(yy+'-'+mm, periods=44640, freq='T')

yy 是文件名中的年份(例如:2001),mm 是月份(例如:01)。 然后我使用上面指出的方法每分钟创建一个时间序列的日期。 我看到这个系列的变量类型是:

<class 'pandas.core.indexes.datetimes.DatetimeIndex'>

如何在 netCDF 文件中写入这个时间序列?

提前致谢

python pandas time netcdf4
1个回答
0
投票

我做了一个这样的脚本,我使用 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" ; 
}
© www.soinside.com 2019 - 2024. All rights reserved.