Ascii 到 NetCDF 转换

Question

我是 netCDF 库的新用户，所以对于这个基本问题我深表歉意。我正在尝试将具有 4 列（lon、lat、depth、velocity）的 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'

Answer 1

这里有几个问题

不要使用内置函数名作为变量名，例如file
不要将坐标数据放入变量中，而是放入单独的变量中
如果将数据添加到 NetCDF 变量，切片大小必须相等
你当前的错误只是从 numpy 数组中询问错误的维度，在你的示例文件中是 1000,4

这是一个有效的例子，但我没有提取正确的变量数据，因为我不知道实际的结构并且在重塑时也要小心。

from netCDF4 import Dataset
import numpy as np
import sys


input_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('longitude', 128)
mydata.createDimension('latitude', 4)
lat = mydata.createVariable('latitude', 'f4','longitude')
lat.long_name = "latitude"
lat.units = "degrees_north"
lat.standard_name = "latitude"
lon = mydata.createVariable('longitude', 'f4','latitude')
lon.long_name = "longitude"
lon.units = "degrees_east"
lon.standard_name = "longitude"

## mydata.createVariable('depth', 1)

Vp = mydata.createVariable('Vp', 'f4', ('longitude','latitude'), fill_value=0)

Vp[:128,:4] = input_file[:512,3].reshape(128,4)
## lon[:] = ???
## lat[:] = ???
Vp.units = 'km/s'

Ascii 到 NetCDF 转换

Ascii to NetCDF conversion

python

netcdf