如何使用python获取由不规则数据点给定的曲面下的体积？

Question

给定一组z>；0（大约10^4个点）的数据（x，y，z），我想用python找到它描述的曲面下的体积。 所有这些数据点都是随机生成和未排序的，但是从三维绘图中可以看出，曲面是平滑的并且消失的。请看这里plotted data

到目前为止，我已经检查了以下相关主题
Using Convex Hull：除了其他错误之外，这个错误：QhullError:QH6114 qhull精度错误：初始单纯形不是凸的。
Trapezoidal：虽然它在运行，但我不相信结果，因为即使生成了半球体壳，这个体积也是错误的。
Other Convex Hull：和以前一样的错误。在

我也认为我可以做一个网格，插值，然后计算二重积分

# Given independents numpy.ndarray x2, y2, and z2
# Set up a regular grid of interpolation points
x1i = np.linspace(x2.min(), x2.max(), 1000)
y1i = np.linspace(y2.min(), y2.max(), 1000)
x1i, y1i = np.meshgrid(x1i, y1i)

# Interpolate; 
z1i = scipy.interpolate.griddata((x2, y2), z2, (x1i, y1i), method='linear')

但是我不知道如何实现给定griddata的dblquad

^{pr2}$

不起作用。在

我该怎么做？在