ガボールフィルタ

2次元ガボールフィルタを漢字に作用させた例。左上から時計回りに原画像、搬送波の方向を45°づつ順に変えた4つのフィルタの作用結果、及びそれら4つのフィルタ結果を重ね合わせた図を表す。

ガボールフィルタ（英: Gabor filter）とは、画像処理のテクスチャー解析等に用いられる線型フィルタの一種。（2次元のガボールフィルタでは）画像の各点周りの局所領域において、方向毎に特定の周波数成分を抽出することができる。 虹彩認識や指紋認証にも応用されている他、哺乳類の脳の一次視覚野にある単純型細胞の活動をモデル化できることが示されている。名称はガーボル・デーネシュに因む。

定義

2次元ガボールフィルタのインパルス応答の例

ガボールフィルタ$g$は、ガウス関数$w$（ガウシアンエンベロープとも呼ばれる）と三角関数$s$（搬送波とも呼ばれる）を乗じたものとして定義される：

$g({\boldsymbol {x}};A,{\boldsymbol {u}},P)=K\exp(iP)w({\boldsymbol {x}};A)s({\boldsymbol {x}};{\boldsymbol {u}})$

ここで$K$ はガウシアンエンベロープの強度のスケール、$A$ はガウシアンエンベロープのスケールや方向を指定する線形変換、${\boldsymbol {u}}$は搬送波の（空間）周波数と方向を表す波数ベクトル、$P$ は位相オフセットを表す。 $w$と$s$の具体的な形は、1次元の場合、通常スケールパラメータ$\sigma$により

$w(x;\sigma )=\exp \left(-{\frac {x^{2}}{2\sigma ^{2}}}\right)$,

$s(x;u)=\exp \left(2\pi iux\right)$,

2次元の場合は、ガウシアンエンベロープの方向$\theta$と軸毎のスケールパラメータ$\sigma _{x},\sigma _{y}$により

$w(x,y;A(\theta ,\sigma _{x},\sigma _{y}))=\exp \left(-{\frac {x'^{2}+y'^{2}}{2}}\right)$,

$s(x,y;{\boldsymbol {u}})=\exp \left(2\pi i(x,y)\cdot {\boldsymbol {u}}\right)$

と表す。ここで$(x',y')=(x,y)A(\theta ,\sigma _{x},\sigma _{y})$で、

$A(\theta ,\sigma _{x},\sigma _{y})={\begin{pmatrix}\cos \theta &-\sin \theta \\\sin \theta &\cos \theta \end{pmatrix}}{\begin{pmatrix}{\frac {1}{\sigma _{x}}}&0\\0&{\frac {1}{\sigma _{y}}}\end{pmatrix}}$

である。${\boldsymbol {u}}=(f,0)$として、実部と虚部に分けた形も用いられる：

$\operatorname {Re} g(x,y)=K\exp \left(-{\frac {x'^{2}+y'^{2}}{2}}\right)\cos(2\pi fx+P)$,

$\operatorname {Im} g(x,y)=K\exp \left(-{\frac {x'^{2}+y'^{2}}{2}}\right)\sin(2\pi fx+P)$,

応用

画像処理では通常、異なる周波数、検出方向のガボールフィルタで構成されるフィルタバンクの形でよく利用される（これを画像に作用させて得られる出力をガボール空間という）。具体的事例は以下：

• スキャンした文書の中のテキスト部分の抽出（一般にテキスト部分が、画像部分に比べ相対的に高周波数成分が多いことを利用している）。

• 顔の表情の認識や髪型のモデリング。

• 脊椎の骨梁（trabecular bone）における方向性の分布の分析。

• 哺乳類の脳の一次視覚野にある単純型細胞の活動のモデル 。

• シーンを記述する低次元の大域特徴。

ガボールフィルタはガボールウェーブレットと直接関係しており、ガボールフィルタに回転や膨張を加えることで設計されるが、拡大は非常に時間がかかるため、ガボールウェーブレットでは適用されない。

実装例

• scikit-imageにおける実装

• MATLAB

function gb=gabor_fn(sigma_x,sigma_y,theta,lambda,psi,gamma)

sz_x=fix(6*sigma_x);
if mod(sz_x,2)==0, sz_x=sz_x+1;end
	

sz_y=fix(6*sigma_y);
if mod(sz_y,2)==0, sz_y=sz_y+1;end

[x y]=meshgrid(-fix(sz_x/2):fix(sz_x/2),fix(-sz_y/2):fix(sz_y/2));

% Rotation 
x_theta=x*cos(theta)+y*sin(theta);
y_theta=-x*sin(theta)+y*cos(theta);

gb=exp(-.5*(x_theta.^2/sigma_x^2+gamma^2*y_theta.^2/sigma_y^2)).*cos(2*pi/lambda*x_theta+psi);

参考文献

• Leveseque, Vincent (2000年12月6日). “Texture Segmentation Using Gabor Filters”. 2008年5月14日閲覧。

関連項目

• ガーボル・デーネシュ

外部リンク

• 初期視覚機能をモデル化したガボールフィルタLSI 九州工業大学大学院生命体工学研究科森江研究室
• Web-enabled image processing operators 2次元のガボールフィルタを画像に適用した結果を見ることができる。