[論文レビュー] Localized deconvolution on the sphere

Dominique Picard|arXiv (Cornell University)|Aug 13, 2009

Numerical methods in inverse problems参考文献 15被引用数 6

ひとこと要約

本稿では、球面上での局所的デコンボリューション手法を提案する。この手法は、特異値分解（SVD）逆問題と2次世代ウェーブレット基底におけるしきい値処理を組み合わせ、L^p損失における適応的ミニマックスレートを達成する。このアプローチは、球面調和関数のグローバルなサポートの制限を克服し、滑らかさやスパarsityの変動する条件下でも最適な性能を発揮する。

ABSTRACT

AbstractWe provide a new algorithm for the treatment of deconvolutionon the sphere which combines the traditional SVD inversion withan appropriate thresholding technique in a well chosen new basis.We establish upper bounds for the behaviour of our procedure forany L p loss. It is important to emphasize the adaptation propertiesof our procedures with respect to the regularity (sparsity) of theobject to recover as well as to inhomogeneous smoothness. We alsoperform a numerical study which proves that the procedure showsvery promising properties in practice as well. Key words and phrases: statistical inverse problems, minimax estima-tion, second- generation waveletsMSC 2000 Subject Classi cation 62G05 62G08 62G20 62C10 1 Introduction The spherical deconvolution problem was rst proposed by Rooij and Ruym-gaart (1991) [14] and subsequently solved in Healy et al. (1998) [3]. Kimand Koo (2002) [7] established minimaxity for the L 2 -rate of convergence.The optimal procedures obtained there are using orthogonal series meth-ods associated with spherical harmonics. One important problem arisingwith these procedures is their poor local performances due to the fact thatspherical harmonics are spread all over the sphere. This explains for in-stance the fact that although they are optimal in the L

研究の動機と目的

• 従来の球面調和関数に基づくデコンボリューション手法がグローバルサポートを持つため、局所的性能が劣ることを是正すること。
• 球面上の潜在関数の正則性およびスパarsityに適応するデコンボリューション手順を開発すること。
• 非一様滑らかさ下でのL^p損失におけるミニマックス最適収束レートを達成すること。
• 球面上の逆問題に対して数値的に安定かつ実用的に効果的な手法を提供すること。

提案手法

• 本手法は、球面上での再構築プロセスを局所化するために2次世代ウェーブレット基底を用いる。
• SVD逆問題とウェーブレット領域におけるしきい値処理を組み合わせ、解の安定化とノイズ増幅の低減を図る。
• 目的関数の滑らかさとスパarsityの変動に適応できるように手順を設計する。
• ウェーブレットベースの関数表現としきい値処理ルールを用いて、L^pリスクの理論的上界を導出する。
• 2次世代ウェーブレットのマルチスケール構造を活用し、計算効率と局所化を実現する。
• 数値実験を通じて、手法の実用的性能とロバストネスを検証する。

実験結果

リサーチクエスチョン

• RQ1 球面上のデコンボリューション手法は、L^p損失下でミニマックス最適レートを達成しつつ、局所化を維持できるか？
• RQ2 提案手法は、回復すべき関数の非一様滑らかさおよびスパarsityにどのように適応するか？
• RQ3 古典的な球面調和関数に基づく手法と比較して、本手法の実用的性能はいかがなものか？
• RQ4 局所的ウェーブレット基底におけるしきい値処理は、球面上のデコンボリューションにおける安定性と精度を向上させられるか？

主な発見

• 提案手法は、L^p損失におけるミニマックス最適収束レートを達成しており、理論的下界と一致する。
• 本手法は、潜在関数のスパarsityおよび非一様滑らかさの両面において、優れた適応性を示す。
• 2次世代ウェーブレットの使用により、球面調和関数のグローバルサポートの制限を克服した局所的再構築が可能になる。
• 数値的検証から、有限標本設定下でも本手法の実用的有効性とロバストネスが確認された。
• しきい値処理により、不適切な逆問題におけるノイズ増幅が低減され、安定性が維持される。