itkAdvancedTranslationTransform.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkAdvancedTranslationTransform.h,v $
  Language:  C++
  Date:      $Date: 2007-07-15 16:38:25 $
  Version:   $Revision: 1.36 $

  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 __itkAdvancedTranslationTransform_h
#define __itkAdvancedTranslationTransform_h

#include <iostream>
#include "itkAdvancedTransform.h"
#include "itkExceptionObject.h"
#include "itkMatrix.h"




namespace itk
{

template <
    class TScalarType=double,          // Data type for scalars (float or double)
    unsigned int NDimensions=3>        // Number of dimensions
class ITK_EXPORT AdvancedTranslationTransform :
          public AdvancedTransform< TScalarType, NDimensions, NDimensions >
{
public:
  typedef AdvancedTranslationTransform Self;
  typedef AdvancedTransform< TScalarType, NDimensions, NDimensions > Superclass;
  typedef SmartPointer<Self>        Pointer;
  typedef SmartPointer<const Self>  ConstPointer;

  itkNewMacro( Self );

  itkTypeMacro( AdvancedTranslationTransform, AdvancedTransform );

  itkStaticConstMacro(SpaceDimension, unsigned int, NDimensions);
  itkStaticConstMacro(ParametersDimension, unsigned int, NDimensions);

  typedef typename Superclass::ScalarType ScalarType;

  typedef typename Superclass::ParametersType ParametersType;

  typedef typename Superclass::JacobianType JacobianType;

  typedef Vector<TScalarType, itkGetStaticConstMacro(SpaceDimension)> InputVectorType;
  typedef Vector<TScalarType, itkGetStaticConstMacro(SpaceDimension)> OutputVectorType;

  typedef CovariantVector<TScalarType, itkGetStaticConstMacro(SpaceDimension)> InputCovariantVectorType;
  typedef CovariantVector<TScalarType, itkGetStaticConstMacro(SpaceDimension)> OutputCovariantVectorType;

  typedef vnl_vector_fixed<TScalarType, itkGetStaticConstMacro(SpaceDimension)> InputVnlVectorType;
  typedef vnl_vector_fixed<TScalarType, itkGetStaticConstMacro(SpaceDimension)> OutputVnlVectorType;

  typedef Point<TScalarType, itkGetStaticConstMacro(SpaceDimension)> InputPointType;
  typedef Point<TScalarType, itkGetStaticConstMacro(SpaceDimension)> OutputPointType;

  typedef typename Superclass
    ::NonZeroJacobianIndicesType                    NonZeroJacobianIndicesType;
  typedef typename Superclass::SpatialJacobianType  SpatialJacobianType;
  typedef typename Superclass
    ::JacobianOfSpatialJacobianType                 JacobianOfSpatialJacobianType;
  typedef typename Superclass::SpatialHessianType   SpatialHessianType;
  typedef typename Superclass
    ::JacobianOfSpatialHessianType                  JacobianOfSpatialHessianType;
  typedef typename Superclass::InternalMatrixType   InternalMatrixType;

  virtual bool GetHasNonZeroSpatialHessian( void ) const
  {
    return false;
  }
  virtual bool HasNonZeroJacobianOfSpatialHessian( void ) const
  {
    return false;
  }

  const OutputVectorType & GetOffset(void) const
    { return m_Offset; }

  void SetParameters(const ParametersType & parameters);

  virtual const ParametersType& GetParameters(void) const;

  void SetOffset(const OutputVectorType &offset)
    { m_Offset = offset; return; }

  void Compose(const Self * other, bool pre=0);

  void Translate(const OutputVectorType &offset, bool pre=0);

  OutputPointType     TransformPoint(const InputPointType  &point ) const;
  OutputVectorType    TransformVector(const InputVectorType &vector) const;
  OutputVnlVectorType TransformVector(const InputVnlVectorType &vector) const;
  OutputCovariantVectorType TransformCovariantVector(
    const InputCovariantVectorType &vector) 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;

  bool GetInverse(Self* inverse) const;

  virtual const JacobianType & GetJacobian(const InputPointType  &point ) const;

  virtual void GetJacobian(
    const InputPointType &,
    JacobianType &,
    NonZeroJacobianIndicesType & ) const;

  virtual void GetSpatialJacobian(
    const InputPointType &,
    SpatialJacobianType & ) const;

  virtual void GetSpatialHessian(
    const InputPointType &,
    SpatialHessianType & ) const;

  virtual void GetJacobianOfSpatialJacobian(
    const InputPointType &,
    JacobianOfSpatialJacobianType &,
    NonZeroJacobianIndicesType & ) const;

  virtual void GetJacobianOfSpatialJacobian(
    const InputPointType &,
    SpatialJacobianType &,
    JacobianOfSpatialJacobianType &,
    NonZeroJacobianIndicesType & ) const;

  virtual void GetJacobianOfSpatialHessian(
    const InputPointType &,
    JacobianOfSpatialHessianType &,
    NonZeroJacobianIndicesType & ) const;

  virtual void GetJacobianOfSpatialHessian(
    const InputPointType & ipp,
    SpatialHessianType & sh,
    JacobianOfSpatialHessianType & jsh,
    NonZeroJacobianIndicesType & nonZeroJacobianIndices ) const;

  void SetIdentity(void);

  virtual unsigned int GetNumberOfParameters(void) const
                      { return NDimensions; }

  virtual bool IsLinear() const { return true; }

  virtual void SetFixedParameters( const ParametersType & )
    { /* purposely blank */ };

  virtual const ParametersType& GetFixedParameters(void) const
    {
    this->m_FixedParameters.SetSize(0);
    return this->m_FixedParameters;
    };

protected:
  AdvancedTranslationTransform();
  ~AdvancedTranslationTransform();
  void PrintSelf(std::ostream &os, Indent indent) const;

private:
  AdvancedTranslationTransform(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

  OutputVectorType   m_Offset;       // Offset of the transformation

  SpatialJacobianType m_SpatialJacobian;
  SpatialHessianType m_SpatialHessian;
  NonZeroJacobianIndicesType m_NonZeroJacobianIndices;
  JacobianOfSpatialJacobianType m_JacobianOfSpatialJacobian;
  JacobianOfSpatialHessianType m_JacobianOfSpatialHessian;

}; //class AdvancedTranslationTransform


// Back transform a point
template<class TScalarType, unsigned int NDimensions>
inline
typename AdvancedTranslationTransform<TScalarType, NDimensions>::InputPointType
AdvancedTranslationTransform<TScalarType, NDimensions>::
BackTransform(const OutputPointType &point) const {
  return point - m_Offset;
}




// Back transform a vector
template<class TScalarType, unsigned int NDimensions>
inline
typename AdvancedTranslationTransform<TScalarType, NDimensions>::InputVectorType
AdvancedTranslationTransform<TScalarType, NDimensions>::
BackTransform(const OutputVectorType &vect ) const
{
  return  vect;
}




// Back transform a vnl_vector
template<class TScalarType, unsigned int NDimensions>
inline
typename AdvancedTranslationTransform<TScalarType, NDimensions>::InputVnlVectorType
AdvancedTranslationTransform<TScalarType, NDimensions>::
BackTransform(const OutputVnlVectorType &vect ) const
{
  return  vect;
}


// Back Transform a CovariantVector
template<class TScalarType, unsigned int NDimensions>
inline
typename AdvancedTranslationTransform<TScalarType, NDimensions>::InputCovariantVectorType
AdvancedTranslationTransform<TScalarType, NDimensions>::
BackTransform(const OutputCovariantVectorType &vect) const
{
  return vect;
}

}  // namespace itk

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

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

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

#endif /* __itkAdvancedTranslationTransform_h */