itkAffineTransform2D3D.h

Go to the documentation of this file.

/*=============================================================================

  NifTK: A software platform for medical image computing.

  Copyright (c) University College London (UCL). All rights reserved.

  This software is distributed WITHOUT ANY WARRANTY; without even
  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  PURPOSE.

  See LICENSE.txt in the top level directory for details.

=============================================================================*/

#ifndef itkAffineTransform2D3D_h
#define itkAffineTransform2D3D_h

#include <iostream>

#include <itkMatrix.h>
#include <itkMatrixOffsetTransformBase.h>
#include <itkExceptionObject.h>
#include <itkMacro.h>

namespace itk
{

template <
 class TScalarType=double,         // Data type for scalars 
                                   //    (e.g. float or double)
 unsigned int NDimensions=3>       // Number of dimensions in the input space
class ITK_EXPORT AffineTransform2D3D
: public MatrixOffsetTransformBase< TScalarType, NDimensions, NDimensions >
{
public:
  typedef AffineTransform2D3D                           Self;
  typedef MatrixOffsetTransformBase< TScalarType,
                                     NDimensions,
                                     NDimensions >  Superclass;
  typedef SmartPointer<Self>                        Pointer;
  typedef SmartPointer<const Self>                  ConstPointer;
  
  itkTypeMacro( AffineTransform2D3D, MatrixOffsetTransformBase );

  itkNewMacro( Self );

  itkStaticConstMacro(InputSpaceDimension, unsigned int, NDimensions);
  itkStaticConstMacro(OutputSpaceDimension, unsigned int, NDimensions);
  itkStaticConstMacro(SpaceDimension, unsigned int, NDimensions);
  itkStaticConstMacro(ParametersDimension, unsigned int,
                                           NDimensions*(NDimensions+1));

  
  typedef typename Superclass::ParametersType         ParametersType;
  typedef typename Superclass::JacobianType           JacobianType;
  typedef typename Superclass::ScalarType             ScalarType;
  typedef typename Superclass::InputPointType         InputPointType;
  typedef typename Superclass::OutputPointType        OutputPointType;
  typedef typename Superclass::InputVectorType        InputVectorType;
  typedef typename Superclass::OutputVectorType       OutputVectorType;
  typedef typename Superclass::InputVnlVectorType     InputVnlVectorType;
  typedef typename Superclass::OutputVnlVectorType    OutputVnlVectorType;
  typedef typename Superclass::InputCovariantVectorType 
                                                      InputCovariantVectorType;
  typedef typename Superclass::OutputCovariantVectorType
                                                      OutputCovariantVectorType;
  typedef typename Superclass::MatrixType             MatrixType;
  typedef typename Superclass::InverseMatrixType      InverseMatrixType;
  typedef typename Superclass::CenterType             CenterType;
  typedef typename Superclass::OffsetType             OffsetType;
  typedef typename Superclass::TranslationType        TranslationType;

  typedef typename Superclass::InverseTransformBaseType InverseTransformBaseType;
  typedef typename InverseTransformBaseType::Pointer    InverseTransformBasePointer;

  virtual void SetParameters( const ParametersType & parameters );

  virtual const ParametersType& GetParameters(void) const;

  void Translate(const OutputVectorType &trans);

  void Scale(const OutputVectorType &factor);

  void Rotate(const OutputVectorType &angle);

  void Shear(const OutputVectorType &coef);

  bool GetInverse(Self* inverse) const;

  virtual InverseTransformBasePointer GetInverseTransform() const;

  inline InputPointType   BackTransform(const OutputPointType  &point ) const;
  inline InputVectorType  BackTransform(const OutputVectorType &vector) const;
  inline InputVnlVectorType BackTransform(
                                     const OutputVnlVectorType &vector) const;
  inline InputCovariantVectorType BackTransform(
                              const OutputCovariantVectorType &vector) const;

  inline InputPointType BackTransformPoint(const OutputPointType  &point) const;

  ScalarType Metric(const Self * other) const;

  ScalarType Metric(void) const;

protected:
  AffineTransform2D3D(const MatrixType &matrix,
                  const OutputVectorType &offset);
  AffineTransform2D3D(unsigned int outputDims,
                  unsigned int paramDims);
  AffineTransform2D3D();

  // implementation of the method in the parent class
  virtual void ComputeMatrixParameters(void);

  // implementation of the method in the parent class
  virtual void ComputeMatrix(void);
  
  virtual ~AffineTransform2D3D();

  void PrintSelf(std::ostream &s, Indent indent) const;

  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Rotations;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Scales;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Shears; 
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Translations;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_TranslateToCentre;  
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_BackTranslateCentre;

private:

  AffineTransform2D3D(const Self & other);
  const Self & operator=( const Self & );

}; //class AffineTransform

template<class TScalarType, unsigned int NDimensions>
inline
typename AffineTransform2D3D<TScalarType, NDimensions>::InputVectorType
AffineTransform2D3D<TScalarType, NDimensions>::
BackTransform(const OutputVectorType &vect ) const 
{
  return this->GetInverseMatrix() * vect;
}


template<class TScalarType, unsigned int NDimensions>
inline
typename AffineTransform2D3D<TScalarType, NDimensions>::InputVnlVectorType
AffineTransform2D3D<TScalarType, NDimensions>::
BackTransform(const OutputVnlVectorType &vect ) const 
{
  return this->GetInverseMatrix() * vect;
}


template<class TScalarType, unsigned int NDimensions>
inline
typename AffineTransform2D3D<TScalarType, NDimensions>::InputCovariantVectorType
AffineTransform2D3D<TScalarType, NDimensions>::
BackTransform(const OutputCovariantVectorType &vec) const 
{
  InputCovariantVectorType result;    // Converted vector

  for (unsigned int i = 0; i < NDimensions; i++) 
    {
    result[i] = NumericTraits<ScalarType>::Zero;
    for (unsigned int j = 0; j < NDimensions; j++) 
      {
      result[i] += this->GetMatrix()[j][i]*vec[j]; // Direct matrix transposed
      }
    }
  return result;
}


template<class TScalarType, unsigned int NDimensions>
inline
typename AffineTransform2D3D<TScalarType, NDimensions>::InputPointType
AffineTransform2D3D<TScalarType, NDimensions>::
BackTransformPoint(const OutputPointType &point) const
{
  return this->BackTransform(point);
}

template<class TScalarType, unsigned int NDimensions>
inline
typename AffineTransform2D3D<TScalarType, NDimensions>::InputPointType
AffineTransform2D3D<TScalarType, NDimensions>::
BackTransform(const OutputPointType &point) const 
{
  InputPointType result;       // Converted point
  ScalarType temp[NDimensions];
  unsigned int i, j;

  for (j = 0; j < NDimensions; j++) 
    {
    temp[j] = point[j] - this->GetOffset()[j];
    }

  for (i = 0; i < NDimensions; i++) 
    {
    result[i] = 0.0;
    for (j = 0; j < NDimensions; j++) 
      {
      result[i] += this->GetInverseMatrix()[i][j]*temp[j];
      }
    }
  return result;
}

}  // namespace itk

// Define instantiation macro for this template.
#define ITK_TEMPLATE_AffineTransform2D3D(_, EXPORT, x, y) namespace itk { \
  _(2(class EXPORT AffineTransform2D3D< ITK_TEMPLATE_2 x >)) \
  namespace Templates { typedef AffineTransform2D3D< ITK_TEMPLATE_2 x > \
                                                  AffineTransform2D3D##y; } \
  }

#if ITK_TEMPLATE_EXPLICIT
# include "Templates/itkAffineTransform2D3D+-.h"
#endif

#if ITK_TEMPLATE_TXX
# include "itkAffineTransform2D3D.txx"
#endif

#endif /* __itkAffineTransform2D3D_h */