Estoy tratando de colorear un diagrama de Voronoi creado usandoscipy.spatial.Voronoi. Aquí está mi código:

import numpy as np
import matplotlib.pyplot as plt
from scipy.spatial import Voronoi, voronoi_plot_2d

# make up data points
points = np.random.rand(15,2)

# compute Voronoi tesselation
vor = Voronoi(points)

# plot
voronoi_plot_2d(vor)

# colorize
for region in vor.regions:
    if not -1 in region:
        polygon = [vor.vertices[i] for i in region]
        plt.fill(*zip(*polygon))

plt.show()

La imagen resultante:

Como puede ver, algunas de las regiones Voronoi en el borde de la imagen no están coloreadas. Esto se debe a que algunos índices de los vértices de Voronoi para estas regiones se establecen en-1, es decir, para esos vértices fuera del diagrama de Voronoi. Según los documentos:

regiones: (lista de ints, forma (nregiones, *)) Índices de los vértices Voronoi que forman cada región Voronoi.-1 indica vértice fuera del diagrama de Voronoi.

Para colorear también estas regiones, he intentado eliminar estos vértices "externos" del polígono, pero eso no funcionó. Creo que necesito rellenar algunos puntos en el borde de la región de la imagen, pero parece que no puedo encontrar la forma de lograr esto de manera razonable.

¿Alguien puede ayudar?