Python中具有可变限制的2D集成

Question

我有一个看起来像

的积分

其中f（x，y）是数据点的数组。如何在Python中执行此积分？我知道解决变量限制的方法，但是这些方法仅在被积是解析函数（例如f（x，y）= x + y）而不是数据点数组时才适用，至少就我而言。我知道。