Exact reconstruction formulas for a Radon transform over cones

技術要約：円錐上のラドン変換に対する厳密な再構成公式

問題設定
本論文は、円錐ラドン変換として知られる $R^{(p)}$ と呼ばれる特定の変種に関する数学的な逆変換を取り扱う。この変換は、上半空間 $\mathbb{R}^{d-1} \times (0, \infty)$ 内に支持される滑らかな関数 $f$ を、片側円錐面上にわたって積分するものである。これらの円錐は、超平面 $\mathbb{R}^{d-1} \times \{0\}$ 上に頂点を持ち、対称軸がこの平面に対して垂直（$y$ 軸に平行）であり、半開角 $\theta \in (0, \pi/2)$ を持つものとして特徴付けられる。積分には、円錐の生成線に沿った距離 $s$ に関する重み $1/s^p$ が含まれる。

この変換を研究する主な動機は、コンプトン散乱を考慮したコンピュータ断層撮影（CT）にある。このような撮像法では、散乱光子の検出により、光源の位置を線や平面ではなく、円錐の表面へと遡って追跡することが可能になる。3次元における特定のケース（$p \in \{0, 2\}$）についてはフーリエ再構成公式が存在し、2次元においては逆投影公式が存在するが、任意の空間次元 $d \geq 2$ および任意の重みパラメータ $p \in \mathbb{R}$ に対する一般的なケースについては、これまで研究されていなかった。

手法
著者の Markus Haltmeier は、フーリエ変換とハンケル変換を活用した調和解析の手法を用いて、逆転公式を導出している。手法の核となるのは、古典的なラドン変換における有名なフーリエスライス定理と同様の、「フーリエスライス恒等式」を円錐ラドン変換に対して確立することである。

導出は以下の3つの論理的ステップで進められる：

1. 非加重ケースへの還元： 任意の重み $p$ を持つ変換 $R^{(p)}$ が、$\cos(\theta)$ と $y^{-p}$ を含む単純な乗法作用素を用いて、非加重の変換 $R^{(0)}$ に関連付けられるという補題を確立する。これにより、著者らは $p=0$ のケースに焦点を当てて導出を行い、後に結果を一般化することが可能となる。
2. フーリエスライス恒等式の導出： 空間変数 $u \in \mathbb{R}^{d-1}$ に関するフーリエ変換を適用し、単位球面 $S^{d-2}$ 上のベッセル関数の積分表示を利用することで、元の関数のフーリエ・ハンケル変換と、円錐ラドン変換のデータのフーリエ変換を結びつける恒等式を導出する。
3. 逆投影による反転： 導出されたフーリエスライス恒等式を用い、逆フーリエ変換および逆ハンケル変換を適用することで、関数を反転させる。このプロセスにより、フィルタリング付き逆投影型の明示的な再構成公式が得られる。「フィルタリング」のステップは、第1引数に関するリースのポテンシャル $I^{(1-d)}$ の適用として特定される。

主な貢献と結果
本論文は、$R^{(p)}f$ から $f(x, y)$ を回復するための2つの主要な明示的再構成公式を提供している：

1. 角度パラメータ化 (Theorem 1.1, Eq. 1.3):
   $$f(x, y) = \frac{y^p}{(2\pi)^{d-1}} \int_0^{\pi/2} \frac{1}{\cos(\theta)^{1+p}} \int_{S^{d-2}} \left( I^{(1-d)} R^{(p)} f \right)(x + y \tan(\theta) n, \theta) \, dn \, d\theta$$
   この公式は、再構成点 $(x, y)$ を通過するすべての円錐を、開き角 $\theta$ と方向 $n$ によってパラメータ化して積分する。

2. 空間パラメータ化 (Theorem 1.1, Eq. 1.4):
   $$f(x, y) = \frac{1}{(2\pi)^{d-1}} \int_{\mathbb{R}^{d-1}} \frac{(|u-x|^2 + y^2)^{\frac{p-1}{2}}}{|u-x|^{d-2}} \left( I^{(1-d)} R^{(p)} f \right)\left(u, \arctan\frac{|u-x|}{y}\right) \, du$$
   この形式は、底面平面 $\mathbb{R}^{d-1}$ 上で積分を行い、頂点が $u$ に位置する円錐の寄与を合算する。

特定の次元について:

• $d=2$ の場合： リースのポテンシャル $I^{-1}$ は、微分とヒルベルト変換の合成（$\partial_u H_u$）に対応する。
• $d=3$ の場合： リースのポテンシャル $I^{-2}$ は、負のラプラシアン（$-\Delta_u$）に対応する。したがって、3次元の再構成公式は、データにラプラシアンを適用した後に円錐逆投影を行うことを含む。

意義と主張
本論文は、これらが任意の次元 $d \geq 2$ および任意の実数重み $p$ における、円錐ラドン変換に対する最初のフィルタリング付き逆投影型の明示的な再構成公式であると主張している。

• 局所性： 著者は、次元 $d$ のパリティ（奇偶）に基づく区別に注目している。奇数次元の場合、再構成は局所的であり、これは $f$ の値が、その点の近傍を通過する円錐の積分のみに依存することを意味する。偶数次元の場合、再構成は非局所的であり、すべての円錐面上での積分を必要とする。この挙動は、古典的なラドン変換と同様である。
• 新規性： $d=2$ における同等の公式は他の著作ですでに導出されているが、$d=3$ およびそれ以上の高次元における公式は、新たな貢献として提示されている。
• 実用性： 著者は、3次元の再構成公式が特に単純であり、ラプラシアンの適用と円錐逆投影から構成されることを強調している。また、数値的な実装は、古典的または球面ラドン変換で確立された手法に従えば、容易であると示唆している。ただし、数値的研究については今後の課題としている。

本研究は、コンプトンカメラを用いた放出断層撮影のための理論的基礎として位置づけられており、そのような装置によって生成されるデータを反転させるために必要な数学的ツールを提供している。