跳转到内容
Jason's Study

Loop 1987 — 三角形网格的递归光滑细分

是什么

Loop 细分（Loop scheme）是一套反复把三角形网格切细，最终收敛成一张光滑曲面的方法。日常类比：木雕粗坯打磨——你拿到一块棱角分明的木块，每打一遍砂纸，棱角就钝一点，打三五遍之后摸上去就像鹅卵石。

Loop 不像 Catmull-Clark 那样处理四边形网格，它只吃三角形网格——这一点让它天然契合：游戏引擎里所有渲染最终都被切成三角形、3D 扫描出来的网格也是三角形、Mesh shader / GPU tessellation 硬件原生处理的也是三角形。

数学角度看，Loop 把规则三角网格上的 4 次盒样条（quartic box spline） 推广到了任意拓扑。一次细分把每个三角形切成 4 个小三角形（不是切两半，是切成 4 个），然后按一套固定权重重算所有顶点位置。

为什么重要

不知道 Loop 这件事，下面这些解释不通：

  • 为什么 Catmull-Clark 是四边形派、Loop 是三角形派——同一个时代两个世界两套算法
  • 为什么游戏引擎做 LOD（远处低模、近处高模）能”实时”提升精度——Loop 是其中一种主流方案
  • 为什么 OpenSubdiv（Pixar 2012 开源）里两套并存——CC 给电影艺术家、Loop 给游戏管线
  • 为什么 DirectX 11 引入硬件 tessellation 时直接选了类似 Loop 的三角分裂模式

核心要点

Loop 把一次细分拆成 两类新点：

1. 边上新点（odd vertex）

每条边中点会生成一个新顶点。它不是简单的两端点平均，而是按 3/8 + 3/8 + 1/8 + 1/8 的权重：

        v3
        /\
       /  \
      /    \
   v1 ------ v2     新点 e = 3/8*v1 + 3/8*v2 + 1/8*v3 + 1/8*v4
      \    /
       \  /
        \/
        v4

也就是边的两个端点各占 3/8，边两侧的两个对顶点各占 1/8。直觉：新点离边重，离对顶点轻，但完全不忽略对顶点——这正是它能光滑的关键。

2. 老顶点的新位置（even vertex）

每个原顶点也要被”重算”成一个新位置（不会原地不动）。设这个老顶点的邻居数为 n（图论里叫度数 / valence），新位置的公式：

v_new = (1 - n*β) * v_old + β * (邻居之和)

其中 β 是按 n 算出来的权重：

β = (1/n) * [ 5/8 − (3/8 + (1/4)*cos(2π/n))² ]

公式看着吓人，但记一个常用情形就够：正则点 n = 6 时，β = 1/16。也就是新位置 = 5/8 老位置 + 1/16 × 6 个邻居（每个邻居权重 1/16）。

为什么 n = 6 是”正则”？因为平面上等边三角形铺开时，每个顶点恰好接 6 个三角形——这是三角网格的”理想态”。其他度数（3、4、5、7…）叫奇异点（extraordinary vertex），它们的处理只能 C¹ 连续。

3. 切完之后

一次细分后，每个老三角形 ABC 变成 4 个小三角形：中间一个由三个新边点组成，剩下三个分别在 A、B、C 三个角上。三角形数 ×4，顶点数大约 ×4。

为什么会光滑（一句话直觉）

Loop 在他的论文里证明：不断重复这两条规则，曲面在正则点处收敛到盒样条曲面（一种已知光滑的连续曲面），所以是 C² 连续（二阶导数处处连续）。

奇异点处证明用了特征值分析——把”细分”看成一个矩阵作用，要求矩阵的最大特征值是 1，第二第三大相等且小于 1，这样反复作用之后曲面会收敛到一个切平面。Loop 的两条权重就是为了凑出这个特征值结构反推出来的。

实践案例

案例 1：Blender 里你能直接看到

Blender 选中一个三角网格 → Modifier → Subdivision Surface → 把”Type”切到 “Simple”（线性细分，不光滑）和 “Catmull-Clark”（默认，四边形）你都能选，但实际三角网格被分时，引擎里跑的就是 Loop 的近亲。Pixar OpenSubdiv 库里 Far::TopologyRefiner 可以直接选 Loop scheme。

案例 2：游戏 LOD 怎么用

近距离：跑 2 级 Loop 细分（×16 三角形数），表面光滑、看不到棱 中距离：跑 1 级（×4） 远距离：用原始低模

这套切换由 GPU 实时做，每帧根据距离决定细分级别。DirectX 11 / OpenGL 4 的 tessellation shader 硬件原生支持。

案例 3：3D 扫描后处理

LiDAR / 摄影测量出来的点云三角化后，表面有阶梯感。跑 1-2 次 Loop 细分能把阶梯变成平滑斜面，但奇异点处会留下小平面（C¹ 不够 C²），这是在反光强的金属件上偶尔能看到的瑕疵。

踩过的坑

  1. Loop 只吃三角网格：拿到混合 n 边形（含四边形、五边形），必须先三角化再细分。三角化方式不同 → 细分结果不同（不是天然唯一的）。
  2. 奇异点只 C¹：n ≠ 6 的顶点处只能保证一阶光滑，二阶导数会跳。在反射强的渲染下能看到一个小亮斑。游戏里通常无所谓，影视渲染要在意。
  3. 没有锐边/折痕：原始 Loop 假设网格处处光滑。要做”硬边”（比如方块的角不能被磨圆），得用 Hoppe 1994 的扩展（“piecewise smooth”）—— 给边打标签，标了的边在细分时按特殊规则处理。
  4. 内存爆炸：每级 ×4 三角形。3 级 = 64×。要做 LOD 必须懒求值，不能预先全展开。
  5. β 的公式有多个版本：原版 Loop 1987 的 β 和 Warren 1995 的简化版（β = 3/(8n) for n > 3, β = 3/16 for n = 3）数值上接近但不一样。OpenSubdiv 用 Warren 版。

适用 vs 不适用

适用：

  • 三角网格的光滑化（游戏 / 实时渲染 / GPU tessellation）
  • 角色建模、有机形状（人脸、动物、布料）
  • 3D 扫描点云后处理

不适用：

  • 四边形主导的网格 → 用 Catmull-Clark（catmull-clark-1978）
  • 需要精确锐边的硬表面（CAD / 工业件）→ 用 NURBS 或带 crease 的扩展
  • 拓扑会变（比如撕裂、合并） → 细分假设拓扑稳定
  • 极高精度科学可视化（要 C∞）→ 用解析曲面

历史小故事（可跳过）

  • 1974 年：Catmull 在博士论文里提出”递归把曲面切细”的想法，但只对规则四边形可用
  • 1978 年：Catmull-Clark 把它推广到任意拓扑四边形网格 → catmull-clark-1978
  • 1987 年：Charles Loop 在犹他大学（Tony DeRose 指导）的硕士论文 44 页，把三角形版本搞定
  • 1994 年：Hoppe 等扩展支持锐边、角点（“Piecewise Smooth Surface Reconstruction”, SIGGRAPH 94）
  • 1996-2000：Reif、Schweitzer 等用特征值分析正式证明 Loop 在奇异点的连续性
  • 2012 年：Pixar 开源 OpenSubdiv，CC + Loop 双方案同时进入工业级 GPU 实现

四十年后，三角网格细分仍然是 Loop 的形状。

学到什么

  1. 细分曲面 = 极限过程：定义不是”曲面方程”，而是”反复施加一条规则，看它收敛到哪”。这是一种全新的几何定义方式
  2. 三角 vs 四边形是两个并行宇宙：CC 和 Loop 同时代发展，互不替代——选哪个看你的网格类型
  3. 特征值反推权重：先想”我希望曲面光滑（一阶/二阶连续）“，再用特征值条件倒推出权重公式。这是从需求到算法的经典推理路径
  4. 正则 vs 奇异：连续性在”理想态”和”边缘态”是两种不同的保证。工程上要承认边缘态的妥协

延伸阅读

关联

反向链接