如何在Python中找到非闭合网格的体积？

import numpy as np
from scipy.spatial import Delaunay
import matplotlib.pyplot as plt

vertices = np.array([[0, 0, 0], [1, 0, 0], [0, 1, 0], [0, 0, 1]])

tri = Delaunay(vertices)

volumes = []
for simplex in tri.simplices:
    p1, p2, p3 = vertices[simplex]
    volume = np.abs(np.dot(p1 - p3, np.cross(p2 - p3, p1 - p3))) / 6
    volumes.append(volume)
total_volume = sum(volumes)

print("Non-closed mesh volume:", total_volume)