Estoy usando Scipy 0.13.0 en Python 2.7 para calcular un conjunto de celdas Voronoi en 3d. Necesito obtener el volumen de cada celda para (des) ponderar la salida de una simulación propietaria. ¿Hay alguna forma sencilla de hacerlo? Seguramente es un problema común o un uso común de las células Voronoi, pero no puedo encontrar nada. El siguiente código se ejecuta, y vuelca todo lo que elscipy.spatial.Voronoi manual sabe sobre

from scipy.spatial import Voronoi
x=[0,1,0,1,0,1,0,1,0,1]
y=[0,0,1,1,2,2,3,3.5,4,4.5]
z=[0,0,0,0,0,1,1,1,1,1]
points=zip(x,y,z)
print points
vor=Voronoi(points)
print vor.regions
print vor.vertices
print vor.ridge_points
print vor.ridge_vertices
print vor.points
print vor.point_region