1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
|
/*=========================================================================
*
* Copyright NumFOCUS
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* https://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
#ifndef itkPointSetToPointSetMetric_h
#define itkPointSetToPointSetMetric_h
#include "itkImageBase.h"
#include "itkTransform.h"
#include "itkMultipleValuedCostFunction.h"
#include "itkMacro.h"
#include "itkGradientRecursiveGaussianImageFilter.h"
namespace itk
{
/** \class PointSetToPointSetMetric
* \brief Computes similarity between two point sets.
*
* This Class is templated over the type of the two point-sets. It
* expects a Transform to be plugged in. This particular
* class is the base class for a hierarchy of point-set to point-set metrics.
*
* This class computes a value that measures the similarity between the fixed point-set
* and the transformed moving point-set.
*
* \ingroup RegistrationMetrics
*
* \ingroup ITKRegistrationCommon
*/
template <typename TFixedPointSet, typename TMovingPointSet>
class ITK_TEMPLATE_EXPORT PointSetToPointSetMetric : public MultipleValuedCostFunction
{
public:
ITK_DISALLOW_COPY_AND_MOVE(PointSetToPointSetMetric);
/** Standard class type aliases. */
using Self = PointSetToPointSetMetric;
using Superclass = MultipleValuedCostFunction;
using Pointer = SmartPointer<Self>;
using ConstPointer = SmartPointer<const Self>;
/** Type used for representing point components */
using CoordinateRepresentationType = Superclass::ParametersValueType;
/** \see LightObject::GetNameOfClass() */
itkOverrideGetNameOfClassMacro(PointSetToPointSetMetric);
/** Type of the moving Pointset. */
using MovingPointSetType = TMovingPointSet;
using MovingPointSetPixelType = typename TMovingPointSet::PixelType;
using MovingPointSetConstPointer = typename MovingPointSetType::ConstPointer;
/** Type of the fixed Pointset. */
using FixedPointSetType = TFixedPointSet;
using FixedPointSetConstPointer = typename FixedPointSetType::ConstPointer;
/** Constants for the pointset dimensions */
static constexpr unsigned int MovingPointSetDimension = TMovingPointSet::PointDimension;
static constexpr unsigned int FixedPointSetDimension = TFixedPointSet::PointDimension;
using FixedPointIterator = typename FixedPointSetType::PointsContainer::ConstIterator;
using FixedPointDataIterator = typename FixedPointSetType::PointDataContainer::ConstIterator;
using MovingPointIterator = typename MovingPointSetType::PointsContainer::ConstIterator;
using MovingPointDataIterator = typename MovingPointSetType::PointDataContainer::ConstIterator;
/** Type of the Transform Base class */
using TransformType =
Transform<CoordinateRepresentationType, Self::MovingPointSetDimension, Self::FixedPointSetDimension>;
using TransformPointer = typename TransformType::Pointer;
using InputPointType = typename TransformType::InputPointType;
using OutputPointType = typename TransformType::OutputPointType;
using TransformParametersType = typename TransformType::ParametersType;
using TransformJacobianType = typename TransformType::JacobianType;
/** Type of the measure. */
using MeasureType = Superclass::MeasureType;
/** Type of the derivative. */
using DerivativeType = Superclass::DerivativeType;
/** Type of the parameters. */
using ParametersType = Superclass::ParametersType;
/** Get/Set the Fixed Pointset. */
itkSetConstObjectMacro(FixedPointSet, FixedPointSetType);
itkGetConstObjectMacro(FixedPointSet, FixedPointSetType);
/** Get/Set the Moving Pointset. */
itkSetConstObjectMacro(MovingPointSet, MovingPointSetType);
itkGetConstObjectMacro(MovingPointSet, MovingPointSetType);
/** Connect the Transform. */
itkSetObjectMacro(Transform, TransformType);
/** Get a pointer to the Transform. */
itkGetModifiableObjectMacro(Transform, TransformType);
/** Set the parameters defining the Transform. */
void
SetTransformParameters(const ParametersType & parameters) const;
/** Return the number of parameters required by the Transform */
unsigned int
GetNumberOfParameters() const override
{
return m_Transform->GetNumberOfParameters();
}
/** Initialize the Metric by making sure that all the components are present and plugged together correctly. */
virtual void
Initialize();
protected:
PointSetToPointSetMetric();
~PointSetToPointSetMetric() override = default;
void
PrintSelf(std::ostream & os, Indent indent) const override;
FixedPointSetConstPointer m_FixedPointSet{};
MovingPointSetConstPointer m_MovingPointSet{};
mutable TransformPointer m_Transform{};
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
# include "itkPointSetToPointSetMetric.hxx"
#endif
#endif
|
