[論文レビュー] Discrete Laplace-Beltrami Operator Determines Discrete Riemannian Metric

Xianfeng Gu, Ren Guo+2 more|arXiv (Cornell University)|Oct 19, 2010

3D Shape Modeling and Analysis参考文献 13被引用数 13

ひとこと要約

本稿では、三角メッシュとしてのユークリッド多面的表面における離散ラプラシアン＝ベルトラミー作用素が、凸エネルギー最小化に基づく変分的アプローチを用いて、スケーリングを除いて離散リーマン計量を一意に決定することを証明している。この手法は、勾配がコタンジェント辺重みと一致する凸エネルギー関数を構築し、作用素行列から計量を安定的かつ計算的に効率的に回復するためのアルゴリズムを可能にする。

ABSTRACT

The Laplace-Beltrami operator of a smooth Riemannian manifold is determined by the Riemannian metric. Conversely, the heat kernel constructed from its eigenvalues and eigenfunctions determines the Riemannian metric. This work proves the analogy on Euclidean polyhedral surfaces (triangle meshes), that the discrete Laplace-Beltrami operator and the discrete Riemannian metric (unique up to a scaling) are mutually determined by each other. Given an Euclidean polyhedral surface, its Riemannian metric is represented as edge lengths, satisfying triangle inequalities on all faces. The Laplace-Beltrami operator is formulated using the cotangent formula, where the edge weight is defined as the sum of the cotangent of angles against the edge. We prove that the edge lengths can be determined by the edge weights unique up to a scaling using the variational approach. First, we show that the space of all possible metrics of a polyhedral surface is convex. Then, we construct a special energy defined on the metric space, such that the gradient of the energy equals to the edge weights. Third, we show the Hessian matrix of the energy is positive definite, restricted on the tangent space of the metric space, therefore the energy is convex. Finally, by the fact that the parameter on a convex domain and the gradient of a convex function defined on the domain have one-to-one correspondence, we show the edge weights determines the polyhedral metric unique up to a scaling. The constructive proof leads to a computational algorithm that finds the unique metric on a topological triangle mesh from a discrete Laplace-Beltrami operator matrix.

研究の動機と目的

• 多面的表面における離散ラプラシアン＝ベルトラミー作用素と離散リーマン計量の間の一対一対応を確立すること。
• コタンジェント公式から導かれる辺重みが、スケーリングを除いて元の辺長（計量）を一意に決定することを証明すること。
• 与えられた離散ラプラシアン＝ベルトラミー作用素行列から計量を回復するための構成的で計算的に実行可能なアルゴリズムを開発すること。
• 凸最適化による計量回復の理論的理解を拡張することにより、離散微分幾何学の理論的枠組みを発展させること。

提案手法

• 辺重みが対角の余割の和に等しいコタンジェントスキームを用いて、離散ラプラシアン＝ベルトラミー作用素を定式化する。
• 三角不等式を満たすすべての許容可能な辺長の空間Ωᵤを定義し、それが凸であることを示す。
• 計量空間上での線積分を用いて、勾配が辺重みwᵢに一致する凸エネルギー関数E(u)を構築する。
• エネルギー関数Eのヘッセ行列がΩᵤの接空間上で正定値であることを証明し、厳密な凸性および最小値の一意性を保証する。
• 凸領域におけるパラメータと勾配の間の一対一対応を用いて、wがuをスケーリングを除いて一意に決定することを示す。
• エネルギー最小化に基づく数値アルゴリズムを実装し、所定の辺重みから計量を回復する。

実験結果

リサーチクエスチョン

• RQ1 多面的表面の離散リーマン計量は、その離散ラプラシアン＝ベルトラミー作用素から一意に回復可能か？
• RQ2 離散ラプラシアンの辺重みとその背後にある計量を結ぶ変分原理は存在するか？
• RQ3 計量空間とエネルギー関数の凸性が、計量回復のための一意解を保証するか？
• RQ4 離散ラプラシアン＝ベルトラミー作用素は、スケーリングを除いて計量と同等の形状記述子として機能可能か？

主な発見

• 多面的表面におけるすべての許容可能な離散計量の空間は凸であるため、最適化が安定的に行える。
• 勾配がコタンジェント辺重みと一致する凸エネルギー関数が存在し、安定的な回復が可能になる。
• エネルギー関数のヘッセ行列は計量空間の接空間上で正定値であり、一意のグローバル最小値が保証される。
• 計量から辺重みへの写像は滑らかな埋め込みであるため、辺重みが計量をスケーリングを除いて一意に決定することを証明する。
• 数値実験により、さまざまなトポロジーと genus を持つメッシュにおいて、所定の辺重みから元の計量（スケーリングを除いて）が回復されることを確認した。
• 二重被覆技術を用いて、境界を持つ表面を含む閉じた表面および開いた表面に対しても、計量の再構築が成功した。