Объем под «плоскостью», определяемой точками данных - python

У меня есть большая сетка точек данных, которые я создал из моделирования, и с каждой точкой в ​​плоскости xy связано значение z (результат моделирования) .

У меня есть значения x, y, z, выгруженные в обычный текстовый файл, и я бы хотел измерить объем между плоскостью xy (т.е. z = 0) и «плоскостью», определяемой точки данных. Точки данных в настоящее время не распределены равномерно, хотя они ДОЛЖНЫ располагаться после завершения моделирования.

Я просматривал документацию scipy, и я не уверен, предоставляет ли scipy.integrate необходимую мне функциональность - похоже, что есть возможность делать это только в 2d, а не в 3d, как мне нужно.

Начнем с того, что, если это не необходимо, я могу обойтись без интерполяции, интегрирование, основанное исключительно на «правиле трапеции» или аналогичной аппроксимации, является хорошей основой для начала.

Любая помощь приветствуется.

спасибо

РЕДАКТИРОВАТЬ: оба решения, описанные ниже, работают хорошо. В моем случае оказалось, что использование сплайна может вызвать "рябь" вокруг резких максимумов на плоскости, поэтому метод Делоне работает лучше, но я бы посоветовал людям проверить и то, и другое.