template<class TFixedImage, class TMovingImage>
class itk::GradientDifferenceImageToImageMetric< TFixedImage, TMovingImage >
Computes similarity between two objects to be registered.
This Class is templated over the type of the Images to be compared and over the type of transformation and Iterpolator to be used.
This metric computes the sum of squared differences between pixels in the derivatives of the moving and fixed images after passing the squared difference through a function of type .
Spatial correspondance between both images is established through a Transform. Pixel values are taken from the Moving image. Their positions are mapped to the Fixed image and result in general in non-grid position on it. Values at these non-grid position of the Fixed image are interpolated using a user-selected Interpolator.
Implementation of this class is based on: Hipwell, J. H., et. al. (2003), "Intensity-Based 2-D-3D Registration of
Cerebral Angiograms,", IEEE Transactions on Medical Imaging, 22(11):1417-1426.
Definition at line 74 of file itkGradientDifferenceImageToImageMetric2.h .
using CastFixedImageFilterPointer = typename CastFixedImageFilterType::Pointer
using CastFixedImageFilterType = itk::CastImageFilter< FixedImageType, FixedGradientImageType >
using CastMovedImageFilterPointer = typename CastMovedImageFilterType::Pointer
using CastMovedImageFilterType = itk::CastImageFilter< TransformedMovingImageType , MovedGradientImageType >
using CombinationTransformPointer = typename CombinationTransformType::Pointer
using CombinationTransformType = typename itk::AdvancedCombinationTransform < ScalarType , FixedImageDimension >
using ConstPointer = SmartPointer< const Self >
using FixedGradientImageType = itk::Image< RealType, Self::FixedImageDimension >
using FixedGradientPixelType = typename FixedGradientImageType::PixelType
using FixedImagePixelType = typename TFixedImage::PixelType
using InterpolatorPointer = typename InterpolatorType::Pointer
using MovedGradientImageType = itk::Image< RealType, Self::MovedImageDimension >
using MovedGradientPixelType = typename MovedGradientImageType::PixelType
using MovedImagePixelType = typename TMovingImage::PixelType
using MovingImageRegionType = typename MovingImageType::RegionType
using Pointer = SmartPointer< Self >
using RayCastInterpolatorPointer = typename RayCastInterpolatorType::Pointer
using RayCastInterpolatorType = typename itk::AdvancedRayCastInterpolateImageFunction < MovingImageType, ScalarType >
using ScalarType = typename TransformType::ScalarType
using ScalesType = typename Optimizer::ScalesType
using Self = GradientDifferenceImageToImageMetric
using Superclass = AdvancedImageToImageMetric < TFixedImage, TMovingImage >
using TransformedMovingImageType = itk::Image< FixedImagePixelType , Self::FixedImageDimension >
using TransformMovingImageFilterType = itk::ResampleImageFilter< MovingImageType, TransformedMovingImageType >
using AdvancedTransformType = AdvancedTransform < ScalarType , FixedImageDimension, MovingImageDimension >
using BSplineOrder1TransformPointer = typename BSplineOrder1TransformType::Pointer
using BSplineOrder1TransformType = AdvancedBSplineDeformableTransform < ScalarType , FixedImageDimension, 1 >
using BSplineOrder2TransformPointer = typename BSplineOrder2TransformType::Pointer
using BSplineOrder2TransformType = AdvancedBSplineDeformableTransform < ScalarType , FixedImageDimension, 2 >
using BSplineOrder3TransformPointer = typename BSplineOrder3TransformType::Pointer
using BSplineOrder3TransformType = AdvancedBSplineDeformableTransform < ScalarType , FixedImageDimension, 3 >
using CombinationTransformType = AdvancedCombinationTransform < ScalarType , FixedImageDimension >
using ConstPointer = SmartPointer< const Self >
using DerivativeValueType = typename DerivativeType::ValueType
using FixedImageLimiterOutputType = typename FixedImageLimiterType::OutputType
using FixedImageLimiterPointer = typename FixedImageLimiterType::Pointer
using FixedImageLimiterType = LimiterFunctionBase < RealType, FixedImageDimension >
using FixedImageMaskSpatialObject2Type = ImageMaskSpatialObject< Self::FixedImageDimension >
using FixedImagePixelType = typename FixedImageType::PixelType
using FixedImagePointer = typename FixedImageType::Pointer
using HessianType = vnl_sparse_matrix< HessianValueType >
using HessianValueType = typename DerivativeType::ValueType
using ImageSampleContainerPointer = typename ImageSamplerType::OutputVectorContainerPointer
using ImageSampleContainerType = typename ImageSamplerType::OutputVectorContainerType
using ImageSamplerPointer = typename ImageSamplerType::Pointer
using ImageSamplerType = ImageSamplerBase < FixedImageType >
using MovingImageDerivativeScalesType = FixedArray< double , Self::MovingImageDimension >
using MovingImageLimiterOutputType = typename MovingImageLimiterType::OutputType
using MovingImageLimiterPointer = typename MovingImageLimiterType::Pointer
using MovingImageLimiterType = LimiterFunctionBase < RealType, MovingImageDimension >
using MovingImageMaskSpatialObject2Type = ImageMaskSpatialObject< Self::MovingImageDimension >
using MovingImagePointer = typename MovingImageType::Pointer
using MovingImageRegionType = typename MovingImageType::RegionType
using NumberOfParametersType = typename AdvancedTransformType::NumberOfParametersType
using Pointer = SmartPointer< Self >
using ScalarType = typename TransformType::ScalarType
using Self = AdvancedImageToImageMetric
using Superclass = ImageToImageMetric< TFixedImage, TMovingImage >
using ThreaderType = itk::PlatformMultiThreader
using ThreadInfoType = typename ThreaderType::WorkUnitInfo
virtual const char * GetClassName () const
void GetDerivative (const TransformParametersType ¶meters, DerivativeType &derivative) const override
virtual const double & GetDerivativeDelta ()
virtual const ScalesType & GetScales ()
MeasureType GetValue (const TransformParametersType ¶meters) const override
void GetValueAndDerivative (const TransformParametersType ¶meters, MeasureType &Value, DerivativeType &derivative) const override
void Initialize () override
ITK_DISALLOW_COPY_AND_MOVE (GradientDifferenceImageToImageMetric )
itkStaticConstMacro (FixedImageDimension, unsigned int , FixedImageType::ImageDimension)
itkStaticConstMacro (MovedImageDimension, unsigned int , MovingImageType::ImageDimension)
virtual void SetDerivativeDelta (double _arg)
virtual void SetScales (ScalesType _arg)
void WriteGradientImagesToFiles () const
virtual void BeforeThreadedGetValueAndDerivative (const TransformParametersType ¶meters) const
virtual const char * GetClassName () const
virtual const FixedImageLimiterType * GetFixedImageLimiter ()
virtual double GetFixedLimitRangeRatio () const
virtual ImageSamplerType * GetImageSampler () const
virtual const MovingImageDerivativeScalesType & GetMovingImageDerivativeScales ()
virtual const MovingImageLimiterType * GetMovingImageLimiter ()
virtual double GetMovingLimitRangeRatio () const
virtual double GetRequiredRatioOfValidSamples () const
virtual bool GetScaleGradientWithRespectToMovingImageOrientation () const
virtual void GetSelfHessian (const TransformParametersType ¶meters, HessianType &H) const
const AdvancedTransformType * GetTransform () const override
virtual bool GetUseFixedImageLimiter () const
virtual bool GetUseImageSampler () const
virtual const bool & GetUseMetricSingleThreaded ()
virtual bool GetUseMovingImageDerivativeScales () const
virtual bool GetUseMovingImageLimiter () const
virtual const bool & GetUseMultiThread ()
void Initialize () override
ITK_DISALLOW_COPY_AND_MOVE (AdvancedImageToImageMetric )
itkStaticConstMacro (FixedImageDimension, unsigned int , TFixedImage::ImageDimension)
itkStaticConstMacro (MovingImageDimension, unsigned int , TMovingImage::ImageDimension)
virtual void SetFixedImageLimiter (FixedImageLimiterType *_arg)
virtual void SetFixedLimitRangeRatio (double _arg)
virtual void SetImageSampler (ImageSamplerType *_arg)
virtual void SetMovingImageDerivativeScales (MovingImageDerivativeScalesType _arg)
virtual void SetMovingImageLimiter (MovingImageLimiterType *_arg)
virtual void SetMovingLimitRangeRatio (double _arg)
virtual void SetNumberOfWorkUnits (ThreadIdType numberOfThreads)
virtual void SetRequiredRatioOfValidSamples (double _arg)
virtual void SetScaleGradientWithRespectToMovingImageOrientation (bool _arg)
virtual void SetTransform (AdvancedTransformType *arg)
virtual void SetUseMetricSingleThreaded (bool _arg)
virtual void SetUseMovingImageDerivativeScales (bool _arg)
virtual void SetUseMultiThread (bool _arg)
virtual void UseMetricSingleThreadedOff ()
virtual void UseMetricSingleThreadedOn ()
virtual void UseMultiThreadOff ()
virtual void UseMultiThreadOn ()
using FixedSobelFilter = NeighborhoodOperatorImageFilter< FixedGradientImageType , FixedGradientImageType >
using MovedSobelFilter = NeighborhoodOperatorImageFilter< MovedGradientImageType , MovedGradientImageType >
using BSplineInterpolatorFloatPointer = typename BSplineInterpolatorFloatType::Pointer
using BSplineInterpolatorFloatType = BSplineInterpolateImageFunction< MovingImageType, CoordinateRepresentationType, float >
using BSplineInterpolatorPointer = typename BSplineInterpolatorType::Pointer
using BSplineInterpolatorType = BSplineInterpolateImageFunction< MovingImageType, CoordinateRepresentationType, double >
using CentralDifferenceGradientFilterPointer = typename CentralDifferenceGradientFilterType::Pointer
using CentralDifferenceGradientFilterType = GradientImageFilter< MovingImageType, RealType, RealType >
using FixedImageIndexType = typename FixedImageType::IndexType
using FixedImageIndexValueType = typename FixedImageIndexType::IndexValueType
using FixedImagePointType = typename TransformType::InputPointType
using LinearInterpolatorPointer = typename LinearInterpolatorType::Pointer
using LinearInterpolatorType = AdvancedLinearInterpolateImageFunction < MovingImageType, CoordinateRepresentationType >
using MovingImageContinuousIndexType = typename InterpolatorType::ContinuousIndexType
using MovingImageDerivativeType = typename BSplineInterpolatorType::CovariantVectorType
using MovingImageIndexType = typename MovingImageType::IndexType
using MovingImagePointType = typename TransformType::OutputPointType
using NonZeroJacobianIndicesType = typename AdvancedTransformType::NonZeroJacobianIndicesType
using ReducedBSplineInterpolatorPointer = typename ReducedBSplineInterpolatorType::Pointer
using ReducedBSplineInterpolatorType = ReducedDimensionBSplineInterpolateImageFunction < MovingImageType, CoordinateRepresentationType, double >
MeasureType ComputeMeasure (const TransformParametersType ¶meters, const double *subtractionFactor) const
void ComputeMovedGradientRange () const
void ComputeVariance () const
GradientDifferenceImageToImageMetric ()
void PrintSelf (std::ostream &os, Indent indent) const override
~GradientDifferenceImageToImageMetric () override=default
AdvancedImageToImageMetric ()
virtual void AfterThreadedGetValue (MeasureType &value) const
virtual void AfterThreadedGetValueAndDerivative (MeasureType &value, DerivativeType &derivative) const
virtual void CheckForAdvancedTransform ()
virtual void CheckForBSplineInterpolator ()
virtual void CheckForBSplineTransform () const
virtual void CheckNumberOfSamples (unsigned long wanted, unsigned long found) const
virtual bool EvaluateMovingImageValueAndDerivative (const MovingImagePointType &mappedPoint, RealType &movingImageValue, MovingImageDerivativeType *gradient) const
virtual bool EvaluateTransformJacobian (const FixedImagePointType &fixedImagePoint, TransformJacobianType &jacobian, NonZeroJacobianIndicesType &nzji) const
virtual void EvaluateTransformJacobianInnerProduct (const TransformJacobianType &jacobian, const MovingImageDerivativeType &movingImageDerivative, DerivativeType &imageJacobian) const
bool FastEvaluateMovingImageValueAndDerivative (const MovingImagePointType &mappedPoint, RealType &movingImageValue, MovingImageDerivativeType *gradient, const ThreadIdType threadId) const
virtual void InitializeImageSampler ()
virtual void InitializeLimiters ()
virtual void InitializeThreadingParameters () const
virtual bool IsInsideMovingMask (const MovingImagePointType &point) const
itkAlignedTypedef (ITK_CACHE_LINE_ALIGNMENT, PaddedGetValueAndDerivativePerThreadStruct, AlignedGetValueAndDerivativePerThreadStruct)
itkAlignedTypedef (ITK_CACHE_LINE_ALIGNMENT, PaddedGetValuePerThreadStruct, AlignedGetValuePerThreadStruct)
itkPadStruct (ITK_CACHE_LINE_ALIGNMENT, GetValueAndDerivativePerThreadStruct , PaddedGetValueAndDerivativePerThreadStruct)
itkPadStruct (ITK_CACHE_LINE_ALIGNMENT, GetValuePerThreadStruct , PaddedGetValuePerThreadStruct)
void LaunchGetValueAndDerivativeThreaderCallback () const
void LaunchGetValueThreaderCallback () const
void PrintSelf (std::ostream &os, Indent indent) const override
virtual void SetUseFixedImageLimiter (bool _arg)
virtual void SetUseImageSampler (bool _arg)
virtual void SetUseMovingImageLimiter (bool _arg)
virtual void ThreadedGetValue (ThreadIdType threadID)
virtual void ThreadedGetValueAndDerivative (ThreadIdType threadID)
MovingImagePointType TransformPoint (const FixedImagePointType &fixedImagePoint) const
~AdvancedImageToImageMetric () override=default