An Optimal Non-Orthogonal Separation of the Anisotropic Gaussian Convolution Filter

  • We give an analytical and geometrical treatment of what it means to sepa rate a Gaussian kernel along arbitrary axes in Rn, and we present a separation scheme that allows to efficiently implement anisotropic Gaussian convolution filters in arbitrary dimension. Based on our previous analysis we show that this scheme is optimal with regard to the number of memory accesses and nterpolation operations needed. Our method relies on non-orthogonal convolution axes and works com- pletely in image space. Thus, it avoids the need for an FFT-subroutine. Depending on the accuracy and speed requirements, different interpolation schemes and methods to implement the one-dimensional Gaussian (FIR, IIR) can be integrated. The algorithm is also feasible for hardware that does not contain a floating-point unit. Special emphasis is laid on analyzing the performance and accuracy of our method. In particular, we show that withot any special optimization of the source code, our method can perform anisotropic Gaussian filtering faster than methods relyin on the Fast Fourier Transform.

Download full text files

Export metadata

Additional Services

Share in Twitter Search Google Scholar
Author:C.H. Lampert, O. Wirjadi
Serie (Series number):Berichte des Fraunhofer-Instituts für Techno- und Wirtschaftsmathematik (ITWM Report) (82)
Document Type:Report
Language of publication:German
Year of Completion:2005
Year of Publication:2005
Publishing Institute:Fraunhofer-Institut für Techno- und Wirtschaftsmathematik
Date of the Publication (Server):2005/11/23
Tag:Anisotropic Gaussian filter; linear filtering; nD image processing; orientation space; separable filters
Faculties / Organisational entities:Fraunhofer (ITWM)
DDC-Cassification:5 Naturwissenschaften und Mathematik / 510 Mathematik
Licence (German):Standard gemäß KLUEDO-Leitlinien vor dem 27.05.2011