ガウシアンフィルタのフーリエ変換を導出する
本記事は2021年5月6日にアップロードしたGist ⧉を基にしています。
問題
以下の式で示される空間領域のガウシアンフィルタのフーリエ変換を求めなさい。
空間領域のガウシアンフィルタ
hg(x,y)=2πσ21e−2σ2x2+y2
出典は以下の投稿です。
導出過程
略解の行間を埋めるように導出の過程を示す。j2=−1とする。
まず、xに関する積分とyに関する積分とを分離する。
=∫−∞∞∫−∞∞e−2σ2x2+y2e−j2π(ux+vy)dxdy∫−∞∞exp(−2σ2x2−j2πux)dx∫−∞∞exp(−2σ2y2−j2πvy)dy.
左の被積分関数の指数をf(x,u)と置いて平方完成すると以下のようになる。
f(x,u):===−2σ2x2−j2πux−2σ21(x2+j4πσ2ux)−2σ21((x+j2σ2πu)2+(2πσ2)2u2).
これを用いて被積分関数exp(f(x,u))は以下のように整理できる。
exp(f(x,u))==exp(−2σ21((x+j2σ2πu)2+(2πσ2)2u2))exp(−2σ21(x+j2σ2πu)2)exp(−2π2σ2u2).
右の指数関数はその指数に積分変数xを含まないから定数とみなせることに注意して、ガウス積分の結果を使えば以下のように積分できる。
∫−∞∞exp(f(x,u))dx==exp(−2π2σ2u2)∫−∞∞exp(−2σ21(x+j2πσ2u)2)dxexp(−2π2σ2u2)2πσ2.
したがって以下の通り解が得られる。
F(hg(x,y))===2πσ21∫−∞∞exp(f(x,u))dx∫−∞∞exp(f(y,v))dy2πσ21⋅exp(−2π2σ2u2)2πσ2⋅exp(−2π2σ2v2)2πσ2exp(−2π2σ2(u2+v2)).