0
我是netCDF库的新用户,所以我为这个基本问题道了歉。 我试图将具有4列(lon,lat,depth,velocity)的ascii文件转换为netCDF。我收到以下错误:从ASCII到NetCDF的转换
VP [:] =文件[:,:,:,:] IndexError:数组太多指数
能否请你帮我找的原因错误?我已将变量的维数设置为行数。是对的吗?
在此先感谢。
该文件是大,你可以在这里找到一个样本:https://www.dropbox.com/s/qq9529vcvq8t1av/test.txt?dl=0
# -*- coding: utf-8 -*-
from netCDF4 import Dataset
import numpy as np
import sys
file = np.loadtxt('test.txt', delimiter=' ')
# NC file setup
mydata = Dataset('Vp.nc', 'w', format='NETCDF4')
mydata.description = '3D VP velocity model'
# dimensions
mydata.createDimension('latitude', 1000)
mydata.createDimension('longitude', 1000)
mydata.createDimension('depth', 1000)
Vp = mydata.createVariable('Vp', 'f4', ('latitude', 'longitude', 'depth'), fill_value=0)
Vp[:] = file[:,:,:]
Vp.units = 'km/s'