matlab二维三角网格,matlab – 给定多边形和固定点,找到三角网格

您要使用的功能是

DelaunayTri,您可以按照以下步骤操作：

>创建多边形中边缘点的列表.

>获取多边形的所有顶点,并将它们与要包含在多边形内的其他固定点组合.

>创建一个约束三角剖分(正如我在其他答案here和here中所说明的那样).

>如您所述,这将创建凸包的三角剖分(即使您有凹多边形),因此您必须使用方法inOutStatus删除受约束边缘之外的三角形(也在上面链接的答案中说明).

这是一些示例代码：

polygonVertices = [0 0;... %# Concave polygon vertices

0 1;...

1 1;...

0.5 0.5;...

1 0];

polygonEdges = [1 2;... %# Polygon edges (indices of connected vertices)

2 3;...

3 4;...

4 5;...

5 1];

otherVertices = [0.5.*rand(5,1) rand(5,1)]; %# Additional vertices to be added

%# inside the polygon

vertices = [polygonVertices; otherVertices]; %# Collect all the vertices

dt = DelaunayTri(vertices,polygonEdges); %# Create a constrained triangulation

isInside = inOutStatus(dt); %# Find the indices of inside triangles

faces = dt(isInside,:); %# Get the face indices of the inside triangles

使用旧版本的MATLAB …

查看archived version documentation(注意：需要MathWorks帐户),可以看到DelaunayTri首次出现在版本7.8.0(2009a)中.在此之前,可用于执行2-D Delaunay三角测量的唯一内置功能是delaunay,其基于Qhull,因此无法支持受约束的三角测量或非凸面的三角测量.

较新的DelaunayTri使用CGAL.因此,对于早于7.8.0的版本的用户,一个选项是create MEX-files到MATLAB中的CGAL例程接口.例如,如果您要对凹面多边形进行三角剖分,则可以创建一个MEX文件来连接其中一个convex partitioning routines in CGAL,以便将凹多边形分解为一组凸多边形.然后可以使用delaunay对每个凸多边形进行三角测量,并将最终的三角测量组分组为凹多边形的一个较大的三角剖分.