logo

n次元ラドン変換 📂トモグラフィ

n次元ラドン変換

定義1

$s \in \mathbb{R}^{1}$、$\boldsymbol{\theta} \in S^{n-1}$に対して、ラドン変換 $\mathcal{R} : L^{2}(\mathbb{R}^{n}) \to L^{2}(Z_{n})$を次のように定義する。

$$ \mathcal{R} f (s, \boldsymbol{\theta}) = \int\limits_{\mathbf{x} \cdot \boldsymbol{\theta} = s} f(\mathbf{x}) d \mathbf{x} $$

ここで、$Z_{n} := \mathbb{R}^{1} \times S^{n-1}$は$n+1$次元のユニットシリンダーである。

説明

$\mathcal{R} f (s, \boldsymbol{\theta})$の幾何学的意味は、$f$を原点から$s$だけ離れ、$\boldsymbol{\theta}$に垂直な全ての点で積分することである。

$$ \mathbf{x} \cdot (-\boldsymbol{\theta}) = -s \iff \mathbf{x} \cdot \boldsymbol{\theta} = s\quad \forall \mathbf{x} \in \mathbb{R}^{n} $$

上記の式が成り立つので、ラドン変換は偶関数である。

$$ \mathcal{R}f(-s, -\boldsymbol{\theta}) = \mathcal{R}f(s, \boldsymbol{\theta}) $$

別の表現

固定された$\boldsymbol{\theta}$に対して、

$$ \mathcal{R} f (s, \boldsymbol{\theta}) = \mathcal{R}_{\boldsymbol{\theta}}f (s) $$

$\boldsymbol{\theta} ^{\perp} := \left\{ \mathbf{u} : \mathbf{u} \cdot \boldsymbol{\theta} = 0 \right\}$とする。それならば、

$$ \mathcal{R} f (s, \boldsymbol{\theta}) = \int\limits_{\boldsymbol{\theta}^{\perp}} f( s \boldsymbol{\theta} + \mathbf{u}) d \mathbf{u} $$

ディラックのデルタ関数 $\delta$について、

$$ \mathcal{R} f (s, \boldsymbol{\theta}) = \int\limits_{\mathbb{R}^{n}} f( \mathbf{x} ) \delta ( \mathbf{x} \cdot \boldsymbol{\theta} - s) d \mathbf{x} $$

定理

$f \in L^{1}(\mathbb{R})$とする。すると$\mathcal{R}_{\boldsymbol{\theta}}f(s) \in L^{1}(\mathbb{R})$であり、次が成り立つ。

$$ \int\limits_{-\infty}^{\infty} \mathcal{R}_{\boldsymbol{\theta}}f(s) ds = \int\limits_{\mathbb{R}^{n}}f(\mathbf{x})d \mathbf{x} $$

証明

$\mathcal{R}_{\boldsymbol{\theta}}f(s)$は原点から$t$だけ離れ、$\boldsymbol{\theta}$に垂直な平面で$f$を積分したものである。これを全ての$s \in \mathbb{R}$に対して積分すると、$f$を$\mathbb{R}^{n}$の全ての点で積分したものと同じである。したがって、

$$ \int\limits_{-\infty}^{\infty} \mathcal{R}_{\boldsymbol{\theta}}f(s) ds = \int\limits_{\mathbb{R}^{n}}f(\mathbf{x})d \mathbf{x} \lt \infty $$


  1. Boris Rubin, Introduction to Radon Transforms With Elements of Fractional Calculus and Harmonic Analysis (2015), p127-131 ↩︎