精读笔记
Problem Setting
这篇论文真正处理的是 Laguerre tessellation 的 inverse barycenter problem：给定连续参考密度 ρ、每个 cell 的质量 v_i 和 barycenter b_i，恢复一个具有这些 cell moments 的 Laguerre tessellation。困难点在于 barycenter 约束不是 generators/weights 的简单约束；weights 对 fixed generators 可由 semi-discrete OT 调整出 volumes，但 barycenters 反过来约束 generators 时是非凸且退化的。已有路线如直接优化 generators/weights 或 stochastic fitting 可以工作，但没有说明哪些 B 可由 Laguerre tessellation 实现，也无法区分真实不可行、退化可行、数值优化失败。关键矛盾是：tessellation 由 volumes+barycenters 唯一确定，但其 Laguerre 参数不唯一，甚至在退化/层级极限下普通 Laguerre 参数会坍缩。
Motivation
已有方法缺的是 barycenter 数据的“可行域几何”。如果只在 generator space 搜索，问题看起来像一个 ill-conditioned inverse optimization；但从 OT 角度看，Laguerre cells 是把 ρ transport 到离散 measure ν_N(v,Y) 的 Monge cells，而 cell barycenters 是 ρ 在这些 cells 上的一阶矩。作者的关键观察是：所有可能的 barycenter vectors 不应先由 generators 参数化，而应由 convex order 直接定义。这样可以绕开 generator non-uniqueness，把反演问题降到一个凸体 C_N(v,ρ) 的边界结构上。真正缺口是 semi-discrete convex order 的 exposed/extreme point 结构此前没有被显式用于 Laguerre fitting。
Core Idea
核心思想是把“找 Laguerre tessellation”改写成“在 convex-order convex body 上找 exposed point”。定义 C_N(v,ρ) 为所有满足 ν_N(v,B) ⪯_C ρ 的离散 barycenter positions。这个集合包含所有能由任意 ρ 分解产生的 barycenters；而论文证明其中最几何化、最刚性的点——exposed points——恰好对应 Laguerre tessellation 的 cell barycenters。
这改变了建模方式：prior 把 generators 作为 primary variables，这篇把 barycenters 的 convex-order feasibility 作为 primary object，generators 只是 exposing normal。新的 inductive bias 是凸序/一阶矩守恒结构，而不是直接拟合 cell boundaries。直觉上它有效，是因为 Laguerre tessellations 正是 semi-discrete OT support function 的 differentiability points；当 Y 不碰撞时，F_N 的梯度唯一，给出的就是 cells 的 barycenters，因此 exposed geometry 和 OT cells 对齐。
Method
1. Convex-order feasible set C_N：解决“哪些 barycenters 至少可能由某种 ρ 分解产生”的问题。它把非凸几何反演变成凸体分析；核心变化是从 cell boundary 参数转向 measure decomposition。
2. Support function F_N：解决“如何从 generators 读出 C_N 的边界”的问题。F_N 是 semi-discrete OT objective，也是 C_N 的 support function；在非碰撞 generators 处可微，梯度就是 Laguerre cell barycenters。这是把 OT 计算和凸几何连接起来的桥。
3. Exposed/extreme point characterization：解决“Laguerre barycenters 在 C_N 中处于什么位置”的问题。普通 Laguerre barycenters = exposed points；hierarchical Laguerre barycenters = extreme points。这解释了退化 generator collapse 为什么对应层级 tessellation，而不是简单失败。
4. Projection reconstruction：解决“如何数值恢复 generators”的问题。对 λB 投影到 C_N，projection residual Y*=λB−P_C_N(λB) 是 supporting normal；若 Y* 无碰撞，它就是近似 generators。这个机制用强凸 projection 替代不稳定的 Legendre maximization。
5. Subgradient / Frank-Wolfe：解决 projection 的可计算性。它们依赖反复求 semi-discrete OT/Laguerre cells；从贡献归因看是必要工程层，而非最核心数学创新。
Key Insight / Why It Works
最核心的贡献是识别出 barycenters of Laguerre cells are exposed points of a convex-order body。这个 insight 很强，因为它同时解释了可恢复性、不可恢复/退化、hierarchical limits 和 projection-based fitting。方法有效不是因为更强的 optimizer，而是因为变量空间选对了：C_N 把所有一阶矩可行性压进一个凸集合，而 Laguerre structure 出现在其 support geometry 中。
本质上这是 better inductive bias + latent geometric structure，不是 scaling，也不是 retrieval。它把原问题的非凸性转移为“在凸体边界上找 exposed normal”；虽然 projection 数值上仍需 semi-discrete OT，但理论对象变清楚了。最可能的核心贡献是 Theorem 3.6/3.9 这一组 exposed/extreme characterization；projection algorithm 是利用该结构的自然后果。
辅助部分包括 subgradient 和 Frank-Wolfe 的收敛界，比较标准，增益来源不在算法新颖性。实验中 Frank-Wolfe 近似线性收敛文中未充分说明，可能来自具体几何实例较好而非一般理论。EBSD fitting 的成功也更多说明 convex-order projection 是合理 inductive bias，不足以证明真实材料重构中边界几何被准确恢复。
Relation To Prior Work
最接近的三条线是 semi-discrete optimal transport、convex order/Wasserstein projection、以及 balanced clustering/gravity bodies。和 Bourne–Pearce–Roper 2025 的直接反演相比，本质差异是这篇不把 Legendre transform maximization 当作唯一入口，而是先刻画 barycenter feasible set 的凸几何，再从 projection normal 恢复 generators。和 Ekeland–Schachermayer 的 convex order/OT 几何相比，新意在 semi-discrete setting 下把 exposed/extreme points 显式落到 Laguerre/hierarchical Laguerre tessellations。和 fully-discrete gravity body 相比，这篇是把 polytope vertex picture 推到 absolutely continuous ρ，此时 C_N 通常不再是多面体，边界结构更微妙。
看似新的 numerical fitting 其实很大程度是已有 Wasserstein projection、semi-discrete OT 和 Frank-Wolfe/subgradient 的重组；实质创新在于证明这些对象正好对应 Laguerre inverse problem 的几何本质。它属于“OT-based computational geometry / convex-order optimization”谱系，而不是传统材料图像重构算法。
Dataset / Evaluation
evaluation 主要分两类：合成 Laguerre tessellation reconstruction 和真实 EBSD steel image fitting。合成实验能验证 projection normal 可以恢复接近的 barycenters，并展示退化情形下 generators 不必接近真实 generators；这反而支持论文的理论判断，即数据唯一确定 tessellation 但不唯一确定参数。真实 EBSD 例子说明方法能处理非理想 barycenters，也即 fitting 而非 exact reconstruction。
但 evaluation 覆盖仍偏窄：主要是二维、有限 N、ρ 为均匀区域密度；没有系统噪声分析、三维 3DXRD 重构、真实 incomplete measurements 或大规模复杂 microstructure。benchmark 支持“方法可用于 fitting/reconstruction”这个 claim，但不足以支持强 scalability 或工程上优于 stochastic optimization 的泛化结论。实验数字不是核心证据，理论结构才是。
Limitation
最大限制是恢复保证并不覆盖所有实际会遇到的 Laguerre data。Proposition 5.2 需要 B 来自 irreducible Laguerre tessellation，并且 λ>1 足够小；但实际退化、近退化、层级结构很常见，此时 projection residual 可能落在 Δ_N，generators 碰撞，或者恢复出与真实生成参数完全不同的 tessellation。文中承认一般 B∈B_Lag 或 B∈B_hLag 时无法保证 Y*∉Δ_N。
第二个限制是误差控制不足。虽然 B*=P_C_N(λB)→B，但 cell convergence 与 ∥B−B*∥_v 之间没有清晰定量关系；这意味着 barycenter 拟合好不必然意味着 cell geometry 好。对材料科学应用而言，这很关键，因为 grain boundary geometry 往往才是下游关心对象。
第三，算法把困难转移到 repeated semi-discrete OT/Laguerre computation。二维例子可行不代表三维大 N、复杂 density、noisy input 下稳定。scalability 上限文中未充分说明。所谓 fitting 的 generality 也应谨慎理解：它能投影任意 B，但 projection 得到的 Laguerre tessellation 是否物理合理，取决于 convex-order inductive bias 是否匹配数据生成过程。
Takeaway
- 1. 最值得记住的是：Laguerre inverse problem 的自然对象不是 generators，而是 convex-order body C_N(v,ρ)；generators 是 exposed normals。
- 2. hierarchical Laguerre tessellations 不是边缘技术细节，而是普通 Laguerre tessellations 的闭包/退化极限，对理解 non-uniqueness 和 generator collision 很关键。
- 3. Wasserstein projection 在这里不是普通 regularization，而是把不可行或 noisy barycenters 拉回 convex-order feasible geometry，并从 projection residual 提取 tessellation direction。
- 4. 未来真正值得做的是定量稳定性：从 barycenter error 到 cell boundary error、从 λ 到 generator separation、从 noisy 3D measurements 到 physically plausible tessellations 的误差链条。
一句话总结
这篇论文把 Laguerre tessellation fitting 从一个病态的 generator 反演问题，提升为 semi-discrete convex-order convex body 的边界几何问题，核心贡献是证明 Laguerre/hierarchical Laguerre barycenters 分别对应 exposed/extreme points，并据此给出 projection-based reconstruction。