Voronoi Diagram Edges：如何从 scipy.spatial.Voronoi 对象中获取 (point1, point2) 形式的边？

Question

我现在花了很多时间试图从 scipy.spatial.Voronoi 图获取边缘但无济于事。这是主要文档： http://docs.scipy.org/doc/scipy-dev/reference/generated/scipy.spatial.Voronoi.html

如果您像这样创建 Voronoi 图：

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

那么您可以访问以下对象属性：

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

但不清楚如何将这些组合起来以获得 2d voronoi 图的 (point1, point2) 形式的边？我知道边的存在是因为您可以绘制 voronoi 图及其边和顶点，因为您可以执行以下操作：

voronoi_plot_2d(vor)
plt.show()

清楚地描述了 voronoi edgres - 如何获取它们的列表及其起点和终点？如果我只得到实心边缘（而不是不受情节限制的虚线边缘），那也没关系

Answer 1

看看ridge_vertices属性：

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

该列表中的每个元素都是一对整数。每个整数都是一个索引 进入vertices 列表。所以每个元素都定义了一条要绘制的线 沃罗诺伊图。 -1 的索引表示“无限远”的点。

这是绘制 Voronoi 图线条的脚本：

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()

# 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()

它创建：