itkCumulativeGaussianOptimizer.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkCumulativeGaussianOptimizer.h,v $
  Language:  C++
  Date:      $Date: 2009-06-24 12:02:50 $
  Version:   $Revision: 1.7 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/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 __itkCumulativeGaussianOptimizer_h
#define __itkCumulativeGaussianOptimizer_h

#include "itkMultipleValuedNonLinearOptimizer.h"
#include "itkCumulativeGaussianCostFunction.h"

namespace itk
{

class ITK_EXPORT CumulativeGaussianOptimizer : 
    public MultipleValuedNonLinearOptimizer
{

public:

  typedef CumulativeGaussianOptimizer           Self;
  typedef MultipleValuedNonLinearOptimizer      Superclass;
  typedef SmartPointer<Self>                    Pointer;
  typedef SmartPointer<const Self>              ConstPointer;
  
  typedef CumulativeGaussianCostFunction CostFunctionType;
  
  typedef CostFunctionType::MeasureType MeasureType;

  itkNewMacro(Self);

  itkTypeMacro(CumulativeGaussianOptimizer, MultipleValuedNonLinearOptimizer);

  itkSetMacro(DifferenceTolerance, double);
  itkSetMacro(Verbose, bool);
  itkGetMacro(ComputedMean, double);
  itkGetMacro(ComputedStandardDeviation, double);
  itkGetMacro(UpperAsymptote, double);
  itkGetMacro(LowerAsymptote, double);
  itkGetMacro(FinalSampledArray, MeasureType*);
  itkGetMacro(FitError, double);

  void SetDataArray(MeasureType * dataArray);

  void StartOptimization();
  
  void PrintArray(MeasureType * array);

  const std::string GetStopConditionDescription() const;

protected:
  CumulativeGaussianOptimizer();
  virtual ~CumulativeGaussianOptimizer();
  void PrintSelf(std::ostream &os, Indent indent) const;

private:

  double  m_DifferenceTolerance;
  
  double  m_ComputedMean; 
  
  double  m_ComputedStandardDeviation; 
  
  double  m_ComputedAmplitude;  
  
  double  m_ComputedTransitionHeight;
  
  double  m_UpperAsymptote;
  
  double  m_LowerAsymptote;
  
  double  m_OffsetForMean;

  bool m_Verbose;

  double m_FitError;

  MeasureType * m_FinalSampledArray;

  MeasureType * m_CumulativeGaussianArray;

  MeasureType * ExtendGaussian(MeasureType * originalArray, MeasureType * extendedArray, int startingPointForInsertion);
  
  MeasureType * RecalculateExtendedArrayFromGaussianParameters(MeasureType * originalArray,
                                                               MeasureType * extendedArray,
                                                               int startingPointForInsertion);

  double FindAverageSumOfSquaredDifferences(MeasureType * array1, MeasureType * array2);

  void FindParametersOfGaussian(MeasureType * sampledGaussianArray);

  void MeasureGaussianParameters(MeasureType * array);
  
  void PrintComputedParameterHeader();
  
  void PrintComputedParameters();
  
  double VerticalBestShift(MeasureType * originalArray, MeasureType * newArray);

  OStringStream m_StopConditionDescription;

};

} // end namespace itk

#endif