python - 如何从 scipy.spatial.Voronoi 获取边界会合点的坐标
问题描述
我正在使用 scipy.spatial.Voronoi 来计算 Voronoi 图:
import numpy as np
from scipy.spatial import Voronoi
points = np.array([[51.129378, 17.02925 ],
[51.086225, 17.012689],
[50.913433, 15.765608],
[53.121764, 17.987906],
[53.134083, 17.995708],
[51.75805 , 19.529786]])
vor = Voronoi(points)
我收到以下图表:
如何获取橙色点的坐标?似乎没有(至少记录在案)它的属性。
动机:我的目标是计算没有蓝点的最大圆圈。所以我打算计算每个橙色和蓝色点之间的距离。然后我将检查每个橙色点并选择一个具有最高最小距离值的点。它将是我的目标圈子的中心。也许没有Voronoi还有其他方法可以做到这一点?
解决方案
这些点在vertices
属性中可用。Voronoi的文档将属性描述vertices
为“Voronoi 顶点的坐标”。(在链接的网页中向下滚动以查找属性的描述。)
这是您的示例:
In [5]: import numpy as np
In [6]: from scipy.spatial import Voronoi
In [7]: points = np.array([[51.129378, 17.02925 ],
...: [51.086225, 17.012689],
...: [50.913433, 15.765608],
...: [53.121764, 17.987906],
...: [53.134083, 17.995708],
...: [51.75805 , 19.529786]])
...:
In [8]: vor = Voronoi(points)
这些是绘制为橙色点的点:
In [9]: vor.vertices
Out[9]:
array([[52.56952748, 18.87348869],
[51.7974129 , 18.19059283],
[56.91850562, 12.00665177],
[52.80703622, 16.09228084],
[50.53735155, 18.50739102],
[51.36995954, 16.33786426]])
请注意,voronoi_plot_2d(vor)
在绘图上选择的 x 和 y 限制太小而无法看到所有顶点。这是一个显示所有这些的情节:
In [15]: import matplotlib.pyplot as plt
In [16]: from scipy.spatial import voronoi_plot_2d
In [17]: voronoi_plot_2d(vor)
Out[17]: <Figure size 1280x960 with 1 Axes>
In [18]: plt.xlim(50, 58)
Out[18]: (50, 58)
In [19]: plt.ylim(11, 20)
Out[19]: (11, 20)
推荐阅读
- r - 数据框上的摘要显示了字符字段的专业值。当我对其应用 substr 时:我收到“无用”信息
- hive - 将数据加载到单个列中时的 Hive Serde
- react-native - React-native Fb(错误:没有找到适合 AccessToken 的构造函数)
- android - 为什么手机gps有时会在不同的位置提供相同的坐标
- c# - 两个具有常见布尔问题的 If 语句
- docker - Zeppelin Docker 解释器配置
- javascript - api调用后更新组件的正确方法是什么
- python - 数组排序、切片和反向返回空
- python - 在字符串中查找重复模式
- javascript - 替代 eval() 从服务器执行自动生成的 JS 代码