ITK
4.3.0
Insight Segmentation and Registration Toolkit
|
#include <itkRecursiveGaussianImageFilter.h>
Base class for computing IIR convolution with an approximation of a Gaussian kernel.
RecursiveGaussianImageFilter is the base class for recursive filters that approximate convolution with the Gaussian kernel. This class implements the recursive filtering method proposed by R.Deriche in IEEE-PAMI Vol.12, No.1, January 1990, pp 78-87, "Fast Algorithms for Low-Level Vision"
Details of the implementation are described in the technical report: R. Deriche, "Recursively Implementing The Gaussian and Its Derivatives", INRIA, 1993, ftp://ftp.inria.fr/INRIA/tech-reports/RR/RR-1893.ps.gz
Further improvements of the algorithm are described in: G. Farneback & C.-F. Westin, "On Implementation of Recursive Gaussian Filters", so far unpublished.
As compared to itk::DiscreteGaussianImageFilter, this filter tends to be faster for large kernels, and it can take the derivative of the blurred image in one step. Also, note that we have itk::RecursiveGaussianImageFilter::SetSigma(), but itk::DiscreteGaussianImageFilter::SetVariance().
Definition at line 62 of file itkRecursiveGaussianImageFilter.h.
Static Public Member Functions | |
static Pointer | New () |
Private Member Functions | |
void | ComputeDCoefficients (ScalarRealType sigmad, ScalarRealType W1, ScalarRealType L1, ScalarRealType W2, ScalarRealType L2, ScalarRealType &SD, ScalarRealType &DD, ScalarRealType &ED) |
void | ComputeNCoefficients (ScalarRealType sigmad, ScalarRealType A1, ScalarRealType B1, ScalarRealType W1, ScalarRealType L1, ScalarRealType A2, ScalarRealType B2, ScalarRealType W2, ScalarRealType L2, ScalarRealType &N0, ScalarRealType &N1, ScalarRealType &N2, ScalarRealType &N3, ScalarRealType &SN, ScalarRealType &DN, ScalarRealType &EN) |
void | ComputeRemainingCoefficients (bool symmetric) |
void | operator= (const Self &) |
RecursiveGaussianImageFilter (const Self &) | |
Private Attributes | |
bool | m_NormalizeAcrossScale |
OrderEnumType | m_Order |
ScalarRealType | m_Sigma |
typedef SmartPointer< const Self > itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::ConstPointer |
Definition at line 70 of file itkRecursiveGaussianImageFilter.h.
typedef TOutputImage itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::OutputImageType |
Type of the output image
Definition at line 95 of file itkRecursiveGaussianImageFilter.h.
typedef SmartPointer< Self > itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::Pointer |
Definition at line 69 of file itkRecursiveGaussianImageFilter.h.
typedef Superclass::RealType itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::RealType |
Definition at line 72 of file itkRecursiveGaussianImageFilter.h.
typedef Superclass::ScalarRealType itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::ScalarRealType |
Definition at line 73 of file itkRecursiveGaussianImageFilter.h.
typedef RecursiveGaussianImageFilter itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::Self |
Standard class typedefs.
Definition at line 67 of file itkRecursiveGaussianImageFilter.h.
typedef RecursiveSeparableImageFilter< TInputImage, TOutputImage > itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::Superclass |
Definition at line 68 of file itkRecursiveGaussianImageFilter.h.
enum itk::RecursiveGaussianImageFilter::OrderEnumType |
Enum type that indicates if the filter applies the equivalent operation of convolving with a gaussian, first derivative of a gaussian or the second derivative of a gaussian.
Definition at line 92 of file itkRecursiveGaussianImageFilter.h.
|
protected |
|
inlineprotectedvirtual |
Definition at line 160 of file itkRecursiveGaussianImageFilter.h.
|
private |
|
private |
Compute the D coefficients in the recursive filter.
|
private |
Compute the N coefficients in the recursive filter.
|
private |
Compute the M coefficients and the boundary coefficients in the recursive filter.
|
virtual |
Create an object from an instance, potentially deferring to a factory. This method allows you to create an instance of an object that is exactly the same type as the referring object. This is useful in cases where an object has been cast back to a base class.
Reimplemented from itk::Object.
|
virtual |
Type macro that defines a name for this class
Reimplemented from itk::RecursiveSeparableImageFilter< TInputImage, TOutputImage >.
|
virtual |
|
virtual |
Set/Get the Order of the Gaussian to convolve with.
|
virtual |
Set/Get the Sigma, measured in world coordinates, of the Gaussian kernel. The default is 1.0. An exception will be generated if the Sigma value is less than or equal to zero.
|
static |
Method for creation through the object factory.
|
private |
|
protectedvirtual |
Methods invoked by Print() to print information about the object including superclasses. Typically not called by the user (use Print() instead) but used in the hierarchical print process to combine the output of several classes.
Reimplemented from itk::RecursiveSeparableImageFilter< TInputImage, TOutputImage >.
void itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::SetFirstOrder | ( | ) |
Explicitly set a first order derivative.
|
virtual |
Set/Get the flag for normalizing the gaussian over scale-space.
This flag enables the analysis of the differential shape of features independent of their size ( both pixels and physical size ). Following the notation of Tony Lindeberg:
Let
be the scale-space representation of image
where
is the Gaussian function and
denotes convolution. This is a change from above with
.
Then the normalized derivative operator for normalized coordinates across scale is:
The resulting scaling factor is
where N is the order of the derivative.
When this flag is ON the filter will be normalized in such a way that the values of derivatives are not biased by the size of the object. That is to say the maximum value a feature reaches across scale is independent of the scale of the object.
For analyzing an image across scale-space you want to enable this flag. It is disabled by default.
|
virtual |
Set/Get the Order of the Gaussian to convolve with.
void itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::SetSecondOrder | ( | ) |
Explicitly set a second order derivative.
|
virtual |
Set/Get the Sigma, measured in world coordinates, of the Gaussian kernel. The default is 1.0. An exception will be generated if the Sigma value is less than or equal to zero.
|
protectedvirtual |
Set up the coefficients of the filter to approximate a specific kernel. Here it is used to approximate a Gaussian or one of its derivatives. Parameter is the spacing along the dimension to filter.
Implements itk::RecursiveSeparableImageFilter< TInputImage, TOutputImage >.
void itk::RecursiveGaussianImageFilter< TInputImage, TOutputImage >::SetZeroOrder | ( | ) |
Explicitly set a zeroth order derivative.
|
protectedvirtual |
Verifies that the process object has been configured correctly, that all required inputs are set, and needed parameters are set appropriately. If not valid an exceptions will be thrown.
This method is called before UpdateOutputInformation is propagated to the inputs.
The ProcessObject's implementation verifies that the NumberOfRequiredInputs are set and not null.
Reimplemented from itk::ProcessObject.
|
private |
Normalize the image across scale space
Definition at line 198 of file itkRecursiveGaussianImageFilter.h.
|
private |
Definition at line 200 of file itkRecursiveGaussianImageFilter.h.
|
private |
Sigma of the gaussian kernel.
Definition at line 195 of file itkRecursiveGaussianImageFilter.h.