Initial implementation of Seb Ourselin's block matching algorithm. More...

Inheritance diagram for itk::BlockMatchingMethod< TImageType, TScalarType >:

Collaboration diagram for itk::BlockMatchingMethod< TImageType, TScalarType >:

Classes
class	ResidualHeapDataType
	So we can have an ordered list of points, based on the residual between transformed fixed and moving. More...

class	VarianceHeapDataType
	So we can have an ordered list of indexes, based on variance of a block. More...

Public Types
typedef BlockMatchingMethod	Self

typedef MaskedImageRegistrationMethod < TImageType >	Superclass

typedef SmartPointer< Self >	Pointer

typedef SmartPointer< const Self >	ConstPointer

typedef TImageType	ImageType

typedef ImageType::SizeType	ImageSizeType

typedef ImageType::PixelType	ImagePixelType

typedef ImageType::SpacingType	ImageSpacingType

typedef ImageType::RegionType	ImageRegionType

typedef ImageType::IndexType	ImageIndexType

typedef ResampleImageFilter < TImageType, TImageType >	ResampleFilterType

typedef ResampleFilterType::Pointer	ResampleFilterPointer

typedef RegionOfInterestImageFilter < TImageType, TImageType >	RegionOfInterestFilterType

typedef RegionOfInterestFilterType::Pointer	RegionOfInterestFilterPointer

typedef itk::Statistics::ImageToListSampleAdaptor < ImageType >	ImageTypeListAdaptorType

typedef ImageTypeListAdaptorType::Pointer	ImageTypeListAdaptorPointer

typedef itk::Statistics::CovarianceSampleFilter < ImageTypeListAdaptorType >	ImageTypeCovarianceSampleFilterType

typedef ImageTypeCovarianceSampleFilterType::Pointer	ImageTypeCovarianceSampleFilterPointer

typedef Image< float, TImageType::ImageDimension >	GradientImageType

typedef GradientMagnitudeImageFilter < TImageType, GradientImageType >	GradientMagnitudeFilterType

typedef GradientMagnitudeFilterType::Pointer	GradientMagnitudeFilterPointer

typedef itk::Statistics::ImageToListSampleAdaptor < GradientImageType >	GradientImageTypeListAdaptorType

typedef GradientImageTypeListAdaptorType::Pointer	GradientImageTypeListAdaptorPointer

typedef itk::Statistics::CovarianceSampleFilter < GradientImageTypeListAdaptorType >	GradientImageTypeCovarianceSampleFilterType

typedef GradientImageTypeCovarianceSampleFilterType::Pointer	GradientImageTypeCovarianceSampleFilterPointer

typedef SimilarityMeasure < TImageType, TImageType >	SimilarityMeasureType

typedef SimilarityMeasureType *	SimilarityMeasurePointer

typedef NearestNeighborInterpolateImageFunction < TImageType, TScalarType >	DummyInterpolatorType

typedef DummyInterpolatorType::Pointer	DummyInterpolatorPointer

typedef Superclass::ParametersType	ParametersType

typedef PointSet< TScalarType, TImageType::ImageDimension, DefaultStaticMeshTraits < TScalarType, TImageType::ImageDimension, TImageType::ImageDimension, TScalarType, TScalarType > >	PointSetType

typedef PointSetType::PointType	PointType

typedef PointSetType::PointsContainer::ConstIterator	PointIterator

typedef PointSetType::Pointer	PointSetPointer

typedef PointSetType::PointsContainer	PointsContainerType

typedef PointsContainerType::Pointer	PointsContainerPointer

typedef PointSetToPointSetSingleValuedMetric < PointSetType, PointSetType >	PointSetMetricType

typedef PointSetMetricType::Pointer	PointSetMetricPointer

typedef Superclass::TransformType	TransformType

typedef TransformType *	TransformPointer

typedef MinimumMaximumImageCalculator < TImageType >	MinimumMaximumImageCalculatorType

typedef MinimumMaximumImageCalculatorType::Pointer	MinimumMaximumImageCalculatorPointer

Public Types inherited from itk::MaskedImageRegistrationMethod< TImageType >
typedef MaskedImageRegistrationMethod	Self

typedef SingleResolutionImageRegistrationMethod < TImageType, TImageType >	Superclass

typedef SmartPointer< Self >	Pointer

typedef SmartPointer< const Self >	ConstPointer

typedef TImageType	InputImageType

typedef InputImageType::PixelType	InputImagePixelType

typedef InputImageType::Pointer	InputImagePointer

typedef InputImageType::ConstPointer	InputImageConstPointer

typedef InputImageType::RegionType	InputImageRegionType

typedef InputImageType::SizeType	InputImageSizeType

typedef InputImageType::IndexType	InputImageIndexType

typedef InputImageType::SpacingType	InputImageSpacingType

typedef InputImageType::PointType	InputImageOriginType

typedef InputImageType::DirectionType	InputImageDirectionType

typedef BoundaryValueRescaleIntensityImageFilter < InputImageType >	RescaleFilterType

typedef RescaleFilterType::Pointer	RescaleFilterPointer

typedef SmoothingRecursiveGaussianImageFilter < InputImageType, InputImageType >	SmoothingFilterType

typedef SmoothingFilterType::Pointer	SmoothingFilterPointer

typedef BinaryThresholdImageFilter < InputImageType, InputImageType >	ThresholdFilterType

typedef ThresholdFilterType::Pointer	ThresholdFilterPointer

typedef BinaryCrossStructuringElement < InputImagePixelType, InputImageType::ImageDimension >	StructuringType

typedef BinaryDilateImageFilter < InputImageType, InputImageType, StructuringType >	DilateMaskFilterType

typedef DilateMaskFilterType::Pointer	DilateMaskFilterPointer

typedef MultiplyImageFilter < InputImageType, InputImageType >	MultiplyFilterType

typedef MultiplyFilterType::Pointer	MultiplyFilterPointer

typedef unsigned char	MaskPixelType

typedef Image< MaskPixelType, InputImageType::ImageDimension >	MaskImageType

typedef CastImageFilter < InputImageType, MaskImageType >	CastToMaskImageTypeFilterType

typedef CastToMaskImageTypeFilterType::Pointer	CastToMaskImageTypeFilterPointer

typedef ImageMaskSpatialObject < InputImageType::ImageDimension >	MaskFilterType

typedef MaskFilterType::Pointer	MaskFilterPointer

typedef Superclass::TransformType	TransformType

Public Types inherited from itk::SingleResolutionImageRegistrationMethod< TImageType, TImageType >
typedef SingleResolutionImageRegistrationMethod	Self

typedef ImageRegistrationMethod < TImageType, TImageType >	Superclass

typedef SmartPointer< Self >	Pointer

typedef SmartPointer< const Self >	ConstPointer

typedef itk::IterationUpdateCommand	IterationUpdateCommandType

typedef IterationUpdateCommandType::Pointer	IterationUpdateCommandPointer

typedef Superclass::ParametersType	ParametersType

typedef Superclass::TransformType	TransformType

typedef Superclass::InterpolatorType	InterpolatorType

Public Member Functions
virtual const char *	GetClassName () const

	itkStaticConstMacro (Dimension, unsigned int, TImageType::ImageDimension)

virtual void	SetPointSetMetric (PointSetMetricType *_arg)

virtual PointSetMetricType *	GetPointSetMetric ()

virtual void	SetMaximumNumberOfIterationsRoundMainLoop (unsigned int _arg)

virtual unsigned int	GetMaximumNumberOfIterationsRoundMainLoop ()

void	SetBlockParameters (double blockSize, double blockHalfWidth, double blockSpacing, double blockSubSampling)

virtual double	GetBlockSize ()

virtual double	GetBlockHalfWidth ()

virtual double	GetBlockSpacing ()

virtual double	GetBlockSubSampling ()

virtual void	SetTransformedMovingImageFileName (std::string _arg)

virtual std::string	GetTransformedMovingImageFileName ()

virtual void	SetTransformedMovingImageFileExt (std::string _arg)

virtual std::string	GetTransformedMovingImageFileExt ()

virtual void	SetWriteTransformedMovingImage (bool _arg)

virtual bool	GetWriteTransformedMovingImage ()

virtual void	SetEpsilon (double _arg)

virtual double	GetEpsilon ()

virtual void	SetParameterReductionFactor (double _arg)

virtual double	GetParameterReductionFactor ()

virtual void	SetMinimumBlockSize (double _arg)

virtual double	GetMinimumBlockSize ()

virtual void	SetPercentageOfPointsToKeep (int _arg)

virtual int	GetPercentageOfPointsToKeep ()

virtual void	SetPercentageOfPointsInLeastTrimmedSquares (int _arg)

virtual int	GetPercentageOfPointsInLeastTrimmedSquares ()

virtual void	SetUseGradientMagnitudeVariance (bool _arg)

virtual bool	GetUseGradientMagnitudeVariance ()

virtual void	SetScaleByMillimetres (bool _arg)

virtual bool	GetScaleByMillimetres ()

virtual void	SetPointSetFileNameWithoutExtension (std::string _arg)

virtual std::string	GetPointSetFileNameWithoutExtension ()

virtual void	SetWritePointSet (bool _arg)

virtual bool	GetWritePointSet ()

virtual void	SetNoZero (bool _arg)

virtual bool	GetNoZero ()

virtual void	SetTransformedMovingImagePadValue (ImagePixelType _arg)

virtual ImagePixelType	GetTransformedMovingImagePadValue ()

Public Member Functions inherited from itk::MaskedImageRegistrationMethod< TImageType >
void	SetFixedMask (const InputImageType *fixedMask)

virtual const InputImageType *	GetFixedMask ()

void	SetMovingMask (const InputImageType *movingMask)

virtual const InputImageType *	GetMovingMask ()

virtual void	SetUseFixedMask (bool _arg)

virtual bool	GetUseFixedMask () const

virtual void	SetUseMovingMask (bool _arg)

virtual bool	GetUseMovingMask () const

virtual void	SetNumberOfDilations (unsigned int _arg)

virtual unsigned int	GetNumberOfDilations ()

virtual void	SetMaskImageDirectly (bool _arg)

virtual bool	GetMaskImageDirectly ()

virtual void	SetRescaleFixedImage (bool _arg)

virtual bool	GetRescaleFixedImage () const

virtual void	SetRescaleFixedMinimum (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleFixedMinimum () const

virtual void	SetRescaleFixedMaximum (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleFixedMaximum () const

virtual void	SetRescaleFixedBoundaryValue (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleFixedBoundaryValue () const

virtual void	SetRescaleFixedLowerThreshold (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleFixedLowerThreshold () const

virtual void	SetRescaleFixedUpperThreshold (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleFixedUpperThreshold () const

virtual void	SetRescaleMovingImage (bool _arg)

virtual bool	GetRescaleMovingImage () const

virtual void	SetRescaleMovingMinimum (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleMovingMinimum () const

virtual void	SetRescaleMovingMaximum (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleMovingMaximum () const

virtual void	SetRescaleMovingBoundaryValue (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleMovingBoundaryValue () const

virtual void	SetRescaleMovingLowerThreshold (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleMovingLowerThreshold () const

virtual void	SetRescaleMovingUpperThreshold (InputImagePixelType _arg)

virtual InputImagePixelType	GetRescaleMovingUpperThreshold () const

virtual void	SetThresholdFixedMask (bool _arg)

virtual bool	GetThresholdFixedMask () const

virtual void	SetFixedMaskMinimum (InputImagePixelType _arg)

virtual InputImagePixelType	GetFixedMaskMinimum () const

virtual void	SetFixedMaskMaximum (InputImagePixelType _arg)

virtual InputImagePixelType	GetFixedMaskMaximum () const

virtual void	SetThresholdMovingMask (bool _arg)

virtual bool	GetThresholdMovingMask () const

virtual void	SetMovingMaskMinimum (InputImagePixelType _arg)

virtual InputImagePixelType	GetMovingMaskMinimum () const

virtual void	SetMovingMaskMaximum (InputImagePixelType _arg)

virtual InputImagePixelType	GetMovingMaskMaximum () const

virtual void	SetSigma (float _arg)

virtual float	GetSigma ()

InputImageType *	GetFixedImageCopy ()

InputImageType *	GetMovingImageCopy ()

Public Member Functions inherited from itk::SingleResolutionImageRegistrationMethod< TImageType, TImageType >
virtual void	SetIterationUpdateCommand (IterationUpdateCommandType *_arg)

virtual IterationUpdateCommandType *	GetIterationUpdateCommand ()

virtual void	SetFixedImageInterpolator (InterpolatorType *_arg)

virtual void	SetMovingImageInterpolator (InterpolatorType *_arg)

virtual InterpolatorType *	GetFixedImageInterpolator ()

virtual InterpolatorType *	GetMovingImageInterpolator ()

Static Public Member Functions
static Pointer	New ()

Static Public Member Functions inherited from itk::MaskedImageRegistrationMethod< TImageType >
static Pointer	New ()

Static Public Member Functions inherited from itk::SingleResolutionImageRegistrationMethod< TImageType, TImageType >
static Pointer	New ()

Protected Types
typedef std::priority_queue < VarianceHeapDataType >	VarianceHeap

typedef std::priority_queue < ResidualHeapDataType >	ResidualHeap

Protected Member Functions
	BlockMatchingMethod ()

virtual	~BlockMatchingMethod ()

void	PrintSelf (std::ostream &os, Indent indent) const

virtual void	Initialize () throw (ExceptionObject)

virtual void	DoRegistration () throw (ExceptionObject)

virtual void	WritePointSet (const PointsContainerPointer &fixedPointContainer, const PointsContainerPointer &movingPointContainer)

virtual void	GetPointCorrespondencies2D (ImageSizeType &size, ImageSizeType &bigN, ImageSizeType &bigOmega, ImageSizeType &bigDeltaOne, ImageSizeType &bigDeltaTwo, PointsContainerPointer &fixedPointContainer, PointsContainerPointer &movingPointContainer)

virtual void	GetPointCorrespondencies3D (ImageSizeType &size, ImageSizeType &bigN, ImageSizeType &bigOmega, ImageSizeType &bigDeltaOne, ImageSizeType &bigDeltaTwo, PointsContainerPointer &fixedPointContainer, PointsContainerPointer &movingPointContainer)

virtual bool	CheckEpsilon (ParametersType &previousParameters, ParametersType &currentParameters)

virtual double	CheckSinglePoint (ParametersType &previousParameters, ParametersType &currentParameters, ImageIndexType &index)

virtual void	TrimPoints (const TransformType transform, const PointsContainerType fixedPoints, const PointsContainerType movingPoints, PointsContainerType trimmedFixedPoints, PointsContainerType *trimmedMovingPoints)

Protected Member Functions inherited from itk::MaskedImageRegistrationMethod< TImageType >
	MaskedImageRegistrationMethod ()

virtual	~MaskedImageRegistrationMethod ()

void	PrintSelf (std::ostream &os, Indent indent) const

void	CopyImage (const InputImageType source, InputImageType target)

Protected Member Functions inherited from itk::SingleResolutionImageRegistrationMethod< TImageType, TImageType >
	SingleResolutionImageRegistrationMethod ()

virtual	~SingleResolutionImageRegistrationMethod ()

void	PrintSelf (std::ostream &os, Indent indent) const

void	GenerateData ()

Additional Inherited Members
Protected Attributes inherited from itk::SingleResolutionImageRegistrationMethod< TImageType, TImageType >
InterpolatorType::Pointer	m_FixedImageInterpolator

InterpolatorType::Pointer	m_MovingImageInterpolator

Detailed Description

template<typename TImageType, class TScalarType>
class itk::BlockMatchingMethod< TImageType, TScalarType >

Initial implementation of Seb Ourselin's block matching algorithm.

This class implements the main algorithm of Ourselin et. al. Image and Vision Computing 19 (2000) 25-31. It's a subclass of the itkSingleResolutionImageRegistrationMethod so that we can easily plug it into our usual multi-resolution stuff in itkMultiResolutionImageRegistrationWrapper, and hence put it into itkImageRegistrationFilter.

This implementation uses the same components as the other methods in the RegistrationToolbox, but it doesn't use the same plumbing at all. This class overrides a lot of the base class functionality, so you are advised to look at the implementation.

The algorithm proceeds by:

Take initial transform parameters, and resample the moving image.
Subdivide the fixed image into blocks, and find the corresponding block in the resampled moving image that is most similar.
This results in a list of pairs of point correspondences between the centre of the fixed image block, and the centre of the moving image block.
Take these point correspondences, and perform a trimmed least squares fit, optimising your choice of transformation.
Repeat from 1.

For further details read Ourselin et. al. Image and Vision Computing 19 (2000) 25-31.

See also: MultiResolutionImageRegistrationWrapper; ImageRegistrationFilter