[Insight-developers] Gaussian filter performance

[Joshua Cates]

Hi,

Note that the kernel size of the discrete Gaussian is very much dependent
on the MaximumError parameter setting.  You can constrain the kernel to a
certain size with the MaximumKernelWidth parameter.

The discrete Gaussian kernel is based on work by Tony Lindeberg (Discrete
Scale-Space Theory and the Scale-Space Primal Sketch.  Dissertation. Royal
Institute of Technology, Stockholm, Sweden. May 1991.), which has some
nice properties over traditional sampled discretizations of the Gaussian.  

Josh.

______________________________
 Josh Cates			
 School of Computer Science	
 University of Utah
 Email: cates@sci.utah.edu
 Phone: (801) 587-7697
 URL:   http://www.sci.utah.edu/~cates


On Fri, 9 May 2003, Paul Yushkevich wrote:

> I have been trying to choose an appropriate Gaussian blurring filter to 
> use in the SnAP project.  While doing so, I tried to compare the 
> performance of three candidates from ITK and VTK on a sample image.  The 
> candidates were itk::DiscreteGaussianImageFilter, 
> itk::RecursiveGaussianImageFilter and vtkImageGaussianSmooth .  I ran 
> the filters with default parameters, only changing the standard deviation.
> 
> Perhaps the result is useful to others, so I post it here.  For one 
> thing, it shows that the recursive gaussian filter seems to perform 
> better on even very small sigmas.  What is the error of the recursive 
> filter w.r.t to the true solution for subpixel values of sigma?
> 
> The times listed below are for VC6 Release program running on a 
> two-month old machine.
> 
> Loaded image     : image01.mha
> Image size       : 7109137 bytes
> Image dimensions : 181 x 217 x 181
> Component count  : 1
> File type        : Binary
> File byte order  : Little Endian
> 
> Standard Deviation: 0.5
>    ITK Discrete Gaussian   : 3688ms.
>    ITK Recursive Gaussian  : 515ms.
>    VTK Gaussian            : 1907ms.
> 
> Standard Deviation: 1.5
>    ITK Discrete Gaussian   : 4578ms.
>    ITK Recursive Gaussian  : 484ms.
>    VTK Gaussian            : 2281ms.
> 
> Standard Deviation: 2.5
>    ITK Discrete Gaussian   : 5672ms.
>    ITK Recursive Gaussian  : 485ms.
>    VTK Gaussian            : 2406ms.
> 
> Standard Deviation: 3.5
>    ITK Discrete Gaussian   : 7079ms.
>    ITK Recursive Gaussian  : 500ms.
>    VTK Gaussian            : 2781ms.
> 
> Standard Deviation: 4.5
>    ITK Discrete Gaussian   : 8578ms.
>    ITK Recursive Gaussian  : 484ms.
>    VTK Gaussian            : 3078ms.
> 
> Paul.
> 
> 
> _______________________________________________
> Insight-developers mailing list
> Insight-developers@public.kitware.com
> http://public.kitware.com/mailman/listinfo/insight-developers
>