itkCenteredTransformInitializer2.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright UMC Utrecht and contributors
 *
 *  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
 *
 *        http://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.
 *
 *=========================================================================*/
/*=========================================================================
 
  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkCenteredTransformInitializer2.h,v $
  Date:      $Date: 2010-07-04 10:30:49 $
  Version:   $Revision: 1.11 $
 
  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
 
     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notices for more information.
 
=========================================================================*/
#ifndef itkCenteredTransformInitializer2_h
#define itkCenteredTransformInitializer2_h
 
#include "itkObject.h"
#include "itkObjectFactory.h"
#include "itkImageMaskSpatialObject.h"
#include "itkAdvancedImageMomentsCalculator.h"
 
#include <iostream>
 
namespace itk
{

template <class TTransform, class TFixedImage, class TMovingImage>
class ITK_TEMPLATE_EXPORT CenteredTransformInitializer2 : public Object
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(CenteredTransformInitializer2);

  using Self = CenteredTransformInitializer2;
  using Superclass = Object;
  using Pointer = SmartPointer<Self>;
  using ConstPointer = SmartPointer<const Self>;

  itkNewMacro(Self);

  itkTypeMacro(CenteredTransformInitializer2, Object);

  using TransformType = TTransform;
  using TransformPointer = typename TransformType::Pointer;

  itkStaticConstMacro(InputSpaceDimension, unsigned int, TransformType::InputSpaceDimension);
  itkStaticConstMacro(OutputSpaceDimension, unsigned int, TransformType::OutputSpaceDimension);

  using FixedImageType = TFixedImage;
  using MovingImageType = TMovingImage;
 
  using FixedImagePointer = typename FixedImageType::ConstPointer;
  using MovingImagePointer = typename MovingImageType::ConstPointer;
 
  using FixedImageMaskType = Image<unsigned char, InputSpaceDimension>;
  using MovingImageMaskType = Image<unsigned char, OutputSpaceDimension>;
  using FixedImageMaskPointer = typename FixedImageMaskType::ConstPointer;
  using MovingImageMaskPointer = typename MovingImageMaskType::ConstPointer;
  using FixedImageMaskConstPointer = typename FixedImageMaskType::ConstPointer;
  using MovingImageMaskConstPointer = typename MovingImageMaskType::ConstPointer;

  using FixedImageCalculatorType = AdvancedImageMomentsCalculator<FixedImageType>;
  using MovingImageCalculatorType = AdvancedImageMomentsCalculator<MovingImageType>;
 
  using FixedImageCalculatorPointer = typename FixedImageCalculatorType::Pointer;
  using MovingImageCalculatorPointer = typename MovingImageCalculatorType::Pointer;

  using OffsetType = typename TransformType::OffsetType;

  using InputPointType = typename TransformType::InputPointType;

  using OutputVectorType = typename TransformType::OutputVectorType;
 
  using InputPixelType = typename FixedImageType::PixelType;

  itkSetObjectMacro(Transform, TransformType);

  itkSetConstObjectMacro(FixedImage, FixedImageType);

  itkSetConstObjectMacro(MovingImage, MovingImageType);

  itkSetConstObjectMacro(FixedImageMask, FixedImageMaskType);
  itkSetConstObjectMacro(MovingImageMask, MovingImageMaskType);

  itkSetMacro(NumberOfSamplesForCenteredTransformInitialization, SizeValueType);
  itkSetMacro(LowerThresholdForCenterGravity, InputPixelType);
  itkSetMacro(CenterOfGravityUsesLowerThreshold, bool);

  virtual void
  InitializeTransform();

  void
  GeometryOn()
  {
    m_UseMoments = false;
    m_UseOrigins = false;
    m_UseTop = false;
  }
  void
  MomentsOn()
  {
    m_UseMoments = true;
    m_UseOrigins = false;
    m_UseTop = false;
  }
  void
  OriginsOn()
  {
    m_UseMoments = false;
    m_UseOrigins = true;
    m_UseTop = false;
  }
  void
  GeometryTopOn()
  {
    m_UseMoments = false;
    m_UseOrigins = false;
    m_UseTop = true;
  }

  itkGetConstObjectMacro(FixedCalculator, FixedImageCalculatorType);
  itkGetConstObjectMacro(MovingCalculator, MovingImageCalculatorType);
 
protected:
  CenteredTransformInitializer2();
  ~CenteredTransformInitializer2() override = default;
 
  void
  PrintSelf(std::ostream & os, Indent indent) const override;
 
  itkGetModifiableObjectMacro(Transform, TransformType);

  SizeValueType  m_NumberOfSamplesForCenteredTransformInitialization{};
  InputPixelType m_LowerThresholdForCenterGravity{};
  bool           m_CenterOfGravityUsesLowerThreshold{};
 
private:
  TransformPointer m_Transform{};
 
  FixedImagePointer      m_FixedImage{};
  MovingImagePointer     m_MovingImage{};
  FixedImageMaskPointer  m_FixedImageMask{};
  MovingImageMaskPointer m_MovingImageMask{};
 
  bool m_UseMoments{};
  bool m_UseOrigins{};
  bool m_UseTop{};
 
  FixedImageCalculatorPointer  m_FixedCalculator{};
  MovingImageCalculatorPointer m_MovingCalculator{};
};
 
} // namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkCenteredTransformInitializer2.hxx"
#endif
 
#endif /* itkCenteredTransformInitializer2_h */