itkExtensionVelocitiesImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkExtensionVelocitiesImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2008-12-17 18:52:03 $
  Version:   $Revision: 1.26 $

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

#include "itkLevelSetVelocityNeighborhoodExtractor.h"
#include "itkFastMarchingExtensionImageFilter.h"
#include "itkReinitializeLevelSetImageFilter.h"

namespace itk
{
template <
  class TLevelSet,
  class TAuxValue = float,
  unsigned int VAuxDimension = 1
  >
class ITK_EXPORT ExtensionVelocitiesImageFilter :
    public ReinitializeLevelSetImageFilter<TLevelSet>
{
public:

  typedef ExtensionVelocitiesImageFilter             Self;
  typedef ReinitializeLevelSetImageFilter<TLevelSet> Superclass;
  typedef SmartPointer<Self>                         Pointer;
  typedef SmartPointer<const Self>                   ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(ExtensionVelocitiesImageFilter, ReinitializeLevelSetImageFilter);

  typedef LevelSetTypeDefault<TLevelSet>               LevelSetType;
  typedef typename LevelSetType::LevelSetPointer       LevelSetPointer;
  typedef typename LevelSetType::LevelSetConstPointer  LevelSetConstPointer;
  typedef typename LevelSetType::PixelType             PixelType;
  typedef typename LevelSetType::NodeType              NodeType;
  typedef typename LevelSetType::NodeContainer         NodeContainer;
  typedef typename LevelSetType::NodeContainerPointer  NodeContainerPointer;

  itkStaticConstMacro(SetDimension, unsigned int,LevelSetType::SetDimension);

  typedef AuxVarTypeDefault<TAuxValue,VAuxDimension,
                            itkGetStaticConstMacro(SetDimension)>
                                                     AuxVarType;
  typedef typename AuxVarType::AuxValueType          AuxValueType;
  typedef typename AuxVarType::AuxValueVectorType    AuxValueVectorType;
  typedef typename AuxVarType::AuxValueContainer     AuxValueContainer;
  typedef typename AuxVarType::AuxImageType          AuxImageType;
  typedef typename AuxVarType::AuxImagePointer       AuxImagePointer;
  typedef typename AuxVarType::AuxImageConstPointer  AuxImageConstPointer;

  itkStaticConstMacro(AuxDimension, unsigned int,VAuxDimension);

  void SetInputVelocityImage(const AuxImageType * ptr, unsigned int idx = 0);
  const AuxImageType * GetInputVelocityImage(unsigned int idx = 0);

  AuxImageType * GetOutputVelocityImage( unsigned int idx = 0 );

#ifdef ITK_USE_CONCEPT_CHECKING

  itkConceptMacro(AuxValueHasNumericTraitsCheck,
    (Concept::HasNumericTraits<TAuxValue>));
  itkConceptMacro(LevelSetOStreamWritableCheck,
    (Concept::OStreamWritable<PixelType>));

#endif

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

  virtual void GenerateDataFull();
  virtual void GenerateDataNarrowBand();
  virtual void AllocateOutput();

  virtual void EnlargeOutputRequestedRegion( DataObject * );

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

  typedef Image<float, itkGetStaticConstMacro(SetDimension) > SpeedImageType;

  typedef LevelSetVelocityNeighborhoodExtractor<TLevelSet,TAuxValue,VAuxDimension> 
                                                  LocatorType;
  typedef FastMarchingExtensionImageFilter<TLevelSet,TAuxValue,VAuxDimension,SpeedImageType> 
                                                  FastMarchingImageFilterType;

  typename LocatorType::Pointer                    m_Locator;
  typename FastMarchingImageFilterType::Pointer    m_Marcher;

};

} // namespace itk

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

#endif