Method Compute

Modules/ModelFit/src/Functors/mitkModelFitFunctorBase.cpp:15–101

#include <mitkModelFitFunctorBase.h>

mitk::ModelFitFunctorBase::OutputPixelArrayType
mitk::ModelFitFunctorBase::
Compute(const InputPixelArrayType& value, const ModelBase* model,
        const ModelBase::ParametersType& initialParameters) const
{
  if (!model)
  {
    itkExceptionMacro("Cannot compute fit. Passed model is not defined.");
  }

  if (model->GetNumberOfParameters() != initialParameters.Size())
  {
    itkExceptionMacro("Cannot compute fit. Parameter count of passed model and passed initial parameters differ. Model parameter count: "
                      << model->GetNumberOfParameters() << "; Initial parameters: " << initialParameters);
  }

  SignalType sample(value.size());

  for (SignalType::SizeValueType i = 0; i < sample.Size(); ++i)
  {
    sample[i] = value [i];
  }

  DebugParameterMapType debugParams;
  ParameterNamesType debugNames;
  if (this->m_DebugParameterMaps)
  {
    debugNames = this->GetDebugParameterNames();
  }

  ParametersType fittedParameters = DoModelFit(sample, model, initialParameters, debugParams);

  OutputPixelArrayType derivedParameters = this->GetDerivedParameters(model, fittedParameters);

  OutputPixelArrayType criteria = this->GetCriteria(model, fittedParameters, sample);

  OutputPixelArrayType evaluationParameters = this->GetEvaluationParameters(model, fittedParameters,
      sample);

  if (criteria.size() != this->GetCriterionNames().size())
  {
    itkExceptionMacro("ModelFitInfo implementation seems to be inconsistent. Number of criterion values is not equal to number of criterion names.");
  }

  OutputPixelArrayType result(fittedParameters.Size() + derivedParameters.size() + criteria.size() +
                              evaluationParameters.size() + debugNames.size());

  for (ParametersType::SizeValueType i = 0; i < fittedParameters.Size(); ++i)
  {
    result[i] = fittedParameters[i];
  }

  OutputPixelArrayType::size_type offset = fittedParameters.Size();

  for (OutputPixelArrayType::size_type j = 0; j < derivedParameters.size(); ++j)
  {
    result[offset + j] = derivedParameters[j];
  }