[論文レビュー] The Higher-Dimensional Skeletonization Problem

Sara Kališnik, Vitaliy Kurlin+1 more|arXiv (Cornell University)|Jan 29, 2017

Topological and Geometric Data Analysis参考文献 20被引用数 6

ひとこと要約

この論文は単体複体における高次元のスケルトン化問題を導入し分析し、k-スパンニング部分複体の基礎的性質とそのホモロジー的挙動を確立する。k-スパンニング部分複体が次元k−2までホモロジーを保存し、次元k−1で全射となることを証明し、面の追加によるベッチ数の変化を用いた最小スパンニングツリーおよびフォレストの構成基準を提示する。

ABSTRACT

Real data is often given as a point cloud, i.e. a finite set of points with pairwise distances between them. An important problem is to detect the topological shape of data --- for example, to approximate a point cloud by a low-dimensional non-linear subspace such as an embedded graph or a simplicial complex. Classical clustering methods and principal component analysis work well when given data points split into well-separated clusters or lie near linear subspaces of a Euclidean space. Methods from topological data analysis in general metric spaces detect more complicated patterns such as holes and voids that persist for a long time in a 1-parameter family of shapes associated to a cloud. These features can be visualized in the form of a 1-dimensional homologically persistent skeleton, which optimally extends a minimal spanning tree of a point cloud to a graph with cycles. We generalize this skeleton to higher dimensions and prove its optimality among all complexes that preserve topological features of data at any scale.

研究の動機と目的

• 高次元単体複体におけるk-スパンニング部分複体の概念を形式化し分析すること。
• 特にホモロジー同型と全射に関して、そのような部分複体の基礎的ホモロジー的性質を確立すること。
• 面の追加に伴うベッチ数の変化を用いて、最小スパンニング構造（k-ツリーおよびk-フォレスト）を特徴づけること。
• 抽象理論を、特にČechおよびVietoris-Ripsフィルトレーションを通じたトポロジカル・データ解析への実用的応用に結びつけること。
• 標準単体および点群複体における最大k-フォレストおよびスパンニングk-ツリーを特定する構成的基準を提供すること。

提案手法

• 次元i ≤ k−1について包含写像$H_i(\text{S} \to \text{Q})$を単体ホモロジー理論を用いて分析し、k−2まで同型であり、k−1で全射であることを示す。
• オイラー特性の公式を用いて、$\mathcal{S}$と$\mathcal{Q}$間の面数とベッチ数のバランス式を導出する：$\#k\text{-faces in } \mathcal{Q} + \beta_{k-1}(\mathcal{Q}^{(k)}) - \beta_k(\mathcal{Q}^{(k)}) = \#k\text{-faces in } \mathcal{S} + \beta_{k-1}(\mathcal{S}) - \beta_k(\mathcal{S})$。
• $\mathcal{Q} \setminus \mathcal{S}$内のk-面の集合に対して最小性の議論を適用し、$\beta_{k-1}(\mathcal{S}) > \beta_{k-1}(\mathcal{Q})$ならば、ある単一のk-面$F$が$\beta_{k-1}$を減少させつつ$\beta_k$を変化させないことを証明する。
• k-サイクルの基底を$\mathcal{S}$から取り除いて$F \subseteq \mathcal{S}$を構成し、$F$が$(k-1)$-フィッティングかつk-スパンニングのままであることを保証する。
• 相対ホモロジーの長完全系列を用いて、包含写像$H_k((\mathcal{S},F) \to (\mathcal{Q},F))$が同型であるための必要十分条件が$H_k(\mathcal{S} \to \mathcal{Q})$が同型であることであることを示す。
• 神経補題を適用して、点群データにおけるオフセット和集合のホモトピー近似としてČech複体の使用を正当化し、Vietoris-Ripsフィルトレーションと比較する。

実験結果

リサーチクエスチョン

• RQ1 単体複体におけるk-部分複体がk-スパンニングであるための必要十分条件は何か？
• RQ2 k-フォレストに1つのk-面を追加したとき、ベッチ数はどのように変化するのか？これはホモロジー保存に何を意味するか？
• RQ3 単体複体$\mathcal{Q}$のk-部分複体$\mathcal{S}$が$H_{k-1}$で全射であり、$i < k-1$について$H_i$で同型であるための条件は何か？
• RQ4 単体複体における最大k-フォレストはどのように特徴づけられるか？また、これはk-ツリーとどのように関係するか？
• RQ5 k-スパンニング部分複体の理論を、ČechまたはVietoris-Ripsフィルトレーションを用いた点群からの形状再構築にどのように応用できるか？

主な発見

• 複体$\mathcal{Q}$のk-スパンニング部分複体$\mathcal{S}$は、すべての$i \leq k-2$について$H_i$で同型であり、$H_{k-1}$で全射となる。これは強いホモロジー的制御を保証する。
• 面数とベッチ数のバランス式が成り立つ：$\#k\text{-faces in } \mathcal{Q} + \beta_{k-1}(\mathcal{Q}^{(k)}) - \beta_k(\mathcal{Q}^{(k)}) = \#k\text{-faces in } \mathcal{S} + \beta_{k-1}(\mathcal{S}) - \beta_k(\mathcal{S})$。
• $\beta_{k-1}(\mathcal{S}) > \beta_{k-1}(\mathcal{Q})$ならば、$\mathcal{Q} \setminus \mathcal{S}$に属するある単一のk-面$F$が、$\beta_{k-1}(\mathcal{S} \cup \{F\}) = \beta_{k-1}(\mathcal{S}) - 1$かつ$\beta_k(\mathcal{S} \cup \{F\}) = \beta_k(\mathcal{S})$を満たす。
• k-スパンニングk-フォレスト$F \subseteq \mathcal{S}$が存在し、$\mathcal{S}$および$\mathcal{Q}$の両方で$(k-1)$-フィッティングかつk-スパンニングである。ただし$\mathcal{S}$が$(k-1)$-フィッティングであることが前提。
• $\Delta_n$におけるk-フォレスト$F$について、$F$がk-ツリーであるための必要十分条件は、任意の$k$-面$E \notin F$に対して$H_k(F \cup \{E\}) \neq 0$、すなわち$\beta_k(F \cup \{E\}) = 1$であることである。
• 標準n-単体$\Delta_n$において、すべてのスパンニングk-ツリーは$k \geq 1$に対してちょうど$\binom{n}{k}$個のk-面を含み、$k = 0$の場合は0個である。