如何使用Gabor滤波器进行特征提取？

Question

我有一组灰度图像（96,96），我想逐像素地对整个图像应用gabor滤镜。我听说gabor滤波器在数据上有最大的变化，对方向和闪电不敏感。

这是我使用的一段代码。我将滤镜应用于一个图像，输出数据是应用的Image + Gabor滤镜。

我尝试将此gabor过滤器应用的图像提供给机器学习算法，并且分类准确性没有变化或改进。为什么呢？

import cv2
import numpy as np
import pandas as pd
import pylab as pl

df = pd.read_csv('test.csv',header=0)

df['Image'] = df['Image'].apply(lambda im: np.fromstring(im, sep=' ') )
X = np.vstack (df['Image'].values) 

X = X.astype(np.uint8)
X = X.reshape(-1,96,96)


def image_histogram_equalization(image, number_bins=256):
    # from http://www.janeriksolem.net/2009/06/histogram-equalization-with-python-and.html

    # get image histogram
    image_histogram, bins = np.histogram(image.flatten(), number_bins, normed=True)
    cdf = image_histogram.cumsum() # cumulative distribution function
    cdf = 255 * cdf / cdf[-1] # normalize

    # use linear interpolation of cdf to find new pixel values
    image_equalized = np.interp(image.flatten(), bins[:-1], cdf)

    return image_equalized.reshape(image.shape), cdf

def build_filters():
 filters = []
 ksize = 31
 for theta in np.arange(0, np.pi, np.pi / 16):
  kern = cv2.getGaborKernel((ksize, ksize), 4.0, theta, 10.0, 0.5, 0, ktype=cv2.CV_32F)
  kern /= 1.5*kern.sum()
  filters.append(kern)
 return filters

def process(img, filters):
 accum = np.zeros_like(img)
 for kern in filters:
  fimg = cv2.filter2D(img, cv2.CV_8UC3, kern)
  np.maximum(accum, fimg, accum)
 return accum

filters = []

for k in xrange(5):
 img = X[k]

 X[k, :, :] = image_histogram_equalization(X[k, :,:])[0]

 filters = build_filters()
 filters = np.asarray(filters)
 print filters.shape
 res1 = process(img, filters)
 cv2.imshow('result', res1)
 cv2.waitKey(0)
 cv2.destroyAllWindows()

f = np.asarray(filters)
print 'Filters', f.shape

它输出具有16个方向的滤波器（对于一个特定频率）

(16, 31, 31) - Image 1
(16, 31, 31) - Image 2
(16, 31, 31) - Image 3
..

抱歉，我的问题可能为时过早，但我很想知道更多。

1. 我应该向机器学习分类器提供什么？只有gabor过滤器数据或过滤器本身应用于所有图像。
2. 如何调整上述代码以获得10个不同频率的滤波器响应？我想为每张图片提取16x10 = 160个滤镜。
3. 我为96x96像素图像设置了gabor内核大小为31x31。它是否正确？