Denoise an image using min/max curvature flow. More...

#include <itkMinMaxCurvatureFlowImageFilter.h>

Inheritance diagram for itk::MinMaxCurvatureFlowImageFilter< TInputImage, TOutputImage >:

Collaboration diagram for itk::MinMaxCurvatureFlowImageFilter< TInputImage, TOutputImage >:

Public Types
typedef SmartPointer< const Self >	ConstPointer
typedef Superclass::FiniteDifferenceFunctionType	FiniteDifferenceFunctionType
typedef MinMaxCurvatureFlowFunction < OutputImageType >	MinMaxCurvatureFlowFunctionType
typedef Superclass::OutputImageType	OutputImageType
typedef SmartPointer< Self >	Pointer
typedef FiniteDifferenceFunctionType::RadiusType	RadiusType
typedef RadiusType::SizeValueType	RadiusValueType
typedef MinMaxCurvatureFlowImageFilter	Self
typedef CurvatureFlowImageFilter < TInputImage, TOutputImage >	Superclass
Public Member Functions
virtual ::itk::LightObject::Pointer	CreateAnother (void) const
virtual const char *	GetNameOfClass () const
	typedef (Concept::GreaterThanComparable< typename TOutputImage::PixelType, unsigned long >) OutputGreaterThanUnsignedLongCheck
	typedef (Concept::MultiplyAndAssignOperator< typename TOutputImage::PixelType, double >) OutputDoubleMultiplyAndAssignOperatorCheck
	typedef (Concept::Comparable< typename TOutputImage::PixelType, double >) OutputDoubleComparableCheck
	typedef (Concept::Convertible< long, typename TOutputImage::PixelType >) LongConvertibleToOutputCheck
	typedef (Concept::LessThanComparable< typename TOutputImage::PixelType >) OutputLessThanComparableCheck
	typedef (Concept::Convertible< unsigned long, typename TOutputImage::PixelType >) UnsignedLongConvertibleToOutputCheck
	typedef (Concept::AdditiveOperators< unsigned long, typename TOutputImage::PixelType >) UnsignedLongOutputAditiveOperatorsCheck

virtual void	SetStencilRadius (RadiusValueType _arg)
virtual RadiusValueType	GetStencilRadius () const
Static Public Member Functions
static Pointer	New ()
Static Public Attributes
static const unsigned int	ImageDimension = Superclass::ImageDimension
Protected Member Functions
virtual void	InitializeIteration ()

	MinMaxCurvatureFlowImageFilter ()
	~MinMaxCurvatureFlowImageFilter ()
void	PrintSelf (std::ostream &os, Indent indent) const
Private Member Functions
	MinMaxCurvatureFlowImageFilter (const Self &)
void	operator= (const Self &)
Private Attributes
RadiusValueType	m_StencilRadius

Detailed Description

template<class TInputImage, class TOutputImage>
class itk::MinMaxCurvatureFlowImageFilter< TInputImage, TOutputImage >

Denoise an image using min/max curvature flow.

MinMaxCurvatureFlowImageFilter implements a curvature driven image denoising algorithm. Iso-brightness contours in the grayscale input image are viewed as a level set. The level set is then evolved using a curvature-based speed function:

$I_t = F_{\mbox{minmax}} |\nabla I|$

where $F_{\mbox{minmax}} = \max(\kappa,0)$ if $\mbox{Avg}_{\mbox{stencil}}(x)$ is less than or equal to $T_{thresold}$ and $\min(\kappa,0)$ , otherwise. $\kappa$ is the mean curvature of the iso-brightness contour at point $ x $ .

In min/max curvature flow, movement is turned on or off depending on the scale of the noise one wants to remove. Switching depends on the average image value of a region of radius $ R $ around each point. The choice of $ R $ , the stencil radius, governs the scale of the noise to be removed.

The threshold value $T_{threshold}$ is the average intensity obtained in the direction perpendicular to the gradient at point $ x $ at the extrema of the local neighborhood.

This filter make use of the multi-threaded finite difference solver hierarchy. Updates are computed using a MinMaxCurvatureFlowFunction object. A zero flux Neumann boundary condition is used when computing derivatives near the data boundary.

Warning:: This filter assumes that the input and output types have the same dimensions. This filter also requires that the output image pixels are of a real type. This filter works for any dimensional images, however for dimensions greater than 3D, an expensive brute-force search is used to compute the local threshold.

Reference: "Level Set Methods and Fast Marching Methods", J.A. Sethian, Cambridge Press, Chapter 16, Second edition, 1999.

See also:: MinMaxCurvatureFlowFunction; CurvatureFlowImageFilter; BinaryMinMaxCurvatureFlowImageFilter

Definition at line 77 of file itkMinMaxCurvatureFlowImageFilter.h.