itkDiffusionTensor3D.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkDiffusionTensor3D.h,v $
  Language:  C++
  Date:      $Date: 2010-04-09 21:20:50 $
  Version:   $Revision: 1.14 $

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

// Undefine an eventual DiffusionTensor3D macro
#ifdef DiffusionTensor3D
#undef DiffusionTensor3D
#endif

#include "itkSymmetricSecondRankTensor.h"


namespace itk
{

template < typename TComponent >
class DiffusionTensor3D: public SymmetricSecondRankTensor<TComponent,3>
{
public:
  typedef DiffusionTensor3D                           Self;
  typedef SymmetricSecondRankTensor<TComponent, 3>    Superclass;
  
  typedef typename Superclass::ValueType             ValueType;
  typedef typename Superclass::ComponentType         ComponentType;
#if defined(__GNUC__) && !defined(__INTEL_COMPILER) && (__GNUC__ == 3)
  typedef ComponentType ComponentArrayType[6];
#else
  typedef typename Superclass::ComponentArrayType    ComponentArrayType;
#endif
  typedef typename Superclass::AccumulateValueType   AccumulateValueType;
  typedef typename Superclass::RealValueType         RealValueType;

  typedef typename Superclass::EigenValuesArrayType    EigenValuesArrayType;
  typedef typename Superclass::EigenVectorsMatrixType  EigenVectorsMatrixType;

  DiffusionTensor3D();

  DiffusionTensor3D(const Superclass& r);
  DiffusionTensor3D(const ComponentType& r);
  DiffusionTensor3D(const ComponentArrayType r);
  
  template < typename TCoordRepB >
  DiffusionTensor3D( const DiffusionTensor3D<TCoordRepB> & pa )
    :SymmetricSecondRankTensor<TComponent,3>(pa) { }

  Self& operator= (const Superclass & r);
  Self& operator= (const ComponentType& r);
  Self& operator= (const ComponentArrayType r);

  template < typename TCoordRepB >
  Self& operator= ( const DiffusionTensor3D<TCoordRepB> & pa )
  {
    SymmetricSecondRankTensor<TComponent,3>::operator=(pa);
    return *this;
  }

  AccumulateValueType GetTrace() const;

  RealValueType GetFractionalAnisotropy() const;

  RealValueType GetRelativeAnisotropy() const;

  RealValueType GetInnerScalarProduct() const;

};

} // end namespace itk
#include "itkNumericTraitsDiffusionTensor3DPixel.h"

// Define instantiation macro for this template.
#define ITK_TEMPLATE_DiffusionTensor3D(_, EXPORT, x, y) namespace itk { \
  _(1(class EXPORT DiffusionTensor3D< ITK_TEMPLATE_1 x >)) \
  namespace Templates { typedef DiffusionTensor3D< ITK_TEMPLATE_1 x > \
                                            DiffusionTensor3D##y; } \
  }

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

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

#endif