企业官网建设流程全解析

简介：VTK与CGAL联合实现的3D模型高级处理代码包，面向从事三维可视化、网格处理与逆向工程的开发者，重点解决布尔运算、三角剖分、补洞、模型偏置（膨胀）、降采样和骨架提取等常见需求。压缩包共2个文件，含1个头文件与1个实现文件，头文件声明接口与类，实现文件承载算法逻辑，整体仅6KB，代码量小但封装完整，适合已配置VTK与CGAL5.5环境、希望按需抽取核心算法的中高级开发者阅读。已有2041人学习下载，文件结构简洁，接口划分清晰，可快速定位到对应功能模块。代码基于CGAL与VTK接口封装，将底层几何算法组织成便于调用的函数与类，使用者无需从零编写复杂的计算几何逻辑，即可在已有工程中完成两个模型求交、求并、求差，将散乱点云或破损网格重建为高质量三角网格，并进一步实现孔洞修复、厚度调整及结构骨架抽取。整体为模型修复、简化与结构分析提供了可运行、易扩展的参考实现，可直接集成到自有3D处理流程中。

1. 从VTK到CGAL：为什么网格处理要两套库协作

做逆向设计、医疗建模或三维打印预处理时，很多人和我一样，第一个拿到的工具是VTK。VTK读写STL、显示模型、做简单滤波都顺手，但真到布尔运算、三角剖分、补洞、模型偏置（膨胀）和骨架提取这五个环节，VTK自带组件往往不够用——vtkBooleanOperationPolyDataFilter遇到共面或退化三角形就翻车，vtkFillHolesFilter只能补个平面呢，膨胀和骨架更是没有现成接口。某开发者做植入体设计时，卡在“给骨头模型布尔减孔后再补洞膨胀”这一步整整两周，换CGAL之后一个下午跑通。CGAL是计算几何算法库，布尔用Nef多面体或corefinement，补洞有专门的重建算法，骨架提取有基于平均曲率流的标准实现。本文把这套“VTK做I/O和显示、CGAL做计算、转回VTK渲染”的流程拆开讲，适合想解决实际网格处理问题、又不想被库内部细节劝退的从业者。

2. 布尔运算与三角剖分：VTK读入、CGAL算、再回VTK显示

2.1 数据流设计：内存在里的VTK PolyData怎么交给CGAL

VTK用vtkPolyData存储点、线、面，底层是隐式共享的拓扑结构；CGAL的Surface_mesh则基于半边数据结构，能显式表示邻接关系。两套模型不能互相替代，只能转换。

我一般会写一个转换函数，把vtkPolyData拷贝成CGAL::Surface_mesh。关键点是：VTK的面片是三角形或四边形，CGAL的Face是有方向的顶点序列，所以转换时要保持顶点索引一致。

#include <vtkPolyData.h> #include <vtkCellArray.h> #include <CGAL/Surface_mesh.h> using K = CGAL::Exact_predicates_inexact_constructions_kernel; using Mesh = CGAL::Surface_mesh<K::Point_3>; using Vertex_index = Mesh::Vertex_index; using Face_index = Mesh::Face_index; bool vtkPolyDataToSurfaceMesh(vtkPolyData* pd, Mesh& mesh) { vtkPoints* points = pd->GetPoints(); if (!points || points->GetNumberOfPoints() == 0) return false; // 1. 先加入所有顶点，记录VTK点索引到CGAL顶点索引的映射 std::vector<Vertex_index> vtkToCgal; vtkToCgal.resize(points->GetNumberOfPoints()); for (vtkIdType i = 0; i < points->GetNumberOfPoints(); ++i) { double p[3]; points->GetPoint(i, p); Vertex_index vi = mesh.add_vertex(K::Point_3(p[0], p[1], p[2])); vtkToCgal[i] = vi; } // 2. 遍历所有多边形单元，构建面 vtkCellArray* polys = pd->GetPolys(); vtkIdType npts = 0; const vtkIdType* pts = nullptr; polys->InitTraversal(); while (polys->GetNextCell(npts, pts)) { std::vector<Vertex_index> face_vids; for (vtkIdType j = 0; j < npts; ++j) { face_vids.push_back(vtkToCgal[pts[j]]); } mesh.add_face(face_vids); } return true; }

逻辑说明：vtkPolyData的GetPolys()返回的是三角形和四边形单元，但实际可能是任意凸多边形。这里按面单元逐个add_face，CGAL会检查面的方向，如果相邻面的方向不一致，add_face可能返回nullptr，表示网格非法。

参数说明：vtkToCgal映射数组是必须的，因为VTK的点索引是连续的，但CGAL的顶点索引在add_vertex时自动分配，不能假设一致。更稳妥的做法是先用std::unordered_map<double*, Vertex_index>做哈希，避免重复点，但简单场景下线性映射够了。转换后要检查mesh.is_valid()，如果返回false，多半是面方向杂乱或者有重复边。

2.2 布尔运算流程与参数说明

CGAL做布尔运算有两条路线：老牌的是Nef_polyhedron_3，它基于Nef多面体表示，能精确表达布尔结果；新推荐的是Polygon_mesh_processing包里的corefine_and_compute_union/difference/intersection，这套操作直接在Surface_mesh上工作，速度快且支持非凸模型。

我用的是corefinement方案，它要求输入的两个网格都是闭合、无边界的三角网格。运算前一定要先验证：

#include <CGAL/Polygon_mesh_processing/corefinement.h> #include <CGAL/Polygon_mesh_processing/repair.h> namespace PMP = CGAL::Polygon_mesh_processing; // 假设mesh1和mesh2都已从VTK转换而来 bool is_closed1 = PMP::is_closed(mesh1); bool is_closed2 = PMP::is_closed(mesh2); if (!is_closed1 || !is_closed2) { // 先补洞，见下一章 return false; } // 求并集 bool success = PMP::corefine_and_compute_union(mesh1, mesh2, result); if (!success) { // 求差集 // success = PMP::corefine_and_compute_difference(mesh1, mesh2, result); }

逻辑说明：corefine_and_compute_union会原地修改mesh1和mesh2，把它们相交的边切碎、交换连接性，然后生成result。这个操作不是复制输入，所以原始网格会被改变；如果还需要保留原始数据，必须先复制一份。

参数说明：内核Exact_predicates_inexact_constructions_kernel是这里最常用组合。它保证几何判断（哪个点在哪个面内、边是否相交）使用精确谓词，避免浮点误差导致拓扑错误；构造点的坐标仍然是近似值，但布尔运算结果不会出现播放器那种“裂缝”。如果模型非常大或自相交，可能需要换Exact_constructions_kernel，代价是速度慢很多。

另一个容易踩的点：两个网格的空间位置必须预先对齐。CGAL布尔运算不判断“接触”，它只处理面之间的交叉，如果两个网格完全分离，并集就是简单的把mesh2的面加到mesh1上，但corefinement要求输入最好有交集，否则结果可能包含退化面。

2.3 三角剖分：约束Delaunay与边界保留

布尔运算返回的result可能包含多边形面（四边形、五边形），因为corefinement会在相交处形成新的平面边界。VTK渲染不喜欢非三角形面，很多算法（如骨架提取）也明确要求三角网格，所以需要把多边形面三角剖分。

CGAL的Polygon_mesh_processing::triangulate_faces可以直接完成：

#include <CGAL/Polygon_mesh_processing/triangulate_faces.h> std::size_t n_fails = PMP::triangulate_faces(result); std::cout << "Failed to triangulate faces: " << n_fails << std::endl;

逻辑说明：这个函数遍历所有非三角形面，以面的顶点为输入构造一个简单多边形，然后使用保形Delaunay三角化内部三角剖分。它不改变面的边界顶点，所以保留原多边形的轮廓——这正是布尔运算后想要的效果。

参数说明：triangulate_faces可以接受可选参数，比如CGAL::parameters::angle_max和aspect_ratio_max限制生成三角形的质量。默认值是180度（无限制），但在工程中我发现设置angle_max = 60可以避免长出极扁的钝角三角形，这些钝角三角形在后续偏置或骨架计算中会引发数值问题。如果剖分失败，n_fails返回失败面数，常见原因是面顶点不共面、自相交或者顶点索引不连续。这时先用PMP::repair清理一遍网格再剖分。

3. 补洞与模型偏置（膨胀）：从修复到生成厚度

3.1 补洞：孔洞识别与CGAL补洞算法

扫描得到的模型或布尔运算产生的结果，常常有一些开放边界。CGAL把网格中的“洞”定义为一组形成环的边界半边。补洞的目标是生成一个曲面片，使环闭合，并且新生成的三角面与周围网格平滑过渡。

我自己写了一个识别洞的函数：遍历所有半边，如果一条半边没有对面（halfedge的opposite是空），那它就属于边界环。一个环通常可以沿next指针走一圈收集所有边界极小半边。

#include <CGAL/Polygon_mesh_processing/repair.h> #include <CGAL/Polygon_mesh_processing/orientation.h> int fill_holes(Mesh& mesh) { std::vector<Mesh::Halfedge_index> border_cycles; PMP::extract_boundary_cycles(mesh, std::back_inserter(border_cycles)); int filled = 0; for (auto h : border_cycles) { // 记录洞的边界顶点数 int boundary_size = 0; auto hh = h; do { ++boundary_size; hh = mesh.next(hh); } while (hh != h); // 如果环太小，跳过 if (boundary_size < 3) continue; // 对单个洞进行补洞 bool ok = PMP::triangulate_hole(mesh, h, PMP::parameters::vertex_point_map(mesh.points())); if (ok) ++filled; } return filled; }

逻辑说明：PMP::extract_boundary_cycles返回所有边界环的半边索引。每个环以半边为起点，用mesh.next(hh)绕过一周。triangulate_hole以这个环为边界，在内部生成一个三角剖分补丁。注意它只补单个洞，不自动识别嵌套关系。

参数说明：triangulate_hole还支持use_delaunay_triangulation（是否用Delaunay提高三角形质量）和max_edge_length（补洞块最长边长）。我会设置max_edge_length为模型平均边长的2到3倍，否则遇到很大的洞，生成的补丁会有超长边，显示时像被扯开的膜。补洞前最好先看看洞的边界环方向是否与周围面法向一致。如果洞的方向反了，补出来的面法向会朝内，渲染时前面看起来是黑的。可以用PMP::is_outward_oriented(mesh)检查整体方向，反转后再补。

3.2 模型偏置（膨胀）：让模型长出厚度

模型偏置也叫膨胀，目的让模型表面沿法向外扩，比如给内骨骼模型生成植入体外壳。VTK没有直接提供偏置滤波器。CGAL也没有一个函数叫“offset”，但工程上最通用的做法是基于顶点法向偏移：

#include <CGAL/Polygon_mesh_processing/compute_normal.h> void inflate_mesh(Mesh& mesh, double offset, int iterations) { for (int iter = 0; iter < iterations; ++iter) { // 1. 每个顶点法向 std::vector<K::Vector_3> normals(mesh.num_vertices(), K::Vector_3(0,0,0)); PMP::compute_vertex_normals(mesh, CGAL::parameters::vertex_normal_map(boost::make_assoc_property_map(normals))); // 2. 记录偏移后位置 std::vector<K::Point_3> new_positions; new_positions.reserve(mesh.num_vertices()); for (auto v : mesh.vertices()) { K::Vector_3 n = normals[v]; double len = sqrt(n.squared_length()); if (len == 0) { new_positions.push_back(mesh.point(v)); continue; } // 归一化 K::Vector_3 unit_n = n / len; // 偏移 new_positions.push_back(mesh.point(v) + unit_n * offset); } // 3. 更新顶点位置 auto v_it = mesh.vertices().begin(); auto p_it = new_positions.begin(); for (; v_it != mesh.vertices().end(); ++v_it, ++p_it) { mesh.point(*v_it) = *p_it; } } }

逻辑说明：顶点法向是周围三角形的面法向按面积加权得到的。对每个顶点沿法向移动offset单位，迭代多次会让模型整体膨胀。每次迭代后顶点位置改变，法向也随之变化，所以必须每轮重新计算法向。

参数说明：offset表示单次偏移距离，iterations表示重复次数。两种用法：一种是一步到位iterations=1，另一种是小步多次（比如offset=0.1, iterations=10）模拟平滑偏移。关键陷阱是：凹区域的顶点法置信度低，沿法向偏移容易让相邻面互相穿插导致自相交。我通常会在每次迭代后用PMP::remove_self_intersections(mesh)清理交替面，这个函数把自相交区域切开、重建面片。代价是速度慢，但比后面渲染出错好。

更保守的方案是Minkowski和：用一个半径很小的球对模型做球扫掠。CGAL的Mesh_3包可以生成这个结果，但搭建很麻烦，普通制作周期内不建议碰。

3.3 参数调优与失败时的现象

一个真实失败案例：某开发者对一块薄板模型膨胀，offset=0.5、iterations=1，结果中间薄壁区域出现一个大空洞，原因是该区域法向在两侧高度变化，顶点偏移后变成Z字形，自相交清理把整条窄带删掉了。后来改成iterations=5, offset=0.1，并在每步之间做一次各向同性重网格化，问题才消失。所以“膨胀”不是一次性整形，要当作迭代过程来观察。

4. VTK常见问题排查：坐标、法向、多边形走向

这一章整理我自己在处理CGAL到VTK往返时踩过的四个主要坑，每条都是现象、原因、解决三步走。

4.1 现象一：布尔运算结果转回VTK后，模型表面有大片黑色

现象：在VTK窗口里渲染并集或差集，结果模型表面有奇怪的黑色区域，旋转视角时黑斑在亮暗间闪烁。这是典型的法向不一致渲染问题。

原因：CGAL的Surface_mesh有统一的面方向约定，但VTK的vtkPolyData默认使用的vtkPolyDataNormals滤波按单元顺序计算法向。如果从CGAL转回VTK时，多边形走向没有按右手定则统一，某些面法向会相反，导致渲染器认为这些面是背面。

解决：转换后用vtkPolyDataNormals强制计算并统一法向。注意要设置ConsistencyOn和SplittingOff，否则相邻面会被错误拆分。

vtkSmartPointer<vtkPolyDataNormals> normals = vtkSmartPointer<vtkPolyDataNormals>::New(); normals->SetInputData(vtkOutputMesh); normals->ConsistencyOn(); normals->SplittingOff(); normals->Update();

逻辑说明：ConsistencyOn会沿着网格拓扑传播法向方向，确保相邻面的法向朝向同一侧。SplittingOff防止因为法向不一致而在重合顶点上产生分裂点。这是VTK渲染闭环里的必做项。

解决后如果还有黑面，多半是原始网格本身包含退化面，需要用PMP::remove_degenerate_faces提前清理。

4.2 现象二：补洞结果与周围表面有可见折痕

现象：补出来的洞片，颜色和周围一致，但从侧视图看，补片比周围表面高出或低一块，像鼓起一个包。

原因：triangulate_hole只保证边界拓扑闭合，不保证新生成三角形与邻接面共平面。特别当洞边界是一个空间多边形时，三角剖分选出的点是边界上的，但连接方式不同，曲面补片的弯曲方向就可能与周围不一致。

解决：不用triangulate_hole直接生成一个patch，而是把洞先用PMP::triangulate_hole铺上三角形，再用PMP::smooth_shape对补片的内部顶点做平滑迭代。准确说，对洞边界上的相邻顶点做双调和变形。我用的步骤是：

1. 调用triangulate_hole得到初步补丁。
2. 对补丁内部顶点执行10～20次拉普拉斯平滑，边界顶点固定不动。

PMP::smooth_shape(mesh, CGAL::parameters::number_of_iterations(20) .vertex_is_constrained_map(vertex_constrained_map));

参数说明：vertex_constrained_map是一个std::map，把洞边界顶点标记为’true’不参与平滑，内部顶点标记为’false’。number_of_iterations太多会让补片塌缩，太少折痕还明显，20次是我常用的中位数。

4.3 现象三：膨胀后模型体积变小了

现象：明明沿法向向外偏移，结果模型的总体积反而比原始模型小，尤其对凹的球壳和管道模型。

原因：顶点法向是局部表面的平均朝向，对于凹区域，法向指向包围盒内部（因为开口形状的边缘处法向指向中心线的垂线），导致顶点被推向内部。

解决：在膨胀前先补好所有大开孔，让模型具备“整体外凸”的流形结构。如果只是薄壁件，改用中面偏移——先把模型抽成中面再用侧向厚度方向扩展。更简单的方法是：膨胀后检查模型是否满足PMP::volume(mesh) > 0。如果体积变小，把offset变号向内收缩一次，再与原始模型做并集。

4.4 现象四：降采样后尖角消失或拓扑断裂

现象：用vtkDecimatePro把面数从50万降到5万，结果原本锋利的棱边被圆角化，有时还出现破洞。

原因：vtkDecimatePro默认的PreserveEdges选项对非流形边的处理不够保守，它按照几何误差来合并顶点，但不会检查是否消除特征边。CGAL的edge_collapse也有类似问题，但可以通过设置停止条件来限制。

解决：改在CGAL里做edge_collapse，把特征边标记为不可折叠。具体做法：先用PMP::sharp_edges_segmentation检测二面角大于某阈值的边，生成边标记，再传给edge_collapse作为约束。

PMP::edge_collapse(mesh, stop_condition, CGAL::parameters::edge_is_constrained_map(edge_is_constrained_map));

参数说明：stop_condition我用CGAL::Surface_mesh_simplification::Edge_length_ratio，当边长度与包围盒对角线之比小于0.01时停止。edge_is_constrained_map标记尖锐边，防止折叠跨越特征线。这样降采样结果既能减面，又保留棱边特征，回到VTK时尖角依然清晰。

5. 降采样与骨架提取：从点云到中线结构

5.1 降采样：在VTK做还是CGAL做？

VTK的vtkDecimatePro速度极快，适合交互式预览，但它不给你精细控制。CGAL的Surface_mesh_simplification::edge_collapse则有丰富的停止条件和约束映射，更适合对结果有要求的场景（骨架提取、有限元）。

我通常根据模型面数分流：小于10万面直接用VTK，大于10万面或者需要保留特征边时用CGAL。核心调用如下：

namespace SMS = CGAL::Surface_mesh_simplification; SMS::Edge_collapse<Mesh> collapse(mesh); SMS::Count_ratio_stop_ratio stop(0.2); // 保留20%的边 bool ok = collapse.stop_when(stop);

逻辑说明：Edge_collapse是迭代式折叠，每次选择折叠后误差最小的边（Garland-Heckbert二次误差度量），Count_ratio_stop_ratio设置折叠到原始边数20%时停止。过程中不允许改变网格拓扑（不合并非流形边或设标签）。

参数说明：除了Count_ratio_stop_ratio，常用Cost_stop_ratio（当折叠代价超过初始平均代价的倍数时停止）和Edge_length_ratio（边长度比）。如果点云场景，还能用Constrained_vertices_map固定某些关键顶点。

5.2 骨架提取：CGAL的骨架化与VTK可视化

CGAL的Mean_curvature_flow_skeleton把表面模型看成弹性体，不断向内收缩迭代，同时在表面顶点间形成骨架线。这个算法输入必须是三角网格，最好是封闭的。

#include <CGAL/Mean_curvature_flow_skeleton.h> using Skeleton = CGAL::Mean_curvature_flow_skeleton<Mesh>; Skeleton skeleton(mesh); Skeleton::Skeleton skeleton_data; double quality = 0.15; // 骨架质量参数 bool converged = skeleton.contract_until(quality); skeleton.convert_to_skeleton(mesh);

逻辑说明：contract_until是迭代收缩过程，quality控制最终骨架的复杂度。数值越大骨架越稀疏，0.15是比较均衡的值；收敛返回false表示失败，通常因为网格质量太差。convert_to_skeleton把收缩后的网格转换成Skeleton对象，里面包含骨架顶点和边。

输出后，需要转成VTK可显示的线段数据：

vtkSmartPointer<vtkPolyData> skelPoly = vtkSmartPointer<vtkPolyData>::New(); vtkSmartPointer<vtkPoints> skelPoints = vtkSmartPointer<vtkPoints>::New(); vtkSmartPointer<vtkCellArray> skelLines = vtkSmartPointer<vtkCellArray>::New(); for (auto v : skeleton_data.vertices()) { auto p = skeleton_data.vertex_point(v); skelPoints->InsertNextPoint(p.x(), p.y(), p.z()); } // 边转换为线段并插入 skelLines

逻辑说明：骨架顶点位置是浮点坐标，与原始模型坐标一致，所以可以直接放进vtkPoints。线段索引在原骨架中是一个邻接列表，映射到VTK时，把每个顶点索引加进对应线段cell。

参数说明：quality不是线性比例。0.3会让骨架只剩几个大分支，0.05则产生密集的毛发状细枝。对树状血管模型，我推荐从0.2开始逐步调小，观察骨架分支点是否落入模型内部。

5.3 验证：骨架与原始模型的拟合度

骨架提取后必须检查骨架是否真的在模型内部，而不是飘在表面。最简单的验证方法：用vtkSelectEnclosedPoints，把骨架点作为输入，原始模型作为封闭面，判断每个点是否在模型内。

vtkSmartPointer<vtkSelectEnclosedPoints> select = vtkSmartPointer<vtkSelectEnclosedPoints>::New(); select->SetSurface(originalMesh); select->SetInputData(skelPoly); select->Update(); int insideCount = 0; for (int i = 0; i < skelPoly->GetNumberOfPoints(); ++i) { if (select->IsInside(i)) insideCount++; } double ratio = static_cast<double>(insideCount) / skelPoly->GetNumberOfPoints();

逻辑说明：vtkSelectEnclosedPoints基于射线求交法，判断每个点相对封闭曲面的位置。如果ratio低于0.9，说明骨架有相当一部分伸出模型，常见原因是模型有开口或凹陷。这时需要先补洞、再膨胀一点点（让骨架“芯”离表面更远）、重新做骨架提取。

我还会算骨架每个点到模型曲面的最近距离，用vtkImplicitPolyDataDistance，如果距离大于局部半径的某个倍数，就是过大的悬空分支。经验上，距离超过模型平均厚度一半的骨架线都值得警惕。

6. 复合流程：把布尔、补洞、偏置、降采样、骨架提取串成一条管线

6.1 管线编排与内存管理

在实际解决一个复杂模型时，顺序不是随便写出来的。我的标准顺序是：

1. VTK读STL → 转CGAL Surface_mesh
2. 补洞（保证闭合）
3. 布尔运算（在闭合网格上进行）
4. 再次补洞（布尔可能撕开网格）
5. 三角剖分（把多边形面切碎）
6. 偏置（膨胀）
7. 自相交清理
8. 降采样（提高后续骨架效率）
9. 骨架提取

第6步偏置前，网格必须是流形；第8步降采样会关闭一部分顶点，但骨架提取对顶点密度不敏感，只关心拓扑。每两步之间，我会调用一次mesh.is_valid()和PMP::is_closed(mesh)，避免错误沿管线传播。

内存方面，VTK的vtkPolyData用智能指针管理，但CGAL的Surface_mesh是STL容器，没有自动清理机制。当中间结果不再需要时，我直接让它离开作用域；如果保留在堆里，记得在循环末尾mesh.clear()。布尔运算特别消耗内存，两个100万面网格做相交，结果可能有600万面，建议先降采样再布尔。

6.2 一个可复现的最小示例

下面这个函数展示了完整流程的骨架，省略了每个步骤内部的细节函数，但顺序和调用都是实际可跑的。

vtkSmartPointer<vtkPolyData> runPipeline(vtkPolyData* input) { // 1. 转CGAL网格 Mesh mesh; if (!vtkPolyDataToSurfaceMesh(input, mesh)) return nullptr; // 2. 补洞（确保闭合） fill_holes(mesh); // 3. 布尔运算：和另一个网格做差集（假设cutMesh源） Mesh cutMesh; // ... 从cutSource转换 ... if (PMP::is_closed(mesh) && PMP::is_closed(cutMesh)) { Mesh result; if (!PMP::corefine_and_compute_difference(mesh, cutMesh, result)) return nullptr; mesh = std::move(result); } // 4. 再次补洞 fill_holes(mesh); // 5. 三角剖分 PMP::triangulate_faces(mesh); // 6. 偏置（膨胀） inflate_mesh(mesh, 0.02, 5); PMP::remove_self_intersections(mesh); // 7. 降采样 Mesh simplified = mesh; SMS::Count_ratio_stop_ratio stop(0.3); SMS::edge_collapse(simplified, stop); // 8. 骨架提取 CGAL::Mean_curvature_flow_skeleton<Mesh> skeleton(simplified); // ... // 9. 转回VTK vtkSmartPointer<vtkPolyData> resultPoly = vtkSmartPointer<vtkPolyData>::New(); // ... 转换代码 ... return resultPoly; }

逻辑说明：第3步中如果布尔失败，返回nullptr是合理的，因为后续步骤都依赖一个闭合流形。第6步的inflate_mesh和remove_self_intersections是成对出现的，必须放在一起。第7步降采样用edge_collapse直接修改副本，不污染原始网格。

参数说明：偏置的0.02是相对模型包围盒对角线的比例，我在这里用绝对长度单位，读者应替换为自己的模型尺度。降采样保留30%的边，对骨架提取来说足够了，因为骨架不依赖高密顶点。

6.3 性能观察与取舍

对某模拟项目X的一个6万面骨盆模型，布尔运算大约耗时数秒，偏置加自相交清理约占一半时间，骨架提取在降采样后只需要不到一秒。如果把降采样放到最前面，布尔运算会快很多，但布尔结果精度也会下降——所以我的经验是：只在布尔运算前做“轻降采样”（保留70%），骨架前做“重降采样”（保留20%）。

另一个取舍是法向一致性的检测频率。每步都做PMP::is_outward_oriented会花费额外遍历，但能救命。我在一次做差集时，因为输入网格法向颠倒，差集结果变成内部壳体，后续膨胀全部朝向内部，最后骨架追着外壳跑，丢了三个分支。从那以后我每次在转换和修改后都强制走一遍法向统一和自相交检查，哪怕多花十几秒，也比最后视觉翻车强。希望帮到你。

本文还有配套的精品资源，点击获取