Python使用libpysal/scipy生成Voronoi镶嵌时触发KeyError:5问题

核心触发原因 #

这个KeyError: 5属于libpysal旧版本的已知实现漏洞，触发逻辑如下：
libpysal的voronoi_regions函数会先读取scipy返回的ridge_points结果构建all_ridges字典，存储每个输入点对应的所有脊线信息。当你加入问题点后，输入点的点位排列刚好触发以下两种异常场景之一，导致索引为5的输入点没有被写入all_ridges字典，调用时直接报错：

• 浮点精度触发的脊线合并逻辑异常：你的输入点位都是EPSG:4326地理坐标，点位间距仅为1e-4度级别（约10米间隔），问题点加入后，和周围点的位置刚好处于scipy Voronoi计算的浮点精度阈值边缘，scipy判定索引5的点和相邻点的脊线重合/无效，不会写入ridge_points结果，自然libpysal就不会为它生成all_ridges对应的键。
• 边缘点处理逻辑漏洞：旧版本libpysal处理Voronoi区域延伸到无穷远的边缘输入点时，存在漏写all_ridges键的逻辑漏洞，问题点加入后刚好让索引为5的点成为所有脊线都带无穷顶点的边缘点，触发漏洞。

可行解决方案 #

• 升级libpysal到最新稳定版，该索引漏洞已经在后续版本的迭代中修复。
• 计算Voronoi前先做坐标投影：EPSG:4326是球面地理坐标系，scipy的Voronoi基于欧氏距离计算，直接使用会有距离误差，也容易触发浮点异常。你可以把坐标转换为埃及区域对应的平面投影坐标系（如EPSG:22992 UTM 36N），用米为单位的平面坐标计算，既提升结果准确性也能规避精度问题。
• 若无法升级依赖，可以给所有输入点添加量级为1e-8的微小随机扰动，避开刚好触发精度阈值的点位排列。

内容的提问来源于stack exchange，提问作者Kaizzen