itkRigidPCADeformationModelTransform.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 itkRigidPCADeformationModelTransform_h
#define itkRigidPCADeformationModelTransform_h

#include <iostream>
#include <itkImage.h>
#include <itkTransform.h>
#include <itkExceptionObject.h>
#include <itkMatrix.h>

#include <itkVectorLinearInterpolateNearestNeighborExtrapolateImageFunction.h>

namespace itk
{

template <
    class TScalarType=float, // Type for coordinate representation type (float or double)
        unsigned int NDimensions = 3 >  // Number of dimensions
        class ITK_EXPORT RigidPCADeformationModelTransform : public PCADeformationModelTransform< TScalarType, NDimensions>
{
public:
  typedef RigidPCADeformationModelTransform Self;
  typedef PCADeformationModelTransform< TScalarType, NDimensions >  Superclass;
  typedef SmartPointer<Self>        Pointer;
  typedef SmartPointer<const Self>  ConstPointer;
  
  typedef typename Superclass::VectorPixelType  VectorPixelType; 
  typedef typename Superclass::FieldType        FieldType;
  typedef typename Superclass::DisplacementType DisplacementType;

  typedef typename Superclass::FieldPointer        FieldPointer;
  typedef typename Superclass::FieldConstPointer   FieldConstPointer;
  typedef typename Superclass::FieldIterator       FieldIterator;
  typedef typename Superclass::FieldConstIterator  FieldConstIterator;

  typedef typename Superclass::FieldConstPointerArray FieldConstPointerArray;

  typedef typename Superclass::FieldPointerArray   FieldPointerArray;
  typedef typename Superclass::FieldIteratorArray  FieldIteratorArray;

  typedef typename Superclass::ContinuousIndexType ContinuousIndexType;

  typedef typename Superclass::FieldIndexType      FieldIndexType;
  typedef typename Superclass::FieldRegionType     FieldRegionType;
  typedef typename Superclass::FieldSizeType       FieldSizeType;

  typedef typename Superclass::FieldInterpolatorType         FieldInterpolatorType;
  typedef typename Superclass::FieldInterpolatorPointer      FieldInterpolatorPointer;
  typedef typename Superclass::FieldInterpolatorPointerArray FieldInterpolatorPointerArray;
  typedef itk::Matrix<double, NDimensions, NDimensions>      MatrixType;

  itkNewMacro( Self );

  itkTypeMacro( RigidPCADeformationModelTransform, PCADeformationModelTransform );

  virtual void SetNumberOfComponents(unsigned int numberOfComponents);

  virtual FieldPointer GetSingleDeformationField();

  itkStaticConstMacro(SpaceDimension, unsigned int, NDimensions);

  typedef typename Superclass::ScalarType  ScalarType;

  typedef typename Superclass::ParametersType  ParametersType;

  typedef typename Superclass::JacobianType  JacobianType;

  typedef typename Superclass::InputVectorType InputVectorType;
  typedef typename Superclass::OutputVectorType OutputVectorType;
  
  typedef typename Superclass::InputCovariantVectorType InputCovariantVectorType;
  typedef typename Superclass::OutputCovariantVectorType OutputCovariantVectorType;
  
  typedef typename Superclass::InputVnlVectorType InputVnlVectorType;
  typedef typename Superclass::OutputVnlVectorType OutputVnlVectorType;
  
  typedef typename Superclass::InputPointType InputPointType;
  typedef typename Superclass::OutputPointType OutputPointType;
  
  // Removed and to be filled in the future
  //virtual const JacobianType & GetJacobian( const InputPointType & point ) const;

  virtual OutputPointType TransformPoint(const InputPointType  &point ) const;

  virtual void Initialize() throw ( ExceptionObject );

  void SetParameters(const ParametersType & parameters);

  void SetCentre( InputPointType & );

  const InputPointType & GetCentre( void );

protected:
  RigidPCADeformationModelTransform();

  virtual ~RigidPCADeformationModelTransform();

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

  void ComputeMatrix(void);

  //Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Matrix;
  //Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Rotations;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_InPlateMatrix;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_RotationsZ;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_RollingMatrix;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_RotationsY;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_Translations;
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_TranslateToCentre;  
  Matrix<TScalarType,NDimensions+1,NDimensions+1> m_BackTranslateCentre;
 
  InputPointType m_centre;

private:
  RigidPCADeformationModelTransform(const Self & other); //purposely not implemented
  const Self & operator=( const Self & ); //purposely not implemented

}; //class RigidPCADeformationModelTransform

}  // namespace itk


#ifndef ITK_MANUAL_INSTANTIATION
#include "itkRigidPCADeformationModelTransform.txx"
#endif


#endif /* __itkRigidPCADeformationModelTransform_h */