如何使用pykrige插值对降水数据进行插值？

Question

我想对许多日期的台站降水数据进行插值。我在整个数组中有些日期带有零（降水= 0.0），因为它没有下雨。当我使用Pykrige OrdinaryKriging时，出现以下错误：

  

ValueError：“每个下限必须小于每个上限。”

我猜想它与零值有关，因为我对温度值没有问题（它没有所有值都为零的日期）。

在API页面中，它表示OK对象使用L1范数自动计算变异函数变量。如果所有值都相同（在降水情况下为零），则距离为零，则无法计算L1范数。

import numpy as np
from pykrige.ok import OrdinaryKriging

OK = OrdinaryKriging(x, y, vals, variogram_model='gaussian',
     verbose=False, enable_plotting=False)