[論文レビュー] Fast Finite Shearlet Transform

S. Häuser, Gabriele Steidl|arXiv (Cornell University)|Feb 8, 2012

Mathematical Analysis and Transform Methods参考文献 21被引用数 57

ひとこと要約

本稿では、高速な有限シアーレット変換（FFST）を提案する。これは、FFTに基づく並進不変な離散シアーレット変換であり、画像処理におけるシアーレット係数の効率的計算を可能にする。有限で離散的なシアーレットを用いてパーセバルフレームを構築し、フーリエスペクトルに対称性を保つ変更を加えることで、偶数サイズの画像に対しても実数値の係数が得られることを保証する。これにより、カートゥーン風画像の最適なスパース近似が達成され、$ olimits_{L_2}^2 \leq C N^{-2} (\log N)^3$ を満たす。実装はGPLライセンスのもとで自由に利用可能であり、FFTおよび適切なインデキシングによりパフォーマンス最適化がなされている。

ABSTRACT

In recent years it has turned out that shearlets have the potential to retrieve directional information so that they became interesting for many applications. Moreover the continuous shearlet transform has the outstanding property to stem from a square integrable group representation. However, to use shearlets and the shearlet transform for reasonable applications one needs fast algorithms to compute a discrete shearlet transform. In this tutorial we present the steps towards an implementation of a fast and finite shearlet transform that is only based on the FFT. Using band-limited shearlets we construct a Parseval frame that provides a simple and straightforward inverse shearlet transform. We provide all proofs and discuss several aspects of our implementation.

研究の動機と目的

• 画像処理のための高速で有限的かつ並進不変な離散シアーレット変換の開発。
• 偶数サイズの画像に対して実数値のシアーレット係数を保証するため、フーリエスペクトルを変更して共役対称性を維持する。
• FFTを用いたパフォーマンス最適化を施した実用的でオープンソースのMATLAB実装（FFST）の提供。
• 有限次元のユークリッド空間におけるパーセバルフレームを形成する離散シアーレットフレームの確立により、随伴変換による単純な逆変換が可能になる。
• 位相に基づく解析に適した実数および複素数の両方のシアーレット係数をサポートし、明確な構築ガイドラインを提供する。

提案手法

• 数学的整合性を保つために、一様で可積分な群表現を用いて、全グリッド上に離散シアーレットを構築する。
• すべてのスケールと方向におけるシアーレット係数の計算を高速化するために、高速フーリエ変換（FFT）を用いる。
• 偶数サイズの画像におけるフーリエスペクトルの共役対称性を保つために、シアーレットスペクトルの最初の列および行をミラー画像化し、$1/\sqrt{2}$ でスケーリングする。これにより、逆FFT後に実数値の係数が得られる。
• 離散シアーレットは有限次元の$L_2$空間においてパーセバルフレームを形成し、随伴変換により正確な再構成が可能になる。
• 片側周波数サポートを用いることで複素数シアーレットを構築し、位相解析を可能にしながらフレーム性質を維持する。
• 最適化されたインデキシング、スペクトル計算、およびモジュラー化されたドキュメンテーションを実装に含め、インストールと利用を容易にする。

実験結果

リサーチクエスチョン

• RQ1 FFTを用いて効率的に計算され、シアーレットの数学的性質を保つ有限的で並進不変なシアーレット変換をどのように実現できるか？
• RQ2 偶数サイズの画像に対して実数値のシアーレット係数を保証するために、フーリエスペクトルに必要な変更は何か？
• RQ3 提案された離散シアーレットシステムは、フレーム性質をどのように維持し、随伴変換による単純な逆変換を可能にするか？
• RQ4 スペクトルの対称性は、実用的な画像処理におけるシアーレット係数の数値的安定性と精度にどのような影響を与えるか？
• RQ5 本手法は、位相解析や幾何的分離などの多様な応用に適した、実数および複素数の両方のシアーレット係数をサポートできるか？

主な発見

• FFSTは、$N \to \infty$ のとき、誤差の減少が $\|f - f_N\|_{L_2}^2 \leq C N^{-2} (\log N)^3$ に従い、シアーレットシステムの理論的限界と一致する最適なスパース近似を達成する。
• シアーレットスペクトルの最初の列および行をミラー化し、スケーリングすることで、偶数サイズの画像に対しても実数値のシアーレット係数が保証され、虚数部が無視できるほど小さくなる。
• 離散シアーレットは有限次元の$L_2$空間においてパーセバルフレームを形成し、随伴変換を用いた正確な再構成が可能で、計算オーバーヘッドは最小限に抑えられる。
• パフォーマンスベンチマークの結果、FFTベースの実装は非FFTアプローチよりも顕著に高速であり、効率的なインデキシングとメモリアクセスパターンが実現されている。
• 片側周波数サポートを用いることで、複素数シアーレットが実装され、フレーム性質と計算効率を維持したまま位相に基づく解析が可能になる。
• オープンソースのMATLABパッケージはGPLライセンスのもとで自由に利用可能であり、研究者および実務家が容易にインストール・利用できる明確なドキュメンテーションとインストール手順が整えられている。