freTransformBasedRegistrationProcessorInterface.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: freTransformBasedRegistrationProcessorInterface.h,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 __freTransformBasedRegistrationProcessorInterface_txx
#define __freTransformBasedRegistrationProcessorInterface_txx

#include "freTransformBasedRegistrationProcessorInterface.h"

namespace FREE
{


template <class TTransformBase,unsigned int VImageDimension>
void
TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::
SetTransform(TransformType* pTransform)
{
  if (pTransform != m_Transform.GetPointer())
  {
    m_Transform = pTransform;
  }
};

template <class TTransformBase,unsigned int VImageDimension>
TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::
TransformBasedRegistrationProcessorInterface()
{
    m_Transform = NULL;
}

template <class TTransformBase,unsigned int VImageDimension>
typename TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::TransformFieldPointer
TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::
ComputeTransformationField(const PointType& origin, const RegionType& region, const SpacingType& spacing)
{
    NotifyProgress(RSTFinalizing, "Computing transformation field");

    TransformFunctionMediaPointer transformFunction = GetTransformationFunction(origin,region,spacing);

    TransformFieldPointer smpField = transformFunction->GenerateField();

    return smpField;
};

template <class TTransformBase,unsigned int VImageDimension>
typename TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::TransformFunctionMediaPointer
TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::
GetTransformationFunction(const PointType& origin, const RegionType& region, const SpacingType& spacing)
{

        if (!TransformationIsValid()) throwExceptionMacro("Cannot return transformation function before registration is processed.");
        if (m_Transform.IsNull()) throwExceptionMacro("Cannot return transformation function when transform is not set.");

  NotifyProgress(RSTFinalizing, "Computing transformation function");

  TransformPointer newTransform = dynamic_cast<TransformType*>((m_Transform->CreateAnother()).GetPointer());
  if (newTransform.IsNull()) throwExceptionMacro("Transform is not based on the TransformType. Class of transform: "<<m_Transform->GetNameOfClass());
  
  try
  {
    newTransform->SetFixedParameters(m_Transform->GetFixedParameters());
  }
  catch(...)
  {
    //not every transform has fixed parameters and the default implementation of itk
    //throws an exception if you try to get or set exception in these cases.
    //catch the exception an go on, with just setting the normal parameters
  }

  newTransform->SetParametersByValue(m_CurrentTransformParameters);

  TransformFunctionMediaPointer transformFunction = TransformFunctionMediaType::New();
  transformFunction->SetTransform(newTransform);
  transformFunction->SetFieldOrigin(origin);
  transformFunction->SetFieldSpacing(spacing);
  transformFunction->SetFieldRegion(region);

  return transformFunction;
};

template <class TTransformBase,unsigned int VImageDimension>
typename TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::TransformFunctionMediaPointer
TransformBasedRegistrationProcessorInterface<TTransformBase, VImageDimension>::
GetInverseTransformationFunction(const PointType& origin, const RegionType& region, const SpacingType& spacing)
{
        if (!TransformationIsValid()) throwExceptionMacro("Cannot return transformation function before registration is processed.");
        if (m_Transform.IsNull()) throwExceptionMacro("Cannot return transformation function when transform is not set.");

  NotifyProgress(RSTFinalizing, "Computing transformation function");

  TransformPointer newTransform = dynamic_cast<TransformType*>((m_Transform->CreateAnother()).GetPointer());
  if (newTransform.IsNull()) throwExceptionMacro("Transform is not based on the TransformType. Class of transform: "<<m_Transform->GetNameOfClass());

    // get last actual transform parameter and set the transform component
  m_Transform->SetParametersByValue(m_CurrentTransformParameters);
  if (! m_Transform->GetInverse(newTransform)) throwExceptionMacro("Cannot generate inverse of the used transform. Transform: "<<m_Transform->GetNameOfClass());

  TransformFunctionMediaPointer transformFunction = TransformFunctionMediaType::New();
  transformFunction->SetTransform(newTransform);
  transformFunction->SetFieldOrigin(origin);
  transformFunction->SetFieldSpacing(spacing);
  transformFunction->SetFieldRegion(region);

  return transformFunction;
};

}//End of Namespace free

#endif

---