python - 坐标处的采样值

Question

我有某种格式的坐标列表，例如：

my_lats = [12.233, -21.123, 85.112]
my_lons = [-12.543, -51.231, 5.562]

该列表实际上是大约 10 000 个位置。

我还有一个带有坐标的 xarray 数据集。例如 3600x1800 大：

data = xr.open_dataarray('my_global_data.nc')

如何最好地提取列表中每个坐标处的插值？

我可以做这个：

for my_lat, my_lon in zip(my_lats, my_lons):
    data_value = data.interp({'LAT': my_lat, 'LON':my_lon})

但是，使用迭代很慢。

我可以给出一个列表或数组，但是，这会生成一个新的 2D 网格。

data_value = data.interp({'LAT': my_lats, 'LON':my_lons})

标签： pythonpython-xarray