python - 如何在 voronoi_finite_polygons_2d 函数中获取与给定点关联的区域?
问题描述
我正在尝试修改我在 stackoverflow 上找到的这段代码,以创建一个具有有限边界的 voronoi 单元。
但是我的问题是我不知道如何将区域关联到给定点。这是在普通的 Voronoi 中使用该point_region
方法完成的,但这在这里不起作用,因为区域已经改变。
我使用的数据点是:
points = array([[289255.176 , 921667.461 ],
[289296.31699, 921687.13826],
[289463.30305, 921770.12504],
[289725.08002, 921905.75745],
[289960.48198, 922099.46056],
[290106.98928, 922361.79529],
[289255.184 , 921646.244 ],
[289307.48677, 921627.05485],
[289500.80493, 921555.50067],
[289825.14532, 921435.65147],
[290141.79326, 921322.77935],
[290454.91721, 921211.09355],
[289255.187 , 921635.627 ],
[289327.07776, 921558.85263],
[289565.21795, 921298.17707],
[289875.40176, 920978.013 ],
[290192.86361, 920656.82017],
[289255.185 , 921630.386 ],
[289318.54181, 921453.18492],
[289421.06861, 921167.57934],
[289565.42462, 920770.1386 ],
[289701.83141, 920376.28627],
[289833.6501 , 919990.66467]])
vor = Voronoi(points)
min_x = vor.min_bound[0] - 100
max_x = vor.max_bound[0] + 100
min_y = vor.min_bound[1] - 100
max_y = vor.max_bound[1] + 100
regions, vertices, pts = voronoi_finite_polygons_2d(vor)
解决方案
此代码使用多边形分析中的点为您提供与给定有限Voronoi 区域关联的点。让我们假设我们希望找到与 vor_regions 列表中第二位的 Voronoi 区域相关联的点。
vor = Voronoi(points)
vor_regions = vor.regions
vor_vertices = vor.vertices
from shapely.geometry import MultiPoint
from shapely.geometry import Point
region = vor_regions[1]
coords = tuple(map(tuple, vor_vertices[region]))
poly = MultiPoint(coords).convex_hull
for i in range(0,len(points)):
pt = Point(tuple(map(tuple, points[i:i+1])))
if poly.contains(pt) == True:
print('The input region is ')
print(region)
print('Index of the associated point in the vor_vertices list is '+str(i))
print('The coordinates of the associated point is ')
print(pt)
输出:
The input region is
[6, 2, 1, 0, 5]
Index of the associated point in the vor_vertices list is 4
The coordinates of the associated point is
POINT (289960.48198 922099.46056)
推荐阅读
- rest - req.FormValue 不映射 req.body 值
- ruby-on-rails - Rails webpacker 为所有 JSX 表达式返回 SyntaxError
- angular - Angular 中 setTimeout 的替代方法
- python - pyproj FutureWarning: '+init=
: ' syntax is deprecated
- java - ANTLR4 只能解析一个类语句,不能解析整个文件
- django - Django REST Framework:自定义序列化程序变量
- c++ - C++ 宏——返回指定格式
- python - 时间转换器不工作,让我“无法将字符串转换为浮点数”
- r - 获取用于注释的行名
- javascript - Gulp watch 导致无限构建循环