使用 R 从另一个 nc 文件中提取基于 netcdf 文件的值的网格单元

Question

我想从 MERRA2 netcdf 文件中提取网格单元，这些文件仅用 R 标记为 sea：

然而，由于土地区域下边界的弯曲性质，部分土地区域在我的选择范围内。 我想选择离陆地最近的海域区域，但是我的选择中有些陆地区域要过滤掉。

我有另一个 nc 文件，其中包含陆地和海洋区域的分数，我知道可以从中提取我需要的网格单元：

所以我必须定义一个遮罩，我可以使用它来仅从我的第一个文件中选择具有特定阈值的海洋点，因为陆地和海洋的覆盖范围很小。 请问有没有办法用R或者在ArcGIS中做到这一点

Answer 1

也许你可以使用 CDO 做一些事情。 如果您从您的 Landsea 掩码文件 landsea.nc 中决定要用于掩码的土地分数阈值“c”（从您的描述中可以看出 c 需要是一个非常小的正实数），那么您可以定义一个海面为 0，陆地为 1 的掩码

cdo lec,c landsea.nc mask.nc 

这将值 < c 的所有点设置为 1，否则为 0（对于土地）

或者，如果您想将所有陆地点设置为“缺失”，那么您可以

cdo setrtomiss,0.5,2 mask.nc mask_miss.nc

现在土地点被设置为“失踪”

然后您可以使用掩码缩放数据文件以将陆地点设置为零

cdo mul merra2.nc mask.nc out1.nc 

或者您可以将它们设置为 missing ：

cdo mul merra2.nc mask_miss.nc out2.nc 

然后您可以对其进行处理。

当然，这里我假设 landsea.nc 掩码文件与您的 merra 文件位于同一网格上。 如果不是，您将需要先重新网格化掩码文件，您也可以使用 CDO 来完成。