Рассчитать среднее значение по выбранной области с помощью cdo

Question

Я пытаюсь вычислить среднее значение SPI для выбранной области с определенным значением lon / lat. Я сделал следующую область, где я хочу рассчитать среднее значение.

ds_mlw1 = ds.sel(lat=slice(-16.74833, -16.75613), lon=slice(35.27023, 35.27915))

Я увидел, что лучше всего работать с пакетом cdo. Я получаю сообщение об ошибке, что мой CD не распознается Кто-нибудь знает, как решить эту проблему.

Это мой сценарий

from cdo import *
import matplotlib.pyplot as plt
file = 'spi3_6_12_1deg_cru_ts_3_21_1949_2012.nc'
cdo = Cdo()
cdo.debug=True


meanspi  = cdo.fldmean(input=file,returnCdf=True).variables['spi3'][:] 

Я получаю следующую ошибку

---------------------------------------------------------------------------
FileNotFoundError                         Traceback (most recent call last)
<ipython-input-77-7eeb6d09369c> in <module>()
      9 import matplotlib.pyplot as plt
      10 file = 'spi3_6_12_1deg_cru_ts_3_21_1949_2012.nc'
  ---> 11 cdo = Cdo()
     12 cdo.debug=True
     13 

FileNotFoundError: [WinError 2] Het systeme kan het opgegeven bestand niet vinden

значение (система не может найти данный файл)

Answer 1

В конце я не использовал пакет cdo, но использовал функцию slice в xarray.

ds_mlw = ds.sel(lon=slice(33,37), lat=slice(-18,-13), time=slice("2002-01-01", "2017- 
11-01"))
spi_avg = ds_mlw.mean(dim=('lon', 'lat'))
plt.figure(figsize=(10,5))
spi_avg['spi3'].plot()

Answer 2

Я не знаком с cdo, но используя ipython (интерактивный питон), я отладил этот раздел кода:

from cdo import *
cdo=Cdo() 

И смог получить вашу ошибку:

FileNotFoundError: [Errno 2] No such file or directory: 'cdo': 'cdo'

Я также вижу в трассировке исключений, что она ищет исполняемый файл 'cdo'. Похоже, что модуль python cdo является оберткой вокруг двоичного файла программы cdo.

Если вы установите программу cdo, она, вероятно, решит вашу проблему. Информация появляется здесь? : https://code.mpimet.mpg.de/projects/cdo/wiki/Cdo#Documentation

Связано ли это изменение климата? Это очень круто!