[論文レビュー] Dynamic shape analysis and comparison of leaf growth

Stephan Huckemann|arXiv (Cornell University)|Feb 3, 2010

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

ひとこと要約

本稿では、明示的な公式が利用できない形状空間における測地線に 沿った平行移動を計算する数値的手法を提案し、動的な葉の成長軌跡の比較を可能にした。ポプラの葉に適用した結果、初期および終端形状が異なる葉がしばしば平行な測地線に 沿って成長することが判明し、植物形態学における「平行成長仮説」を支持する。

ABSTRACT

In the statistical analysis of shape a goal beyond the analysis of static shapes lies in the quantification of `same' deformation of different shapes. Typically, shape spaces are modelled as Riemannian manifolds on which parallel transport along geodesics naturally qualifies as a measure for the `similarity' of deformation. Since these spaces are usually defined as combinations of Riemannian immersions and submersions, only for few well featured spaces such as spheres or complex projective spaces (which are Kendall's spaces for 2D shapes), parallel transport along geodesics can be computed explicitly. In this contribution a general numerical method to compute parallel transport along geodesics when no explicit formula is available is provided. This method is applied to the shape spaces of closed 2D contours based on angular direction and to Kendall's spaces of shapes of arbitrary dimension. In application to the temporal evolution of leaf shape over a growing period, one leaf's shape-growth dynamics can be applied to another leaf. For a specific poplar tree investigated it is found that leaves of initially and terminally different shape evolve rather parallel, i.e. with comparable dynamics.

研究の動機と目的

• 明示的な公式が利用できないリーマン形状空間における平行移動を計算する一般化された数値的手法の開発。
• 特に時間経過に伴う葉の成長といった、生物学的形態における動的形状変化の定量的比較を可能にする。
• 葉の形状進化が測地線に従い、かつ異なる初期および終端形状を持つ個体間でそのダイナミクスが平行であるかどうかの検証。
• 特徴点ベースのケンダール形状空間が、効率的な計算および非パラメトリックな境界表現への再マッピングを可能にするかの評価。
• 「測地線仮説」を拡張して「平行仮説」にすることで、生物学的妥当性を評価する。

提案手法

• リーマン埋め込みおよび準同型から導かれる微分方程式を用いて、形状空間からアンビエントなユークリッド空間またはヒルバート空間への平行移動計算の移行。
• 曲線に沿った平行ベクトル場を定義するための共変微分方程式（例：∇γ̇W = 0）を用い、閉形式の解が得られない場合には数値的に解く。
• 2つの形状空間にこの手法を適用：Σ_ZR（角度方向を用いた2次元閉曲線）およびケンダールの特徴点ベース形状空間。
• Σ_ZRにおける形状表現に、境界の角度方向のフーリエ係数を用い、全葉の境界を非パラメトリックに再構築可能にする。
• 幾何的整合性を評価するため、正多角形の配置を用いてΣ_ZRとケンダール空間における平行移動結果を比較。
• 高次元ノイズを補正するための相関に基づく測定指標μ(v,w)（式20）を導入し、測地線の初期速度ベクトル間の平行性を定量化。

実験結果

リサーチクエスチョン

• RQ1 複雑な生物学的形状において、明示的な公式が利用できない形状空間でも、平行移動を信頼性高く計算できるか？
• RQ2 初期および終端形状が異なる葉が、形状空間における測地線経路を用いて測定した場合、類似した成長ダイナミクスを示すか？
• RQ3 角度方向に基づく境界のΣ_ZRの幾何構造は、ケンダールの形状空間と比較して、曲率および平行移動挙動においてどのように異なるか？
• RQ4 特徴点ベースの形状空間は、平行移動を介して非パラメトリック表現への再マッピングにより、全境界ダイナミクスをどの程度正確に近似できるか？
• RQ5 観察された葉の成長軌跡の類似性は、「測地線仮説」を「平行成長仮説」に拡張する根拠となるか？

主な発見

• カナダの黒ポプラの例では、初期および終端形状が異なる葉が、高いμ(v,w)値（例：be1b8で0.99、be1b3で0.96）を示すことで、非常に平行な測地線的成長ダイナミクスを示した。
• 相関に基づく平行性指標μ(v,w)は、201次元という高次元空間においても、ρ(v,w) ≈ 0.17–0.44 のベクトルが、高次元球面の統計的性質により強く平行であることを示している。
• Σ_ZRにおける閉曲線の部分空間は平坦であるのに対し、ケンダールの形状空間は断面曲率が1以上で有界であるため、根本的な幾何的差異が存在する。
• be1b8の移植された測地線に沿って投影された葉の境界は、後期段階でわずかな自己交差を示し、Σ_ZRの現在の計量が自然な成長をモデル化する上で限界を示している可能性がある。
• 本手法により、ケンダール空間における効率的計算と、その後のΣ_ZRにおける非パラメトリック境界再構築が可能となり、動的形状解析のスケーラブルなフレームワークが提供された。
• 葉be2b1における移植された測地線からの逸脱は、形状欠損への生物学的適応を示唆しており、自然な成長が理想化された測地線経路から逸脱しうることを示している。