mitkMVModelFitCostFunction.cpp

Go to the documentation of this file.

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

The Medical Imaging Interaction Toolkit (MITK)

Copyright (c) German Cancer Research Center (DKFZ)
All rights reserved.

Use of this source code is governed by a 3-clause BSD license that can be
found in the LICENSE file.

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

#include "mitkMVModelFitCostFunction.h"

#include <iostream>


mitk::MVModelFitCostFunction::MeasureType mitk::MVModelFitCostFunction::GetValue(const ParametersType &parameter) const
{
  MeasureType measure;

  SignalType signal = m_Model->GetSignal(parameter);

  if(signal.GetSize() != m_Sample.GetSize()) itkExceptionMacro("Signal size does not matche sample size!");
  if(signal.GetSize() == 0)  itkExceptionMacro("Signal is empty!");

  measure = CalcMeasure(parameter, signal);

  return measure;
}

void mitk::MVModelFitCostFunction::GetDerivative (const ParametersType &parameters, DerivativeType &derivative) const
{
  ParametersType::SizeValueType paramCount = parameters.Size();
  MeasureType::SizeValueType measureCount = GetNumberOfValues();

  derivative.SetSize(paramCount,m_Sample.Size());

  for ( ParametersType::SizeValueType i = 0; i < paramCount; i++ )
  {
    ParametersType newParameters = parameters;
    newParameters[i] -= m_DerivativeStepLength;

    MeasureType e0 = GetValue(newParameters);

    newParameters = parameters;
    newParameters[i] += m_DerivativeStepLength;

    MeasureType e1 = GetValue(newParameters);

    for(MeasureType::SizeValueType j = 0; j<measureCount; ++j)
    {
      derivative[i][j] = (e1[j] - e0[j]) / ( 2 * m_DerivativeStepLength );
    }
  }


};

unsigned int mitk::MVModelFitCostFunction::GetNumberOfParameters() const
{
  return m_Model->GetNumberOfParameters();
}

unsigned int mitk::MVModelFitCostFunction::GetNumberOfValues() const{
  return m_Sample.GetSize();
}

void mitk::MVModelFitCostFunction::SetSample(const SignalType &sampleSet)
{
  itkDebugMacro("setting Sample to " << sampleSet);
  this->m_Sample = sampleSet;
  this->Modified();
}