i3dlib

Documentation

i3d::MCDEFilterFunction< VOXEL > Class Template Reference

This class implements together with class MCDEFilter the Modified Curvature Diffusion Flow anisotropic diffusion filter, which was published by Whitaker and Xue. As the other anisotropic filters, this filter smooths the image in the relatively continuos region and do not smooth the significante edges. Moreover it generates the so called morphological scale space. In the typical output image is suppresed the noise and the edges remain sharp and well localised. The diffused image is the solution of the equation $ u_t = |\nabla u| \cdot \mathrm{div} \left( g(| \nabla u|)\frac{ \nabla u }{ |\nabla u|} \right) $. The particular implementation is based on explicit finite difference solver and therefore there exist some limitations to the time step. There is no need to instantiate this class in the user program. This class is instantiated automatically in the MCDEFilter class. More...

#include <LevelSet.h>

Inheritance diagram for i3d::MCDEFilterFunction< VOXEL >:

i3d::ExplicitSchemeFunction< VOXEL > i3d::PDESolverFunction< VOXEL > List of all members.

Public Member Functions

virtual VOXEL ComputeUpdate (size_t *neighbors, size_t center)
 Compute the update value for one voxel, returns the ideal timestep for the particular voxel.
virtual void UpdateTimeStep (size_t *neighbors, size_t index, VOXEL tempTimeStep)
 Udates the Globaltime step from the value returned by ComputeUpdate() method.

Detailed Description

template<class VOXEL>
class i3d::MCDEFilterFunction< VOXEL >

This class implements together with class MCDEFilter the Modified Curvature Diffusion Flow anisotropic diffusion filter, which was published by Whitaker and Xue. As the other anisotropic filters, this filter smooths the image in the relatively continuos region and do not smooth the significante edges. Moreover it generates the so called morphological scale space. In the typical output image is suppresed the noise and the edges remain sharp and well localised. The diffused image is the solution of the equation $ u_t = |\nabla u| \cdot \mathrm{div} \left( g(| \nabla u|)\frac{ \nabla u }{ |\nabla u|} \right) $. The particular implementation is based on explicit finite difference solver and therefore there exist some limitations to the time step. There is no need to instantiate this class in the user program. This class is instantiated automatically in the MCDEFilter class.

The idea of this equation is to smooth the image in regions with high curvature (i.e. noise) and do not diffuse the image near significant edges. The diffusion is more agresive in the interior of objects and there is almost no diffusion on the objects boundaries due to the conductance function. For the detailed description of the conductance function $ g(|\nabla u|) $ see the documentation if the class PMFilter.

The typical setting of coefficients for biomedical images produced in LOM should be:
A few of iterations (1 to 10) should be sufficient to to smooth the image. due to the limitation of used explicit numerical discretisation should be the time step lower than 0.25 for 2D images and lower than 0.125 for 3D images (See also the example code). If there is no time step set, the possible largest is set automatically.

The numerical scheme which finds the solution of the equation is implemented in ExplicitScheme class. This class inherits the solution.

The implementation of this class is based on the article

R. T. Whitaker, X. W. Xue, Variable Conductance , Level Set Curvature For Image Denoising

The implementation of this class is based on the chapter 16 from book

J. A. Sethian, Level Set Methods and Fast Marching Methods , Cambridge University Press 1999.

Author:
Jan Hubený xhubeny@fi.muni.cz
Date:
2005


The documentation for this class was generated from the following files: