[論文レビュー] An algorithm to reconstruct convex polyhedra from their face normals and areas

Giuseppe Sellaroli|arXiv (Cornell University)|Dec 3, 2017

Point processes and geometric inequalities参考文献 1被引用数 4

ひとこと要約

本稿では、与えられた面法線と面積から3次元凸多面体を再構築するための数値的アルゴリズム、polyhedrec を提示する。この手法は、目標とする面積と一致するように支持パラメータ $ h_i $ を解くことで実現される。Lasserre の体積計算手法を応用して面積を計算し、明示的に導出したヤコビ行列を用いたLevenberg–Marquardtの根を求め方を採用することで、収束の安定性と一意解の保証が達成される。これは、3つの $ h_i $ 値を固定することによって実現される。

ABSTRACT

A well-known result in the study of convex polyhedra, due to Minkowski, is that a convex polyhedron is uniquely determined (up to translation) by the directions and areas of its faces. The theorem guarantees existence of the polyhedron associated to given face normals and areas, but does not provide a constructive way to find it explicitly. This article provides an algorithm to reconstruct 3D convex polyhedra from their face normals and areas, based on an method by Lasserre to compute the volume of a convex polyhedron in $\mathbb{R}^n$. A Python implementation of the algorithm is available at https://github.com/gsellaroli/polyhedrec.

研究の動機と目的

• Minkowskiの定理は凸多面体の存在を保証するが、具体的な構成法を示さないため、面法線と面積から凸多面体を構成的にアルゴリズム的に再構築する解決策を提供すること。
• 収束の信頼性を損なう、多重根や明示的ヤコビ行列計算の欠如といった、先行研究の限界を克服すること。
• ループ量子重力理論など、面法線と面積が基本的な幾何変数を占める理論物理学分野における多面体の実用的再構築を可能にすること。
• 再現可能性と使いやすさを確保するため、明確な疑似コードと公開されたPythonライブラリを提供すること。

提案手法

• 多面体の族は、各面が平面 $ \langle \mathbf{u}_i, \mathbf{x} \rangle = h_i $ 上にあることにより、支持パラメータ $ h_i $ でパrameter化される。ここで $ \mathbf{u}_i $ は外向きの面法線を表す。
• 面積は、$ \mathbb{R}^n $ 内の凸多面体の体積を計算するLasserreのアルゴリズムを応用し、ここでは多角形面の2次元面積を計算するために適用する。
• 計算された面積と目標面積との差を測るベクトル関数 $ \widetilde{g}(\mathbf{h}) $ を定義し、その根を数値的に探索する。
• 関数 $ \widetilde{g}(\mathbf{h}) $ のヤコビ行列を明示的に導出することで、根の探索プロセスの収束性と安定性が向上する。
• 翻訳不変性を回避し、一意解を得るために、多面体の内部に位置する点として原点を固定し、3つの $ h_i $ 値を固定する。
• 最適化の開始時に非退化でかつすべての面積が正である多面体を得るために、初期推定値 $ \widetilde{\mathbf{h}}_0 $ を構築する。この際、すべての面が中心点から固定距離 $ D $ 離れているようにし、$ D $ を調整して平均面積が目標値と一致するようにスケーリングする。

実験結果

リサーチクエスチョン

• RQ1 Minkowskiの定理は存在を保証するが、具体的な構成法を示さない。では、面法線と面積が与えられたとき、3次元凸多面体を明示的に再構築するにはどうすればよいか？
• RQ2 スパurious解や物理的に非現実的な解に収束しないように、目標面積を持つ多面体が得られるように、支持パラメータ $ h_i $ を信頼性高く求める数値的手法とは何か？
• RQ3 この逆問題における根の探索アルゴリズムの収束性を向上させるために、面積関数のヤコビ行列をどのように明示的に計算できるか？
• RQ4 最適化の開始段階で多面体が退化せず、すべての面積が正であるように保証する初期推定値の戦略は何か？

主な発見

• 本稿の数値的例により、面法線と面積から凸多面体を高い精度で再構築できることを示した。
• ヤコビ行列を明示的に計算することで、導出不能または暗黙的ヤコビ行列を用いた手法と比較して、収束速度と信頼性が顕著に向上した。
• 3つの $ h_i $ 値を固定することにより、Minkowskiの定理に内在する曖昧性を解消し、多面体は平行移動を除いて一意に特定される。
• 固定された中心点とスケーラブルな距離 $ D $ を用いた初期推定値戦略により、初期段階ですべての面積が正となることが保証され、退化解への収束が防止された。
• 本手法は、Bianchiら[2]の手法と比較して、複数の根を持つ3次元多様体ではなく、一意の根を持つという点で優れており、実装の詳細も完全に提供されている。
• 公開されたPython実装は https://github.com/gsellaroli/polyhedrec で入手可能であり、再現性と下流の応用への統合を可能にしている。