Estou tentando colorir um diagrama de Voronoi criado usandoscipy.spatial.Voronoi. Aqui está meu 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()

A imagem resultante:

Como você pode ver algumas das regiões de Voronoi na borda da imagem não são coloridas. Isso porque alguns índices para os vértices de Voronoi para essas regiões são definidos como-1, isto é, para aqueles vértices fora do diagrama de Voronoi. Segundo os docs:

regiões: (lista de lista de ints, shape (nregions, *)) Índices dos vórgonos de Voronoi formando cada região de Voronoi.-1 indica vértice fora do diagrama de Voronoi.

Para colorir também essas regiões, tentei apenas remover esses vértices "externos" do polígono, mas isso não funcionou. Eu acho que preciso preencher alguns pontos na fronteira da região da imagem, mas não consigo descobrir como conseguir isso razoavelmente.

Alguém pode ajudar?