logo

ラドン逆変換:フィルタリングバックプロジェクション(FBP) 📂トモグラフィ

ラドン逆変換:フィルタリングバックプロジェクション(FBP)

定理

f:R2Rf : \mathbb{R}^{2} \to \mathbb{R}について、次の式が成り立つ。

説明

フィルタードバックプロジェクション公式とも言う。

ffラドン変換Rf\mathcal{R}fが与えられた時、フーリエ変換バックプロジェクションを使ってffを得られるということだ。つまり、ラドン変換にフーリエ変換を適用し、S|S|を乗じた後、再びフーリエ逆変換を適用し、バックプロジェクションを行うことがラドン逆変換ということ。

証明

フーリエ逆変換の定理により、次が成り立つ。

  • f(x,y)=F21F2f(x,y) f(x,y)={\mathcal{F}_2}^{-1}\mathcal{F}_2 f(x,y)

ここで、F2\mathcal{F}_22次元フーリエ変換だ。フーリエ逆変換の定義により、上記の式の右辺は次のようになる。

  • 14π2F2f(X,Y)ei(xX+yY)dXdY \dfrac{1}{4\pi^2} \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} \mathcal{F}_2f(X,Y)e^{i(xX+yY)}dXdY

直交座標(X,Y)(X,Y)を極座標(S,θ)(S,\theta)で表す。すると、X=ScosθX=S\cos\thetaY=SsinθY=S\sin\thetaになる。そして、次が成り立つ。

  • XSXθYSYθ=S \begin{vmatrix} \frac{\partial X}{\partial S} & \frac{\partial X}{\partial \theta} \\ \frac{\partial Y}{\partial S} & \frac{\partial Y}{\partial \theta} \end{vmatrix} =|S|

従って、dXdY=SdSdθdXdY=|S|dSd\thetaであり、上記の積分を極座標で表すと次のようになる。

  • 14π20πF2f(Scosθ,Ssinθ)eiS(xcosθ+ysinθ)SdSdθ \dfrac{1}{4\pi^2} \int_{0}^{\pi} \int_{-\infty}^{\infty} \mathcal{F}_2f(S\cos\theta,S\sin\theta)e^{iS(x\cos\theta+y\sin\theta)}|S|dSd\theta

フーリエスライス定理

F2f(Scosθ, Ssinθ)=F(Rf)(S, θ) \mathcal{F}_2 f(S \cos\theta,\ S \sin\theta)=\mathcal{F}(\mathcal{R}f)(S ,\ \theta)

それにより、フーリエスライス定理により、上記の式は次のようになる。

  • 12π12π0πF(Rf)(S, θ)eiS(xcosθ+ysinθ)SdSdθ \dfrac{1}{2\pi}{\color{blue}\dfrac{1}{2\pi}} \int_{0}^{\pi} {\color{blue}\int_{-\infty}^{\infty} \mathcal{F}(\mathcal{R}f)(S, \ \theta)e^{iS(x\cos\theta+y\sin\theta)}|S|dS}d\theta

青く塗られた部分は、フーリエ逆変換の定義により、次のようになる。

  • 12π0πF1[SF(Rf)(S, θ)](xcosθ+ysinθ, θ)dθ \dfrac{1}{2\pi} \int_{0}^{\pi} {\color{blue} \mathcal{F}^{-1} \Big[ |S| \mathcal{F}(\mathcal{R}f)(S,\ \theta) \Big] (x\cos\theta+y\sin\theta,\ \theta) } d\theta

バックプロジェクション

Bf(x,y):=0πf(xcosθ+ysinθ, θ)dθ \mathcal{B}f(x,y) := \int_{0}^\pi f(x\cos\theta+y\sin\theta,\ \theta) d\theta

そして、上記の式は、バックプロジェクションの定義により、次のようになる。

  • 12B{F1[SF(Rf)(S, θ)]}(x,y) \dfrac{1}{2} \mathcal{B} \left\{ \mathcal{F}^{-1} \Big[ |S|\mathcal{F}(\mathcal{R}f)(S,\ \theta) \Big]\right\} (x,y)

従って、次を得る。

  • f(x,y)=12B{F1[SF(Rf)(S, θ)]}(x,y) f(x,y)= \dfrac{1}{2} \mathcal{B} \left\{ \mathcal{F}^{-1} \Big[ |S|\mathcal{F}(\mathcal{R}f)(S,\ \theta) \Big]\right\} (x,y)