python - 在 scipy.spatial.Delaunay python 中找到一个点是其一部分的所有单纯形

Question

有没有办法使用 scipy.spatial.Delaunay 得到所有的单纯形/三角形，某个点是 Delaunay 三角剖分的一部分？

我知道有一个 find_simplex() 函数，它只返回一个点所属的 1 个三角形，但我想获取它所属的所有三角形。

所以在示例中，当我对点 6 执行 find_simplex() 时，它只返回三角形 2，但我希望它返回三角形 1、2、3、4、10，和 9，因为点 6 是所有这些三角形的一部分。

如有任何帮助，我们将不胜感激!

Answer 1

您不需要 find_simplex，因为它是几何的，而不是拓扑的。也就是说，它将一个点视为一个位置而不是三角剖分的一个组成部分:几乎所有点都只位于一个单纯形中，所以这就是它报告的内容。

相反，使用顶点编号。简单的答案是使用 simplices 属性:

vert=6
[i for i,s in enumerate(d.simplices) if vert in s]

使用更多代码，可以使用 vertex_to_simplex 更有效地进行搜索和 neighbors 属性。