使用Python将netcdf“land”变量转换为纬度和经度

Question

我有一个全球气象数据集，我想访问某个网格（lat，lon）的数据。但是，数据是压缩的，即感兴趣的参数没有维度（lat，lon），而是“land”land”是一个1D整数数组。在

我在python中使用

import scipy.io.netcdf as netcdf

path = '/path/.../ncfile.nc'

ncfile = netcdf.netcdf_file(path,'r')

然后我检查了有哪些变量，并发现，例如，“rain”变量具有维度（tstep，land）。我在网上查了一下，找到了文件landmask_gswp.nc公司（http://dods.ipsl.jussieu.fr/gswp/Fixed/landmask_gswp.nc），它应该包含我需要的信息，即如何从“land”中提取信息（lat，lon）。 此文件包含变量nav_lat、nav_lon和landmask。据我所知，nav_lat和nav_lon将坐标变量x和y与纬度和经度联系起来。”landmask”是一个二维数组，包含海洋=0或陆地=1的信息。 实际上，着陆点的数量与我的“陆地”1D阵列的长度一致。但是，我不知道如何从中提取（lat，lon）信息。任何帮助都将不胜感激。在

我希望我的问题有点可以理解；我没有编程和/或使用netcdf的经验，所以我希望你能帮忙！提前谢谢！在

Answer 1

其他人可能会觉得这很有帮助，所以下面是问题的答案。一个朋友帮我弄明白了。 “land”是扁平阵列nav_lat、nav_lon的索引，即“land”的第一个条目对应纬度：平坦的[land[0]]和长平头[土地[0]]。在

Answer 2

可以使用从ncfile提取变量

lat = ncfile.variables['nav_lat'][:,:]
lon = ncfile.variables['nav_lon'][:,:] 

这将创建2dnumpy数组lat和lon。在