我有两个 netcdf 格式的空间数据集。它们具有相同的时间、维度、坐标和数据变量。但它们适用于不同的空间坐标。在下面，我尝试通过多边形显示我的两个数据集：

第一个数据集图：

第二个数据集图：

我想将这两个 netcdf 文件与 xarray 合并。

错误：

尝试过：

但情节combined与上面的第一个情节相同。然后尝试：

错误：

然后尝试：

并且情节combined再次与第一个情节相同。根据ThomasNicolas 的建议，我使用了以下代码：

但它返回此错误：

结果中没有任何数据：

第一个数据集的打印（例如）显示：

对这些文件的合并或组合有什么建议吗？

更新

根据@dl.meteo 的建议，我使用satpy库来解决我的问题，它似乎可以合并两个 netcdf 文件但不能完全合并，你可以在合并的图像中看到不正确的部分（红色边界）。能satpy正确做到吗？

我正在按照Pyresample 示例绘制 Sentinel-3 场景和海岸线。但是，似乎pyresample.kd_tree.resample_nearest并pyresample.save_quicklook期望将 2D 数组作为数据，而不是我试图绘制的 3D（3 波段）Sentinel 场景：

伪代码和错误是：

如何绘制这些数据？

当我使用 satpy 读取 Himawari-8 标准数据时，终端总是说“'filenames' is provided but is empty。”

Himawari-8 有 16 个观测波段，每 10 分钟拍摄一张全盘照片，因此一个观测结果文件夹包含 16 个文件夹（一个文件夹为一个波段），每个波段分为 10 个部分，因此在一个波段中文件夹它有10个文件。
Himawari-8 标准数据的格式为 .dat，数据命名为 .dat HS_aaa_yyyymmdd_hhnn_Bbb_cccc_Rjj_Skkll.DAT。

例如HS_H08_20210518_1100_B01_FLDK_R10_S0110.DAT

这些是来自 vscode 的所有消息：

我使用来自 Github 的代码，他上传了代码“ https://github.com/gSasikala/ftp-himawari8-hsd/blob/main/examples/Processing_Satellite_Imagery.ipynb ”

我使用编辑器 vs 代码编写，python 版本是 3.9.7，我已经在我的电脑上安装了 anaconda 和 satpy。

抱歉，我是 Python 和 Himawari-8 标准数据的新手，如果有什么遗漏的，请提醒我。感谢您对我的任何回复或建议。

示例文件可以从https://noaa-goes17.s3.amazonaws.com/ABI-L1b-RadF/2021/213/00/OR_ABI-L1b-RadF-M6C13_G17_s20212130000319_e20212130009396_c20212130009445.nc下载

我正在尝试将投影转换为板 carree。

我使用 Satpy 加载了 netcdf。

然后我想使用https://satpy.readthedocs.io/en/stable/resample.html#create-custom-area-definition将其重新采样到 Plate carree， 但没有示例代码。

我想使用 GOES17、EUMETSAT 和 GK-2A 生成全球气象卫星图像。我想让它成为 Plate carree 坐标。（GOES 17 netcdf 文件转换为 Plate Carree）

首先，使用 Satpy，我制作了盘子 carree 图像。

像这种方法，我想制作其他卫星图像并合并到 1 个图像文件。但是有没有任何简单或最简单的方法来生成全球天气图像？我的最终目标是生成全球卫星图像的 numpy 数组。

-- 我的整个代码 --