[論文レビュー] Fast, exact (but unstable) spin spherical harmonic transforms

Jason D. McEwen|arXiv (Cornell University)|Jul 28, 2008

Geophysics and Gravity Measurements参考文献 27被引用数 7

ひとこと要約

本稿では、問題を2次トーラス上のフーリエ変換に再定式化することにより、球面上のスピン球面調和関数変換の高速かつ理論的正確なアルゴリズムを提示する。シャノン標本化定理を用いることで理論的正確性を達成する。O(L³)の計算量と正確性を達成しているが、L > 32では線形系の悪条件性のため、前方変換が数値的に不安定である。逆変換は安定であると予想されるが、正確な回復のための安定な前方実装がないため、検証されていない。

ABSTRACT

In many applications data are measured or defined on a spherical manifold; spherical harmonic transforms are then required to access the frequency content of the data. We derive algorithms to perform forward and inverse spin spherical harmonic transforms for functions of arbitrary spin number. These algorithms involve recasting the spin transform on the two-sphere S^2 as a Fourier transform on the two-torus T^2. Fast Fourier transforms are then used to compute Fourier coefficients, which are related to spherical harmonic coefficients through a linear transform. By recasting the problem as a Fourier transform on the torus we appeal to the usual Shannon sampling theorem to develop spherical harmonic transforms that are theoretically exact for band-limited functions, thereby providing an alternative sampling theorem on the sphere. The computational complexity of our forward and inverse spin spherical harmonic transforms scale as O(L^3) for any arbitrary spin number, where L is the harmonic band-limit of the spin function on the sphere. Numerical experiments are performed and unfortunately the forward transform is found to be unstable for band-limits above L~32. The instability is due to the poorly conditioned linear system relating Fourier and spherical harmonic coefficients. The inverse transform is expected to be stable, although it is not possible to verify this hypothesis.

研究の動機と目的

• 任意のスピン数に対して、球面上のスピン球面調和関数変換の高速かつ理論的正確なアルゴリズムを開発すること。
• 球面上の新しい標本化定理を用いて、O(L³)の計算量スケーリングを達成しながら正確性を維持すること。
• スピン数に依存せず、バンド制限が2の累乗でなければならない必要がない汎用アルゴリズムを提供すること。
• 提案されたアルゴリズムの数値的安定性、特に高調波バンド制限における安定性を調査すること。

提案手法

• 変数変換を用いて、S²上のスピン球面調和関数変換を2次トーラスT²上のフーリエ変換に再定式化する。
• トーラス上のフーリエ係数を高速フーリエ変換（FFTs）で計算し、線形変換により球面調和関数係数に変換する。
• トーラス上でのシャノン標本化定理を活用し、バンド制限関数に対して理論的正確性を保証する。
• θₜおよびφₚで定義されるグリッドベースのピクセル化を用いて関数を標本化し、FFTの適用を可能にする。
• フーリエ係数からスピン球面調和関数係数への写像に線形系を用い、L > 32で悪条件性が生じる。
• 正規化や幾何代数を含む代替定式化を検討し、数値的安定性の向上を図る。

実験結果

リサーチクエスチョン

• RQ1 任意のスピン数に対して、O(L³)スケーリングを達成する高速かつ正確なスピン球面調和関数変換を導出可能か？
• RQ2 球面変換をトーラス上のフーリエ変換に再定式化することで、球面上で理論的正確な標本化定理が得られるか？
• RQ3 なぜL ≈ 32より大きな調和バンド制限において前方変換が不安定になるのか？
• RQ4 前方変換の不安定性にもかかわらず、逆変換は数値的に安定しているのか？
• RQ5 悪条件な線形系を正則化または再定式化することで、数値的安定性を回復できるか？

主な発見

• 提案されたアルゴリズムは、トーラス上でのシャノン標本化定理への依拠により、バンド制限関数に対して理論的正確性を達成する。
• 前方および逆変換はO(L³)スケーリングを示し、任意のスピン数に対して計算的に効率的である。
• 数値実験により、計算時間のO(L³)スケーリングが確認され、理論的計算量と整合する。
• L > 32では、条件数がL=50、m₀=38で6×10¹⁴に達するなど、深刻な悪条件性が確認され、前方変換が不安定であることが判明した。
• 逆変換は安定すると予想されるが、正確な回復のための安定な前方変換が存在しないため、検証されていない。
• 数値的手法による安定化の試みは現在までに失敗しており、幾何代数などの代替定式化が安定性向上の可能性を模索されている。