Я сделал DataArray из файла NetCDF, содержащего местоположение (широту и долготу), время и среднюю температуру. У меня также есть список координат, которые представляют местоположения станций (список кортежей с широтой, долготой). То, что я хочу сделать, это сохранить местоположения, которые точно совпадают с одним из кортежей в списке местоположений.
Я пробовал много вещей, включая методы .sel и .where, но пока ничего не получалось так, как я хотел.
Это то, что у меня так далеко:
ds = xr.open_dataset (r "путь \ tg_2019_12.n c")
da = ds.to_array ()
da = da.sel (долгота = срез (-7,5, 31))
da = da.sel (широта = срез (41,5, 68,5))
da = da.stack (latlon = ['широта', 'долгота'])
list_of_unique_coordinates = [(62,875, 30,875), (63,875 , 23.875), ...]
Кроме того, я хотел бы дать каждому местоположению site_id на основе индекса списка кортежей. Идеальным идентификатором site_id будет значение индекса соответствующего списка кортежей + 1 (так что первому местоположению будет присвоено значение site_id = 1 вместо значения индекса 0).
Надеюсь, вы, ребята, можете мне помочь!