当我想用python在现有的netcdf文件上添加新的组时,我感到很困难。我们在我的机器上工作。
米维
import os
import numpy as np
from scipy.io import netcdf
from netCDF4 import Dataset
try:
os.remove('test.nc')
except OSError:
pass
mydata = np.linspace(1,30,30)
g=Dataset('test.nc','w',format='NETCDF4')
g.description = 'The data in this file was created for a test'
g.close()
#Do stuff
g=Dataset('test.nc','r+',format='NETCDF4')
g1 = g.createGroup('grp1')
g1.createDimension('dim', len(mydata))
ncdata = g1.createVariable('data1', 'f8', ('dim',)); ncdata[:] = mydata
ncdata = g1.createVariable('data2', 'f8', ('dim',)); ncdata[:] = mydata
g.close()
#Do stuff
g=Dataset('test.nc','r+',format='NETCDF4')
g1 = g.createGroup('grp2')
g1.createDimension('dim', len(mydata))
ncdata = g1.createVariable('data1', 'f8', ('dim',)); ncdata[:] = mydata
ncdata = g1.createVariable('data2', 'f8', ('dim',)); ncdata[:] = mydata
g.close()但是我不能用ncview打开我的文件,相反,我明白了
Ncview 1.93g大卫W.皮尔斯2009年2月24日page.html版权(C) 1993至2009年,大卫W.皮尔斯Ncview提供绝对没有保证;关于细节类型
ncview -w'. This is free software licensed under the Gnu General Public License version 3; typencview -c‘用于重新分配详细信息。 没有找到可显示的变量!
上面写着no displayable variables found!,但通过ncdump -h test.nc,我可以看到以下几行;
netcdf测试{ //全局属性::description =“该文件中的数据是为测试创建的”; 组: grp1 {维数: dim = 30;变量:双data1(dim);双data2(dim);} //组grp1 组: grp2 {维数: dim = 30;变量:双data1(dim);双data2(dim);} //组grp2 }
我是不是在我的Python脚本中遗漏了什么?(我认为ncview是可以的,因为当我不将它们保存在netcdf组中时,我能够看到非常相同的数据)
发布于 2016-01-29 20:23:38
您可以尝试升级ncview版本;1.93g已经相当老了。2.1.6 (2015年秋季发布)可在以下网站获得:page.html
https://stackoverflow.com/questions/35071673
复制相似问题