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

Question

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

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

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

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 等。