[論文レビュー] Inversion of the Spherical Mean Transform with Sources on a Hyperplane

Aleksei Beltukov|ArXiv.org|Oct 8, 2009

Numerical methods in inverse problems参考文献 12被引用数 19

ひとこと要約

本稿では、n次元上半空間における超平面に位置する源を持つ球面平均変換の新しい逆変換公式を提示する。この公式は、球面平均変換と放物面カーネルによる畳み込みを結ぶ相互作用作用素を用いて構築される。主な貢献は、フーリエ解析および古典的ラドン理論に基づく明示的で安定な再構成公式であり、任意の急速減少関数に対して有効であり、任意の次元で適用可能である。

ABSTRACT

The object of this study is an integral operator $\mathcal{S}$ which averages functions in the Euclidean upper half-space $\mathbb{R}_{+}^{n}$ over the half-spheres centered on the topological boundary $\partial \mathbb{R}_{+}^{n}$. By generalizing Norton's approach to the inversion of arc means in the upper half-plane, we intertwine $\mathcal{S}$ with a convolution operator $\mathcal{P}$. The latter integrates functions in $\mathbb{R}^{n}$ over the translates of a paraboloid of revolution. Our main result is a set of inversion formulas for $\mathcal{P}$ and $\mathcal{S}$ derived using a combination of Fourier analysis and classical Radon theory. These formulas appear to be new and are suitable for practical reconstructions.

研究の動機と目的

• 上半空間 ℝⁿ⁺¹ における超平面上に源を持つ球面平均変換の新しい明示的逆変換公式を開発すること。
• 2次元の場合のノートンの手法を、任意の次元 n ≥ 1 に一般化すること。
• 球面平均変換と古典的ラドン変換との間に、畳み込み作用素を介して直接的な関係を確立すること。
• トモグラフィー応用に適した安定なバックプロジェクション型再構成アルゴリズムを提供すること。

提案手法

• 球面平均変換 𝒮 と ℝⁿ⁺¹ 上の関数に作用する畳み込み作用素 𝒫 を結ぶ相互作用作用素のペアを導入する。
• 𝒫 を回転放物面の平行移動に沿った積分として定義する畳み込み作用素とする。
• フーリエ解析および分布論を用いて、𝒫 の逆作用素 𝒫⁻¹ を導出する。
• 𝒮⁻¹ = 𝒮* ∘ 𝒦_y ∘ 𝒮 の形で球面平均変換の逆作用素を確立する。ここで 𝒦_y は古典的 Λ-作用素に関連する1次元作用素である。
• 古典的ラドン逆変換理論を活用して、逆変換公式の構造と有効性を裏付ける。
• 微分幾何的技法（体積形式の引き戻し、ヤコビ行列式の計算など）を用いて公式を検証する。

実験結果

リサーチクエスチョン

• RQ1 球面平均変換が超平面上の源を持つ場合、関数定義域に強い解析的制約を課さずに逆変換可能か？
• RQ2 球面平均変換を畳み込み作用素を介して古典的ラドン変換とどのように関連づけられるか？
• RQ3 任意の次元 n ≥ 1 に適用可能な明示的で安定な球面平均変換の逆変換公式を導出可能か？
• RQ4 逆変換公式を、数値的安定性を向上させるために随伴作用素 𝒮* を含むバックプロジェクション形式で表現可能か？
• RQ5 特に放物面畳み込みカーネルの観点から、逆作用素の正確な幾何学的・解析的構造は何か？

主な発見

• 本稿では、シュワーツ空間 𝒮(ℝⁿ⁺¹) 内のすべての急速減少関数に対して、制限的な解析的仮定を課さずに有効な球面平均変換 𝒮 の新しい逆変換公式を導出する。
• 逆作用素は 𝒮⁻¹(f) = 𝒮* ∘ 𝒦_y ∘ 𝒮(f) の形で表され、𝒦_y は古典的 Λ-作用素の1次元的引き戻しであるため、単純性と安定性を保証する。
• 𝒮 と相互作用する畳み込み作用素 𝒫 は、分布論を用いて可逆であることが示され、𝒫⁻¹ の明示的表現が得られる。
• 体積形式の引き戻しに用いられる変換のヤコビアン行列式は 1/(1+|z|²)⁽ⁿ⁺¹⁾/² に等しく、逆変換に適切なスケーリングが保証されている。
• 本手法により、球面平均変換と古典的ラドン変換との間の直接的で自然な関係が確立され、後者により逆変換の理論的基盤が提供される。
• 構成的構造とフーリエに基づく数値アルゴリズムとの整合性を考慮すると、提案された逆変換公式は実用的トモグラフィー再構成に適している。