1. 先从一块残缺的扫描件说起我做了几年三维扫描和模型修复被问得最多的问题永远绕不开两个这个模型破了个洞怎么补和补完之后表面皱巴巴的怎么弄平。前者对应补洞hole filling后者对应的就是fairing算法。这两个词听着专业拆开看逻辑其实非常朴素——一个负责解决缺一个负责解决糙。补洞说穿了就是在网格或者点云上找到数据缺失的开口然后想办法把这块缺口用合理的曲面填回去。fairing算法则是给已经存在的曲面做美容把噪声、褶皱、锯齿压下去让表面更顺、更光。你可能觉得这是两件独立的事但在真实的逆向工程、3D打印前处理、数字雕刻、医学建模这些场景里它们几乎是绑在一起出现的你补完洞补丁和原曲面之间必然有接缝和形状突变这时候必须上fairing来过渡而如果你只做fairing不补洞那些破口处的褶皱和尖刺根本无从下手。这篇文章我想按一条真实的处理链路来讲怎么发现洞、怎么判断这个洞该走哪条技术路线、补完以后fairing到底在光顺什么又偷走了什么、两者的衔接处有哪些坑。内容偏实操也会把背后的数学直觉尽量说明白适合做点云处理、网格修复、逆向建模的同行也适合刚开始接触这部分的朋友。不追求把每种算法公式推导一遍但追求讲清楚为什么这么选、参数怎么定、翻车了怎么查。1.1 模型上的洞到底是从哪冒出来的很多人以为洞是扫描没扫到这么简单实际来源多得超乎想象。最常见的是遮挡三维扫描仪从一个视角拍过去物体的凹陷、背面、被其他结构挡住的地方就成了数据盲区重建出来自然是个洞。第二种是数据缺失扫描过程中反光、透明、太薄或者太黑的地方传感器直接拿不到有效点孔洞边缘往往还很毛躁。第三种是格式和拓扑转换带来的模型在不同软件之间导来导去三角面片精度不一致原本闭合的表面被撕开一条缝。第四种是手工编辑的副产品删掉一块错误的面、切了一刀边缘就露出来了。理解洞的来源很重要因为它决定了这个洞的性格。遮挡造成的洞通常形状规整、边界平滑周围曲面质量也不错直接用几何方法补就行补出来的效果一般不会太差。而反光或黑色材质造成的洞边界往往参差不齐周围还伴随大量噪声点这种洞你要是直接补补丁会和周围打起来必须先去噪再补。这就是为什么我一再强调动手之前先看边界质量别拿到洞就急着点填充按钮。提示补洞之前先确认模型是不是水密watertight。很多下游流程比如3D打印切片、有限元网格划分对水密性有硬性要求一个洞没补干净可能整条流水线都卡住。2. 洞的识别与分类动手前必须先把边界摸清楚补洞算法再强也得先知道洞在哪、长什么样。这一步的输入是一个半残缺的网格输出是若干条闭合的边界环。听起来简单但里面有不少门道。2.1 从半边结构里把边界边挑出来主流做法是基于半边half-edge数据结构来查边界。半边结构把每条无向边拆成两条方向相反的有向半边每条半边的对偶半边的对应关系记录在案。对于一个正常闭合的流形网格每条边的两条半边都能找到配对而边界边只有一条相邻三角形它的对偶半边是空的。换句话说只要遍历所有半边把那些没有配对的挑出来剩下的就是边界。挑出边界边以后还要把它们连成一条条闭合环。做法是从任意一条边界边出发沿着这条半边的起点等于下一条半边的终点这个方向一路走直到回到起点就得到一个洞的边界环。这里要小心一个常见坑如果模型拓扑不干净比如存在非流形边一条边被三个以上三角面共享边界追踪可能走进死胡同或者把两个本该分开的洞连成一个。我的经验是在追踪之前先跑一遍非流形检测把那些三面共边的位置标出来单独处理。注意边界环追踪失败十有八九不是代码 bug而是模型本身拓扑有问题。先怀疑数据再怀疑算法。2.2 洞的几种形态处理策略完全不同把洞都当成一个圈来处理是最容易踩的坑。实际上洞有很多形态我习惯分成四类洞的类型特征常用处理策略简单洞单条闭合边界环内部无其他元素直接三角化或曲面拟合带岛洞边界环内部还有独立的小网格岛先处理岛再补外环嵌套洞洞里有洞形成环形区域逐层追踪按层级补槽形缺口边界不闭合像个开口的 U 形需要先延伸投影闭合边界简单洞最好对付边界环一抓直接上算法就行。带岛洞最容易被忽略你以为抓住了一个洞其实里面还飘着一小片孤立的网格补的时候要么把它挤掉要么得绕开它。嵌套洞和槽形缺口更麻烦尤其是槽形缺口它的边界根本不是闭合的补洞算法默认你要补一个圈内部遇到开口的槽它会直接报错或者补出奇怪的东西。我的处理习惯是先把槽形缺口的边界向两个方向虚拟延伸人为闭合再交给补洞算法最后把延伸出来的部分裁掉。分类这一步花不了多少时间但它能避免后面一大半的返工。我见过太多人上来就跑补洞结果带岛洞补得一团糟回头还要删掉重来。3. 补洞算法的三条技术路线与各自的适用边界补洞算法发展到现在能用的方法很多但归归类基本就三条路线直接三角剖分、隐式曲面拟合、模板匹配。它们的复杂度和适用场景差别很大选错路线会让后面的fairing变得异常艰难。3.1 直接三角剖分快但它只解决有没有面最朴素的补洞就是把边界环上的顶点直接连成三角面。经典方法有耳切法ear clipping和基于 Delaunay 准则的三角剖分。耳切法思路简单每次从边界多边形上切掉一个耳朵三角形直到多边形被剖完。快是快但缺点是容易切出又细又长的狭长三角形这种三角形在视觉上就是一排尖刺而且会让后续fairing非常难受。所以工程上更常用的是带权重的最优三角剖分比如 Liepa 那套加权思路给每个候选三角形算一个代价代价里同时考虑三角形面积、和邻面的二面角、以及边长然后选总代价最小的剖分方案。它的好处是补出来的三角面分布比较均匀不容易出现狭长三角形而且剖分完还能顺手优化一下内部新增顶点的位置让补丁稍微贴合周围曲率。这条路线最大的优点是快几十个顶点的小洞基本瞬间出结果而且不依赖洞周围的曲面信息。但它的短板也很明显它只保证把面填上不保证补出来的形状是那个物体本该有的形状。如果洞很大直接剖分出来的补丁通常是个平面或者轻微弯曲的片和周围曲面接不自然这时候就必须靠fairing硬揉效果往往一般。我的经验是边界顶点数在几十个以内、洞周围曲率变化不大的时候直接剖分是首选洞一大就别硬上。3.2 隐式曲面拟合大洞和复杂曲面的主力当洞比较大或者周围曲面有明显曲率时我会转向隐式曲面方法代表就是径向基函数RBF和泊松重建。思路是把洞周围一圈点的空间分布拟合成一个隐式函数这个函数的零等值面就是我们要找的补丁曲面然后在洞的范围内对这个等值面做三角化。RBF 的做法是用一组基函数的加权和去插值已知点求解一个线性系统得到权重泊松方法则更进一层它把问题转成求解泊松方程从离散点重建出一个全局的指示函数再提取等值面。这两类方法的威力在于它们不是就洞补洞而是利用了洞周围的整体形状信息补出来的补丁能自然延续周围的曲率趋势。一个弯曲的圆柱壳破了个洞用泊松方法补出来还是一段弯曲的柱面而直接剖分只会给你一块平板。代价是计算量RBF 求解线性系统的规模随点数量增长很快几百个点还行上万点就得考虑稀疏化或者分块求解泊松重建虽然相对高效但也需要控制到合适的八叉树深度不然内存和耗时都会失控。用这类方法有个关键前提洞周围得有一圈质量足够好的点。如果洞边缘本身噪声就大拟合出来的隐式函数会把这些噪声当成真实形状补出来的补丁跟着一起抖。所以我在用隐式方法之前几乎一定先对洞周围一圈点做一次局部去噪。3.3 模板与相似性补洞有重复结构时的捷径还有一类方法很讨巧如果模型上存在大量重复结构比如机械零件上成排的孔、建筑模型上重复的窗格那就可以在模型的其他位置找一个形状相似的完整区域把它对齐、变形、贴到洞上。这就是基于模板或相似性的补洞。它的优点是补出来的形状特别对因为它本来就是从同一个物体上抄来的连细节纹理都能对上。但它的适用面很窄只有当模型确实存在重复结构时才有效而且怎么找相似块、怎么对齐、怎么变形适配边界每一步都要额外投入。我一般把它当成一个备选方案在直接剖分和隐式拟合都补不好的特定场景下才考虑。4. fairing算法光顺的数学直觉以及它偷偷拿走的东西补丁放上去只是第一步接下来轮到fairing算法登场。很多人对fairing的理解停留在把表面磨平但光顺从来不是免费的午餐它在让表面变光滑的同时一定会在别的地方拿走一些东西——最常见的代价就是体积收缩和特征丢失。理解这一点你才能知道参数该往哪个方向调。4.1 拉普拉斯光顺为什么会把模型饿瘦最经典的fairing就是拉普拉斯光顺。它的操作极其简单对每个顶点把它朝它所有邻居顶点的平均位置挪一步。用公式说就是 v_new v λ × (邻居平均位置 − v)λ 是步长一般取 0 到 1 之间。它的直觉是让每个点都往周围点的中心靠一靠反复迭代表面上的锯齿和小褶皱就被抹平了。问题是这个操作本质上等价于在曲面上做热扩散。热传导会把温度差抹平同时也会让整个曲面的面积和体积朝内收缩。你迭代个几十次一个球体会肉眼可见地变小一圈这就是著名的拉普拉斯收缩问题。对于补洞场景收缩的后果就是补丁比原本该有的形状小了一圈和周围曲面之间出现一个凹进去的台阶。我刚开始做的时候踩过这个坑光顺完看着挺光滑一量尺寸发现关键部位的直径少了百分之几对于精密零件来说这是致命的。提示拉普拉斯光顺的迭代次数一定要和模型的尺度挂钩。同样的参数小模型上迭代五十次可能刚好大模型上可能还没光顺到位或者已经缩没了。4.2 Taubin 与曲率流在光顺和保体积之间找平衡为了对付收缩Taubin 提出了一个很聪明的做法不要一直用正的步长往前推而是正负步长交替。他把拉普拉斯光顺拆成两步——先做一次正步长的扩散 λ再做一次负步长的扩散 μ其中 μ 是负数且绝对值略大于 λ经典搭配是 λ 0.5、μ −0.53。正步长负责把高频噪声压下去负步长负责把低频的收缩反弹回来一正一负抵消掉整体体积的变化。这就是 Taubin 光顺也叫 λ|μ 光顺在实际工程里比纯拉普拉斯稳得多。更进阶的一条路线是基于曲率流curvature flow的光顺。它让顶点沿着曲率方向移动目标是降低曲面能量本质上和平均曲率流是一个东西。它的好处是光顺效果里有明确的几何意义能控制曲面面积和曲率的变化但直接用它同样会收缩所以通常要配合体积约束或者额外的罚项。工程上我的取舍是小块补丁用 Taubin 就够了性价比最高大块补丁、曲率变化剧烈的地方才考虑曲率流加约束。4.3 保特征的双边光顺别把该有的棱角也磨没了有一类模型特别怕光顺就是那种有明确棱角、有锐利特征的结构比如机械零件的倒角、建筑的檐口、雕塑的衣褶。单纯的拉普拉斯或 Taubin 光顺是最一视同仁的它不管这里是平面还是棱边一律抹平结果就是棱角被磨圆细节全丢了。这时候要用保特征的光顺代表就是基于双边滤波的思路在决定一个顶点往哪挪的时候不光看邻居的空间距离还看法向量的差异。法向差别大的邻居说明这里可能是棱边权重降低法向接近的邻居权重拉高。这样曲面上平坦的区域被充分光顺而棱边处因为两侧法向差异大几乎不动特征就被保住了。我在做文物数字化扫描件的时候特别依赖这一类方法。那些刻字、浮雕的边缘一旦被磨圆文物的形貌信息就损失了而保特征光顺能在去噪的同时把线脚和字口留住。这里有个参数要重点调法向差异的阈值。阈值设小了本该保的棱边被当成噪声磨掉设大了噪声也会被当成特征留下来。这个阈值没有万能值得对着模型反复试我的习惯是从着色的法向图上看棱边分布再定阈值。5. 补洞和 fairing 的衔接连续性才是真正的难点把补洞和 fairing 各做一遍并不等于问题解决了。真正决定成品质量的是两者之间的衔接处也就是补丁和原曲面的接缝。这条缝处理不好前面做得再漂亮都白搭。5.1 C0、C1、C2 连续性到底要保到哪一级连续性这个概念通俗说就是接缝处顺不顺。C0 连续性只要求两边位置对上也就是补丁边界顶点和原边界顶点重合接缝处没有裂缝。C1 连续性要求切平面连续也就是接缝两侧的法向一致没有折角——视觉上就是摸上去是顺的。C2 连续性要求曲率也连续接缝处连曲率的变化都不跳变这是最高级别通常只有要求极高的场合才会追。实际操作中直接三角剖分补出来的补丁一般只能做到 C0补丁和原曲面之间会有一道明显的折角。隐式曲面拟合如果处理得好能比较自然地到 C1因为它利用了周围的曲率信息。但要注意光顺这一步很可能会破坏已经达到的连续性——如果你只对补丁区域做 fairing补丁内部变平了补丁和原曲面的交界处反而会形成新的台阶。我的建议是fairing 的作用范围一定要跨过接缝把补丁边界往外扩展一圈到两圈顶点一起参与光顺让过渡是渐变的而不是在接缝处硬切一刀。5.2 网格密度不一致带来的过渡难题还有一个特别容易被忽视的问题补丁的网格密度往往和原曲面的密度对不上。隐式曲面拟合出来的补丁要么特别密、要么特别疏而原曲面是扫描出来的密度又是另一种分布。当你把密度差异很大的两块网格缝在一起接缝处的三角形会拉得很长、很畸形光顺的时候这些畸形三角形会被进一步放大出现抖动或者尖刺。我的处理习惯是在补洞之前先做一次局部重网格化把洞周围一圈的网格密度统一到差不多水平再让补洞算法去补。另外一个办法是在补洞之后对补丁区域做一次自适应简化把过密的三角形减下去让密度和周围接近。这两步都会多花点时间但它们能让后面的 fairing 少踩很多坑。说白了光滑的接缝是缝出来的不是磨出来的。6. 一套能直接抄的实操流程和参数经验讲完原理我把平时处理一个带洞模型的完整流程整理出来。这套流程不是唯一解但覆盖了大多数常见场景你可以按自己的模型情况删减。6.1 预处理洞的筛选和局部清洗第一步永远是清洗别急着补。先跑一遍非流形检测把三面共边、重复顶点、退化三角面标出来能自动修的自动修修不了的单独记下来。然后做洞的追踪和分类按前面的表格判断每个洞属于哪一类。这一步的输出是一张清单有几个洞、每个洞有多少个边界顶点、周围有没有噪声、有没有岛。对于边界点数少于十个的碎洞我的习惯是直接用一个简单的封盖处理掉不值得为它上复杂算法大洞才进入后面的精细流程。6.2 补洞参数怎么定从边界点数和曲率来推断补洞阶段最关键的是路线选择和参数。我一般这么判断边界顶点数小于五十、周围曲率平缓的走加权三角剖分边界顶点数超过五十或者洞横跨了明显的曲率变化比如洞里穿过一个圆弧面走隐式拟合。隐式拟合里泊松方法的八叉树深度我通常从八到十开始试深度太浅补丁丢掉细节太深则计算量爆炸还可能过拟合噪声。RBF 则要注意基函数的选择和正则化项的强度正则化太弱会被噪声带偏太强补丁会变得过于平坦。这些参数没有一个放之四海皆准的值我的做法是先在小洞上试出合适的档位再套到大洞上微调。6.3 fairing 的迭代次数和步长宁少勿多光顺阶段我的核心原则是宁少勿多。因为光顺是不可逆的磨过头就回不来了而磨得不够可以再来一轮。具体到参数拉普拉斯和 Taubin 的迭代次数我一般从十次开始观察接缝处的着色法向图不够再往上加。步长方面拉普拉斯 λ 取 0.3 到 0.5 比较稳太大会震荡Taubin 的 λ|μ 用 0.5 和 −0.53 的经典搭配通常不用改。保特征光顺的法向阈值先定一个偏大的值多保留特征再逐步调小边调边看棱边有没有被磨掉。这里再强调一次作用范围我习惯把 fairing 的作用区域设成补丁区域 往外扩两圈顶点这样过渡自然又不会把整个模型都动一遍。对整个模型做全局光顺是新手最常犯的错辛辛苦苦扫出来的细节一轮全局光顺全没了。7. 我踩过的那些坑以及怎么排查最后这部分是我觉得最有价值的地方因为这些问题在教科书里基本不会写但实操中一定会撞上。7.1 补完也光顺了模型却整体变胖或者变瘦这是最隐蔽的一类问题。表面上看接缝也顺了但一对比原始尺寸发现关键部位的直径或者厚度不对。常见原因是fairing 的作用范围或者迭代次数没控制好导致补丁区域整体内缩拉普拉斯收缩或者补丁本身拟合时就偏大隐式曲面外扩。排查思路是先把 fairing 关掉只看补洞结果量一下尺寸再逐步打开 fairing看每一步尺寸怎么变。这样你能立刻定位到底是补洞阶段的偏差还是光顺阶段的收缩。定位清楚以后要么换 Taubin 这类保体积的方法要么给补丁加上边界约束强制它贴合原始边界。7.2 补丁边缘出现一圈抖动或者尖刺这个现象几乎都出在接缝的网格质量上。要么是补丁和原曲面的网格密度差太多接缝处三角形被拉畸要么是洞的边界本身不干净有几个噪声顶点没清掉补洞算法把这些噪声当成真实形状来拟合。排查方法很直接单独放大接缝区域看网格。如果是一片畸形三角形就去做局部重网格化如果是几个孤立的尖刺多半是边界顶点的问题回去清边界。我吃过一次亏一个黑色哑光材质的扫描件洞边缘全是噪声点我偷懒没清直接补补完边缘一圈刺返工重头做。7.3 大洞补出来的形状完全不对大洞补错八成是路线选错了。用一个简单的加权三角剖分去补一个横跨复杂曲率的大洞结果必然是一块平板糊上去怎么 fairing 都救不回来。这时候别在光顺上死磕回到补洞阶段换隐式拟合。还有一种情况是洞周围的数据本身就不够——周围点太稀隐式函数拟合不出正确的形状补出来的东西凭空想象。这种时候要考虑是不是先在洞周围做一次上采样或者干脆接受这个洞补不完美用人工干预的方式手工拉几根引导线。排查这类问题的通用思路是逐层剥离先关掉 fairing 看原始补丁再关掉补洞看原始边界一层层往回找问题总能定位到具体的某一步。别一上来就盯着最后的结果干瞪眼问题从来不在最后一步只是最后一步把它暴露出来了。我个人在实际操作里的体会是补洞和 fairing 这两件事慢就是快。花十分钟认真做边界清洗和洞的分类能省掉后面一两个小时的返工。算法本身其实都成熟了真正拉开差距的是对模型的理解和对参数的拿捏。多做几个不同类型的模型你自然会形成一套自己的判断顺序——什么时候该快、什么时候该稳、什么时候该老老实实手工修几根线。这套直觉比任何一个具体算法都值钱。