寻求Cartopy解决方案：陆地覆盖含数据的栅格区域（zorder失效时）

补充说明 #

我的问题实际出现在工作流的后续环节，是在使用Earthpy叠加两个tiff文件时产生的。相关可复现资源包括test_overlay.py脚本、Earthpy环境yml配置文件及叠加用的两个tiff文件。

问题描述 #

GHRSST-MUR卫星产品生成的海表温度（SST）数据中，部分栅格单元实际覆盖了陆地，这种情况在岛屿区域尤为明显。例如绘制包含美国长岛区域的地图时，用ncview、NOAA ERDDAP或Cartopy都会出现SST数据覆盖长岛，但不覆盖大陆区域的情况。调整zorder让Cartopy陆地图层最后绘制无法解决该问题，推测原因是该卫星产品中岛屿上方的栅格单元包含有效数据。我需要找到Cartopy中的相关设置，能在实际存在陆地的位置绘制陆地，即使这些经纬度对应的SST数据数组包含数值。

已尝试的方法 #

• 将Cartopy陆地要素的zorder设为4，SST的pcolormesh绘图zorder设为1，但未达到预期效果：陆地要素会将海洋和陆地区域全部灰化，遮挡了海洋上的SST数据。

验证情况 #

• 用ncview打开同一netCDF文件，确认SST数据确实覆盖了长岛（图中长岛被SST数据遮挡无法显示）。
• NOAA ERDDAP生成的对应区域图像显示，长岛区域的陆地掩膜失效。

数据源信息 #

数据源为NOAA ERDDAP的jplMURSST41数据集，选取的是2023-11-26T09:00:00Z时间点、纬度范围38.2至42.0、经度范围-74.2至-67.9的analysed_sst数据。

内容的提问来源于stack exchange，提问作者Sam McClatchie