在掩码 netcdf

如何解决在掩码 netcdf

我在 Python 中解析 netcdf 时遇到了困难，似乎有一个简单的解决方案，但我无法解决。

我试图使表面温度接近这个纬度/经度（57.13156000,18.21576000），尽管它在 netcdf 文件中用破折号（--）表示。我明白这是因为该点在陆地上，我如何找到数据存在的最近点？

我可以在寻找最近的纬度/经度之前以某种方式掩盖它吗？

我使用哥白尼的 netcdf 文件，20210720000000-DMI-L4_GHRSST-sstfnd-DMI_OI-NSEABALTIC-v02.0-fv01.0.nc

import os

import netCDF4

import numpy as np


year = datetime.datetime.today().year

month = datetime.datetime.today().strftime('%m')

day = datetime.datetime.today().strftime('%d')

def geo_idx(dd,dd_array):
   geo_idx = (np.abs(dd_array - dd)).argmin()
   return geo_idx

pathtonc = 'copernicus.nc' 
nci = netCDF4.Dataset(pathtonc)

analysed_sst = nci.variables['analysed_sst']

#lat lng of beach
in_lat = 57.13156000
in_lon = 18.21576000

lats = nci.variables['lat'][:]
lons = nci.variables['lon'][:]

lat_idx = geo_idx(in_lat,lats)
lon_idx = geo_idx(in_lon,lons)

analysed_sst_temp = analysed_sst[:,lat_idx,lon_idx]

print (lats[lat_idx])
print (lons[lon_idx])
print(analysed_sst_temp-273.15)

输出：

57.14
18.22
[--]