Computes the mutual information between two images to be registered using the method of Mattes et al. More...

#include <itkMattesMutualInformationImageToImageMetric.h>

Inheritance diagram for itk::MattesMutualInformationImageToImageMetric< TFixedImage, TMovingImage >:

Collaboration diagram for itk::MattesMutualInformationImageToImageMetric< TFixedImage, TMovingImage >:

Public Types
typedef Superclass::BSplineTransformIndexArrayType	BSplineTransformIndexArrayType
typedef Superclass::BSplineTransformWeightsType	BSplineTransformWeightsType
typedef SmartPointer< const Self >	ConstPointer
typedef Superclass::CoordinateRepresentationType	CoordinateRepresentationType
typedef Superclass::DerivativeType	DerivativeType
typedef Superclass::FixedImageConstPointer	FixedImageConstPointer
typedef Superclass::FixedImageSampleContainer	FixedImageSampleContainer
typedef Superclass::FixedImageType	FixedImageType
typedef Superclass::ImageDerivativesType	ImageDerivativesType
typedef Superclass::IndexValueType	IndexValueType
typedef Superclass::InterpolatorType	InterpolatorType
typedef Image< PDFValueType, 3 >	JointPDFDerivativesType
typedef Image< PDFValueType, 2 >	JointPDFType
typedef Superclass::MeasureType	MeasureType
typedef Superclass::MovingImageConstPointer	MovingImageConstPointer
typedef Superclass::MovingImagePointType	MovingImagePointType
typedef Superclass::MovingImageType	MovingImageType
typedef FixedImageType::OffsetValueType	OffsetValueType
typedef Superclass::ParametersType	ParametersType
typedef float	PDFValueType
typedef SmartPointer< Self >	Pointer
typedef MattesMutualInformationImageToImageMetric	Self
typedef ImageToImageMetric < TFixedImage, TMovingImage >	Superclass
typedef Superclass::TransformJacobianType	TransformJacobianType
typedef Superclass::TransformPointer	TransformPointer
typedef Superclass::TransformType	TransformType
typedef Superclass::WeightsValueType	WeightsValueType
Public Member Functions
virtual ::itk::LightObject::Pointer	CreateAnother (void) const
void	GetDerivative (const ParametersType &parameters, DerivativeType &Derivative) const
virtual const char *	GetNameOfClass () const
MeasureType	GetValue (const ParametersType &parameters) const
void	GetValueAndDerivative (const ParametersType &parameters, MeasureType &Value, DerivativeType &Derivative) const
virtual void	Initialize (void) throw ( ExceptionObject )

virtual void	SetNumberOfHistogramBins (SizeValueType _arg)
virtual const SizeValueType &	GetNumberOfHistogramBins ()

virtual void	SetUseExplicitPDFDerivatives (bool _arg)
virtual const bool &	GetUseExplicitPDFDerivatives ()
virtual void	UseExplicitPDFDerivativesOn ()
virtual void	UseExplicitPDFDerivativesOff ()

const JointPDFType::Pointer	GetJointPDF () const

const JointPDFDerivativesType::Pointer	GetJointPDFDerivatives () const
Static Public Member Functions
static Pointer	New ()
Static Public Attributes
static const unsigned int	MovingImageDimension = MovingImageType::ImageDimension
Protected Member Functions
	MattesMutualInformationImageToImageMetric ()
void	PrintSelf (std::ostream &os, Indent indent) const
virtual	~MattesMutualInformationImageToImageMetric ()
Private Types
typedef BSplineDerivativeKernelFunction< 3 >	CubicBSplineDerivativeFunctionType
typedef BSplineKernelFunction< 3 >	CubicBSplineFunctionType
typedef JointPDFDerivativesType::IndexType	JointPDFDerivativesIndexType
typedef JointPDFDerivativesType::RegionType	JointPDFDerivativesRegionType
typedef JointPDFDerivativesType::SizeType	JointPDFDerivativesSizeType
typedef JointPDFDerivativesType::PixelType	JointPDFDerivativesValueType
typedef JointPDFType::IndexType	JointPDFIndexType
typedef JointPDFType::RegionType	JointPDFRegionType
typedef JointPDFType::SizeType	JointPDFSizeType
typedef JointPDFType::PixelType	JointPDFValueType
typedef Array2D< PRatioType >	PRatioArrayType
typedef double	PRatioType
Private Member Functions
void	ComputeFixedImageParzenWindowIndices (FixedImageSampleContainer &samples)
void	ComputePDFDerivatives (ThreadIdType threadID, unsigned int sampleNumber, int movingImageParzenWindowIndex, const ImageDerivativesType &movingImageGradientValue, double cubicBSplineDerivativeValue) const
void	GetValueAndDerivativeThreadPostProcess (ThreadIdType threadID, bool withinSampleThread) const
void	GetValueAndDerivativeThreadPreProcess (ThreadIdType threadID, bool withinSampleThread) const
bool	GetValueAndDerivativeThreadProcessSample (ThreadIdType threadID, SizeValueType fixedImageSample, const MovingImagePointType &mappedPoint, double movingImageValue, const ImageDerivativesType &movingImageGradientValue) const
void	GetValueThreadPostProcess (ThreadIdType threadID, bool withinSampleThread) const
void	GetValueThreadPreProcess (ThreadIdType threadID, bool withinSampleThread) const
bool	GetValueThreadProcessSample (ThreadIdType threadID, SizeValueType fixedImageSample, const MovingImagePointType &mappedPoint, double movingImageValue) const
	MattesMutualInformationImageToImageMetric (const Self &)
void	operator= (const Self &)
Private Attributes
CubicBSplineDerivativeFunctionType::Pointer	m_CubicBSplineDerivativeKernel
CubicBSplineFunctionType::Pointer	m_CubicBSplineKernel
double	m_FixedImageBinSize
double	m_FixedImageNormalizedMin
double	m_FixedImageTrueMax
double	m_FixedImageTrueMin
bool	m_ImplicitDerivativesSecondPass
double	m_MovingImageBinSize
std::vector< PDFValueType >	m_MovingImageMarginalPDF
double	m_MovingImageNormalizedMin
double	m_MovingImageTrueMax
double	m_MovingImageTrueMin
SizeValueType	m_NumberOfHistogramBins
PRatioArrayType	m_PRatioArray
std::vector< std::vector < PDFValueType > >	m_ThreaderFixedImageMarginalPDF
std::vector < JointPDFType::Pointer >	m_ThreaderJointPDF
std::vector < JointPDFDerivativesType::Pointer >	m_ThreaderJointPDFDerivatives
std::vector< int >	m_ThreaderJointPDFEndBin
std::vector< int >	m_ThreaderJointPDFStartBin
std::vector< double >	m_ThreaderJointPDFSum
std::vector< DerivativeType >	m_ThreaderMetricDerivative
bool	m_UseExplicitPDFDerivatives

Detailed Description

template<class TFixedImage, class TMovingImage>
class itk::MattesMutualInformationImageToImageMetric< TFixedImage, TMovingImage >

Computes the mutual information between two images to be registered using the method of Mattes et al.

MattesMutualInformationImageToImageMetric 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().

If a BSplineInterpolationFunction is used, this class obtain image derivatives from the BSpline interpolator. Otherwise, image derivatives are computed using central differencing.

Warning:: This metric assumes that the moving image has already been connected to the interpolator outside of this class.

The method GetValue() computes of the mutual information while method GetValueAndDerivative() computes both the mutual information and its derivatives with respect to the transform parameters.

The calculations are based on the method of Mattes et al [1,2] where the probability density distribution are estimated using Parzen histograms. Since the fixed image PDF does not contribute to the derivatives, it does not need to be smooth. Hence, a zero order (box car) BSpline kernel is used for the fixed image intensity PDF. On the other hand, to ensure smoothness a third order BSpline kernel is used for the moving image intensity PDF.

On Initialize(), the FixedImage is uniformly sampled within the FixedImageRegion. The number of samples used can be set via SetNumberOfSpatialSamples(). Typically, the number of spatial samples used should increase with the image size.

The option UseAllPixelOn() disables the random sampling and uses all the pixels of the FixedImageRegion in order to estimate the joint intensity PDF.

During each call of GetValue(), GetDerivatives(), GetValueAndDerivatives(), marginal and joint intensity PDF's values are estimated at discrete position or bins. The number of bins used can be set via SetNumberOfHistogramBins(). To handle data with arbitray magnitude and dynamic range, the image intensity is scale such that any contribution to the histogram will fall into a valid bin.

One the PDF's have been contructed, the mutual information is obtained by doubling summing over the discrete PDF values.

Notes: 1. This class returns the negative mutual information value. 2. This class in not thread safe due the private data structures used to the store the sampled points and the marginal and joint pdfs.

References: [1] "Nonrigid multimodality image registration" D. Mattes, D. R. Haynor, H. Vesselle, T. Lewellen and W. Eubank Medical Imaging 2001: Image Processing, 2001, pp. 1609-1620. [2] "PET-CT Image Registration in the Chest Using Free-form Deformations" D. Mattes, D. R. Haynor, H. Vesselle, T. Lewellen and W. Eubank IEEE Transactions in Medical Imaging. Vol.22, No.1, January 2003. pp.120-128. [3] "Optimization of Mutual Information for MultiResolution Image Registration" P. Thevenaz and M. Unser IEEE Transactions in Image Processing, 9(12) December 2000.

Definition at line 111 of file itkMattesMutualInformationImageToImageMetric.h.