极坐标投影到规则的纬度/经度网格

Question

这是我上一个问题的后续，因为示例和数据更加具体。 我有一个 304 x 448 的冰数据阵列（在极地立体投影中），其中的 netcdf 可以在这里获得。 每个网格框对应的纬度/经度点可以在这个文件中找到。 我想将此数据转换为常规的纬度/经度网格（例如 180x360），其中每个网格框都是一个度数。 有没有python的方法来做到这一点？ 使用 Basemap 等，还是使用 CDO（气候数据运营商）？

提前谢谢了。 詹姆士

Answer 1

如果您首先将网格添加到数据文件中，例如，

ncks -A -v lat,lon NSIDC_LatLon.nc sample.nc

然后你可以使用 NCO 的ncremap来重新网格，例如，

ncremap -i sample.nc -d 1x1.nc -o out.nc

不幸的是，这将需要 NCO 的当前快照（而不是已发布的版本），因为输入文件包含的维度顺序不同寻常，我刚刚修补了 NCO 以进行处理。 重新网格化输出的示例是here 。

在 OP 的要求下，我重新整理了他的文件，然后上传了它（这里）。 我没有一个 1x1 数据集，所以我首先创建了一个 1x1 网格文件，如手册中所示，然后用它重新网格化（使用 ncremap -g 选项）：

ncremap -i ~/NSIDC.nc -g ${DATA}/grids/180x360_SCRIP.20150901.nc -o ~/NSIDC_1x1.nc

Answer 2

通常在 CDO 中执行此操作的命令是

cdo remapcon,r360x180 in.nc out.nc 

代替 con=conservative 重映射，您可以使用 bil=bilinear, nn=nearest neighbor, con2=2nd order 保守重映射。

根据这篇文章， https ://code.mpimet.mpg.de/boards/1/topics/8302?r = 8326 最新版本的 cdo 应该能够处理极地立体投影。