cmtkAnatomicalOrientation.h

Go to the documentation of this file.

/*
//
//  Copyright 1997-2009 Torsten Rohlfing
//
//  Copyright 2004-2011 SRI International
//
//  This file is part of the Computational Morphometry Toolkit.
//
//  http://www.nitrc.org/projects/cmtk/
//
//  The Computational Morphometry Toolkit is free software: you can
//  redistribute it and/or modify it under the terms of the GNU General Public
//  License as published by the Free Software Foundation, either version 3 of
//  the License, or (at your option) any later version.
//
//  The Computational Morphometry Toolkit is distributed in the hope that it
//  will be useful, but WITHOUT ANY WARRANTY; without even the implied
//  warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License along
//  with the Computational Morphometry Toolkit.  If not, see
//  <http://www.gnu.org/licenses/>.
//
//  $Revision: 2749 $
//
//  $LastChangedDate: 2011-01-14 14:44:48 -0800 (Fri, 14 Jan 2011) $
//
//  $LastChangedBy: torstenrohlfing $
//
*/

#ifndef __cmtkAnatomicalOrientation_h_included_
#define __cmtkAnatomicalOrientation_h_included_

#include <cmtkconfig.h>

#include <Base/cmtkAnatomicalOrientationBase.h>

#include <System/cmtkSmartPtr.h>

#include <Base/cmtkTypes.h>
#include <Base/cmtkFixedVector.h>
#include <Base/cmtkMetaInformationObject.h>
#include <Base/cmtkAffineXform.h>

namespace
cmtk
{


class AnatomicalOrientation
  : public AnatomicalOrientationBase
{
public:
  typedef AnatomicalOrientation Self;

  typedef AnatomicalOrientationBase Superclass;

  static void GetOrientationFromDirections( char* orientation, const AffineXform::MatrixType& directions, const char* spaceAxes );

  static void GetImageToSpaceAxesPermutation( int (&imageToSpaceAxesPermutation)[3][3], const char* orientation, const char* spaceAxes );

  class PermutationMatrix
  {
  public:
    typedef PermutationMatrix Self;

    typedef SmartPointer<Self> SmartPtr;

    PermutationMatrix( const FixedVector<3,int>& sourceDims, const std::string& curOrientation, const char newOrientation[3] );
   
    template<class T>
    const FixedVector<3,T> GetPermutedArray( const FixedVector<3,T>& source ) const
    {
      FixedVector<3,T> target;
      for ( int i = 0; i < 3; i++ )
        {
        target[i] = source[this->m_Axes[i]];
        }
      return target;
    }

    AffineXform::MatrixType GetPermutedMatrix( const AffineXform::MatrixType& inMatrix ) const;
    
    void GetReorientedIndex( const int* origPoint, int* newPoint ) const
    {
      for ( int i = 0; i < 3; ++i )
        newPoint[i] = this->m_Multipliers[i] * origPoint[this->m_Axes[i]] + this->m_Offsets[i];
    }

    size_t NewOffsetFromPoint( const int* origPoint ) const
    {
      return ( this->m_Multipliers[0] * origPoint[this->m_Axes[0]] + this->m_Offsets[0] ) + 
              this->m_NewDims[0] * 
              ( ( this->m_Multipliers[1] * origPoint[this->m_Axes[1]] + this->m_Offsets[1] ) + 
                  this->m_NewDims[1] * 
                  ( this->m_Multipliers[2] * origPoint[this->m_Axes[2]] + this->m_Offsets[2] ) );
    }

  private:
    FixedVector<3,int> m_Axes;

    FixedVector<3,int> m_Multipliers;

    FixedVector<3,int> m_NewDims;
    
    FixedVector<3,int> m_Offsets;
  };
};


} // namespace cmtk

#endif // #ifndef __cmtkAnatomicalOrientation_h_included_