Borders of a Voronoi diagram: how to get edges in the form (point1, point2) from a scipy.spatial.Voronoi object?

Question

I spent a lot of time trying to get the edges with scipy.spatial.Voronoi diagramt to no avail. Here is the main documentation: http://docs.scipy.org/doc/scipy-dev/reference/generated/scipy.spatial.Voronoi.html

If you create a Voronoi diagram like this:

points = np.array([[0, 0], [0, 1], [0, 2], [1, 0], [1, 1], [1, 2],
[2, 0], [2, 1], [2, 2]]) //Or feel free to use any set of points

then you have access to the following properties of the object:

vor.regions
vor.max_bound
vor.ndim
vor.ridge_dict
vor.ridge_points
vor.ridge_vertices
vor.npoints
vor.point_region
vor.points

But it is unclear how to combine them to get edges in the form (point1, point2) for 2d voronoi diagrams? I know that faces exist because you can build a voronoi diagram and its edgres and vertices, because you can do the following:

voronoi_plot_2d(vor)
plt.show()

- , ? , ( , )

Answer 1

ridge_vertices:

    ridge_vertices  (list of list of ints, shape (nridges, *))
        Indices of the Voronoi vertices forming each Voronoi ridge.

. vertices. , , . -1 , .

script, :

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


points = np.array([[0, 0], [0, 1], [0, 2],
                   [1, 0], [1, 1], [1, 2],
                   [2, 0], [2, 1], [2, 2]])

vor = Voronoi(points)


fig = plt.figure()
plt.hold(True)

# Mark the Voronoi vertices.
plt.plot(vor.vertices[:,0], vor.vertices[:, 1], 'ko', ms=8)

for vpair in vor.ridge_vertices:
    if vpair[0] >= 0 and vpair[1] >= 0:
        v0 = vor.vertices[vpair[0]]
        v1 = vor.vertices[vpair[1]]
        # Draw a line from v0 to v1.
        plt.plot([v0[0], v1[0]], [v0[1], v1[1]], 'k', linewidth=2)

plt.show()

: