如何在python中使用库三角形的三角

Question

当我从这些三角形中执行α形状时，我无法得到点集的边界。所以我认为我应该使用约束Delaunay三角测量。我选择三角形库来执行此操作。但麻烦的是，我不知道给函数triangle.triangulate（tri， opts="）提供什么。我把我所有的点集输入字典，但它返回我的点集。所以任何人都可以帮助我使用这个功能或另一种选择来执行轮廓检测？谢啦

Answer 1

我不知道三角形库，但您可以使用替代方案来执行此操作。

这是一个普遍存在的问题，你可以得到一个简单的代码来解决这个问题。在算法中，这个问题也被称为凸包，scipy中有一个类可用。

代码如下，并附有注释。

#imports
from scipy.spatial import ConvexHull
import matplotlib.pyplot as plt

## Init graphics
axis = plt.gca(); axis.cla();

#Get points sample
points = np.random.rand(20, 2)

#get Convex hull instance for the points set
hull = ConvexHull(points)

#plor the points
plt.plot(points[:,0], points[:,1], 'v')

#Get and plot the boundary
for simplex in hull.simplices:
    axis.plot(points[simplex, 0], points[simplex, 1], 'r-')

结果图可以很快看到：

参考文献：

• http://www.qhull.org/
• https://docs.scipy.org/doc/scipy/reference/generated/scipy.spatial.ConvexHull.html
• https://gis.stackexchange.com/questions/5426/finding-boundary-co-ordinates-from-given-set-of-point-co-ordinates

Answer 2

为了我的目的，我使用了另一个图书馆：shapely。下面是如何使用它对凹面对象进行三角剖分：

from shapely.geometry import MultiPoint
from shapely.ops import triangulate
points=MultiPoint(data)
triangles = triangulate(points)