freAccuracySOMetricMonitor.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: freAccuracySOMetricMonitor.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 "freAccuracySOMetricMonitor.h"

#include "itkNumericTraits.h"

namespace FREE
{

  AccuracySOMetricResult::
  AccuracySOMetricResult()
  {
    this->Reset();
  };

  void
  AccuracySOMetricResult::
  Reset()
  {
                m_MinError = itk::NumericTraits< double >::max();
                m_MaxError = 0.0;
                m_Error = 0.0;
                m_Var = 0.0;

    m_UnevaluatedPoints = 0;

                m_Duration = 0;
    m_Samplesize = 0;

    m_PointErrors.clear();
  };

AccuracySOMetricMonitor::
AccuracySOMetricMonitor()
{};

AccuracySOMetricMonitor::
~AccuracySOMetricMonitor()
{
};

void
AccuracySOMetricMonitor::
InsertEvaluationResults(const ResultType& results)
{
  m_Errors.push_back(results.GetError());
  m_Vars.push_back(results.GetVariance());
  m_MinErrors.push_back(results.GetMinError());
  m_MaxErrors.push_back(results.GetMaxError());
  m_PointErrors.push_back(results.GetPointErrors());
  m_Samplesize.push_back(results.GetSamplesize());
  m_Durations.push_back(results.GetDuration());
  m_UnevaluatedPoints.push_back(results.GetUnevaluatedPoints());
};

void
AccuracySOMetricMonitor::
CopyAdaptationResults(const ResultIndexType& index, ResultType& results) const
{
  results.SetError(m_Errors[index]);
  results.SetVariance(m_Vars[index]);
  results.SetMinError(m_MinErrors[index]);
  results.SetMaxError(m_MaxErrors[index]);
  results.SetPointErrors(m_PointErrors[index]);
  results.SetSamplesize(m_Samplesize[index]);
  results.SetDuration(m_Durations[index]);
  results.SetUnevaluatedPoints(m_UnevaluatedPoints[index]);
};

void
AccuracySOMetricMonitor::
EnlistResultsToStatistic(const ResultType& results, StatisticEntry& entry)
{
  entry.AddValue(Convert::ToStr(results.GetError()),2);
  entry.AddValue(Convert::ToStr(sqrt(results.GetVariance())),3);
  entry.AddValue(Convert::ToStr(results.GetMinError()),4);
  entry.AddValue(Convert::ToStr(results.GetMaxError()),5);

  std::string sPointErrors ="";
  for (AccuracySOMetricResult::PointErrorListType::const_iterator pos = results.GetPointErrors().begin(); pos!=results.GetPointErrors().end(); pos++)
  {
    sPointErrors += Convert::ToStr(*pos) + "; ";
  }
  entry.AddValue(sPointErrors,6);

  entry.AddValue(Convert::ToStr(results.GetDuration()),7);
  entry.AddValue(Convert::ToStr(results.GetUnevaluatedPoints()),8);
};

void
AccuracySOMetricMonitor::
EnlistDictionaryEntries(StatisticDictionary& dictionary)
{
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Error", "Mean error; deviation in real world spacing");
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Std dev", "Standard deviation of the error");
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Minimum", "Minimum error");
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Maximum", "Maximum error");
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Point error", "Absolut error of each reference point. Only valid if metric uses points to evaluate accuracy.");
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Duration", "Duration of the registration in 1/10 sec");
  dictionary.AddValueDefinition(this->GetNameOfClass(),"Unevaluated points", "Number of points that were not evaluated. Only valid if accuracy is estimated via point references.");
};


} // end namespace FREE

---