#include <itkMutualInformationImageToImageMetric.h>
Inheritance diagram for itk::MutualInformationImageToImageMetric< TFixedImage, TMovingImage >:


MutualInformationImageToImageMetric computes the mutual information between a fixed and moving image to be registered.
This class is templated over the FixedImage type and the MovingImage type.
The fixed and moving images are set via methods SetFixedImage() and SetMovingImage(). This metric makes use of user specified Transform and Interpolator. The Transform is used to map points from the fixed image to the moving image domain. The Interpolator is used to evaluate the image intensity at user specified geometric points in the moving image. The Transform and Interpolator are set via methods SetTransform() and SetInterpolator().
The calculations are based on the method of Viola and Wells where the probability density distributions are estimated using Parzen windows.
By default a Gaussian kernel is used in the density estimation. Other option include Cauchy and spline-based. A user can specify the kernel passing in a pointer a KernelFunction using the SetKernelFunction() method.
Mutual information is estimated using two sample sets: one to calculate the singular and joint pdf's and one to calculate the entropy integral. By default 50 samples points are used in each set. Other values can be set via the SetNumberOfSpatialSamples() method.
Quality of the density estimate depends on the choice of the kernel's standard deviation. Optimal choice will depend on the images. It is can be shown that around the optimal variance, the mutual information estimate is relatively insensitive to small changes of the standard deviation. In our experiments, we have found that a standard deviation of 0.4 works well for images normalized to have a mean of zero and standard deviation of 1.0. The variance can be set via methods SetFixedImageStandardDeviation() and SetMovingImageStandardDeviation().
Implementaton of this class is based on: Viola, P. and Wells III, W. (1997). "Alignment by Maximization of Mutual Information" International Journal of Computer Vision, 24(2):137-154
Definition at line 90 of file itkMutualInformationImageToImageMetric.h.
Public Types | |
| typedef SmartPointer< const Self > | ConstPointer |
| typedef Superclass::DerivativeType | DerivativeType |
| typedef Superclass::FixedImageConstPointer | FixedImageConstPointer |
| typedef FixedImageType::IndexType | FixedImageIndexType |
| typedef FixedImageIndexType::IndexValueType | FixedImageIndexValueType |
| typedef FixedImageMaskType::Pointer | FixedImageMaskPointer |
| typedef SpatialObject< itkGetStaticConstMacro(FixedImageDimension) | FixedImageMaskType ) |
| typedef TransformType::InputPointType | FixedImagePointType |
| typedef FixedImageType::RegionType | FixedImageRegionType |
| typedef Superclass::FixedImageType | FixedImageType |
| typedef GradientImageFilterType::Pointer | GradientImageFilterPointer |
| typedef GradientRecursiveGaussianImageFilter< MovingImageType, GradientImageType > | GradientImageFilterType |
| typedef SmartPointer< GradientImageType > | GradientImagePointer |
| typedef Image< GradientPixelType, itkGetStaticConstMacro(MovingImageDimension) | GradientImageType ) |
| typedef CovariantVector< RealType, itkGetStaticConstMacro(MovingImageDimension) | GradientPixelType ) |
| typedef TransformType::InputPointType | InputPointType |
| typedef InterpolatorType::Pointer | InterpolatorPointer |
| typedef Superclass::InterpolatorType | InterpolatorType |
| typedef Superclass::MeasureType | MeasureType |
| typedef MovingImageType::ConstPointer | MovingImageConstPointer |
| typedef Superclass::MovingImageConstPointer | MovingImageCosntPointer |
| typedef MovingImageType::IndexType | MovingImageIndexType |
| typedef MovingImageMaskType::Pointer | MovingImageMaskPointer |
| typedef SpatialObject< itkGetStaticConstMacro(MovingImageDimension) | MovingImageMaskType ) |
| typedef TMovingImage::PixelType | MovingImagePixelType |
| typedef TransformType::OutputPointType | MovingImagePointType |
| typedef Superclass::MovingImageType | MovingImageType |
| typedef TransformType::OutputPointType | OutputPointType |
| typedef Superclass::ParametersType | ParametersType |
| typedef double | ParametersValueType |
| typedef SmartPointer< Self > | Pointer |
| typedef MutualInformationImageToImageMetric | Self |
| typedef ImageToImageMetric< TFixedImage, TMovingImage > | Superclass |
| typedef Superclass::TransformJacobianType | TransformJacobianType |
| typedef TransformType::ParametersType | TransformParametersType |
| typedef Superclass::TransformPointer | TransformPointer |
| typedef Superclass::TransformType | TransformType |
Public Member Functions | |
| virtual void | ComputeGradientOff () |
| virtual void | ComputeGradientOn () |
| virtual LightObject::Pointer | CreateAnother () const |
| virtual void | DebugOff () const |
| virtual void | DebugOn () const |
| virtual void | Delete () |
| Command * | GetCommand (unsigned long tag) |
| virtual const bool & | GetComputeGradient () |
| bool | GetDebug () const |
| void | GetDerivative (const ParametersType ¶meters, DerivativeType &Derivative) const |
| virtual const FixedImageType * | GetFixedImage () |
| virtual const FixedImageMaskType * | GetFixedImageMask () |
| virtual const FixedImageRegionType & | GetFixedImageRegion () |
| virtual double | GetFixedImageStandardDeviation () |
| virtual const GradientImageType * | GetGradientImage () |
| virtual const InterpolatorType * | GetInterpolator () |
| virtual KernelFunction * | GetKernelFunction () |
| const MetaDataDictionary & | GetMetaDataDictionary (void) const |
| MetaDataDictionary & | GetMetaDataDictionary (void) |
| virtual const MovingImageType * | GetMovingImage () |
| virtual const MovingImageMaskType * | GetMovingImageMask () |
| virtual const double & | GetMovingImageStandardDeviation () |
| virtual unsigned long | GetMTime () const |
| virtual const char * | GetNameOfClass () const |
| unsigned int | GetNumberOfParameters (void) const |
| virtual const unsigned long & | GetNumberOfPixelsCounted () |
| virtual const unsigned int & | GetNumberOfSpatialSamples () |
| virtual int | GetReferenceCount () const |
| virtual const TransformType * | GetTransform () |
| MeasureType | GetValue (const ParametersType ¶meters) const |
| void | GetValueAndDerivative (const ParametersType ¶meters, MeasureType &Value, DerivativeType &Derivative) const |
| bool | HasObserver (const EventObject &event) const |
| virtual void | Initialize (void) throw ( ExceptionObject ) |
| void | InvokeEvent (const EventObject &) const |
| void | InvokeEvent (const EventObject &) |
| itkStaticConstMacro (FixedImageDimension, unsigned int, TFixedImage::ImageDimension) | |
| Constants for the image dimensions * | itkStaticConstMacro (MovingImageDimension, unsigned int, TMovingImage::ImageDimension) |
| itkStaticConstMacro (MovingImageDimension, unsigned int, MovingImageType::ImageDimension) | |
| virtual void | Modified () const |
| void | Print (std::ostream &os, Indent indent=0) const |
| virtual void | Register () const |
| void | ReinitializeSeed (int) |
| void | ReinitializeSeed () |
| void | RemoveAllObservers () |
| void | RemoveObserver (unsigned long tag) |
| Set Get gradient computation *virtual void | SetComputeGradient (bool _arg) |
| void | SetDebug (bool debugFlag) const |
| virtual void | SetFixedImage (const FixedImageType *_arg) |
| Set Get the fixed image mask *virtual void | SetFixedImageMask (FixedImageMaskType *_arg) |
| virtual void | SetFixedImageRegion (FixedImageRegionType _arg) |
| Set Get the fixed image intensitiy standard deviation This defines *the kernel bandwidth used in the joint probability distribution *calculation Default value is which works well for image intensities *normalized to a mean of and standard deviation of *Value is clamped to be always greater than zero *virtual void | SetFixedImageStandardDeviation (double _arg) |
| virtual void | SetInterpolator (InterpolatorType *_arg) |
| Set Get the kernel function This is used to calculate the joint *probability distribution Default is the GaussianKernelFunction *virtual void | SetKernelFunction (KernelFunction *_arg) |
| void | SetMetaDataDictionary (const MetaDataDictionary &rhs) |
| virtual void | SetMovingImage (const MovingImageType *_arg) |
| Set Get the moving image mask *virtual void | SetMovingImageMask (MovingImageMaskType *_arg) |
| Set Get the moving image intensitiy standard deviation This defines *the kernel bandwidth used in the joint probability distribution *calculation Default value is which works well for image intensities *normalized to a mean of and standard deviation of *Value is clamped to be always greater than zero *virtual void | SetMovingImageStandardDeviation (double _arg) |
| void | SetNumberOfSpatialSamples (unsigned int num) |
| virtual void | SetReferenceCount (int) |
| virtual void | SetTransform (TransformType *_arg) |
| void | SetTransformParameters (const ParametersType ¶meters) const |
| virtual void | UnRegister () const |
Static Public Member Functions | |
| static void | BreakOnError () |
| static bool | GetGlobalWarningDisplay () |
| static void | GlobalWarningDisplayOff () |
| static void | GlobalWarningDisplayOn () |
| static Pointer | New () |
| This is a global flag that controls whether any warning *or error messages are displayed *static void | SetGlobalWarningDisplay (bool flag) |
Public Attributes | |
| Allow people to add remove invoke observers(callbacks) to any ITK *object.This is an implementation of the subject/observer design *pattern.An observer is added by specifying an event to respond to *and an itk unsigned lon | AddObserver )(const EventObject &event, Command *) const |
| This is a global flag that controls whether any | debug |
| Gaussian filter to compute the gradient of the Moving Image *typedef NumericTraits< MovingImagePixelType >::RealType | RealType |
Protected Member Functions | |
| MutualInformationImageToImageMetric () | |
| bool | PrintObservers (std::ostream &os, Indent indent) const |
| void | PrintSelf (std::ostream &os, Indent indent) const |
| virtual void | PrintTrailer (std::ostream &os, Indent indent) const |
| virtual | ~MutualInformationImageToImageMetric () |
Protected Attributes | |
| bool | m_ComputeGradient |
| FixedImageConstPointer | m_FixedImage |
| FixedImageMaskPointer | m_FixedImageMask |
| GradientImagePointer | m_GradientImage |
| InterpolatorPointer | m_Interpolator |
| MovingImageConstPointer | m_MovingImage |
| MovingImageMaskPointer | m_MovingImageMask |
| unsigned long | m_NumberOfPixelsCounted |
| int | m_ReferenceCount |
| SimpleFastMutexLock | m_ReferenceCountLock |
| TransformPointer | m_Transform |
| Methods invoked by virtual 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.*/virtual void PrintSelf(std voi | PrintHeader )(std::ostream &os, Indent indent) const |
Classes | |
| class | SpatialSample |
|
|||||
|
Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 99 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Type of the derivative. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 113 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 117 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Index and Point typedef support. Definition at line 121 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Definition at line 122 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Definition at line 122 of file itkImageToImageMetric.h. |
|
|||||
|
Type for the mask of the fixed image. Only pixels that are "inside" this mask will be considered for the computation of the metric Definition at line 121 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 124 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
|||||
|
Type of the fixed Image. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 115 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Definition at line 110 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 109 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 106 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 105 of file itkImageToImageMetric.h. |
|
|||||
|
|||||
|
Reimplemented in itk::HistogramImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 89 of file itkImageToImageMetric.h. |
|
|||||
|
Reimplemented in itk::CompareHistogramImageToImageMetric< TFixedImage, TMovingImage >, and itk::KullbackLeiblerCompareHistogramImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 114 of file itkImageToImageMetric.h. |
|
|||||
|
Type of the Interpolator Base class Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 111 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Type of the measure. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 112 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
|||||
|
Definition at line 118 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Definition at line 123 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Definition at line 128 of file itkImageToImageMetric.h. |
|
|||||
|
Type for the mask of the moving image. Only pixels that are "inside" this mask will be considered for the computation of the metric Definition at line 127 of file itkImageToImageMetric.h. |
|
|||||
|
Reimplemented in itk::HistogramImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 68 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 125 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Type of the moving Image. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 116 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Reimplemented in itk::HistogramImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 90 of file itkImageToImageMetric.h. |
|
|||||
|
Type of the parameters. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 114 of file itkMutualInformationImageToImageMetric.h. |
|
|
ParametersType typedef. It defines a position in the optimization search space. Definition at line 46 of file itkCostFunction.h. |
|
|||||
|
Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 98 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Standard class typedefs. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 96 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 97 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 110 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
|||||
|
Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 109 of file itkMutualInformationImageToImageMetric.h. |
|
|||||
|
Types inherited from Superclass. Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 105 of file itkMutualInformationImageToImageMetric.h. |
|
|||||||||
|
|
|
|||||||||
|
Definition at line 183 of file itkMutualInformationImageToImageMetric.h. |
|
|
This method is called when itkExceptionMacro executes. It allows the debugger to break on error. |
|
|||||||||
|
|
|
|||||||||
|
|
|
|
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::LightObject. |
|
|
Turn debugging output off. |
|
|
Turn debugging output on. |
|
|
Delete an itk object. This method should always be used to delete an object when the new operator was used to create it. Using the C delete method will not work with reference counting. |
|
|
Get the command associated with the given tag. NOTE: This returns a pointer to a Command, but it is safe to asign this to a Command::Pointer. Since Command inherits from LightObject, at this point in the code, only a pointer or a reference to the Command can be used. |
|
|||||||||
|
|
|
|
Get the value of the debug flag. |
|
||||||||||||||||
|
Get the derivatives of the match measure. Implements itk::SingleValuedCostFunction. |
|
|||||||||
|
Get the Fixed Image. |
|
|||||||||
|
|
|
|||||||||
|
Get the region over which the metric will be computed |
|
|||||||||
|
|
|
|
|
|
|||||||||
|
Get Gradient Image. |
|
|||||||||
|
Get a pointer to the Interpolator. |
|
|||||||||
|
|
|
|
|
|
|
|
|
|||||||||
|
Get the Moving Image. |
|
|||||||||
|
|
|
|||||||||
|
|
|
|
|||||||||
|
Run-time type information (and related methods). Reimplemented from itk::ImageToImageMetric< TFixedImage, TMovingImage >. |
|
||||||||||
|
Return the number of parameters required by the Transform Implements itk::CostFunction. Reimplemented in itk::CompareHistogramImageToImageMetric< TFixedImage, TMovingImage >, and itk::KullbackLeiblerCompareHistogramImageToImageMetric< TFixedImage, TMovingImage >. Definition at line 197 of file itkImageToImageMetric.h. |
|
|||||||||
|
Get the number of pixels considered in the computation. |
|
|||||||||
|
Get the number of spatial samples. |
|
|
Gets the reference count on this object. Definition at line 98 of file itkLightObject.h. |
|
|||||||||
|
Get a pointer to the Transform. |
|
||||||||||
|
Get the value. Implements itk::SingleValuedCostFunction. |
|
||||||||||||||||||||
|
Get the value and derivatives for single valued optimizers. Reimplemented from itk::SingleValuedCostFunction. |
|
|
Definition at line 100 of file itkObject.h. References itk::Object::SetGlobalWarningDisplay(). |
|
|
Definition at line 98 of file itkObject.h. References itk::Object::SetGlobalWarningDisplay(). |
|
|
Return true if an observer is registered for this event. |
|
||||||||||
|
Initialize the Metric by making sure that all the components are present and plugged together correctly Reimplemented in itk::CompareHistogramImageToImageMetric< TFixedImage, TMovingImage >, itk::GradientDifferenceImageToImageMetric< TFixedImage, TMovingImage >, itk::HistogramImageToImageMetric< TFixedImage, TMovingImage >, itk::KullbackLeiblerCompareHistogramImageToImageMetric< TFixedImage, TMovingImage >, and itk::MattesMutualInformationImageToImageMetric< TFixedImage, TMovingImage >. |
|
|
Call Execute on all the Commands observing this event id. The actions triggered by this call doesn't modify this object. |
|
|
Call Execute on all the Commands observing this event id. |
|
||||||||||||||||||||
|
Reimplemented in itk::GradientDifferenceImageToImageMetric< TFixedImage, TMovingImage >. |
|
||||||||||||||||||||
|
|
|
||||||||||||||||||||
|
Enum of the moving image dimension. |
|
|
|||||||||
|
Method for creation through the object factory. Reimplemented from itk::Object. |
|
||||||||||||
|
Cause the object to print itself out. |
|
||||||||||||
|
|
|
||||||||||||||||
|
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::ImageToImageMetric< TFixedImage, TMovingImage >. |
|
||||||||||||
|
|
|
|
Increase the reference count (mark as used by another object). Reimplemented from itk::LightObject. |
|
||||||||||
|
|
|
|||||||||
|
|
|
|
Remove all observers . |
|
|
Remove the observer with this tag value. |
|
||||||||||
|
|
|
|
Set the value of the debug flag. A non-zero value turns debugging on. |
|
||||||||||
|
Connect the Fixed Image. |
|
||||||||||
|
|
|
||||||||||
|
Set the region over which the metric will be computed |
|
||||||||||
|
|
|
|
Referenced by itk::Object::GlobalWarningDisplayOff(), and itk::Object::GlobalWarningDisplayOn(). |
|
||||||||||
|
Connect the Interpolator. |
|
||||||||||
|
|
|
|
|
|
||||||||||
|
Connect the Moving Image. |
|
||||||||||
|
|
|
||||||||||
|
|
|
||||||||||
|
Set the number of spatial samples. This is the number of image samples used to calculate the joint probability distribution. The number of spatial samples is clamped to be a minimum of 1. Default value is 50. |
|
|
Sets the reference count (use with care) Reimplemented from itk::LightObject. |
|
||||||||||
|
Connect the Transform. |
|
||||||||||
|
Set the parameters defining the Transform. |
|
|
Decrease the reference count (release by another object). Reimplemented from itk::LightObject. |
|
|
|
|
|
Definition at line 94 of file itkObject.h. |
|
|||||
|
Definition at line 217 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 211 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 220 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 218 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 215 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 212 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 221 of file itkImageToImageMetric.h. |
|
|||||
|
Definition at line 209 of file itkImageToImageMetric.h. |
|
|
Number of uses of this object by other objects. Definition at line 119 of file itkLightObject.h. |
|
|
Mutex lock to protect modification to the reference count Definition at line 122 of file itkLightObject.h. |
|
|||||
|
Definition at line 214 of file itkImageToImageMetric.h. |
|
|
|
|
|||||
1.4.2 written by Dimitri van Heesch,
© 1997-2000