ITK  5.4.0
Insight Toolkit
itkDemonsRegistrationFunction.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright NumFOCUS
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * https://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 #ifndef itkDemonsRegistrationFunction_h
19 #define itkDemonsRegistrationFunction_h
20 
22 #include "itkPoint.h"
25 #include <mutex>
26 
27 namespace itk
28 {
53 template <typename TFixedImage, typename TMovingImage, typename TDisplacementField>
54 class ITK_TEMPLATE_EXPORT DemonsRegistrationFunction
55  : public PDEDeformableRegistrationFunction<TFixedImage, TMovingImage, TDisplacementField>
56 {
57 public:
58  ITK_DISALLOW_COPY_AND_MOVE(DemonsRegistrationFunction);
59 
65 
67  itkNewMacro(Self);
68 
70  itkOverrideGetNameOfClassMacro(DemonsRegistrationFunction);
71 
73  using typename Superclass::MovingImageType;
74  using typename Superclass::MovingImagePointer;
75 
77  using typename Superclass::FixedImageType;
78  using typename Superclass::FixedImagePointer;
81  using SpacingType = typename FixedImageType::SpacingType;
82 
84  using typename Superclass::DisplacementFieldType;
85  using typename Superclass::DisplacementFieldTypePointer;
86 
88  static constexpr unsigned int ImageDimension = Superclass::ImageDimension;
89 
91  using typename Superclass::PixelType;
92  using typename Superclass::RadiusType;
93  using typename Superclass::NeighborhoodType;
94  using typename Superclass::FloatOffsetType;
95  using typename Superclass::TimeStepType;
96 
98  using CoordRepType = double;
103 
106 
110 
114 
116  void
118  {
119  m_MovingImageInterpolator = ptr;
120  }
121 
123  InterpolatorType *
125  {
126  return m_MovingImageInterpolator;
127  }
128 
130  TimeStepType
131  ComputeGlobalTimeStep(void * itkNotUsed(GlobalData)) const override
132  {
133  return m_TimeStep;
134  }
135 
138  void *
139  GetGlobalDataPointer() const override
140  {
141  auto * global = new GlobalDataStruct();
142 
143  global->m_SumOfSquaredDifference = 0.0;
144  global->m_NumberOfPixelsProcessed = 0L;
145  global->m_SumOfSquaredChange = 0;
146  return global;
147  }
148 
153  void
154  ReleaseGlobalDataPointer(void * gd) const override;
155 
157  void
158  InitializeIteration() override;
159 
164  PixelType
165  ComputeUpdate(const NeighborhoodType & it, void * gd, const FloatOffsetType & offset = FloatOffsetType(0.0)) override;
166 
170  virtual double
171  GetMetric() const
172  {
173  return m_Metric;
174  }
175 
177  virtual double
178  GetRMSChange() const
179  {
180  return m_RMSChange;
181  }
182 
186  virtual void
188  {
189  m_UseMovingImageGradient = flag;
190  }
191  virtual bool
193  {
194  return m_UseMovingImageGradient;
195  }
202  virtual void
203  SetIntensityDifferenceThreshold(double);
204 
205  virtual double
206  GetIntensityDifferenceThreshold() const;
207 
208 protected:
210  ~DemonsRegistrationFunction() override = default;
211  void
212  PrintSelf(std::ostream & os, Indent indent) const override;
213 
216 
220  {
224  };
225 
226 private:
228  // SpacingType m_FixedImageSpacing;
229  // PointType m_FixedImageOrigin;
230  PixelType m_ZeroUpdateReturn{};
231  double m_Normalizer{};
232 
234  GradientCalculatorPointer m_FixedImageGradientCalculator{};
235 
237  MovingImageGradientCalculatorPointer m_MovingImageGradientCalculator{};
238  bool m_UseMovingImageGradient{};
239 
241  InterpolatorPointer m_MovingImageInterpolator{};
242 
244  TimeStepType m_TimeStep{};
245 
247  double m_DenominatorThreshold{};
248 
250  double m_IntensityDifferenceThreshold{};
251 
255  mutable double m_Metric{};
256  mutable double m_SumOfSquaredDifference{};
257  mutable SizeValueType m_NumberOfPixelsProcessed{};
258  mutable double m_RMSChange{};
259  mutable double m_SumOfSquaredChange{};
260 
262  mutable std::mutex m_MetricCalculationMutex{};
263 };
264 } // end namespace itk
265 
266 #ifndef ITK_MANUAL_INSTANTIATION
267 # include "itkDemonsRegistrationFunction.hxx"
268 #endif
269 
270 #endif
Pointer
SmartPointer< Self > Pointer
Definition: itkAddImageFilter.h:93
itk::FiniteDifferenceFunction< TDisplacementField >::RadiusType
typename ConstNeighborhoodIterator< TDisplacementField >::RadiusType RadiusType
Definition: itkFiniteDifferenceFunction.h:97
itk::DemonsRegistrationFunction
Definition: itkDemonsRegistrationFunction.h:54
itk::DemonsRegistrationFunction::SpacingType
typename FixedImageType::SpacingType SpacingType
Definition: itkDemonsRegistrationFunction.h:81
itk::DemonsRegistrationFunction::GlobalDataStruct::m_NumberOfPixelsProcessed
SizeValueType m_NumberOfPixelsProcessed
Definition: itkDemonsRegistrationFunction.h:222
itkLinearInterpolateImageFunction.h
itkPoint.h
itk::DemonsRegistrationFunction::GlobalDataStruct::m_SumOfSquaredChange
double m_SumOfSquaredChange
Definition: itkDemonsRegistrationFunction.h:223
itk::GTest::TypedefsAndConstructors::Dimension2::PointType
ImageBaseType::PointType PointType
Definition: itkGTestTypedefsAndConstructors.h:51
itkPDEDeformableRegistrationFunction.h
itk::CentralDifferenceImageFunction
Calculate the derivative by central differencing.
Definition: itkCentralDifferenceImageFunction.h:76
itk::Vector
A templated class holding a n-Dimensional vector.
Definition: itkVector.h:62
itk::GTest::TypedefsAndConstructors::Dimension2::SizeType
ImageBaseType::SizeType SizeType
Definition: itkGTestTypedefsAndConstructors.h:49
itk::SmartPointer< Self >
itk::Indent
Control indentation during Print() invocation.
Definition: itkIndent.h:49
itk::DemonsRegistrationFunction::GetUseMovingImageGradient
virtual bool GetUseMovingImageGradient() const
Definition: itkDemonsRegistrationFunction.h:192
itk::DemonsRegistrationFunction::PointType
typename InterpolatorType::PointType PointType
Definition: itkDemonsRegistrationFunction.h:101
itk::DemonsRegistrationFunction::SizeType
typename FixedImageType::SizeType SizeType
Definition: itkDemonsRegistrationFunction.h:80
itk::GTest::TypedefsAndConstructors::Dimension2::IndexType
ImageBaseType::IndexType IndexType
Definition: itkGTestTypedefsAndConstructors.h:50
itk::DemonsRegistrationFunction::GetGlobalDataPointer
void * GetGlobalDataPointer() const override
Definition: itkDemonsRegistrationFunction.h:139
itk::LinearInterpolateImageFunction
Linearly interpolate an image at specified positions.
Definition: itkLinearInterpolateImageFunction.h:51
itk::LightObject
Light weight base class for most itk classes.
Definition: itkLightObject.h:55
itk::DemonsRegistrationFunction::IndexType
typename FixedImageType::IndexType IndexType
Definition: itkDemonsRegistrationFunction.h:79
itkCentralDifferenceImageFunction.h
itk::DemonsRegistrationFunction::CoordRepType
double CoordRepType
Definition: itkDemonsRegistrationFunction.h:98
itk::DemonsRegistrationFunction::InterpolatorPointer
typename InterpolatorType::Pointer InterpolatorPointer
Definition: itkDemonsRegistrationFunction.h:100
itk::DemonsRegistrationFunction::GetMetric
virtual double GetMetric() const
Definition: itkDemonsRegistrationFunction.h:171
itk::PDEDeformableRegistrationFunction
Definition: itkPDEDeformableRegistrationFunction.h:41
itk::FiniteDifferenceFunction< TDisplacementField >::PixelType
typename ImageType::PixelType PixelType
Definition: itkFiniteDifferenceFunction.h:83
itk::DemonsRegistrationFunction::GetRMSChange
virtual double GetRMSChange() const
Definition: itkDemonsRegistrationFunction.h:178
itk::DemonsRegistrationFunction::ComputeGlobalTimeStep
TimeStepType ComputeGlobalTimeStep(void *) const override
Definition: itkDemonsRegistrationFunction.h:131
itk::CovariantVector
A templated class holding a n-Dimensional covariant vector.
Definition: itkCovariantVector.h:70
itk::DemonsRegistrationFunction::SetMovingImageInterpolator
void SetMovingImageInterpolator(InterpolatorType *ptr)
Definition: itkDemonsRegistrationFunction.h:117
itk::DemonsRegistrationFunction::SetUseMovingImageGradient
virtual void SetUseMovingImageGradient(bool flag)
Definition: itkDemonsRegistrationFunction.h:187
itk
The "itk" namespace contains all Insight Segmentation and Registration Toolkit (ITK) classes....
Definition: itkAnnulusOperator.h:24
itk::ConstNeighborhoodIterator
Const version of NeighborhoodIterator, defining iteration of a local N-dimensional neighborhood of pi...
Definition: itkConstNeighborhoodIterator.h:51
itk::DemonsRegistrationFunction::MovingImageGradientCalculatorPointer
typename MovingImageGradientCalculatorType::Pointer MovingImageGradientCalculatorPointer
Definition: itkDemonsRegistrationFunction.h:113
itk::DemonsRegistrationFunction::GlobalDataStruct::m_SumOfSquaredDifference
double m_SumOfSquaredDifference
Definition: itkDemonsRegistrationFunction.h:221
itk::DemonsRegistrationFunction::GetMovingImageInterpolator
InterpolatorType * GetMovingImageInterpolator()
Definition: itkDemonsRegistrationFunction.h:124
itk::FiniteDifferenceFunction< TDisplacementField >
itk::DemonsRegistrationFunction::GlobalDataStruct
Definition: itkDemonsRegistrationFunction.h:219
itk::DemonsRegistrationFunction::GradientCalculatorPointer
typename GradientCalculatorType::Pointer GradientCalculatorPointer
Definition: itkDemonsRegistrationFunction.h:109
itk::InterpolateImageFunction
Base class for all image interpolators.
Definition: itkInterpolateImageFunction.h:45
itk::SizeValueType
unsigned long SizeValueType
Definition: itkIntTypes.h:83
itk::FiniteDifferenceFunction< TDisplacementField >::TimeStepType
double TimeStepType
Definition: itkFiniteDifferenceFunction.h:90