如何通过特定变量值从 netcdf 文件中获取 select 数据？

Question

我正在从 NetCDF 文件中以特定变量值搜索 select 数据的选项。 数据集包含时间、纬度和经度坐标以及一系列变量。 其中一个变量是具有陆地/开阔海洋/海冰/湖泊的特定值的掩码。 由于公海由 ds.mask = 1 表示，我只想提取位于 mask = 1 的坐标（时间和空间）的海面温度值。但是，我不希望海面温度值位于其他坐标设置为 NaN，但只保留 ds.mask = 1 的那些坐标和变量值。我知道如何使用 xarray.sel/isel 来 select 和数据，但是，这只适用于通过坐标选择，而不是通过我正在尝试的变量值。 任何帮助将不胜感激。

lati = stormtrack_lat.values
loni = stormtrack_lon.values
timei = stormtrack_datetime.values
tmax = timei.max() + np.timedelta64(10,'D')
tmin = timei.min() - np.timedelta64(10,'D')
SSTskin_subfile = SSTskin_file.sel(time=slice(tmin, tmax))

#HERE I NEED HELP:
#extract data where mask = ocean (1) and use only these data points and keep these only!
SSTskin_subfile_masked = SSTskin_subfile.sel(SSTskin_subfile.mask == 1) #does not work yet (Thrown error: ValueError: the first argument to .isel must be a dictionary)

这是 NetCDF 文件的结构：

Answer 1

您可以使用.where 应用海洋面具：

SSTskin_subfile_masked = SSTskin_subfile.where(SSTskin_subfile.mask)

由于数据是网格化的，因此无法删除所有被屏蔽的点。 例如，如果给定纬度只有一个定义值，则必须保留所有值。 但是，您可以使用以下命令删除所有值均为 NaN 的坐标：

SSTskin_subfile_masked.dropna(dim = ['lat', 'lon'], how = 'all')