The Intersection Euler Characteristic Profile: Euler Calculus and Stability for Topological Interaction of Ball Unions

技術要約：交差オイラー標数プロファイル（The Intersection Euler Characteristic Profile）

問題提起
トポロジカル・データ・アナリシス（TDA）は、複数の点集合（$X_1, \dots, X_k \subset \mathbb{R}^d$）間の空間スケールにわたるトポロジカルな相互作用を定量化するという根本的な課題に直面している。既存の手法には以下のような重大な限界がある：

• 単一クラウドのパーシステンスは、相対的な位置関係に対して盲目である。2つの離れたクラウドは、その分離距離に関わらず同一のダイアグラムを持つ。
• 和集合ベースのパーシステンス（$X \cup Y$）は、固有のトポロジーと相互作用を混同し、クラウドが離れている場合でも相対的な位置に反応してしまう。
• マルチパラメータ・パーシステンス（個別のスケールによるフィルタリング）は、計算量的に困難であり、区間分解（バーコード）が存在せず、NP困難なインターリービング距離を持つ。
• 近年の相互作用不変量（Mixup barcodesやChromatic TDAなど）は、非対称であったり、次元削減が重かったり（単体の数に対して3次の計算量が必要）、安定性やサンプリング理論を欠いていたりすることが多い。

目標は、5つの欲求事項（P1）分離時には消失する、（P2）高密度サンプルにおいて基礎となる形状のトポロジカル不変量に収束する、（P3）対称であり任意の $k$ に対して定義される、（P4）摂動に対して安定である、（P5）計算量的に扱いやすく、ベクトル化を介さずに統計的推論に直接利用できる、という性質を満たす、単一のスケール依存的な整数不変量を定義することである。

手法
本論文では、点集合のタプルおよびスケールベクトル $t = (t_1, \dots, t_k)$ に対して定義される交差オイラー標数プロファイル（Intersection Euler Characteristic Profile: Intersection ECP）、$\Delta\chi(t)$ を導入する。

1. オイラー・カルキュラスの枠組み： この構成は、構成可能関数（constructible functions）上のオイラー積分理論に基づいている。各クラウド $X_i$ に対して、スケールパラメータ化された指示関数 $1_{U(X_i; t_i)}$ が定義される。ここで $U(X_i; t_i)$ は、$X_i$ の各点を中心とする半径 $t_i$ の球体の和集合である。
2. 交差定理： コアとなる恒等式は、積のオイラー積分の公式から導かれる：
   $$\Delta\chi(t) = \int_{\mathbb{R}^d} \prod_{i=1}^k 1_{U(X_i; t_i)} \, d\chi = \chi\left( \bigcap_{i=1}^k U(X_i; t_i) \right)$$
   これは、指示関数の積のオイラー積分が、幾何学的な球体和集合の共通部分（intersection）のオイラー標数に等しいことを確立している。
3. 可換正方形： この構成は、幾何学的な交差（$\cap$）と代数的な乗算（$\prod$）がオイラー積分下で互換可能であるという可換正方形を経由して構築される。
4. 精緻化： ベティ数が相殺されてオイラー標数が消失する場合（例：ドーナツ型の環状構造）に対処するため、本論文では、対称化された「leave-one-out」相対モジュール（$RPH_*$）を用いた相対ホモロジー精緻化を提案している。これにより、$\Delta\chi$ が圧縮してしまう「隠れた」相互作用を捉えることができる。
5. 相互作用スペクトル： 本論文では、少なくとも $j$ 個のクラウドが重なっている領域を $C_j$ としたとき、次数付きの記述子ファミリー $\chi(C_j)$ を定義する。$\Delta\chi$ は最上層（$j=k$）に対応する。

主な貢献

• 普遍性と標準性： 交差ECPは、分離（空の重なりにおいて消失する）および正規化（単一の可縮なレンズにおいて1となる）を満たす、唯一の点値オイラー相互作用プロファイルであることが証明されている。さらに、多値関数論（ハドウィジャーの定理の一般化）により、この要約はオイラー標数そのものであることが公理的に強制される。
• 安定性： プロファイルは、点集合のハウスドルフ距離に関して $L^1$ 安定かつボトルネック安定であることが証明されている。具体的には、基礎となる交差フィルタリングは、クラウドがハウスドルフ距離で $\epsilon$ 近い場合、$\epsilon$-インターリーブである。
• サンプリングの一貫性： 正の到達範囲（positive reach）仮説の下で、プロファイルは基礎となる形状のトポロカルな相互作用を一貫して回復する。それは、正則なスケールにおいて永続的に、またサンプルに対して正則なスケールにおいて点ごとに、基礎となる集団のオイラー標数を正確に回復する。これには明示的なサンプル複雑度の境界も示されている。
• アルゴリズムの効率性： 本論文では、Alpha複体の一回のソート済みスイープを用いて対角プロファイル $\Delta\chi(r)$ を計算するアルゴリズムを提示している。
  • 計算量： $O(n^{\lceil d/2 \rceil} \log n)$ の時間計算量と $O(n^{\lceil d/2 \rceil})$ の空間計算量（$n$ は全点数）。
  • 利点： Mixup barcodesやChromatic TDAとは異なり、この手法はパーシステンス減少（境界行列の計算）を必要としないため、単体の数に対する3次の計算量を回避できる。偶数次元においては、フィルタリングを構成するアルゴリズムの中で最悪計算量として最適である。

結果

• 幾何学的識別： プロファイルは、単一クラウドのパーシステンスでは不可能な、幾何学的構成（例：重なっているか、あるいは囲んでいるか）を識別する。
• フラクタル成長率： 本論文は、スケール $r \to 0$ のときの $|\Delta\chi|$ の成長率に関する境界を確立し、それを相互作用領域のボックス計数次元と結びつけている。
• 粗視化の限界： 本論文は、$\Delta\chi$ が「厳格な粗視化（strict coarsening）」であることを明確に示している。$\chi=0$ となるトポロジー的に非自明な重なり（例：アニュラスやトーラスの殻）において、$\Delta\chi$ は消失し得る。このような場合、相互作用を回復するには相対ホモロジーによる精緻化が必要となるが、それにはより高い計算コストを要する。
• 比較： 交差ECPは、対称性、$k$-一般性、幾何学的正確性、計算効率（減少なし）、および推論への即時性（スカラー曲線）を同時に満たす唯一の不変量であることが示されている。

意義と主張

本論文は、交差ECPが単なる新しい記述子ではなく、TDAの「測定領域」のために特別に設計された**測定器（measurement instrument）**であると主張している。その意義は以下の点にある：

1. 公理的に強制されている： 自然な分離および正規化の公理の下で、対称でスケール依存的な相互作用不変量を定義する唯一の解である。
2. 大規模データに対して計算可能である： パーシステンス減少を回避することで、3次の計算量によって失敗する規模のデータ（例：ニューラルネットワークや力学系における多くのレイヤー、エポック、またはシードにわたる解析）の分析を可能にする。
3. 統計的に即時利用可能である： ダイアグラムではなくスカラー曲線であるため、損失のあるベクトル化ステップを経ることなく、統計的推論や学習に直接投入できる。

著者らは、交差ECPを記述子のスペクトルの「最上層」として位置づけており、次のような原理的なトレードオフを提示している：それは、パーシステンスダイアグラムの階層的なホモロジーの詳細を犠牲にする代わりに、正確な整数の回復性、安定性、および計算の実行可能性を得ている。オイラー標数が相殺される場合には、相対ホモロジーへの「エスカレーション」が原理的に提供されているが、$\Delta\chi$ がスケーラブルな相互作用測定のための主要な不変量であるとしている。