[论文解读] Demystifying Latschev's Theorem: Manifold Reconstruction from Noisy Data

Sushovan Majhi|arXiv (Cornell University)|May 26, 2023

Topological and Geometric Data Analysis被引用 3

一句话总结

本文首次提供了Latschev定理的定量版本，建立了针对噪声数据的Vietoris–Rips复形在何种采样条件下可与闭黎曼流形同伦等价的明确条件。通过利用流形的凸性半径和截面曲率上界，作者推导出一个可计算的邻近尺度β，确保从噪声样本中实现流形重建的拓扑保真性，进一步通过基于可达性的条件将结果推广至欧几里得子流形。

ABSTRACT

For a closed Riemannian manifold $\mathcal{M}$ and a metric space $S$ with a small Gromov$\unicode{x2013}$Hausdorff distance to it, Latschev's theorem guarantees the existence of a sufficiently small scale $β>0$ at which the Vietoris$\unicode{x2013}$Rips complex of $S$ is homotopy equivalent to $\mathcal{M}$. Despite being regarded as a stepping stone to the topological reconstruction of Riemannian manifolds from a noisy data, the result is only a qualitative guarantee. Until now, it had been elusive how to quantitatively choose such a proximity scale $β$ in order to provide sampling conditions for $S$ to be homotopy equivalent to $\mathcal{M}$. In this paper, we prove a stronger and pragmatic version of Latschev's theorem, facilitating a simple description of $β$ using the sectional curvatures and convexity radius of $\mathcal{M}$ as the sampling parameters. Our study also delves into the topological recovery of a closed Euclidean submanifold from the Vietoris$\unicode{x2013}$Rips complexes of a Hausdorff close Euclidean subset. As already known for Čech complexes, we show that Vietoris$\unicode{x2013}$Rips complexes also provide topologically faithful reconstruction guarantees for submanifolds.

研究动机与目标

• 为解决Latschev定理中长期存在的理论空白，提供一个定量、可计算的邻近尺度β，以实现拓扑重建。
• 基于闭黎曼流形的内在几何不变量——特别是凸性半径和截面曲率上界——建立采样条件。
• 通过可达性作为采样参数，将重建保证推广至欧几里得子流形。
• 证明在Hausdorff邻近性条件下，Vietoris–Rips复形可忠实恢复子流形的拓扑结构，与已知的Čech复形结果相一致。
• 弥合拓扑数据分析中定性同伦等价与实际流形重建之间的理论鸿沟。

提出的方法

• 利用流形M的凸性半径ρ(M)和截面曲率上界Kmax，推导出邻近尺度β的定量上界。
• 通过基于Rips复形的测地线延拓构造，引入一个单纯映射eg: sd(K) → Rβ(S)。
• 通过构造涉及R(1−2ζ)β(M)、Rβ(S)和R4/3(1−2ζ)β(M)的同伦交换图，建立Rβ(S)与M之间的同伦等价。
• 利用引理A.2证明复合映射φ ∘ eg与包含映射ι是邻接的，从而推出同伦映射。
• 利用Whitehead定理得出诱导同伦群映射是同构，因此为同伦等价。
• 对于欧几里得子流形，利用可达性τ(M)定义β，使得当dH(M, S) < ζβ时，有Rβ(S) ≃ M，其中0 < ζ < 1/14。

实验结果

研究问题

• RQ1 在何种显式、可计算的邻近尺度β采样条件下，可确保噪声样本S的Vietoris–Rips复形与闭黎曼流形M同伦等价？
• RQ2 如何利用M的凸性半径和截面曲率推导出Latschev定理的定量版本？
• RQ3 在Hausdorff距离邻近性条件下，Vietoris–Rips复形能否实现对欧几里得子流形的拓扑保真重建？
• RQ4 子流形的可达性在确定实现拓扑保真恢复的尺度β时起何种作用？
• RQ5 是否可能仅依赖计算更简单的Rips复形，而非更复杂的Čech复形，实现重建中的同伦等价？

主要发现

• 建立了Latschev定理的定量版本，其中选择β ≤ (1 − 2ζ)β，使得当dGH(M, S) < δ且δ充分小时，有Rβ(S) ≃ M。
• 邻近尺度β被显式地以凸性半径ρ(M)和截面曲率表示，满足β ≤ (1 − 2ζ)β < (1 − 2ζ)ρ(M)。
• 对于闭黎曼流形M，若β ≤ 3(1 + 2ζ)(1 − 14ζ) / [8(1 − 2ζ)²] ⋅ τ(M)，则Vietoris–Rips复形Rβ(S)与M同伦等价，其中τ(M)为可达性。
• 同伦等价通过一个单纯映射φ: R(1−2ζ)β(M) → Rβ(S)实现，且对所有p ∈ M满足∥p − φ(p)∥ < ζβ。
• 复合映射ψ ∘ φ与包含映射ι: R(1−2ζ)β(M) → R4/3(1−2ζ)β(M)是邻接的，意味着同伦映射，因此诱导出同构的同伦群。
• 结果证实，在Hausdorff邻近性条件下，Vietoris–Rips复形可实现对Rd中子流形的拓扑保真重建，且可达性τ(M)是关键的采样参数。