显示Noaa高程数据

Question

我找到了这个例子，并想复制它：

import numpy as np
import matplotlib.pyplot as plt
from netCDF4 import Dataset

data = Dataset("Path/ETOPO1_Bed_g_gmt4.grd",'r')
print(data.variables.keys())

lon_range = data.variables['x_range'][:]
lat_range = data.variables['y_range'][:]
topo_range = data.variables['z_range'][:]
spacing = data.variables['spacing'][:]
dimension = data.variables['dimension'][:]
z = data.variables['z'][:]
lon_num = dimension[0]
lat_num = dimension[1]

lon = np.linspace(lon_range[0],lon_range[1],dimension[0])
lat = np.linspace(lat_range[0],lat_range[1],dimension[1])

topo = np.reshape(z, (lat_num, lon_num))


plt.imshow(topo, vmax=0)

但是，我不知道x_range是从哪里来的，因为当我打印：

    print(data.variables.keys())
dict_keys(['x', 'y', 'z'])

我怎样才能解决这个问题？

来源：https://earthscience.stackexchange.com/questions/23904/etopo1-region-selection-in-python/23924#23924

Answer 1

你能用X阵列吗？这似乎是一个很好的用例:)

如果没有看到您的数据，就很难确定这是否会开箱即用，但它应该是这样简单的：

ds = xr.open_dataset("Path/ETOPO1_Bed_g_gmt4.grd", engine='netcdf4')
print(ds)
ds.z.plot(vmax=0)