Fast Anisotropic Gauss Filtering J. M. Geusebroek, A. W. M. Smeulders, J. van de Weijer In IEEE Transactions on Image Processing 2003. [bibtex] [pdf] [url]

We derive the decomposition of the anisotropic Gaussian in a one dimensional Gauss filter in the x-direction followed by a one dimensional filter in a non-orthogonal direction phi. So also the anisotropic Gaussian can be decomposed by dimension. This appears to be extremely efficient from a computing perspective. An implementation scheme for normal convolution and for recursive filtering is proposed. Also directed derivative filters are demonstrated. For the recursive implementation, filtering an 512x512 image is performed within 40 msec on a current state of the art PC, gaining over 3 times in performance for a typical filter, independent of the standard deviations and orientation of the filter. Accuracy of the filters is still reasonable when compared to truncation error or recursive approximation error. The anisotropic Gaussian filtering method allows fast calculation of edge and ridge maps, with high spatial and angular accuracy. For tracking applications, the normal anisotropic convolution scheme is more advantageous, with applications in the detection of dashed lines in engineering drawings. The recursive implementation is more attractive in feature detection applications, for instance in affine invariant edge and ridge detection in computer vision. The proposed computational filtering method enables the practical applicability of orientation scale-space analysis.



Download the C source code.
Download the Matlab Mex source code and dll (MS Windows) and mexglx (Linux PC).
[Latest update: January 28, 2005]

@Article{GeusebroekTIP2003,
  author       = "Geusebroek, J. M. and Smeulders, A. W. M. and van de Weijer, J.",
  title        = "Fast Anisotropic Gauss Filtering",
  journal      = "IEEE Transactions on Image Processing",
  number       = "8",
  volume       = "12",
  pages        = "938--943",
  year         = "2003",
  url          = "https://ivi.fnwi.uva.nl/isis/publications/2003/GeusebroekTIP2003",
  pdf          = "https://ivi.fnwi.uva.nl/isis/publications/2003/GeusebroekTIP2003/GeusebroekTIP2003.pdf",
  has_image    = 1
}