18#ifndef itkPCAMetric_F_multithreaded_h
19#define itkPCAMetric_F_multithreaded_h
23#include "itkSmoothingRecursiveGaussianImageFilter.h"
25#include "itkNearestNeighborInterpolateImageFunction.h"
26#include "itkExtractImageFilter.h"
31template <
class TFixedImage,
class TMovingImage>
32class ITK_TEMPLATE_EXPORT PCAMetric :
public AdvancedImageToImageMetric<TFixedImage, TMovingImage>
43 using typename Superclass::FixedImageRegionType;
53 itkSetMacro(SubtractMean,
bool);
55 itkSetMacro(TransformIsStackTransform,
bool);
56 itkSetMacro(NumEigenValues,
unsigned int);
59 using typename Superclass::CoordinateRepresentationType;
60 using typename Superclass::MovingImageType;
61 using typename Superclass::MovingImagePixelType;
62 using typename Superclass::MovingImageConstPointer;
63 using typename Superclass::FixedImageType;
64 using typename Superclass::FixedImageConstPointer;
65 using typename Superclass::TransformType;
66 using typename Superclass::TransformPointer;
67 using typename Superclass::InputPointType;
68 using typename Superclass::OutputPointType;
69 using typename Superclass::TransformParametersType;
70 using typename Superclass::TransformJacobianType;
71 using typename Superclass::InterpolatorType;
72 using typename Superclass::InterpolatorPointer;
73 using typename Superclass::RealType;
74 using typename Superclass::GradientPixelType;
75 using typename Superclass::GradientImageType;
76 using typename Superclass::GradientImagePointer;
77 using typename Superclass::GradientImageFilterType;
78 using typename Superclass::GradientImageFilterPointer;
79 using typename Superclass::FixedImageMaskType;
80 using typename Superclass::FixedImageMaskPointer;
81 using typename Superclass::MovingImageMaskType;
82 using typename Superclass::MovingImageMaskPointer;
83 using typename Superclass::MeasureType;
84 using typename Superclass::DerivativeType;
85 using typename Superclass::ParametersType;
86 using typename Superclass::FixedImagePixelType;
87 using typename Superclass::MovingImageRegionType;
88 using typename Superclass::ImageSamplerType;
89 using typename Superclass::ImageSamplerPointer;
90 using typename Superclass::ImageSampleContainerType;
91 using typename Superclass::ImageSampleContainerPointer;
92 using typename Superclass::FixedImageLimiterType;
93 using typename Superclass::MovingImageLimiterType;
94 using typename Superclass::FixedImageLimiterOutputType;
95 using typename Superclass::MovingImageLimiterOutputType;
96 using typename Superclass::MovingImageDerivativeScalesType;
98 using typename Superclass::ThreaderType;
99 using typename Superclass::ThreadInfoType;
115 GetValue(
const TransformParametersType & parameters)
const override;
119 GetDerivative(
const TransformParametersType & parameters, DerivativeType & derivative)
const override;
125 DerivativeType & Derivative)
const;
130 DerivativeType & Derivative)
const override;
143 PrintSelf(std::ostream & os, Indent indent)
const override;
148 using typename Superclass::FixedImageIndexType;
149 using typename Superclass::FixedImageIndexValueType;
150 using typename Superclass::MovingImageIndexType;
151 using typename Superclass::FixedImagePointType;
153 typename itk::ContinuousIndex<CoordinateRepresentationType, FixedImageDimension>;
154 using typename Superclass::MovingImagePointType;
155 using typename Superclass::MovingImageContinuousIndexType;
156 using typename Superclass::BSplineInterpolatorType;
157 using typename Superclass::CentralDifferenceGradientFilterType;
158 using typename Superclass::MovingImageDerivativeType;
159 using typename Superclass::NonZeroJacobianIndicesType;
167 DerivativeType & imageJacobian)
const override;
184 static ITK_THREAD_RETURN_FUNCTION_CALL_CONVENTION
187 static ITK_THREAD_RETURN_FUNCTION_CALL_CONVENTION
220 PaddedPCAMetricGetSamplesPerThreadStruct,
221 AlignedPCAMetricGetSamplesPerThreadStruct);
229 bool m_SubtractMean{
false };
232 FixedImageSizeType m_GridSize;
235 bool m_TransformIsStackTransform{
false };
238 unsigned int m_NumEigenValues{ 6 };
251#ifndef ITK_MANUAL_INSTANTIATION
252# include "itkPCAMetric_F_multithreaded.hxx"
An extension of the ITK ImageToImageMetric. It is the intended base class for all elastix metrics.
ImageToImageMetric< TFixedImage, TMovingImage > Superclass
typename DerivativeType::ValueType DerivativeValueType
typename BSplineInterpolatorType::CovariantVectorType MovingImageDerivativeType
SmartPointer< const Self > ConstPointer
void GetValueAndDerivativeSingleThreaded(const TransformParametersType ¶meters, MeasureType &Value, DerivativeType &Derivative) const
void InitializeThreadingParameters() const override
void AfterThreadedGetSamples(MeasureType &value) const
void PrintSelf(std::ostream &os, Indent indent) const override
ITK_DISALLOW_COPY_AND_MOVE(PCAMetric)
itkStaticConstMacro(FixedImageDimension, unsigned int, FixedImageType::ImageDimension)
std::vector< unsigned int > m_PixelStartIndex
unsigned int m_LastDimIndex
typename itk::ContinuousIndex< CoordinateRepresentationType, FixedImageDimension > FixedImageContinuousIndexType
PCAMetricMultiThreaderParameterType m_PCAMetricThreaderParameters
itkAlignedTypedef(ITK_CACHE_LINE_ALIGNMENT, PaddedPCAMetricGetSamplesPerThreadStruct, AlignedPCAMetricGetSamplesPerThreadStruct)
void GetValueAndDerivative(const TransformParametersType ¶meters, MeasureType &Value, DerivativeType &Derivative) const override
MeasureType GetValue(const TransformParametersType ¶meters) const override
void EvaluateTransformJacobianInnerProduct(const TransformJacobianType &jacobian, const MovingImageDerivativeType &movingImageDerivative, DerivativeType &imageJacobian) const override
DerivativeMatrixType m_Sv
DerivativeMatrixType m_vSAtmm
typename FixedImageRegionType::SizeType FixedImageSizeType
DerivativeMatrixType m_CSv
itkStaticConstMacro(MovingImageDimension, unsigned int, MovingImageType::ImageDimension)
void LaunchComputeDerivativeThreaderCallback() const
void Initialize() override
void ThreadedComputeDerivative(ThreadIdType threadID)
static ITK_THREAD_RETURN_FUNCTION_CALL_CONVENTION GetSamplesThreaderCallback(void *arg)
std::vector< AlignedPCAMetricGetSamplesPerThreadStruct > m_PCAMetricGetSamplesPerThreadVariables
void AfterThreadedComputeDerivative(DerivativeType &derivative) const
void LaunchGetSamplesThreaderCallback() const
static ITK_THREAD_RETURN_FUNCTION_CALL_CONVENTION ComputeDerivativeThreaderCallback(void *arg)
~PCAMetric() override=default
void ThreadedGetSamples(ThreadIdType threadID)
itkPadStruct(ITK_CACHE_LINE_ALIGNMENT, PCAMetricGetSamplesPerThreadStruct, PaddedPCAMetricGetSamplesPerThreadStruct)
DerivativeMatrixType m_vdSdmu_part1
void GetDerivative(const TransformParametersType ¶meters, DerivativeType &derivative) const override
SizeValueType st_NumberOfPixelsCounted
DerivativeType st_Derivative
std::vector< FixedImagePointType > st_ApprovedSamples