#include <itkMattesMutualInformationImageToImageMetricv4.h>

Inheritance diagram for itk::MattesMutualInformationImageToImageMetricv4< TFixedImage, TMovingImage, TVirtualImage >:

Collaboration diagram for itk::MattesMutualInformationImageToImageMetricv4< TFixedImage, TMovingImage, TVirtualImage >:

Public Types
typedef SmartPointer< const Self >	ConstPointer
typedef Superclass::DerivativeType	DerivativeType
typedef DerivativeType::ValueType	DerivativeValueType
typedef Superclass::FixedImageGradientType	FixedImageGradientType
typedef Superclass::FixedImageIndexType	FixedImageIndexType
typedef Superclass::FixedImagePixelType	FixedImagePixelType
typedef Superclass::FixedImagePointType	FixedImagePointType
typedef Superclass::FixedImageType	FixedImageType
typedef Superclass::FixedSampledPointSetPointer	FixedSampledPointSetPointer
typedef Superclass::JacobianType	JacobianType
typedef Image< PDFValueType, 3 >	JointPDFDerivativesType
typedef Image< PDFValueType, 2 >	JointPDFType
typedef Superclass::MeasureType	MeasureType
typedef Superclass::MovingImageGradientType	MovingImageGradientType
typedef Superclass::MovingImagePixelType	MovingImagePixelType
typedef Superclass::MovingImagePointType	MovingImagePointType
typedef Superclass::MovingTransformType	MovingTransformType
typedef double	PDFValueType
typedef SmartPointer< Self >	Pointer
typedef MattesMutualInformationImageToImageMetricv4	Self
typedef ImageToImageMetricv4 < TFixedImage, TMovingImage, TVirtualImage >	Superclass
typedef Superclass::VirtualImageType	VirtualImageType
typedef Superclass::VirtualIndexType	VirtualIndexType
typedef Superclass::VirtualPointType	VirtualPointType
typedef Superclass::VirtualSampledPointSetType	VirtualSampledPointSetType
Public Member Functions
virtual ::itk::LightObject::Pointer	CreateAnother (void) const
virtual const char *	GetNameOfClass () const
virtual void	GetValueAndDerivative (MeasureType &value, DerivativeType &derivative) const
virtual void	Initialize (void) throw ( itk::ExceptionObject )

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

const JointPDFType::Pointer	GetJointPDF () const

const JointPDFDerivativesType::Pointer	GetJointPDFDerivatives () const
Static Public Member Functions
static Pointer	New ()
Static Public Attributes
static const ImageDimensionType	FixedImageDimension = ::itk::GetImageDimension<TFixedImage>::ImageDimension
static const ImageDimensionType	MovingImageDimension = ::itk::GetImageDimension<TMovingImage>::ImageDimension
static const ImageDimensionType	VirtualImageDimension = ::itk::GetImageDimension<TVirtualImage>::ImageDimension
Protected Types
typedef BSplineDerivativeKernelFunction < 3, PDFValueType >	CubicBSplineDerivativeFunctionType
typedef BSplineKernelFunction < 3, PDFValueType >	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 MattesMutualInformationImageToImageMetricv4GetValueAndDerivativeThreader < ThreadedImageRegionPartitioner < Superclass::VirtualImageDimension > , Superclass, Self >	MattesMutualInformationDenseGetValueAndDerivativeThreaderType
typedef MattesMutualInformationImageToImageMetricv4GetValueAndDerivativeThreader < ThreadedIndexedContainerPartitioner, Superclass, Self >	MattesMutualInformationSparseGetValueAndDerivativeThreaderType
typedef std::vector< PRatioType >	PRatioArrayType
typedef PDFValueType	PRatioType
Protected Member Functions
OffsetValueType	ComputeSingleFixedImageParzenWindowIndex (const FixedImagePixelType &value) const
virtual void	GetValueCommonAfterThreadedExecution ()
	MattesMutualInformationImageToImageMetricv4 ()
void	PrintSelf (std::ostream &os, Indent indent) const
virtual	~MattesMutualInformationImageToImageMetricv4 ()
Protected Attributes
CubicBSplineDerivativeFunctionType::Pointer	m_CubicBSplineDerivativeKernel
CubicBSplineFunctionType::Pointer	m_CubicBSplineKernel
PDFValueType	m_FixedImageBinSize
PDFValueType	m_FixedImageNormalizedMin
PDFValueType	m_FixedImageTrueMax
PDFValueType	m_FixedImageTrueMin
std::vector< OffsetValueType >	m_JointPdfIndex1DArray
std::vector< DerivativeType >	m_LocalDerivativeByParzenBin
PDFValueType	m_MovingImageBinSize
std::vector< PDFValueType >	m_MovingImageMarginalPDF
PDFValueType	m_MovingImageNormalizedMin
PDFValueType	m_MovingImageTrueMax
PDFValueType	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< PDFValueType >	m_ThreaderJointPDFSum
Private Member Functions
	MattesMutualInformationImageToImageMetricv4 (const Self &)
void	operator= (const Self &)
Friends
class	MattesMutualInformationImageToImageMetricv4GetValueAndDerivativeThreader< ThreadedImageRegionPartitioner< Superclass::VirtualImageDimension >, Superclass, Self >
class	MattesMutualInformationImageToImageMetricv4GetValueAndDerivativeThreader< ThreadedIndexedContainerPartitioner, Superclass, Self >

Detailed Description

template<class TFixedImage, class TMovingImage, class TVirtualImage = TFixedImage>
class itk::MattesMutualInformationImageToImageMetricv4< TFixedImage, TMovingImage, TVirtualImage >

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 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.

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.

Warning:: Local-support transforms are not yet supported. If used, an exception is thrown during Initialize().

Note:: The per-iteration post-processing code is not multi-threaded, but could be readily be made so for a small performance gain. See GetValueCommonAfterThreadedExecution(), GetValueAndDerivative() and threader::AfterThreadedExecution().

The algorithm and much of the code was copied from the previous Mattes MI metric, i.e. itkMattesMutualInformationImageToImageMetric.

See MattesMutualInformationImageToImageMetricv4GetValueAndDerivativeThreader::ProcessPoint for poritons of the algorithm implementation.

See ImageToImageMetricv4 for details of common metric operation and options.

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.

See also:: itkImageToImageMetricv4

Definition at line 96 of file itkMattesMutualInformationImageToImageMetricv4.h.