itkWarpTransform3D.h

Go to the documentation of this file.

#ifndef itkWarpTransform3D_h
#define itkWarpTransform3D_h
 
#include <itkTransform.h>
#include "dtiprocessFiles/dtitypes.h"
#include <itkImage.h>
#include "itkConstNeighborhoodIterator.h"
 
namespace itk
{
template <class FieldData>
class WarpTransform3D : public Transform<FieldData, 3, 3>
{
public:
  typedef FieldData FieldDataType;
  typedef WarpTransform3D Self;
  typedef Transform<FieldDataType, 3, 3> Superclass;
  typedef typename Superclass::JacobianType JacobianType;
  typedef typename Superclass::JacobianPositionType JacobianPositionType;
  typedef typename Superclass::InputPointType InputPointType;
  typedef typename Superclass::InputVectorType InputVectorType;
  typedef typename Superclass::InputVnlVectorType InputVnlVectorType;
  typedef typename Superclass::InputCovariantVectorType InputCovariantVectorType;
  typedef typename Superclass::OutputPointType OutputPointType;
  typedef typename Superclass::OutputVectorType OutputVectorType;
  typedef typename Superclass::OutputVnlVectorType OutputVnlVectorType;
  typedef typename Superclass::OutputCovariantVectorType OutputCovariantVectorType;
  typedef DeformationImageType::Pointer DeformationImagePointerType;
  typedef SmartPointer<Self> Pointer;
  typedef SmartPointer<const Self> ConstPointer;
  typedef ConstNeighborhoodIterator<DeformationImageType> ConstNeighborhoodIteratorType;
  typedef typename ConstNeighborhoodIteratorType::RadiusType RadiusType;
  typedef typename Superclass::ParametersType ParametersType;

  static Pointer New();

  ::itk::LightObject::Pointer CreateAnother() const override;

  itkTypeMacro(WarpTransform3D, Transform);
 
  OutputPointType TransformPoint(const InputPointType& inputPoint) const override;
 
  void ComputeJacobianWithRespectToParameters(const InputPointType& p, JacobianType& jacobian) const override;
 
  void ComputeJacobianWithRespectToPosition(const InputPointType& itkNotUsed(x), JacobianType& itkNotUsed(j)) const override
  {
    itkExceptionMacro("ComputeJacobianWithRespectToPosition is not implemented for WarpTransform3D");
  }
 
  void ComputeJacobianWithRespectToPosition(const InputPointType& itkNotUsed(x), JacobianPositionType& itkNotUsed(j)) const override
  {
    itkExceptionMacro("ComputeJacobianWithRespectToPosition is not implemented for WarpTransform3D");
  }
 
  void SetDeformationField(DeformationImagePointerType deformationField);
 
  const DeformationImagePointerType GetDeformationField() const { return m_DeformationField; }
 
  using Superclass::TransformVector;
  OutputVectorType TransformVector(const InputVectorType&) const override { itkExceptionMacro("TransformVector(const InputVectorType&) is not implemented for WarpTransform3D"); }

  OutputVnlVectorType TransformVector(const InputVnlVectorType&) const override
  {
    itkExceptionMacro("TransformVector(const InputVnlVectorType&) is not implemented for WarpTransform3D");
  }
 
  using Superclass::TransformCovariantVector;
  OutputCovariantVectorType TransformCovariantVector(const InputCovariantVectorType&) const override
  {
    itkExceptionMacro("TransformCovariantVector(const InputCovariantVectorType& is not implemented for WarpTransform3D");
  }
 
protected:
  itkGetConstReferenceMacro(NeighborhoodRadius, RadiusType);
  itkSetMacro(NeighborhoodRadius, RadiusType);

  void SetParameters(const ParametersType&) override {}
  void SetFixedParameters(const ParametersType&) override {}
 
  WarpTransform3D();
  void operator=(const Self&); // purposely not implemented
 
  RadiusType m_NeighborhoodRadius;
  double m_DerivativeWeights[3];
  DeformationImagePointerType m_DeformationField;
  //  Vector< double , 3 > m_OutputSpacing ;
  Size<3> m_SizeForJacobian;
};
 
} // end namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
# include "itkWarpTransform3D.txx"
#endif
 
#endif