itkCumulativeGaussianCostFunction.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright Insight Software Consortium
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         http://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
#ifndef __itkCumulativeGaussianCostFunction_h
#define __itkCumulativeGaussianCostFunction_h

#include "itkMultipleValuedCostFunction.h"

namespace itk
{
class CumulativeGaussianCostFunction:public MultipleValuedCostFunction
{
public:

  typedef CumulativeGaussianCostFunction Self;
  typedef MultipleValuedCostFunction     Superclass;
  typedef SmartPointer< Self >           Pointer;
  typedef SmartPointer< const Self >     ConstPointer;

  itkTypeMacro(CumulativeGaussianCostFunction, MultipleValuedCostFunction);

  itkNewMacro(Self);

  typedef Superclass::ParametersType ParametersType;
  typedef Superclass::MeasureType    MeasureType;
  typedef Superclass::DerivativeType DerivativeType;

  enum { SpaceDimension = 4 };

  void GetDerivative( const ParametersType & itkNotUsed(parameters),
                      DerivativeType & itkNotUsed(derivative) ) const {}

  MeasureType GetValue(const ParametersType & parameters) const;

  MeasureType * GetValuePointer(ParametersType & parameters);

  double CalculateFitError(MeasureType *setTestArray);

  double EvaluateCumulativeGaussian(double argument) const;

  unsigned int GetNumberOfParameters() const;

  unsigned int GetNumberOfValues() const;

  void Initialize(unsigned int rangeDimension);

  void SetOriginalDataArray(MeasureType *setOriginalDataArray);

protected:
  CumulativeGaussianCostFunction();
  virtual ~CumulativeGaussianCostFunction();

  void PrintSelf(std::ostream & os, Indent indent) const;

private:

  MeasureType *m_OriginalDataArray;

  unsigned int m_RangeDimension;

  mutable MeasureType    m_Measure;
  mutable MeasureType *  m_MeasurePointer;
  mutable ParametersType m_Parameters;
};
} // end namespace itk

#endif