如何用python自动读入NCEP数据的时间

winston90

RT
如何用python中的Scientific.IO.NetCDF模块读入NCEP数据后，自动将该数据中的时间（从1800-1-1日起计算的小时数）转化为正常的时间？
如果有其他的处理这个问题的模块也可以。

po_po1

我用的是netCDF4模块，和Scientific.IO.NetCDF模块差不多。
from netCDF4 import num2date

from netCDF4 import Dataset

比如q是一个nc数据：

q=Dataset(r'D:\Desktop\model\Bats\gao5\gao2000monmean_lonlat.nc')

t=q.variables['time']

time=t[:]

timed=num2date(time,t.units)

timed就是一个时间序列：

In [4]:timed

Out[4]:

array([datetime.datetime(2000, 1, 31, 18, 0),
       datetime.datetime(2000, 2, 29, 18, 0),
       datetime.datetime(2000, 3, 31, 18, 0),
       datetime.datetime(2000, 4, 30, 18, 0),
       datetime.datetime(2000, 5, 31, 18, 0),
       datetime.datetime(2000, 6, 30, 18, 0),
       datetime.datetime(2000, 7, 31, 18, 0),
       datetime.datetime(2000, 8, 1, 0, 0)], dtype=object)

独一吖

谢谢分享，有一个问题请教：如何批量取出timed序列的年份/月份等等呢，我试了一下只能取单个的，比如timed[0].year

阿阿飞飞

顶楼上的，NetCDF4模块要比Scientific.IO.NetCDF好用，更强大一些，但在xp系统下使用会有问题。
需要注意的是num2date(time,t.units)一般转化出来的是世界时，如果自己再在本系统下创建一个时间，建议使用另外一个时间模块：dateutil ，利用该模块可以方便的实现时区转换，比如想在数据的时间序列中寻找某一年：

1. from dateutil import tz
   
2. year,m,d=2000,8,1
   
3. d1=dt.datetime(year,m,d,tzinfo=tz.tzutc())
   
4. index=np.where(timed>d1)

复制代码

kdy379405950

这个，我只是来水经验的

Arsenal

python自带datetime模块，用datetime模块里的datetime类和timedelta类做个简单加法就行了

Swallow

我也需要这个问题的解决方法！

tiantianyixin

po_po1 发表于 2015-5-19 16:42
我用的是netCDF4模块，和Scientific.IO.NetCDF模块差不多。
from netCDF4 import num2date

正好用到，get

liuxiunuli

今天刚好读时间，又刚好看到这个，，，谢谢分享

杨杨dadadada

谢谢楼上的分享，试试看