The plastimatch filter command is a convenient method for performing common filters on images. This documentation is intended to explain in more detail the mathematics and implementation of these filters.
The Gabor filter is a complex filter, and has real and imaginary components. Plastimatch implements only a uniform and symmetric form of the Gabor kernel, which is described as:
Where x is the voxel position (relative to the center of the kernel), :sigma: is the filter width, and S is an arbitrary scaling factor. The vector k is the direction of modulation, and the frequency of modulation is given by :omega:. Note that the vector k is internally normalized to be a unit vector by plastimatch, so that larger k vectors do not increase the modulation.
The value of k can be set manually, or can be selected according to lie on a Fibanocci spiral. The Fibanocci spiral method is convenient for choosing a set of k vectors that are approximately evenly spaced on the unit sphere.
The values for :omega: and :sigma: can be interlinked. If :omega: is set, but :sigma: is not set, or vice versa, the other will be set automatically according to the following formula, which seems to work well:
At the current time, only the real component of the Gabor filter is generated.
The cannonical form for a Gaussian kernel is given as
Note: because plastimatch uses a fixed width convolution kernel, the above expression may not sum to unity. Therefore, an empirical factor S is used instead.
To be written.