做图形学相关的开发不管是搞网格生成、有限元前处理还是做地形建模、流体模拟几乎都绕不开一个基础问题怎么把二维平面上的点集、多边形区域剖分成合格的三角形网格。早年用C的时候大家基本都靠Triangle库也就是Jonathan Richard Shewchuk写的那一套经典实现。后来转.NET环境做项目我一直在找一套能在C#里直接调用的替代方案Triangle.Net就是我用了很长一段时间、也实实在在踩过不少坑之后觉得值得拿出来聊聊的库。这篇内容就从下载装配讲起一直到具体的三角剖分调用、参数调优和常见坑点排查尽量把我在实际项目里验证过的用法整理清楚。适合刚接触网格生成的初学者也适合已经做过一些剖分工作、想换用Triangle.Net提升效率的开发者参考。我不打算把它包装成什么神兵利器它就是一套务实、稳定、可复现的网格生成工具但用好的话确实能让很多脏活累活变得不再折腾。1. 内容整体设计与思路拆解1.1 为什么网格生成总是绕不开Delaunay三角剖分先聊一个基础但重要的问题为什么大家都在用Delaunay三角剖分而不是随便把点连成三角形就完事。刚开始做网格生成的时候我也尝试过自己写简单的三角化逻辑结果在边界复杂、点数变多的时候生成出来的三角形又扁又长放进求解器里直接导致数值不稳定。后来才明白Delaunay三角剖分的核心特性是“最大化所有三角形的最小角”在某种意义上能避免生成极端细长的三角形这对后续的有限元计算、渲染和几何处理都至关重要。Triangle.Net这一套库本质上是把经典的Triangle算法移植到了.NET生态里。它支持无约束Delaunay三角剖分也支持带边约束的约束Delaunay三角剖分CDT还能做网格细化控制三角形的大小和形状。我在做多边形区域剖分时用得最多的就是带约束的版本因为现实中几乎不会遇到裸点集剖分大多数输入都是带边界的多边形、带孔洞的复杂区域或者带内部加密线的地质剖面这些都需要约束条件来保证边界完整。这里多说一句物理世界里的网格需求和纯数学意义上的三角剖分是有差别的。数学上可能只关心拓扑正确但工程上还要求三角形不过于畸形、面积分布均匀、边界不能丢点。Triangle.Net把这些工程需求的接口都暴露出来了比如最大角度限制、最小角度限制、最大面积限制这几个参数就是它在工程界被广泛使用的重要原因。1.2 为什么选Triangle.Net而不是自己造轮子在自己写剖分算法和用别人现成库这两条路之间我几乎没有任何犹豫就选了后者。原因很简单网格生成是非常成熟的领域Delaunay三角剖分、约束边处理、网格细化每一步都有几十年的研究积累自己造轮子从零推导实现投入产出比实在太低。当然也不是所有现成库都能满足工程需求。我在选型的时候重点考察了几个方面算法质量生成网格是否满足Delaunay性质边界保持是否精确网格细化是否稳定。接口易用性是否能在C#里直接调用能否自定义数据源而不是绑定死特定的几何容器。许可证和依赖商业项目里能不能用有没有一堆莫名其妙的第三方依赖。性能表现几千个点还好说几十万个点的时候如果算法效率不行整个项目就会卡到让人崩溃。Triangle.Net在这几项上的表现都比较均衡算法基于Shewchuk的经典实现质量有保障接口是纯C#的数据结构也清晰许可证对大多数项目友好性能方面虽然不是极致优化但至少在我的项目里从没成为瓶颈。直观一点对比假如自己要写一个带约束边的Delaunay三角化光是处理“插入约束边时需要删除相交三角形再重新剖分”这个逻辑就足够写几百上千行代码而且边界条件稍微变一变就容易出bug。用现成库的话只需要把边界条件和点集交给它内部这些复杂的几何计算和拓扑更新都封装好了。我个人的态度一直是在一个成熟的生态里能用别人写好的稳定模块就尽量不重复造轮子把精力放在业务逻辑上。2. 下载装配与环境准备2.1 从NuGet安装Triangle.Net关于Triangle.Net的下载现在最省事的方式就是通过NuGet这也是我推荐所有.NET开发者采用的方案。打开Visual Studio在解决方案资源管理器里右键项目选择“管理NuGet程序包”然后搜索Triangle.NET直接安装即可。或者在包管理控制台里执行Install-Package Triangle.NET如果你是.NET Core或者.NET 5的项目也可以用dotnet命令dotnet add package Triangle.NET安装完成后在代码里引入命名空间using TriangleNet; using TriangleNet.Geometry; using TriangleNet.Meshing;这套库的依赖非常少不像有些几何库会拖进一堆底层依赖这一点我很喜欢。装完之后直接就能跑不需要额外配置任何环境变量或者原生DLL。提示如果你在NuGet上看到多个相似名称的包认准Triangle.NET这个官方包名。有个别第三方封装的版本更新并不及时接口可能会和当前版本对不上最好以官方发布的为准。2.2 版本选择与项目配置建议我使用过程中接触到的版本主要是0.8.x系列这个系列已经相当稳定API设计也基本定型了。如果你打开一个老项目里面用的可能是0.7.x之类的旧版本那接口上会有一些差异。例如旧版本里通过Vertex类直接构造点新版本则在Geometry命名空间下用Vertex整体结构调整不大但最好还是以你实际引用的版本为准。从目标框架的角度讲Triangle.Net支持.NET Standard 2.0这意味着net472、netcoreapp3.1、net5.0及以上的项目都能直接引用兼容性做得相当好。我自己有一种情况是在Unity里用这个库做程序化网格生成通过.NET Standard 2.0的兼容层也可以正常调用不过Unity的IL2CPP模式下需要稍微注意AOT兼容性这个后面讲坑点的时候会展开说。配置上我只有一个建议把网格生成的逻辑单独封装成一个服务类不要散落在业务代码里。原因很直接Triangle.Net的Mesh对象有比较复杂的内部状态如果你在多个地方直接操作它很容易出现“一个地方改了顶点另一个地方没同步”的问题。我在项目里一般会封装成类似TriangleMeshGenerator这样的静态工具类只暴露输入几何数据和输出网格参数的接口内部细节全部隐藏起来。2.3 快速验证环境是否可用装完包之后第一件事就是跑一个最简单的剖分Demo确认环境可用。我自己习惯用矩形区域加几条对角线来测试因为预期的剖分结果比较直观。using TriangleNet; using TriangleNet.Geometry; var polygon new Polygon(); // 添加四个顶点 polygon.Add(new Vertex(0.0, 0.0)); polygon.Add(new Vertex(10.0, 0.0)); polygon.Add(new Vertex(10.0, 10.0)); polygon.Add(new Vertex(0.0, 10.0)); // 添加外边界环逆时针方向 polygon.Add(new Contour(new[] { new Vertex(0.0, 0.0), new Vertex(10.0, 0.0), new Vertex(10.0, 10.0), new Vertex(0.0, 10.0) })); var mesh polygon.Triangulate(); Console.WriteLine($三角形数量: {mesh.Triangles.Count}); Console.WriteLine($顶点数量: {mesh.Vertices.Count});如果你能正常看到输出的三角形数量和顶点数量说明这套环境已经可以正式开工了。这里三角形数量通常是2个因为四个点构成两个三角形这也是最简单的Delaunay剖分结果。注意Polygon的Contour添加顺序会影响边界方向和孔洞识别。外边界如果用顺时针Triangle.Net也能处理但后续在判断孔洞、计算面积时会因为符号问题出各种莫名其妙的结果。我在实际使用中踩过这个坑后来就固定用逆时针表示外边界孔洞则用顺时针这样语义清晰且符合库的默认处理逻辑。3. 核心细节解析与实操要点3.1 核心数据结构Vertex、Polygon、Contour与Mesh在深入了解Triangle.Net之前需要对它的核心数据结构有一个清晰的认识。很多初学者容易混淆的是Polygon和Contour这两个概念其实它们的定位完全不同。Vertex顶点类包含X、Y坐标和一个可选的ID。在Triangle.Net中顶点还可以附带标记Label方便在后期把不同区域的三角形区分开来。Contour轮廓线也就是一条由多个顶点组成的闭合边界。外边界和孔洞都是通过Contour来表达的。Polygon多边形区域由若干条Contour组成同时还可以包含孔洞Hole和区域标记Region。Mesh剖分结果包含三角形、边和顶点的集合。这是你最终拿到的网格数据。把它们组合起来的方式是先把所有输入条件加到Polygon里然后调用Triangulate或者自定义配置进行三角剖分。var polygon new Polygon(); // 把边界轮廓加进去 var boundary new Contour(vertices, 0); polygon.Add(boundary); // 如果有孔洞就再添加一个顺时针的Contour var hole new Contour(holeVertices, 1); polygon.Add(hole, true); // 第二个参数标记这个轮廓是孔洞这里第二个参数为true表示这个Contour是孔洞区域Triangle.Net会保证剖分结果不覆盖孔洞内部。这一点在带孔机械零件、地形图中带湖泊区域等场景非常实用。3.2 约束Delaunay三角剖分的工作原理简述我之前提到Triangle.Net默认做的事情是约束Delaunay三角剖分CDT。这里稍微展开说下CDT和无约束Delaunay的区别因为这决定了你会拿到什么样的网格。无约束Delaunay三角剖分只针对点集不考虑任何边界线段。它的剖分结果会生成覆盖点集凸包的三角形也就是说如果你的点集构成的是一个L形区域无约束剖分会把L形凹进去的部分也填上三角形这显然不是我们想要的。约束Delaunay三角剖分则会把边界线段作为约束边强行保留在网格中三角形不能跨越这些约束边这样就可以生成贴合真实边界的网格。Triangle.Net内部会在Delaunay剖分之后把约束边逐条插入并删除与约束边相交的三角形再局部重建网格直到所有约束边都存在为止。这也是为什么我强调输入边界一定要闭合且方向正确。如果Contour没有真正闭合或者边界之间出现交叉CDT阶段处理约束边插入时就会出现异常或者生成不规则网格。3.3 网格细化的核心参数解析Triangle.Net的网格细化功能非常强大它允许你通过QualityOptions来控制三角形的大小、形状和整体质量。我用得最多的三个参数是MaxArea最大三角形面积。面积超过这个值的三角形会被细分。MinAngle最小角阈值。任何小于这个角度的三角形都会被分解重建。MaxAngle最大角阈值。用来避免钝角三角形过多但实际中使用频率低于MinAngle。使用方式如下var quality new QualityOptions { MinimumAngle 20.0, MaximumArea 0.5 }; var mesh polygon.Triangulate(quality);我在一个河道网格生成项目中就同时用了这两个参数。河道形状狭长如果只用MaximumArea会导致细长三角形大量出现加上MinimumAngle之后网格质量肉眼可见地提升有限元计算收敛性也好了很多。关于MaximumArea的取值我自己习惯先用输入区域的总面积除以预期目标三角形数量得到一个粗略的初始值然后再根据结果微调。比如区域总面积是100你希望最终有大约1000个三角形那么初始MaximumArea可以设为0.1左右然后看实际生成的三角形数量再往目标方向调整。3.4 如何从Mesh中提取三角形数据三角剖分完成后关键的产出就是从Mesh中读取三角形数据。Triangle.Net的Triangle对象提供了GetVertex(int index)方法可以获取三个顶点的坐标。最常见的做法是遍历所有三角形输出每个三角形的三个顶点坐标。foreach (var triangle in mesh.Triangles) { var v0 triangle.GetVertex(0); var v1 triangle.GetVertex(1); var v2 triangle.GetVertex(2); Console.WriteLine(${v0.X},{v0.Y} | {v1.X},{v1.Y} | {v2.X},{v2.Y}); }如果你要把网格传给渲染管线或者有限元求解器通常需要把顶点缓冲区和索引缓冲区分别组织起来。一个高效的做法是使用网格本身的顶点列表和三角形顶点IDvar vertices mesh.Vertices.ToArray(); var indices new Listint(); foreach (var triangle in mesh.Triangles) { indices.Add(triangle.GetVertexID(0)); indices.Add(triangle.GetVertexID(1)); indices.Add(triangle.GetVertexID(2)); }用VertexID而不是坐标拼接可以避免大量重复顶点数据让网格数据体积更小也便于后续做顶点法线、UV等属性计算时保持索引稳定。4. 实操过程与核心环节实现4.1 一个完整的带孔多边形剖分案例为了把上面的内容串起来我提供一个完整的案例代码示意如何对一个带孔的多边形进行约束Delaunay三角剖分。这类输入形态在地图导航、几何建模、CAD数据预处理中非常常见。假设我们要剖分的形状是一个正方形边框中间挖去一个圆形孔这里用多边形近似表示圆。using System; using System.Collections.Generic; using System.Linq; using TriangleNet; using TriangleNet.Geometry; using TriangleNet.Meshing; class Program { static void Main(string[] args) { var polygon new Polygon(); // 外边界正方形 var outerVertices new ListVertex { new Vertex(0, 0), new Vertex(100, 0), new Vertex(100, 100), new Vertex(0, 100) }; polygon.Add(new Contour(outerVertices), false); // 内孔用一个八边形模拟圆形孔 var holeCenterX 50.0; var holeCenterY 50.0; var radius 20.0; var holeVertices new ListVertex(); for (int i 0; i 8; i) { var angle i * Math.PI / 4.0; holeVertices.Add(new Vertex( holeCenterX radius * Math.Cos(angle), holeCenterY radius * Math.Sin(angle) )); } polygon.Add(new Contour(holeVertices), true); // 设置网格质量参数 var quality new QualityOptions { MinimumAngle 20.0, MaximumArea 5.0 }; // 执行剖分 var mesh polygon.Triangulate(quality); // 输出结果统计 Console.WriteLine($顶点数: {mesh.Vertices.Count}); Console.WriteLine($三角形数: {mesh.Triangles.Count}); // 将网格数据导出示例只做遍历打印 foreach (var tri in mesh.Triangles) { var v0 tri.GetVertex(0); var v1 tri.GetVertex(1); var v2 tri.GetVertex(2); Console.WriteLine(${v0.X:F4},{v0.Y:F4};{v1.X:F4},{v1.Y:F4};{v2.X:F4},{v2.Y:F4}); } } }这里比较关键的点是外边界和孔洞的Contour方向。我使用的规则是外边界逆时针孔洞顺时针但Triangle.Net在Add方法中并不强制你按特定方向传它更多是依赖你标记的布尔值来区分。不过方向正确能避免很多后续几何计算上的异常。执行之后你可以看到输出数据中的三角形都集中在正方形除去圆形孔之外的区域孔洞内部没有任何顶点或三角形穿过。这正是约束Delaunay剖分相比无约束剖分的优势所在。4.2 平面点集的Delaunay三角剖分如果你的输入只是一个点集没有明确的边界约束那么直接用PointSet或者直接把顶点加到Polygon里不加Contour即可。比如我们对100个随机点做Delaunay剖分var random new Random(42); var polygon new Polygon(); for (int i 0; i 100; i) { var x random.NextDouble() * 100.0; var y random.NextDouble() * 100.0; polygon.Add(new Vertex(x, y)); } var mesh polygon.Triangulate(); Console.WriteLine($点数: {mesh.Vertices.Count}, 三角形数: {mesh.Triangles.Count});这种情况下输出网格覆盖的是点集的凸包边界就是凸包边界。如果你需要生成的是凹包或者更贴近实际地形的边界就需要额外提供边界约束否则凹进去的区域一定会被多余的三角形填充。注意无约束点集剖分虽然简单但在做插值、邻域搜索等任务时三角形连接关系可能不够稳定。输入点一旦有微小扰动三角形拓扑可能会发生很大变化。所以在需要稳定网格结构的场景里我会优先考虑带约束边或者带边界区域的剖分。4.3 利用Region实现区域标记与加密Triangle.Net有一个比较隐蔽但很实用的功能——区域标记Region。在有限元或者工程计算中经常需要对不同材料属性的区域分别标记或者对不同区域采用不同的网格密度。使用Region的方式是给Polygon添加RegionPointpolygon.Add(new RegionPoint(25, 25, 1)); // 标记区域1 polygon.Add(new RegionPoint(75, 75, 2)); // 标记区域2剖分完成后每个三角形都会有一个Region属性表示它属于哪个区域。这个功能最大的价值在于不需要自己去判断三角形中心落在哪个多边形内Triangle.Net在剖分时就帮你把区域归属算好了。结合网格加密需求我还可以根据区域ID设置不同的MaximumArea。常见的做法是先生成基础网格再对特定区域的三角形进行局部细化但这个操作Triangle.Net没有直接的Per-Region细化接口通常的做法是在配置里使用统一的MaximumArea然后在后处理时对特定区域做二次细分或者直接用更细的全局参数然后合并小区域网格。实测中先全局粗剖分、局部再加密的思路虽然多一步操作但可控性更强。4.4 与其他几何库的对接技巧在实际项目中你不大可能从零构建Vertex列表几何数据往往是来自CAD导出的DXF、GIS的Shapefile或者是自己引擎里的模型数据。这时就需要把外部几何数据转换成Triangle.Net可以接收的格式。我从Shapefile导入多边形的时候通常会把外环和內环拆开由于Shapefile中多边形的方向约定和Triangle.Net可能不一致我会统一做一次方向矫正再传给Polygon。这一步很关键否则边界约束识别会出现问题。与Unity对接时我一般会把Triangle.Net生成的三角形数据转成Mesh的顶点数组和三角形索引数组var unityVertices new Vector3[mesh.Vertices.Count]; var unityTriangles new int[mesh.Triangles.Count * 3]; int vertexIndex 0; foreach (var vertex in mesh.Vertices) { unityVertices[vertexIndex] new Vector3((float)vertex.X, 0f, (float)vertex.Y); } int triIndex 0; foreach (var tri in mesh.Triangles) { unityTriangles[triIndex] tri.GetVertexID(0); unityTriangles[triIndex] tri.GetVertexID(1); unityTriangles[triIndex] tri.GetVertexID(2); }有一点要特别注意Triangle.Net的顶点坐标是二维的如果你要在三维空间里做平面剖分需要先把三维顶点投影到二维平面剖分完成后再映射回三维。这一步在建筑楼板、地形剖面等场景很容易忽略一旦忘了投影生成的网格在三维空间里就会出现折叠和穿插。5. 常见问题与排查技巧实录5.1 自定义高程点导致边界错乱我有一个项目是做地形三角网格生成输入是散乱的高程点外加一个约定好的边界多边形。刚开始我直接把高程点全部作为Vertex丢给Polygon然后再添加边界Contour。按理说边界Contour会被保留为约束边结果实际输出时发现边界附近的三角形经常出现“越界”现象也就是三角形的边穿过了边界多边形。排查了很久最终定位到问题所在那些落在边界多边形内部的顶点虽然参与了剖分但因为它们在约束边外侧形成了“孤岛”或者凸包延伸导致约束边插入时发生了局部拓扑冲突。解决方案很简单——在添加顶点进Polygon之前过滤一遍只保留落在边界多边形内部的点。几何判断用Point-in-Polygon算法即可。加了这层过滤之后边界错乱的问题就再也没出现过。5.2 顶点点位非常靠近时导致剖分异常还有一种情况比较隐蔽就是输入数据里存在几乎重合但又不完全相同的顶点比如两个顶点之间的距离只有1e-9量级。Triangle.Net对这些极小间距的处理并不总是那么“宽容”有时候会生成退化三角形或者报出“ConformingMesh failed”之类的异常。针对这种情况我在做GIS数据转换时会对顶点坐标做一次去重和聚类尺度太近的点直接合并成一个。具体的阈值根据你的数据精度来定比如坐标单位是米时我通常用1e-6米作为去重阈值。这一步看起来是数据预处理实际上能省掉大量后续排查的麻烦。5.3 大点数模型的性能优化方向如果你要对数十万甚至上百万个点做剖分性能就需要认真对待了。Triangle.Net本身的算法复杂度在大多数情况下表现不错但实际项目里瓶颈往往出在输入数据的处理和输出网格的存储上。我使用的优化手段主要有几条输入顶点用数组而不是List能减少内存碎片和扩容开销。输出数据直接复用预先分配的数组避免不断new新对象。如果场景允许优先使用不带网格细化的快速剖分等做完区域筛选后再对局部做细化避免一次性生成海量微小三角形。实测下来一个10万点级别的点集在普通的开发机上使用无约束剖分可以在几秒内完成但加上严格的网格细化后时间可能飙升至几十秒。这个差距很多时候不是算法的锅而是三角形质量约束本身带来的计算成本。5.4 与Unity等引擎集成时的兼容性问题之前提到过Unity集成的问题。Triangle.Net的核心代码是纯C#的IL2CPP模式下理论上没有问题但如果你在编辑器模式下正常运行打包以后却出现莫名其妙的行为崩溃多半和AOT裁剪有关。三角剖分里用到了不少基于接口的调用和反射遍历如果裁剪策略过于激进部分类型可能被剥离。我的应对方式是在代码里用Link.xml或Unity的Preserve特性把TriangleNet命名空间下的关键类型保留下来确保AOT编译时不会丢失。具体保留哪些类型没有固定清单建议先从Mesh、Triangle、Vertex这几个核心类开始跑一次打包测试看是否报错再逐步补充。此外如果在Unity的Job System或Burst环境下使用Triangle.Net需要特别注意它内部使用了List、Dictionary这类托管容器不能直接在Job中调用。我的做法是先把三角剖分放到主线程做生成Mesh数据后把顶点和索引复制到NativeArray再交给Job System处理后续的计算。这样既保住了剖分稳定性又不影响并行性能。5.5 常见问题速查表为了便于日常查找我把常见的几类问题整理成了一个速查表基本覆盖了我实际开发过程中遇到的高频问题异常现象典型原因解决方案三角剖分后边界被穿透边界Contour未正确闭合或方向混乱检查边界轮廓点序列确保首尾闭合统一外边界逆时针、孔洞顺时针生成的三角形过于细长缺少最小角约束或MinimumAngle设置过小设置QualityOptions.MinimumAngle为20~30度必要时提高剖分结果覆盖了孔洞区域孔洞Contour未标记为孔洞或输入的孔洞轮廓方向错误在polygon.Add(contour)时第二参数传true并检查孔洞方向报错ConformingMesh failed输入存在退化点、重复点或极端细长约束边预处理顶点去重检查约束边是否相交或长度过小性能卡顿严重点数过大且开启严格网格细化适当放宽MaximumArea或者先粗剖分后局部细化在Unity打包后崩溃IL2CPP裁剪掉了TriangleNet类型保留关键类型避免在Job/Burst环境直接调用除了表格里的这些还有一个主观性很强的建议拿到任何新的几何数据先做一次“可视化检查”。哪怕只是把剖分结果输出成OBJ导入MeshLab看一眼都比盯着坐标数据猜问题高效得多。我自己有段时间排查一个拓扑错误怎么都对不上号最后把网格导出来一看才发现是CAD数据里隐藏了一条肉眼几乎看不见的重复边。结尾一些使用上的个人体会最后再聊一点实在的东西。Triangle.Net的文档不多社区资料也远不如一些大型几何库丰富所以遇到问题的时候大部分时间我只能去看源码和跑测试用例来确认行为。但也正因为这样我对这个库的内部机制理解反而更深了。这里想补充一个经验如果你打算在比较复杂的业务项目里用Triangle.Net建议花点时间读一读它的源码至少把Mesh和Triangulate的调用链走一遍。这个库的代码风格清晰没有太多花哨的技巧读完以后你对剖分过程中什么时候会发生顶点插入、什么时候网格会做局部重建会有更具体的把握排查问题时心里就有底了。另外关于网格生成这件事选库只是第一步。真正决定最终网格质量的往往是你对业务数据的预处理是否到位边界是否干净、孔洞是否闭合、极小间距是否合并、区域标记是否合理。这些脏活累活做不到位再好的剖分库也救不回来。希望这篇文章能帮你少走一些弯路至少在下载下来、跑通第一个案例之后不觉得Triangle.Net是个黑盒。
