[論文レビュー] Simplicial Flat Norm with Scale

Sharif Ibrahim, Bala Krishnamoorthy+1 more|arXiv (Cornell University)|May 25, 2011

Topological and Geometric Data Analysis参考文献 16被引用数 4

ひとこと要約

本稿では、任意次元の電流のフラットノルムを計算する手法として、有限の単体複体の部分複体としてモデル化される、マルチスケール単体的フラットノルム（MSFN）を導入する。MSFNが整数上ではNP完全であることを証明するが、複体に相対的ねじれが存在しない場合には、線形計画法の緩和により多項式時間で解けることを示す。また、近似過程における質量拡張を制限する単体的変形定理を確立し、高次元データセットの効率的かつスケーラブルなノイズ除去とスケールに敏感な解析を可能にする。

ABSTRACT

We study the multiscale simplicial flat norm (MSFN) problem, which computes flat norm at various scales of sets defined as oriented subcomplexes of finite simplicial complexes in arbitrary dimensions. We show that the multiscale simplicial flat norm is NP-complete when homology is defined over integers. We cast the multiscale simplicial flat norm as an instance of integer linear optimization. Following recent results on related problems, the multiscale simplicial flat norm integer program can be solved in polynomial time by solving its linear programming relaxation, when the simplicial complex satisfies a simple topological condition (absence of relative torsion). Our most significant contribution is the simplicial deformation theorem, which states that one may approximate a general current with a simplicial current while bounding the expansion of its mass. We present explicit bounds on the quality of this approximation, which indicate that the simplicial current gets closer to the original current as we make the simplicial complex finer. The multiscale simplicial flat norm opens up the possibilities of using flat norm to denoise or extract scale information of large data sets in arbitrary dimensions. On the other hand, it allows one to employ the large body of algorithmic results on simplicial complexes to address more general problems related to currents.

研究の動機と目的

• 任意次元の電流のフラットノルムを計算可能な方法で求めること、特に大規模で複雑なデータセットに対応すること。
• フラットノルムにスケール（λ）を拡張し、単体的複体上でマルチスケール解析を可能にする。
• 一般の電流を単体的電流で近似する際の精度に関する理論的保証を、新規の単体的変形定理を用いて確立すること。
• 相対的ねじれのないトポロジー的条件下で、整数線形計画法によるフラットノルムの効率的計算を可能にし、多項式時間で解けることを示すこと。
• フラットノルムに基づくシグネチャーやスケールに敏感な表現を高次元で可能にすることで、幾何解析とデータサイエンス分野に新たな応用を開くこと。

提案手法

• 連続的な幾何的対象の離散的表現を可能にするために、有限の単体的複体内に電流を単体的電流としてモデル化する。
• マルチスケール単体的フラットノルムを、電流の質量から境界の質量を引いたものにスケールされた表面の質量を加算する形で整数線形計画問題（ILP）として定式化する。
• 単体的複体に相対的ねじれが存在しない場合、ILPの緩和が整数解をもたらす（したがって最適なフラットノルムが得られる）ことを証明し、多項式時間での計算が可能であることを示す。
• 一般の電流を単体的電流で近似する際の質量拡張を制限する単体的変形定理を導入し、複体の細分化に伴い誤差が小さくなることを示す。
• 複雑な幾何構造（例：ノイズと特徴を含むピラミッド）の埋め込みに適した高品質な単体的複体を生成するために、制約付きデローニ四面体分割を用いる。
• 標準的な数値ツール（例：MATLAB）を用いて実装し、3次元データ上で検証を実施。ピラミッドにピーク、谷、ノイズを追加した例を、さまざまなスケールパラメータλで評価した。

実験結果

リサーチクエスチョン

• RQ1 任意次元の電流のフラットノルムは、単体的複体を用いて効率的に計算可能であり、幾何的・位相的忠実性を保持できるか？
• RQ2 単体的複体にどのような位相的条件（相対的ねじれの不在）が満たされている場合、線形計画法の緩和によりマルチスケール単体的フラットノルムが多項式時間で解けるか？
• RQ3 一般の電流は単体的電流でどの程度正確に近似可能か？また、メッシュの細分化によってこの近似誤差を境界づけ・最小化できるか？
• RQ4 スケールλを伴うフラットノルムを用いて、高次元データセットからマルチスケールの幾何的特徴（例：ノイズ除去や特徴検出）を抽出できるか？
• RQ5 さまざまなλ値におけるフラットノルム分解の系列から、形状シグネチャーやコンパクトな表現（例：バーコード）を定義できるか？

主な発見

• ホモロジーが整数上に定義される場合、マルチスケール単体的フラットノルムはNP完全であることが示され、一般設定下での計算の困難性が裏付けられる。
• 単体的複体に相対的ねじれが存在しない場合、MSFNの整数線形計画法定式化は、線形計画法の緩和により多項式時間で解ける。
• 単体的変形定理により、電流近似時の質量拡張に対して明示的かつ定量的な境界が得られ、複体の細分化に伴い誤差を任意に小さくできることが示された。
• 標準PC上で1つのλ値あたり数分の計算時間で、ノイズと特徴を含む3次元ピラミッドの計算が完了し、大規模データへの適用可能性が示された。
• 高λ（例：λ=6）では細かい特徴が保持され、低λ（例：λ=0.01）ではフラットなシートが得られ、効果的なスケールに敏感なノイズ除去と特徴抽出が可能であることが確認された。
• 本手法により、2次元画像や余次元1の境界に限られていた従来の研究を越えて、任意次元におけるフラットノルムを用いたノイズ除去とスケール解析が可能になった。