[論文レビュー] Left-Invariant Diffusion on the Motion Group in terms of the Irreducible Representations of SO(3)

Marco Reisert, Henrik Skibbe|arXiv (Cornell University)|Feb 24, 2012

Advanced Neuroimaging Techniques and Applications参考文献 9被引用数 9

ひとこと要約

本稿では、$SO(3)$ の既約表現、特に球面調和関数とウィグナー D 行列を用いて、3次元運動群 $SE(3)$ 上の左不変な拡散および対流 PDE を解くための新規フレームワークを提示する。左不変なベクトル場を移動座標における微分作用素と回転座標における代数的作用素として表現することで、$SO(3)$ や $S^2$ の明示的離散化を回避し、メモリ効率の良い計算を可能にするとともに、高角分解能拡散MRIや球面ハフ変換への応用を可能にし、回転不変性を完全に保ちながら解析的効率を実現する。

ABSTRACT

In this work we study the formulation of convection/diffusion equations on the 3D motion group SE(3) in terms of the irreducible representations of SO(3). Therefore, the left-invariant vector-fields on SE(3) are expressed as linear operators, that are differential forms in the translation coordinate and algebraic in the rotation. In the context of 3D image processing this approach avoids the explicit discretization of SO(3) or $S_2$, respectively. This is particular important for SO(3), where a direct discretization is infeasible due to the enormous memory consumption. We show two applications of the framework: one in the context of diffusion-weighted magnetic resonance imaging and one in the context of object detection.

研究の動機と目的

• SO(3) や $S^2$ を離散化せずに $SE(3)$ 上の左不変 PDE を解くフレームワークの開発。これは、高メモリ要件のため計算上不可能である。
• SO(3) の調和構造を活用することで、3次元方向データ（特に拡散強調 MRI）の効率的かつメモリスケーラブルな処理を実現する。
• 微分作用素を $SO(3)$ の既約表現に基づいて表現し、回転不変性を保持するとともに、解析的演算を可能にする。
• 高角分解能拡散画像法（HARDI）および球面ハフ変換への実用的応用を示し、等積離散化との比較で優位性を示す。
• 空間的成分と回転の群論的構造を分離することで、$SE(3)$ 上の PDE における空間的および回転的成分の統一的取り扱いを実現する。

提案手法

• SE(3) 上の左不変ベクトル場を、SE(3) のリー代数を用いて、移動座標における微分作用素と回転座標における代数的作用素として表現する。
• 関数を球面調和関数とウィグナー D 行列を用いて表現し、回転部を $SO(3)$ の既約表現に分解する。
• クリービッチ＝ゴルダン係数を用いて $SO(3)$ の表現を結合し、作用素の積の分解と高次微分作用素の構築を可能にする。
• 球面微分作用素とウィグナー D 行列を用いて、$\mathcal{T}_\mp \mathcal{J}_{\pm 1}$ や $\mathcal{J}_{\pm 1}\mathcal{T}_\mp$ などの微分作用素を明示的な角運動量結合の和を含めて実装する。
• ウィグナー D 行列と球面調和関数の直交性および変換性質を活用し、回転不変性を保ちつつ計算を効率化する。
• SO(3) の既約表現による $SE(3)$ 上のフーリエ変換を用いて、回転群を離散化せずに畳み込みや PDE の表現・計算を可能にする。

実験結果

リサーチクエスチョン

• RQ1 SO(3) や $S^2$ を離散化せずに、$SE(3)$ 上の左不変な拡散および対流 PDE を定式化できるか。これにより、$S^2$ や $SO(3)$ における等積サンプリングに起因する高コストなメモリ消費を回避できるか？
• RQ2 ${\mathbb{R}^3 \times SO(3)}$ 上の拡散の生成子を、球面調和関数とウィグナー D 行列を用いて $SO(3)$ の既約表現の観点からどのように表現できるか？
• RQ3 提案フレームワークは、等積離散化と比較して、高角分解能拡散画像法（HARDI）における計算効率と精度をどの程度向上させるか？
• RQ4 球面ハフ変換を、角方向のサンプリングなしに、SO(3) の調和表現を用いて効率的に実装できるか？
• RQ5 クリービッチ＝ゴルダン係数は、微分作用素の結合と $SE(3)$ 上の高次 PDE の分解を可能にする役割を果たすか？

主な発見

• 本フレームワークにより、$SO(3)$ を離散化せずに $SE(3)$ 上の左不変 PDE を解くことが可能となり、$S^2$ や $SO(3)$ における等積サンプリングに起因する膨大なメモリコストを回避できる。
• SO(3) の既約表現の使用により、$S^2$ での畳み込みや回転のコンactかつメモリ効率の良い保存と解析的計算が可能になる。
• 本手法は完全な回転不変性を実現し、調和表現を用いた球面ハフ変換の効率的実装を可能にする。
• 本稿では、空間正則化と高分解能方向推定を角グリッドなしに可能にするという点で、等積離散化よりも優れた性能を示す、HARDI における提案手法の有効性を実証している。
• クリービッチ＝ゴルダン係数とウィグナー D 行列を用いて、$\mathcal{T}_\mp \mathcal{J}_{\pm 1}$ のような混合2次項の明示的公式が導出され、高次微分作用素の構築が可能になる。
• ウィグナー D 行列の直交性と回転下での変換則を用いて、$\overline{\mathcal{T}_{m_1}} \mathcal{T}_{m_2}$ の閉形式表現が導出され、これは $\mathbf{D}^2$ 変換を施した球面微分作用素への分解を示している。