如何使用python获取由不规则数据点给出的表面下的音量？

Question

给定一组数据（x，y，z），其中z> 0（大约10 ^ 4个点），我想找到它使用python描述的表面下的体积。 所有这些数据点都是随机生成和未分类的，但是从3D绘图中，表面平滑且消失。见plotted data

到目前为止，我已经检查了以下相关主题 Using Convex Hull：_其他错误： QhullError：QH6114 qhull精度错误：初始单纯形不凸。
Trapezoidal：虽然它运行，但我不相信结果，因为即使生成半球壳，音量也是错误的。
Other Convex Hull：与以前相同的错误。

我也可以制作网格并进行插值，然后计算双积分

protected $table = 'your table'

但是，鉴于griddata

，我不知道如何实现dblquad

# 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')

不起作用。

我该怎么做？