freTransformFunctionMedia.txx

Go to the documentation of this file.

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

  Program:   F.R.E.E. - flexible registration evaluation engine
  Version:   v.1.0.0
  Date:      $Date: 2006/09/01 12:00:00 $
  Module:    $RCSfile: freTransformFunctionMedia.txx,v $
  Language:  C++



  Copyright (c) 2007 Ralf o Floca (Department of Medical Informatics,
  Institute for Medical Biometry and Informatics, University of Heidelberg,
  Germany). All rights reserved.
  See FREECopyright.txt or http://www.mi.med.uni-hd.de/free/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 __freTransformFunctionMedia_txx
#define __freTransformFunctionMedia_txx

#include "freTransformFunctionMedia.h"

namespace FREE
{



template <class TScalarType, unsigned int NInputDimensions, unsigned int NOutputDimensions>
void
TransformFunctionMedia<TScalarType, NInputDimensions, NOutputDimensions>::
SetFieldRegion(const RegionType &region)
{
  m_FieldRegion = region;
};

template <class TScalarType, unsigned int NInputDimensions, unsigned int NOutputDimensions>
void
TransformFunctionMedia<TScalarType, NInputDimensions, NOutputDimensions>::
Reset()
{
  m_Transform = 0;
  
  m_FieldSpacing.Fill(1.0);
  m_FieldOrigin.Fill(0.0);

  typename RegionType::IndexType index;
  index.Fill(0);
  typename RegionType::SizeType size;
  size.Fill(1);
  m_FieldRegion = RegionType(index, size);
};

template <class TScalarType, unsigned int NInputDimensions, unsigned int NOutputDimensions>
TransformFunctionMedia<TScalarType, NInputDimensions, NOutputDimensions>::
~TransformFunctionMedia()
{
  Reset();
};

template <class TScalarType, unsigned int NInputDimensions, unsigned int NOutputDimensions>
TransformFunctionMedia<TScalarType, NInputDimensions, NOutputDimensions>::
TransformFunctionMedia()
{
  Reset();
};

template <class TScalarType, unsigned int NInputDimensions, unsigned int NOutputDimensions>
typename TransformFunctionMedia<TScalarType, NInputDimensions, NOutputDimensions>::TransformFieldPointer
TransformFunctionMedia<TScalarType, NInputDimensions, NOutputDimensions>::
GenerateField()
{
  if (m_Transform.IsNull()) throwExceptionMacro("Error cannot generate field. Transform function has no transformation set.");

  TransformFieldPointer smpField = TransformFieldType::New();

  try
  {
    //Create field
    smpField->SetOrigin(m_FieldOrigin);
    smpField->SetSpacing(m_FieldSpacing);
    smpField->SetRegions(m_FieldRegion);
    smpField->Allocate();

    typename TransformFieldType::ValueType zeroVector;
    for (int i=0; i<TransformFieldType::ImageDimension; i++)
        zeroVector.SetNthComponent(i,0);
    smpField->FillBuffer(zeroVector);

    // get Iterator to roam the field
    typedef itk::ImageRegionIterator<TransformFieldType> IteratorType;
    IteratorType iterator(smpField,m_FieldRegion);

    //process the field
    typedef typename TransformFieldType::PointType PointType;
    PointType actPoint;
    for (iterator.GoToBegin(); !(iterator.IsAtEnd()); ++iterator)
    {
        smpField->TransformIndexToPhysicalPoint(iterator.GetIndex(),actPoint);
        PointType transPoint = m_Transform->TransformPoint(actPoint);
        iterator.Set(transPoint-actPoint);
    }
  }
  catchAllNPassMacro("Unknown error when computing transformation field.");

  return smpField;
}

} //end of namespace free

#endif

---