finite_voronoi.py

import trimesh
import numpy as np

from scipy.spatial import Voronoi

if __name__ == '__main__':
    points = np.random.random((20,3))
    
    v = Voronoi(points)
    
    meshes = []
    for idx in v.regions:
        # only mesh finite regions
        if -1 in idx: continue
        
        if len(idx) < 4: continue
        
        m = trimesh.convex.convex_hull(v.vertices[idx])
        
        # color BS
        color = trimesh.visual.random_color()
        color[3] = 50
        m.visual.face_colors = color
        
        
        meshes.append(m)
        
    # show with points
    trimesh.Scene(np.append(meshes, trimesh.points.PointCloud(points))).show(smooth=False)