网格生成文献卷 03:四边形网格、方向场与全局参数化

先用六层管线整理文献

四边形网格研究最容易被算法名淹没。把论文放回管线,就能看出它究竟解决哪一层:

flowchart LR
    A[几何/特征] --> B[方向场]
    B --> C[奇异点与拓扑]
    C --> D[可积场/无缝参数化]
    D --> E[整数周期与量化]
    E --> F[网格提取]
    F --> G[布局和质量优化]

方向场只规定局部轴向;参数化需要它是坐标函数的梯度;整数格点图还要切缝跳跃满足90∘90^\circ 旋转与整数平移;最终网格提取必须处理翻转、退化和切缝识别。任一层成功都不自动推出后一层成功。

1991–2005:前沿、调和函数与几何驱动

级别工作贡献今天仍值得读的原因
ABlacker–Stephenson, Paving, 1991从边界向内部铺设四边形并用模板闭合理解工业 quad 的局部规则、前沿冲突与不可预测奇异点
BAll-Quadrilateral Meshing via Paving, 1998 survey context将 paving、mapping、medial decomposition 放入工程分类区分结构化映射法和非结构 paving
AHarmonic Functions for Quadrilateral Remeshing, 2005以调和函数及其 Morse 结构生成拓扑矩形区域标量场临界点如何决定 patch 拓扑

早期工程路线的优点是边界直观、容易和 CAD 交互;缺点是内部拓扑由局部推进事件产生,很难全局控制奇异点。后来的方向场和参数化路线正是在解决这种“全局不可预见性”。

2006–2011:全局参数化成为主线

从周期坐标到整数变量

设切开曲面后的坐标为ϕ=(u,v)\phi=(u,v)。切缝两侧允许

ϕj=Rkπ/2ϕi+t,t∈Z2.\phi_j=R_{k\pi/2}\phi_i+t, \qquad t\in\mathbb Z^2.

旋转编码四向对称,整数平移保证等参数线在切缝处仍落在同一张整数格上。混合整数优化不是附加技巧,而是离散拓扑信息的载体。

级别论文研究转折关键限制
APeriodic Global Parameterization, 2006用周期参数坐标将方向、尺度与无缝拼接统一切图、奇异点与翻转仍影响结果
BDesigning Quadrangulations with Discrete Harmonic Forms, 2006用调和 1-形式控制全局周期与网格线用户约束和同调基选择影响布局
AQuadCover, 2007在分支覆盖上把交叉场提升为向量场,再用 Hodge 分解获得可积结构覆盖空间处理并不自动给出理想整数采样
AMixed-Integer Quadrangulation, 2009将方向对齐、奇异点与整数跳跃纳入 mixed-integer 求解大规模整数优化与局部翻转
BSimple Quad Domains for Field-Aligned Parametrization, 2011用较简单 quad domains 组织参数域和特征分区质量与自动化选择

2012–2017:从“求一个场”到“保证可积并可靠提取”

级别论文它补上的缺口不应夸大的地方
AGlobally Optimal Direction Fields, 2013用谱/线性结构在奇异点配置上求全局平滑解“场的全局最优”不等于网格存在或参数无翻转
AInteger-Grid Maps for Reliable Quad Meshing, 2013明确整数格点图的统一条件并增强可靠性仍需求解离散整数选择和处理退化
CQEx, 2013专门解决从不完美 integer-grid map 提取显式 quad 的工程陷阱它修复提取层,不替代上游场和参数设计
AIntegrable PolyVector Fields, 2015直接优化无旋/可积的多向量场,使其对应参数梯度连续可积仍不自动满足整数周期与目标布局
BInstant Field-Aligned Meshes, 2015联合平滑方向场和位置场,实现交互级大模型处理追求速度与实用性,理论/拓扑控制弱于全局整数法
ADirectional Field Synthesis, Design, and Processing, 2016系统整理NN-方向场表示、离散联络、奇异点与设计读完要能区分表示层、拓扑层和积分层
APartitioning Surfaces Into Quadrilateral Patches, 2017梳理 patch layout、separatrix 与参数化的关系patch 数少、几何对齐与网格质量常互相冲突

2018–2025:鲁棒性、拓扑可实现性与自动化

级别工作贡献研究判断
BQuadriFlow, 2018将方向、位置和奇异结构协调为自动 manifold quad 流程面向图形重拓扑,不等同于严格 CAE 网格认证
BLearning Direction Fields for Quad Mesh Generation, 2021从数据学习方向场候选训练分布、旋转对称表示和下游可积性需单独验证
AWhich Cross Fields Can Be Quadrangulated?, 2022研究 prescribed holonomy 下局部单射无缝参数化的存在与构造说明满足指数约束仍不等于可四边形化
AWhy Cross Fields Are Not Equivalent to Quadrilateral Meshes, 2023从纤维丛/示性类角度解释场与结构化网格的全局差距不是否定交叉场,而是指出缺失的全局约束
BAnisotropy and Cross Fields, 2024重新分析各向异性下交叉场设计与一致性关注度量变化对场和平滑能量的影响
BNeurCross, 2025用神经表示拟合曲面交叉场并处理约束学习得到的场仍需拓扑、可积和提取验证

四种“成功”必须分开报告

层级最低检查仍不能说明
场成功边界对齐、奇异指数总和、平滑能量存在对应参数化
参数成功局部 Jacobian 正、切缝变换一致、周期残差整数等值线能稳定提取
layout 成功分区为拓扑圆盘/矩形、separatrix 无非法交叉细分后单元质量好
mesh 成功流形性、纯 quad、几何误差、scaled Jacobian对目标 PDE 更优

开源实现与复现基线

复现时要保存中间产物:输入三角网格、场、奇异点、切图、参数坐标、整数跳跃和最终 quad。只保存最终 OBJ,无法定位失败发生在哪一层。

真正连续的理论问题

  1. 拓扑必要条件:Poincaré–Hopf 约束奇异指数总和;
  2. holonomy:方向沿闭环平移后如何旋转;
  3. 可积性:目标轴是否来自局部坐标函数的梯度;
  4. 周期量化:闭环积分能否落到整数格;
  5. 局部单射:参数 Jacobian 是否处处保持正定向;
  6. 有限布局:积分轨线是否形成有限 patch,而非螺旋或稠密轨道。

这六问分别对应不同数学工具,不能用一个“场很光滑”的指标代替。

推荐阅读次序

  1. 先读 2013 quad survey 和 2016 direction-field survey;
  2. 用 PGP → QuadCover → MIQ 看三种全局化方式;
  3. 读 Integer-Grid Maps 与 QEx,把“隐式网格”接到显式网格;
  4. 读 Integrable PolyVector Fields,分清光滑与可积;
  5. 读 2022/2023 两篇存在性工作,理解为什么指数条件不够;
  6. 最后比较 Instant Meshes、QuadriFlow、QuadWild 和学习方法的工程取舍。

返回文献地图 · 对应教学正文