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

Question

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

当我使用 Xarray 和 plot 数据打开文件时，我有以下内容：

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

Answer 1

您可以使用我的 package 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 等。