itkActiveShapeModelGradientSearchMethod.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkActiveShapeModelGradientSearchMethod.h,v $
  Language:  C++
  Date:      $Date: 2009-02-01 13:08:40 $
  Version:   $Revision: 1.2 $

  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 __itkActiveShapeModelGradientSearchMethod_h
#define __itkActiveShapeModelGradientSearchMethod_h

#include <time.h>
#include <math.h>

#include "vnl/vnl_vector.h"
#include "vnl/vnl_matrix.h"
#include "vnl/vnl_math.h"

#include "itkImageToImageFilter.h"
#include "itkImage.h"
#include "itkNumericTraits.h"


#include "itkImageRegionConstIterator.h"

#include "itkGradientMagnitudeImageFilter.h"


#include "itkVector.h"
#include "itkListSample.h"
#include "itkExceptionObject.h"

#include "itkConceptChecking.h"

#include <vector>
#include <list>

namespace itk {

template < class TImage >
class ITK_EXPORT ActiveShapeModelGradientSearchMethod : public Object
{
public:

  typedef ActiveShapeModelGradientSearchMethod<TImage>                   Self;
  typedef Object                                                         Superclass;
  typedef SmartPointer<Self>                                             Pointer;
  typedef SmartPointer<const Self>                                       ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(ActiveShapeModelGradientSearchMethod, Object);
  
  typedef TImage                                                         InputImageType;
  typedef Image< double, 2 >                                             Image2DType;

  typedef typename InputImageType::RegionType   RegionType;

  typedef typename RegionType::IndexType  IndexType;

  typedef typename InputImageType::PixelType PixelType;

  typedef typename RegionType::SizeType   SizeType;

  typedef typename InputImageType::ConstPointer InputImagePointer;

  typedef GradientMagnitudeImageFilter< InputImageType, Image2DType>      GradientFilterType;

  typedef Vector< unsigned int, 2 >                                       MeasurementVectorUIntType;
  typedef Vector< double, 2 >                                             MeasurementVectorDoubleType;
  typedef std::vector<unsigned int>                                       VectorType;
  typedef vnl_vector<double>                                              VectorOfDoubleType;
  typedef vnl_matrix<double>                                              MatrixOfDoubleType;

  typedef Statistics::ListSample< MeasurementVectorUIntType >             SampleType;

  typedef ImageRegionConstIterator< Image2DType  >                        ConstIteratorType;

  virtual void SetImage( const InputImageType * image )
    {
    if ( m_Image != image )
      {
      m_Image = image;
      this->Modified();
      m_Valid = false;
      }
    }

  itkSetMacro( LenghtOfProfile, unsigned int );
  itkGetMacro( LenghtOfProfile, unsigned int );

  itkSetMacro( NumberOfIteration, unsigned int );
  itkGetMacro( NumberOfIteration, unsigned int );

  itkSetMacro( MeanShape, VectorOfDoubleType );
  itkGetMacro( MeanShape, VectorOfDoubleType );

  itkSetMacro( EigenValues, VectorOfDoubleType );
  itkGetMacro( EigenValues, VectorOfDoubleType );

  itkSetMacro( EigenVectors, MatrixOfDoubleType );
  itkGetMacro( EigenVectors, MatrixOfDoubleType );


  void GenerateData( );

  VectorOfDoubleType  GetNewShape();

protected:
  ActiveShapeModelGradientSearchMethod();
 ~ActiveShapeModelGradientSearchMethod() {};

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


private:
  ActiveShapeModelGradientSearchMethod(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

  bool                  m_Valid;                      // Have parameters been computed yet?
  VectorOfDoubleType    m_MeanShape; 
  MatrixOfDoubleType    m_EigenVectors;
  VectorOfDoubleType    m_EigenValues;
  unsigned int          m_LenghtOfProfile;
  unsigned int          m_NumberOfIteration;
  VectorOfDoubleType    m_DiffVector;
  VectorOfDoubleType    m_Db;
  VectorOfDoubleType    m_Blimit;
  VectorOfDoubleType    m_NewShape;
  InputImagePointer     m_Image;


}; // class ActiveShapeModelGradientSearchMethod

} // end of namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkActiveShapeModelGradientSearchMethod.txx"
#endif

#endif