python - 使用numpy和scipy将wrf数据的网格点插值到python中的lat lon

Question

我是 python 新手，我正在尝试编写一个函数，该函数将从尺寸为 [时间、高度、经度、纬度] 的 netcdf 文件中获取 numpy 数组，并将函数插值到指定的纬度和经度。我已经研究过 scipy.interpolate 但仍然不确定从哪里开始。有什么帮助吗？

Answer 1

我无法从问题中看出您是否想要(a)将整个数据集插值到一组新的纬度/经度坐标或(b)获取单个纬度/经度位置的值。

无论哪种方式，如果您不习惯在 python 中执行此操作，那么在 Climate Data Operators 中使用 remapgrid 会容易得多(确保安装时支持 netCDF)。这确实要求 netCDF 文件遵循 CF 约定(请参阅此 post 有关使 WRF 文件符合 CF 的信息)。

找到了 remapgrid 的文档 here 。有关网格选项的信息，请特别参阅第 1.3 节(第 9-12 页)，有关插值技术的信息，请参阅第 2.10 节(从第 106 页开始)。对于如何重新映射有多种选项，您使用的命令取决于您的应用程序。

(a) 这是使用双线性插值重新映射到高斯 (128x64) 网格的示例

cdo remapbil,n32 wrfout_d01_1999-01-01-01_00:00:00.nc out.nc

(b) 以下是使用最近邻映射到一个纬度、经度点的重新映射示例:

cdo remapnn,lon=-107.0_lat=34.0 wrfout_d01_1999-01-01-01_00:00:00.nc out.nc