python - 使用 cdo 计算选定区域经纬度的平均值

Question

我正在尝试计算具有特定经度/纬度的选定区域的平均 SPI 值。我制作了以下区域，我想在其中计算平均值。

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

我发现使用 cdo 包效果最好。我收到一条错误消息，指出我的 cdo 无法识别。有人知道如何解决这个问题。

这是我的脚本

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 systeem kan het opgegeven bestand niet vinden

含义(系统找不到给定的文件)

Answer 1

我对cdo不熟悉，但是使用ipython(交互式python)我调试了这一段代码:

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

这与气候变化有关吗？太酷了!