Xarray - 在不同空间区域上定义的合并集

Question

我有来自地球的空间数据（冰分布），分为两个文件（一个用于北极地区，一个用于南极地区）。时间相同，数据变量相同，只是空间暗淡不同（经度和纬度）。

我正在努力使用合并或连接来分组，然后在一个独特的数据集上。

我尝试了不同的组合来合并和具体化。比如

ds = xarray.merge(objects=arrays, join="left", compat="override")
ds = xarray.concat(arrays, dim='time')

我也试过这个解决方案：combining spatial netcdf files using xarray python

但最后，我仍然无法得到我正在寻找的结果：一个独特的数据集，可以在地球上显示它（我为此使用了 cartopy）。

目前最好的选择是什么？

Answer 1

如果您只有两个Datasets 或DataArrays，则使用combine_first 方法可以很好地工作。你可以只做da3 = da1.combine_first(da2)，即使有多个维度需要组合，这也应该有效。如果您在最终的Dataset 中没有完整的经纬度笛卡尔积的数据，您最终可能会得到NaNs。

如果你有两个以上，这种方法也可以工作，你只需要重复使用combine_first来逐个构建完整的结果。