[論文レビュー] Smoothing splines on Riemannian manifolds, with applications to 3D shape space

Kwangrae Kim, Ian L. Dryden+2 more|arXiv (Cornell University)|Jan 15, 2018

Morphological variations and asymmetry参考文献 35被引用数 7

ひとこと要約

本稿では、リーマン多様体上の3次元形状データに滑らかさスプラインをフィッティングするための新規手法を提案する。この手法は、アンロール、アンラップ、ラップの技術を用い、JuppとKent（1987）の球面スプライン手法を拡張したものである。主な貢献は、ケンダールの形状空間における平行移動を1階微分方程式系の数値解法により解き、大ばらつきを示す3次元形状の軌道に対する正確なスプラインフィッティングを可能にしたことである。分子動力学的データを用いた実験的評価により、接空間近似に比べてより優れた測地線の忠実度が得られた。

ABSTRACT

There has been increasing interest in statistical analysis of data lying in manifolds. This paper generalizes a smoothing spline fitting method to Riemannian manifold data based on the technique of unrolling and unwrapping originally proposed in Jupp and Kent (1987) for spherical data. In particular we develop such a fitting procedure for shapes of configurations in general $m$-dimensional Euclidean space, extending our previous work for two dimensional shapes. We show that parallel transport along a geodesic on Kendall shape space is linked to the solution of a homogeneous first-order differential equation, some of whose coefficients are implicitly defined functions. This finding enables us to approximate the procedure of unrolling and unwrapping by simultaneously solving such equations numerically, and so to find numerical solutions for smoothing splines fitted to higher dimensional shape data. This fitting method is applied to the analysis of some dynamic 3D peptide data.

研究の動機と目的

• 曲がったリーマン多様体上での3次元形状データに対する堅牢な非パラメトリック回帰手法の欠如に対処すること。
• JuppとKent（1987）のアンロール・アンラップ・ラップ手法を球面から一般のリーマン多様体、特に3次元配置のケンダール形状空間へ拡張すること。
• 3次元形状空間における平行移動を数値的に安定して計算する手法を開発すること。これはアンロール処理に不可欠である。
• 大ばらつきを示す時間的3次元形状系列の正確な補間およびスムージングを可能にし、接空間近似による歪みを回避すること。
• ペプチドの実際の分子動力学的データを用いて手法を検証し、Procrustes接空間手法に比べて優れた性能を示すこと。

提案手法

• アンロール・アンラップ・ラップフレームワークをリーマン多様体に適応させるために、アンロール用の基準測地線路を定義する。
• 形状空間内の測地線に沿った平行移動を、同次的一階微分方程式系の数値解法により近似する。
• 微分方程式系の係数を形状空間の幾何構造によって暗黙的に定義し、曲率テンソルの明示的知識が不要な形で計算を可能にする。
• アンラップされた接空間でスムージングスプラインフィッティングを実行し、その後、そのパスを元の多様体に戻してラップする。
• 時間的形状系列に最適化された knot 配置を用いた3次および線形スプラインを採用し、主成分スコアプロットと距離比較による妥当性を検証する。
• Procrustes接空間距離と真のリーマン距離の比較により性能を評価し、極端な場合に最大75%の過大評価が生じることを示した。

実験結果

リサーチクエスチョン

• RQ1 アンロール・アンラップ・ラップ手法は、球面から任意のリーマン多様体、特に3次元形状データに一般化可能か？
• RQ2 ケンダールの形状空間における平行移動は、スプラインフィッティングに用いるためにどのように数値的に近似可能か？
• RQ3 提案手法は、高ばらつきを示す3次元形状軌道に対して、接空間ベースの手法よりもより正確な補間およびスムージングを提供するか？
• RQ4 この手法は、分子動力学的データにおける明確なコンformational状態を信頼性高く同定・モデル化できるか？
• RQ5 標準的なProcrustes接空間近似に比べて、この手法は測地線距離をどの程度正確に保持するか？

主な発見

• 本手法は3次ペプチド形状データに3次および線形スプラインをフィットさせることに成功し、先行研究と整合する4つの明確なコンformational状態を明らかにした。
• 3次スプラインモデルおよび4-knot線形スプラインが、複雑な遷移を捉えられなかったgeodesicおよび3-knot線形モデルに比べて優れたフィットを示した。
• Procrustes接空間距離は、特定の時刻で真のリーマン距離に比べ最大75%大きくなった。特にt=26からt=27の間では誤差が73%に達した。
• t=26とt=27の間のリーマン距離は0.374であったが、Procrustes接空間距離は0.648であった。これは線形近似による顕著な歪みを示している。
• 3次スプラインを用いたデータポイント間の滑らかな補間により、中間形状が正確に予測された。主成分スコアプロットでも一貫した遷移経路が観察された。
• 本手法は曲率テンソルの明示的計算を回避しており、曲率を導出するのが困難な多様体に対しても広く適用可能である。