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

编程问答 2022-06-02

问题描述

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

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

enter image description here

有什么方法可以将其重新映射到 Python 或 CDO 中的经纬度网格？ 谢谢！

解决方法

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

cdo-climategdalnetcdfpyprojpython