在三维几何处理与数字内容创作中，布尔运算是最常用的建模手段之一。但当两个网格体进行交集、差集或并集计算时，经常会遇到所谓“破面”现象：原本应该闭合的实体表面出现裂缝、孔洞甚至完全碎裂成无意义碎片。这类问题在复杂模型、非流形网格或存在细微重叠的几何体上尤为突出。Carve求解器作为一类专注于鲁棒布尔计算的算法实现，通过对交线重建与重叠区域的显式处理，提供了比传统BSP树更稳定的解决路径。
Carve求解器的底层原理与破面成因
传统布尔算法多依赖二进制空间分割（BSP）对网格进行递归切分，在交线处直接裁剪三角形。如果输入网格的法线不一致、存在共面重叠或者浮点精度误差，裁剪过程会产生退化边和零面积三角形，最终在渲染或导出时表现为破面。Carve求解器则采用扫掠平面（sweep plane）配合符号距离场（SDF）的思路，将参与运算的两个网格分别离散为空间中的有向距离信息，再沿坐标轴方向推进平面，动态计算交线环。
这种方法的优势在于，它不要求原始网格满足严格的流形条件，也能在重叠区域自动识别内外关系。Carve在交线生成后会执行一轮拓扑清理，将小于容差阈值的裂缝合并为有效边。从实践看，破面往往不是算法本身出错，而是输入数据带有自交、重复顶点或缩放后未应用的变换矩阵。理解这一点，才能正确使用Carve而非盲目切换求解器。
对比Manifold等基于精确算术的库，Carve更偏向工程鲁棒性。它允许一定的容差松弛，适合处理从第三方软件导入的“脏”模型。但松弛过大又会导致细节丢失，因此重叠处理前的预处理成为关键。下面我们通过代码看如何调用典型Carve接口并完成基础差集。
#include <carve/csg.hpp>
#include <carve/input.hpp>
carve::input::PolyhedronData* make_box(double s) {
carve::input::PolyhedronData* p = new carve::input::PolyhedronData();
p->addVertex(-s, -s, -s);
p->addVertex( s, -s, -s);
p->addVertex( s, s, -s);
p->addVertex(-s, s, -s);
p->addVertex(-s, -s, s);
p->addVertex( s, -s, s);
p->addVertex( s, s, s);
p->addVertex(-s, s, s);
p->addFace(0,1,2,3);
p->addFace(4,5,6,7);
p->addFace(0,1,5,4);
p->addFace(2,3,7,6);
p->addFace(1,2,6,5);
p->addFace(0,4,7,3);
return p;
}
int main() {
carve::input::PolyhedronData* a = make_box(1.0);
carve::input::PolyhedronData* b = make_box(0.6);
// 将b沿z轴抬升产生重叠
b->transform(carve::math::Matrix::TRANSLATION(carve::geom::VECTOR(0,0,0.4)));
carve::csg::CSG csg;
carve::mesh::MeshSet<3>* result = csg.compute(a, b, carve::csg::CSG::DIFFERENCE);
return 0;
}
重叠几何的预处理与容差策略
所谓重叠处理，是指在布尔运算前对参与网格做空间关系修正。常见情况是两个实体在边界上有极小重叠（例如公差0.001内穿插），或存在完全一致的共面多边形。若直接交给Carve，虽能算出结果，但交线会极度密集，引发后续网格简化时的破面。推荐做法是先用包围盒相交检测，提取重叠区，对该区域做轻微偏移或抽壳，使重叠变为明确的内含或分离。
具体参数上，Carve暴露了carve::csg::CSG::Options中的epsilon字段，用于控制距离比较容差。经验值设为模型最大尺寸的千分之一到万分之一。若你的模型单位是米，epsilon取1e-4到1e-3即可。过小会放大小数误差，过大则合并了不该合并的特征。配合重叠处理，可先调用carve::mesh::MeshSet的cleanup方法移除重复顶点。
下面示例展示如何使用容差合并重叠顶点，再进行布尔差集，以避免破面。注意代码中转义了比较符号以符合规范。
import carve_wrapper as cw
def safe_boolean(mesh_a, mesh_b, eps=1e-4):
# 预处理：清理小于eps的重叠顶点
mesh_a.cleanup(epsilon=eps)
mesh_b.cleanup(epsilon=eps)
# 检测包围盒重叠
if not cw.bbox_overlap(mesh_a, mesh_b):
return None
# 设置求解器选项
opt = cw.CSGOptions()
opt.epsilon = eps
csg = cw.CSG(opt)
result = csg.difference(mesh_a, mesh_b)
return result
a = cw.load_mesh('base.obj')
b = cw.load_mesh('cutter.obj')
out = safe_boolean(a, b, eps=0.0005)
if out:
out.save('result.obj')
在Blender等软件中，其布尔修改器背后若选用Carve，同样建议在修改器栈前加一个“重合点”与“重算法线”的操作。这相当于把重叠处理前置到可视化界面里，能直观减少破面发生。
典型破面场景排查与替代方案对比
实际项目中，破面高频出现在三种场景：一是模型从CAD导出时带有非三角化曲面；二是多个布尔级联，前一步输出成为下一步输入，误差累积；三是缩放未应用，世界矩阵与本地矩阵冲突。针对第一种，需先使用triangulate将面统一为三角形，Carve虽支持多边形，但混合面型易出错。第二种应控制级联层数，或每步后做网格有效性校验。
我们将Carve与Manifold、Blender原生Exact做简要对比。Carve胜在容错，处理脏模型不易崩；Manifold基于浮点膨胀与精确谓词，理论更干净但拒绝非流形输入；Exact则慢但稳定。若你频繁遇到破面且模型来源复杂，Carve加重叠预处理是最优性价比方案。下表列出关键差异：
| 求解器 | 容错性 | 速度 | 适用模型 |
|---|---|---|---|
| Carve | 高 | 中 | 脏模型、带重叠 |
| Manifold | 低 | 快 | 干净流形 |
| Exact | 极高 | 慢 | 精密件 |
当Carve仍报破面时，可检查输出日志里的degenerate face数量。若占比超过百分之一，基本是输入重叠区过于混乱，此时应回到重叠处理步骤，对交叠部分做手动分割。掌握这些排查逻辑，布尔破面不再是阻碍交付的难题。