freSVLimitedAxisVnlSOMetricAdaptor.cxx

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: freSVLimitedAxisVnlSOMetricAdaptor.cxx,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.

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

#include "freSVLimitedAxisVnlSOMetricAdaptor.h"

namespace FREE
{
  
void
SVLimitedAxisVnlSOMetricAdaptor 
::SetOriginPosition(const OriginPositionType& origin)
{
        m_OriginPosition = origin;
};

void
SVLimitedAxisVnlSOMetricAdaptor 
::SetParameterAxis(const ParameterAxisType& axis)
{
        m_ParameterAxis = axis;
};

SVLimitedAxisVnlSOMetricAdaptor 
::SVLimitedAxisVnlSOMetricAdaptor(unsigned int spaceDimension, unsigned int axisCount):
  SingleValuedVnlSOMetricAdaptor(spaceDimension), m_ParameterAxis(spaceDimension),
        m_OriginPosition(spaceDimension)
{
        m_Scales.SetSize(spaceDimension);
        m_Scales.Fill(1);
        m_ParameterAxis.Fill(-1);
        m_OriginPosition.Fill(0.0);
        m_NumberOfAxis = axisCount;
}    

SVLimitedAxisVnlSOMetricAdaptor::ParametersType
SVLimitedAxisVnlSOMetricAdaptor::
ConvertInternalToParameter(const InternalParametersType & inparameters)
{
        // default scaling for any parameter is 1
        ParametersType parameters(this->get_number_of_unknowns());
                for(unsigned int i=0;i<parameters.size();i++)
                {
                        int iAxisIndex = m_ParameterAxis[i];
                        if (iAxisIndex<0) iAxisIndex = i;

                        if (iAxisIndex>=inparameters.size()) throwExceptionMacro("Error: axis index exceeds internal parameter set. Index: " << iAxisIndex << "; size: " << inparameters.size());
                        
                        parameters[i] = m_OriginPosition[i] + (inparameters[iAxisIndex]*m_Scales[i]);
                }

        return parameters;
};

SVLimitedAxisVnlSOMetricAdaptor::InternalParametersType
SVLimitedAxisVnlSOMetricAdaptor::
ConvertParameterToInternal(const ParametersType & parameters)
{
        InternalParametersType inparameters(m_NumberOfAxis);
        inparameters.fill(0.0);

        for(unsigned int i=0;i<parameters.size();i++)
        {
                int iAxisIndex = m_ParameterAxis[i];
                if (iAxisIndex<0) iAxisIndex = i;

                if (iAxisIndex>=inparameters.size()) throwExceptionMacro("Error: axis index exceeds internal parameter set. Index: " << iAxisIndex << "; size: " << inparameters.size());
                        
                //It is possible that an inparameter value will be overwritten
                //if there are more then one parameter on one axis. Then the last one
                //is the one used, hence the difference between the computed versions
                //should be minor.
                inparameters[iAxisIndex] = (parameters[i] - m_OriginPosition[i])/m_Scales[i];
        }

        return inparameters;
};

} // end namespace itk

---