Main Page   Groups   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Namespace Members   Compound Members   File Members   Concepts

itkCumulativeGaussianOptimizer.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkCumulativeGaussianOptimizer.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-06-24 12:02:50 $
00007   Version:   $Revision: 1.7 $
00008 
00009   Copyright (c) Insight Software Consortium. All rights reserved.
00010   See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
00011 
00012      This software is distributed WITHOUT ANY WARRANTY; without even 
00013      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
00014      PURPOSE.  See the above copyright notices for more information.
00015 
00016 =========================================================================*/
00017 
00018 #ifndef __itkCumulativeGaussianOptimizer_h
00019 #define __itkCumulativeGaussianOptimizer_h
00020 
00021 #include "itkMultipleValuedNonLinearOptimizer.h"
00022 #include "itkCumulativeGaussianCostFunction.h"
00023 
00024 namespace itk
00025 {
00026 
00047 class ITK_EXPORT CumulativeGaussianOptimizer : 
00048     public MultipleValuedNonLinearOptimizer
00049 {
00050 
00051 public:
00052 
00054   typedef CumulativeGaussianOptimizer           Self;
00055   typedef MultipleValuedNonLinearOptimizer      Superclass;
00056   typedef SmartPointer<Self>                    Pointer;
00057   typedef SmartPointer<const Self>              ConstPointer;
00058 
00060   typedef CumulativeGaussianCostFunction CostFunctionType;
00061 
00063   typedef CostFunctionType::MeasureType MeasureType;
00064 
00066   itkNewMacro(Self);
00067 
00069   itkTypeMacro(CumulativeGaussianOptimizer, MultipleValuedNonLinearOptimizer);
00070 
00072   itkSetMacro(DifferenceTolerance, double);
00073   itkSetMacro(Verbose, bool);
00074   itkGetMacro(ComputedMean, double);
00075   itkGetMacro(ComputedStandardDeviation, double);
00076   itkGetMacro(UpperAsymptote, double);
00077   itkGetMacro(LowerAsymptote, double);
00078   itkGetMacro(FinalSampledArray, MeasureType*);
00079   itkGetMacro(FitError, double);
00081 
00082   void SetDataArray(MeasureType * dataArray);
00083 
00085   void StartOptimization();
00086 
00088   void PrintArray(MeasureType * array);
00089 
00091   const std::string GetStopConditionDescription() const;
00092 
00093 protected:
00094   CumulativeGaussianOptimizer();
00095   virtual ~CumulativeGaussianOptimizer();
00096   void PrintSelf(std::ostream &os, Indent indent) const;
00097 
00098 private:
00099 
00101   double  m_DifferenceTolerance;
00102 
00104   double  m_ComputedMean; 
00105 
00107   double  m_ComputedStandardDeviation; 
00108 
00110   double  m_ComputedAmplitude;  
00111 
00114   double  m_ComputedTransitionHeight;
00115 
00117   double  m_UpperAsymptote;
00118 
00120   double  m_LowerAsymptote;
00121 
00123   double  m_OffsetForMean;
00124 
00126   bool m_Verbose;
00127 
00129   double m_FitError;
00130 
00133   MeasureType * m_FinalSampledArray;
00134 
00136   MeasureType * m_CumulativeGaussianArray;
00137 
00139   MeasureType * ExtendGaussian(MeasureType * originalArray, MeasureType * extendedArray, int startingPointForInsertion);
00140 
00142   MeasureType * RecalculateExtendedArrayFromGaussianParameters(MeasureType * originalArray,
00143                                                                MeasureType * extendedArray,
00144                                                                int startingPointForInsertion);
00145 
00148   double FindAverageSumOfSquaredDifferences(MeasureType * array1, MeasureType * array2);
00149 
00151   void FindParametersOfGaussian(MeasureType * sampledGaussianArray);
00152 
00154   void MeasureGaussianParameters(MeasureType * array);
00155 
00157   void PrintComputedParameterHeader();
00158 
00160   void PrintComputedParameters();
00161 
00163   double VerticalBestShift(MeasureType * originalArray, MeasureType * newArray);
00164 
00166   OStringStream m_StopConditionDescription;
00167 
00168 };
00169 
00170 } // end namespace itk
00171 
00172 #endif
00173 

Generated at Mon Jul 12 2010 18:10:36 for ITK by doxygen 1.7.1 written by Dimitri van Heesch, © 1997-2000