18#ifndef itkAdvancedMeanSquaresImageToImageMetric_h
19#define itkAdvancedMeanSquaresImageToImageMetric_h
23#include "itkSmoothingRecursiveGaussianImageFilter.h"
25#include "itkNearestNeighborInterpolateImageFunction.h"
55template <
class TFixedImage,
class TMovingImage>
75 using typename Superclass::CoordinateRepresentationType;
76 using typename Superclass::MovingImageType;
77 using typename Superclass::MovingImagePixelType;
78 using typename Superclass::MovingImageConstPointer;
79 using typename Superclass::FixedImageType;
80 using typename Superclass::FixedImageConstPointer;
81 using typename Superclass::FixedImageRegionType;
82 using typename Superclass::TransformType;
83 using typename Superclass::TransformPointer;
84 using typename Superclass::InputPointType;
85 using typename Superclass::OutputPointType;
86 using typename Superclass::TransformParametersType;
87 using typename Superclass::TransformJacobianType;
88 using typename Superclass::NumberOfParametersType;
89 using typename Superclass::InterpolatorType;
90 using typename Superclass::InterpolatorPointer;
91 using typename Superclass::RealType;
92 using typename Superclass::GradientPixelType;
93 using typename Superclass::GradientImageType;
94 using typename Superclass::GradientImagePointer;
95 using typename Superclass::GradientImageFilterType;
96 using typename Superclass::GradientImageFilterPointer;
97 using typename Superclass::FixedImageMaskType;
98 using typename Superclass::FixedImageMaskPointer;
99 using typename Superclass::MovingImageMaskType;
100 using typename Superclass::MovingImageMaskPointer;
101 using typename Superclass::MeasureType;
102 using typename Superclass::DerivativeType;
103 using typename Superclass::DerivativeValueType;
104 using typename Superclass::ParametersType;
105 using typename Superclass::FixedImagePixelType;
106 using typename Superclass::MovingImageRegionType;
107 using typename Superclass::ImageSamplerType;
108 using typename Superclass::ImageSamplerPointer;
109 using typename Superclass::ImageSampleContainerType;
110 using typename Superclass::ImageSampleContainerPointer;
111 using typename Superclass::FixedImageLimiterType;
112 using typename Superclass::MovingImageLimiterType;
113 using typename Superclass::FixedImageLimiterOutputType;
114 using typename Superclass::MovingImageLimiterOutputType;
115 using typename Superclass::MovingImageDerivativeScalesType;
116 using typename Superclass::HessianValueType;
117 using typename Superclass::HessianType;
118 using typename Superclass::ThreaderType;
119 using typename Superclass::ThreadInfoType;
121 using typename Superclass::FixedImageMaskSpatialObject2Type;
122 using typename Superclass::MovingImageMaskSpatialObject2Type;
135 GetValue(
const TransformParametersType & parameters)
const override;
139 GetDerivative(
const TransformParametersType & parameters, DerivativeType & derivative)
const override;
145 DerivativeType & derivative)
const;
150 DerivativeType & derivative)
const override;
157 itkSetMacro(SelfHessianSmoothingSigma,
double);
158 itkGetConstMacro(SelfHessianSmoothingSigma,
double);
161 itkSetMacro(SelfHessianNoiseRange,
double);
162 itkGetConstMacro(SelfHessianNoiseRange,
double);
165 itkSetMacro(NumberOfSamplesForSelfHessian,
unsigned int);
166 itkGetConstMacro(NumberOfSamplesForSelfHessian,
unsigned int);
182 itkSetMacro(UseNormalization,
bool);
183 itkGetConstMacro(UseNormalization,
bool);
190 itkSetMacro(UseOpenMP,
bool);
197 PrintSelf(std::ostream & os, Indent indent)
const override;
202 using typename Superclass::FixedImageIndexType;
203 using typename Superclass::FixedImageIndexValueType;
204 using typename Superclass::MovingImageIndexType;
205 using typename Superclass::FixedImagePointType;
206 using typename Superclass::MovingImagePointType;
207 using typename Superclass::MovingImageContinuousIndexType;
208 using typename Superclass::BSplineInterpolatorType;
209 using typename Superclass::CentralDifferenceGradientFilterType;
210 using typename Superclass::MovingImageDerivativeType;
211 using typename Superclass::NonZeroJacobianIndicesType;
214 using SmootherType = SmoothingRecursiveGaussianImageFilter<FixedImageType, FixedImageType>;
217 NearestNeighborInterpolateImageFunction<FixedImageType, CoordinateRepresentationType>;
226 const RealType movingImageValue,
227 const DerivativeType & imageJacobian,
229 MeasureType & measure,
230 DerivativeType & deriv)
const;
264#ifndef ITK_MANUAL_INSTANTIATION
265# include "itkAdvancedMeanSquaresImageToImageMetric.hxx"
An extension of the ITK ImageToImageMetric. It is the intended base class for all elastix metrics.
ImageToImageMetric< TFixedImage, TMovingImage > Superclass
typename AdvancedTransformType::NonZeroJacobianIndicesType NonZeroJacobianIndicesType
SmartPointer< const Self > ConstPointer
vnl_sparse_matrix< HessianValueType > HessianType
Compute Mean square difference between two images, based on AdvancedImageToImageMetric....
void ThreadedGetValueAndDerivative(ThreadIdType threadID) override
double m_SelfHessianNoiseRange
BSplineInterpolateImageFunction< FixedImageType, CoordinateRepresentationType > FixedImageInterpolatorType
void AfterThreadedGetValue(MeasureType &value) const override
unsigned int m_NumberOfSamplesForSelfHessian
void GetDerivative(const TransformParametersType ¶meters, DerivativeType &derivative) const override
void UpdateSelfHessianTerms(const DerivativeType &imageJacobian, const NonZeroJacobianIndicesType &nzji, HessianType &H) const
SmoothingRecursiveGaussianImageFilter< FixedImageType, FixedImageType > SmootherType
void ThreadedGetValue(ThreadIdType threadID) override
void AfterThreadedGetValueAndDerivative(MeasureType &value, DerivativeType &derivative) const override
~AdvancedMeanSquaresImageToImageMetric() override=default
virtual MeasureType GetValueSingleThreaded(const TransformParametersType ¶meters) const
AdvancedMeanSquaresImageToImageMetric()
void Initialize() override
void GetSelfHessian(const TransformParametersType ¶meters, HessianType &H) const override
NearestNeighborInterpolateImageFunction< FixedImageType, CoordinateRepresentationType > DummyFixedImageInterpolatorType
double m_NormalizationFactor
void GetValueAndDerivative(const TransformParametersType ¶meters, MeasureType &value, DerivativeType &derivative) const override
ITK_DISALLOW_COPY_AND_MOVE(AdvancedMeanSquaresImageToImageMetric)
void PrintSelf(std::ostream &os, Indent indent) const override
itkStaticConstMacro(FixedImageDimension, unsigned int, FixedImageType::ImageDimension)
itkStaticConstMacro(MovingImageDimension, unsigned int, MovingImageType::ImageDimension)
double m_SelfHessianSmoothingSigma
void UpdateValueAndDerivativeTerms(const RealType fixedImageValue, const RealType movingImageValue, const DerivativeType &imageJacobian, const NonZeroJacobianIndicesType &nzji, MeasureType &measure, DerivativeType &deriv) const
void GetValueAndDerivativeSingleThreaded(const TransformParametersType ¶meters, MeasureType &value, DerivativeType &derivative) const
MeasureType GetValue(const TransformParametersType ¶meters) const override
Samples image voxels on a regular grid.