Numpy/Scipy：集成三维点阵列，或将三维点阵列转换为可积函数f:IR x IR x IR>IR

Question

我正在使用一个复杂的方法（即我不能在合理的时间内重写它），它返回一个表示为3D数组的函数。基础函数在三维空间上定义，并为三维空间中的每个点返回一个浮点值。所谓的标量函数。现在我想在一个矩形体积上积分这个函数。我如何可靠地做到这一点？你知道吗

我花了一整天的时间试着去做，但它变得一团糟。所以我问你是否能给我一个建议。我试过：

（1）寻找从samples/3D数组中集成3D函数的代码（不管你想怎么称呼它）。我没有发现任何有用的东西。Numpy函数，比如trapz()和simps()，只在1D中工作，我没有发现3D的

（2）将我的函数从三维数组转换为三维插值函数。我可以使用scipy.interpolate.RegularGridInterpolator()来实现这一点，但是它没有使用scipy.integrate.tplquad()或者我尝试过的任何类似代码来集成。即使在使用包装器使其输出格式正确之后，此操作仍然失败。你知道吗

在这种情况下什么是最好的？你知道吗