python - 将极地立体投影中的数据重新映射到经纬度

Question

我有一个 NetCDF 文件，其中包含南极立体投影上的风暴轨迹密度数据。X 和 Y 坐标的范围从 -1.924274 到 1.88504（我不知道这里的单位）。.nc 文件的链接在这里。

当我用 Xarray 打开文件并绘制数据时，我有以下内容：

有没有办法在 Python 或 CDO 中将其重新映射到经纬度网格？谢谢！

标签： pythongdalnetcdfcdo-climatepyproj

Answer 1

您可以使用我的包 nctoolkit ( https://nctoolkit.readthedocs.io/en/latest/ ) 来执行此操作，它使用 CDO 作为后端。

import nctoolkit as nc
data = nc.open_data("infile.nc")
data.to_latlon(lon = [lon_min, lon_max], lat = [lat_min, lat_max])
data.plot()

您将需要定义 lon_min 等。