Извлечь значения осадков из netcdf ⇐ Python

[Anonymous]

У меня застрял вопрос об извлечении переменной из файла NetCDF. Вот структура моего входного файла netcdf:

comephore_all
Out[37]:

Размеры: (x: 85, y: 99, время: 236664)
Координаты:

• x (x) float64 1.156e+06 1.158e+06 1.158e+06 ... 1.24e+06 1.24e+06
• y (y) float64 5.075e+05 5.065e+05 5.055e+ 05...4.105e+05 4.095e+05
• объект времени (времени) '1997-01-01 00:00' ... '2023-12-31 23:00'
  Переменные данных:
  rr1 (время, y, x) float32 dask.array
  crs (время) int32 -2147483647 -2147483647 ... -2147483647 -2147483647
  Атрибуты:
  title: Réanalyse des lames d'eau COMEPHORE
  Conventions : CF-1.6
  история: Ср, 24 марта, 13:39:53 2021: ncrcat 1997_rr.nc 19...
  nco_openmp_thread_number: 1

Я хочу извлечь значения осадков для каждой станции в отдельные файлы .csv.
Я провел первоначальный тест, чтобы извлечь значения для одной станции, но это очень медленно. Вот фрагмент кода:

Код: Выделить всё

indice_x = 1205500
indice_y = 439500
precip_1205500_439500 = comephore_all.sel(x=indice_x, y=indice_y)
df_precip = precip_1205500_439500.to_dataframe(name='EH')[['EH','date']]
df_precip[['date', 'EH']].to_csv(output_path, index=False)

Проблема не в записи CSV-файла, а в извлечении значений осадков. Строка «df_precip = precip_1205500_439500.to_dataframe(name='EH')[['EH', 'date']]» занимает 40 минут для одной станции.
Я также пробовал с даск :

Код: Выделить всё

df_precip = precip_1205500_439500.to_dask_dataframe()
df_precip=df_precip[['time', 'rr1']]
df_precip[['time', 'rr1']].to_csv(output_path)

Но выполнение строки df_precip[['time', 'rr1']].to_csv(output_path) занимает бесконечно много времени.
Мне кажется, я неправильно подхожу к получению значений осадков.
Спасибо за помощь!

Подробнее здесь: https://stackoverflow.com/questions/792 ... rom-netcfd