[論文レビュー] Surface Comparison with Mass Transportation

Yaron Lipman, Ingrid Daubechies|arXiv (Cornell University)|Dec 17, 2009

Infrastructure Maintenance and Monitoring参考文献 18被引用数 27

ひとこと要約

本論文は、双曲的ディスク上の等角密度間の最適輸送を用いて、モビウス変換に対して不変な、2次元表面（例：脳皮質、小臼歯）を比較するための新規な尺度を提案する。この手法により、表面の内在的整列が可能となり、計算可能で数学的に整合性のある距離尺度を用いて幾何的類似度を定量化できる。これは、Gromov-Hausdorff距離のようなNP困難な代替手法を凌駕する。

ABSTRACT

We use mass-transportation as a tool to compare surfaces (2-manifolds). In particular, we determine the "similarity" of two given surfaces by solving a mass-transportation problem between their conformal densities. This mass transportation problem differs from the standard case in that we require the solution to be invariant under global Möbius transformations. Our approach provides a constructive way of defining a metric in the abstract space of simply-connected smooth surfaces with boundary (i.e. surfaces of disk-type); this metric can also be used to define meaningful intrinsic distances between pairs of "patches" in the two surfaces, which allows automatic alignment of the surfaces. We provide numerical experiments on "real-life" surfaces to demonstrate possible applications in natural sciences.

研究の動機と目的

• 幾何的類似度に基づいて、ディスク型2次元表面を対象とする堅牢で内在的な尺度の開発。
• 脳皮質や哺乳類の小臼歯のような複雑で等長でない表面の整列の課題に対処する。
• Gromov-Hausdorff距離のようなNP困難な尺度の計算可能代替手法を提供する。
• 一般化された質量輸送フレームワークを用いて、自動的かつ特徴点に依存しない対応点検出を可能にする。
• 生物学的分類や比較解剖学への応用を想定し、表面の差異および類似度を定量化する。

提案手法

• 均質化理論を用いて、各表面をその双曲的ディスク上での等角因子（計量密度）に写像する。
• 局所的近傍の点間の不一致尺度をモビウス不変に一般化した最適輸送問題を定式化する。
• 最適モビウス変換後の等角因子間のL1差を測る輸送コストを定義し、表面全体の向きの変化に対して不変性を確保する。
• 整列と距離計算を、双曲的ディスク上での線形計画法として定式化する。
• モビウス変換の離散的サンプリングと積分則を用いて、一般化された輸送コストを近似する。
• 最小輸送コストに基づくペアワイズ距離行列を構築し、クラスタリングや分類を可能にする。

実験結果

リサーチクエスチョン

• RQ1 グローバルな幾何的変換を要因から除きながら、モビウス不変な最適輸送フレームワークを構築できるか？
• RQ2 等角密度と輸送コストを用いて、等長でない表面間の内在的幾何的類似度をどのように定量化できるか？
• RQ3 本手法により、手動で定義された特徴点に依存せずに、生物学的表面間の意味のある点対応を自動で検出できるか？
• RQ4 提案された距離尺度は、例えば小臼歯の種分類のような生物学的に意味のあるクラスタリングをもたらすか？
• RQ5 Gromov-Hausdorff距離のような既存のNP困難な代替手法と比較して、本手法の計算複雑度はどの程度か？

主な発見

• 提案された距離尺度は、対称性、正定性、三角不等式を満たすため、ディスク型表面の空間上で真の距離関数であることが数学的に証明された。
• 本手法は、異なるレムール種の小臼歯間で120組の一貫性のある点対応を正確に同定し、以前は自動化に失敗していた生物学者の驚きを買った。
• 8つの小臼歯表面の距離行列に対する多次元尺度構成により、完全に種のラベルと一致するクラスタリングが得られ、強い生物学的妥当性を示した。
• 150〜200個の表面点に対して、整列のための線形計画法最適化は15〜20イテレーション（2〜3秒）で収束し、計算上の実行可能性が裏付けられた。
• 主な計算ボトルネックは距離行列の計算であり、N点に対してO(L·P·N³)の時間計算量を要した（N=300の場合約2時間）。ただし、今後さらなる最適化が期待される。