[論文レビュー] A convergent point integral method for isotropic elliptic equations on point cloud

Zhen Li, Zuoqiang Shi|arXiv (Cornell University)|Jun 11, 2015

Numerical methods in engineering参考文献 21被引用数 4

ひとこと要約

本稿では、核関数を用いた積分表現による微分作用素の近似を通じて、滑らかな多様体上の点群上で等方的楕円型方程式を解く収束する点積分法を提案する。この手法は強制性を保ち、メッシュや幾何的事前知識なしに多様体上でのPDEの高精度な解法を可能にする。画像処理やデータ解析への応用が可能である。

ABSTRACT

In this paper, we propose a numerical method to solve isotropic elliptic equations on point cloud by generalizing the point integral method. The idea of the point integral method is to approximate the differential operators by integral operators and discretize the corresponding integral equation on point cloud. The key step is to get the integral approximation. In this paper, with proper kernel function, we get an integral approximation for the elliptic operators with isotropic coefficients. Moreover, the integral approximation has been proved to keep the coercivity of the original elliptic operator. The convergence of the point integral method is also proved.

研究の動機と目的

• メッシュや追加の幾何的情報を必要とせず、点群上での等方的楕円型方程式を解く数値解法の開発。
• ポisson方程式に限らない点積分法を、係数が変化する一般の等方的楕円型方程式へと一般化すること。
• 滑らかな多様体からの点群サンプリングにおいて、この手法の理論的収束性を確立すること。
• 離散化において元の楕円型作用素の強制性を保つことで、安定性と収束性を確保すること。
• 画像復元における多様体モデルを用いた応用、例えばパッチ多様体上の全変動最小化など、画像処理およびデータ解析への応用を可能とすること。

提案手法

• 非局所拡散理論から導かれる核関数 $ R_t(oldsymbol{x}, oldsymbol{y}) $ を用いた積分表現により、楕円型微分作用素を近似する。
• 点群上の積分方程式を、局所的な重みと核関数の評価に基づく数値積分則（四則則）で離散化する。
• 核関数を、等方的楕円型作用素の積分近似を保ちつつ、強制性を維持するように選定する。
• 非一様な点群に対する適応的核幅 $ t(oldsymbol{x}_i) = \rho(\boldsymbol{x}_i)^2 $ を用い、$ \rho(\boldsymbol{x}_i) $ は20個の最近傍点の半径とする。
• 多様体上の勾配および発散作用素を、差分 $ u(\boldsymbol{x}) - u(\boldsymbol{y}) $ を含む積分表現により近似し、幾何的構造と整合性を保つ。
• 離散系は反復的に解き、係数行列の退化を避けるために正則化を適用する。

実験結果

リサーチクエスチョン

• RQ1 点積分法は、係数が変化する等方的楕円型方程式に対して、点群上に一般化可能か？
• RQ2 積分近似は、離散設定において元の楕円型作用素の強制性を保っているか？
• RQ3 点群が密になるにつれて、得られる離散スキームは収束するか？
• RQ4 画像処理に現れるPDE、例えばパッチ多様体上の全変動最小化などに、この手法を適用可能か？
• RQ5 点積分法から導かれる離散固有値問題のスペクトル収束性はいかなるものか？

主な発見

• 提案された点積分法は、滑らかな多様体からの点群サンプリングに対して、等方的楕円型方程式で収束性を示す。
• 積分近似は元の楕円型作用素の強制性を保ち、安定性と収束性に不可欠である。
• 境界あり・なしを問わず多様体に適用可能で、メッシュ、レベルセット関数、最近点関数を必要としない。
• 画像復元における数値実験では、$ L_1 $ 拡張を用いた高品質な結果が得られ、実用的有効性が示された。
• 均一なサンプリング下で、離散固有値問題が連続問題にスペクトル収束する。固有値および固有関数の収束が確立された。
• 局所的近傍密度に基づく適応的核幅選択により、非一様な点群に対しても手法は頑健である。