ガウシアンフィルタのフーリエ変換を導出する

本記事は2021年5月6日にアップロードしたGistを基にしています。

問題

以下の式で示される空間領域のガウシアンフィルタのフーリエ変換を求めなさい。
空間領域のガウシアンフィルタ

hg(x,y)=12πσ2ex2+y22σ2h_{\mathrm{g}}(x, y) = \frac{1}{2\pi\sigma^2} e^{-\frac{x^2+y^2}{2\sigma^2}}

出典は以下の投稿です。

導出過程

略解の行間を埋めるように導出の過程を示す。j2=1j^2 = -1とする。

まず、xxに関する積分とyyに関する積分とを分離する。

ex2+y22σ2ej2π(ux+vy)dxdy=exp ⁣(x22σ2j2πux)dxexp ⁣(y22σ2j2πvy)dy.\begin{align*} & \int_{-\infty}^{\infty} \int_{-\infty}^{\infty} e^{-\frac{x^2+y^2}{2\sigma^2}} e^{-j2\pi(ux+vy)} dx dy \\ {}={} & \int_{-\infty}^{\infty} \exp\!{\left(-\frac{x^2}{2\sigma^2}-j2\pi ux\right)} dx \int_{-\infty}^{\infty} \exp\!{\left(-\frac{y^2}{2\sigma^2}-j2\pi vy\right)} dy . \end{align*}

左の被積分関数の指数をf(x,u)f(x, u)と置いて平方完成すると以下のようになる。

f(x,u):=x22σ2j2πux=12σ2(x2+j4πσ2ux)=12σ2 ⁣((x+j2σ2πu)2+(2πσ2)2u2) ⁣.\begin{align*} f(x, u) := {} & -\frac{x^2}{2\sigma^2} - j2\pi ux \\ {}={} & -\frac{1}{2\sigma^2} (x^2 + j4\pi\sigma^2 ux) \\ {}={} & -\frac{1}{2\sigma^2}\!\left( (x+j2\sigma^2\pi u)^2+(2\pi\sigma^2)^2u^2 \right)\!. \end{align*}

これを用いて被積分関数exp(f(x,u))\exp(f(x, u))は以下のように整理できる。

exp(f(x,u))=exp ⁣(12σ2 ⁣((x+j2σ2πu)2+(2πσ2)2u2))=exp ⁣(12σ2(x+j2σ2πu)2)exp(2π2σ2u2).\begin{align*} \exp(f(x, u)) = {} & \exp\!{\left( -\frac{1}{2\sigma^2}\!\left( (x+j2\sigma^2\pi u)^2+(2\pi\sigma^2)^2u^2 \right) \right)} \\ {}={} & \exp\!{\left( -\frac{1}{2\sigma^2} (x+j2\sigma^2\pi u)^2 \right)} \exp{(-2\pi^2\sigma^2 u^2)} . \end{align*}

右の指数関数はその指数に積分変数xxを含まないから定数とみなせることに注意して、ガウス積分の結果を使えば以下のように積分できる。

exp(f(x,u))dx=exp(2π2σ2u2)exp ⁣(12σ2(x+j2πσ2u)2)dx=exp(2π2σ2u2)2πσ2.\begin{align*} \int_{-\infty}^{\infty} \exp(f(x, u)) \, dx ={} & \exp{(-2\pi^2\sigma^2 u^2)} \int_{-\infty}^{\infty} \exp\!{\left(-\frac{1}{2\sigma^2} (x+j2\pi\sigma^2 u)^2 \right)} dx \\ {}={} & \exp{(-2\pi^2\sigma^2 u^2)} \sqrt{2\pi\sigma^2} . \end{align*}

したがって以下の通り解が得られる。

F(hg(x,y))=12πσ2exp ⁣(f(x,u))dxexp ⁣(f(y,v))dy=12πσ2exp(2π2σ2u2)2πσ2exp(2π2σ2v2)2πσ2=exp ⁣(2π2σ2(u2+v2)).\begin{align*} F( h_{\mathrm{g}}(x, y) ) = {} & \frac{1}{2\pi\sigma^2} \int_{-\infty}^{\infty} \exp\!{\left(f(x, u) \right)} \, dx \int_{-\infty}^{\infty} \exp\!{\left(f(y, v) \right)} \, dy \\ {}={} & \frac{1}{2\pi\sigma^2} \cdot \exp{(-2\pi^2\sigma^2 u^2)} \sqrt{2\pi\sigma^2} \cdot \exp{(-2\pi^2\sigma^2 v^2)} \sqrt{2\pi\sigma^2} \\ {}={} & \exp\!{\left( -2 \pi^2 \sigma^2 (u^2 + v^2) \right)} . \end{align*}
  1. 問題
  2. 導出過程