itkScalarChanAndVeseSparseLevelSetImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkScalarChanAndVeseSparseLevelSetImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2010-02-22 23:40:23 $
  Version:   $Revision: 1.11 $

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

#include "itkMultiphaseSparseFiniteDifferenceImageFilter.h"
#include "itkRegionOfInterestImageFilter.h"
#include "itkScalarChanAndVeseLevelSetFunction.h"

namespace itk
{
template < class TInputImage, class TFeatureImage, class TOutputImage,
     class TFunction=ScalarChanAndVeseLevelSetFunction< TInputImage, TFeatureImage >,
     class TSharedData=ITK_TYPENAME TFunction::SharedDataType,
     typename TIdCell = unsigned int >
class ITK_EXPORT ScalarChanAndVeseSparseLevelSetImageFilter :
public MultiphaseSparseFiniteDifferenceImageFilter< TInputImage, TFeatureImage,
  TOutputImage, TFunction, TIdCell >
{
public:
  typedef ScalarChanAndVeseSparseLevelSetImageFilter      Self;
  typedef MultiphaseSparseFiniteDifferenceImageFilter< TInputImage,
    TFeatureImage, TOutputImage, TFunction, TIdCell >     Superclass;
  typedef SmartPointer<Self>                              Pointer;
  typedef SmartPointer<const Self>                        ConstPointer;

  itkNewMacro( Self );

  itkTypeMacro( ScalarChanAndVeseSparseLevelSetImageFilter,
    MultiphaseSparseFiniteDifferenceImageFilter );

  itkStaticConstMacro( ImageDimension, unsigned int, TInputImage::ImageDimension );

  typedef typename Superclass::InputImageType           InputImageType;
  typedef typename Superclass::InputImagePointer        InputImagePointer;
  typedef typename Superclass::InputPointType           InputPointType;
  typedef typename Superclass::ValueType                ValueType;
  typedef typename Superclass::InputSpacingType         InputSpacingType;

  typedef TFeatureImage                                 FeatureImageType;
  typedef typename FeatureImageType::Pointer            FeatureImagePointer;
  typedef typename FeatureImageType::PixelType          FeaturePixelType;
  typedef typename FeatureImageType::IndexType          FeatureIndexType;
  typedef typename FeatureIndexType::IndexValueType     FeatureIndexValueType;
  typedef typename FeatureImageType::RegionType         FeatureRegionType;

  typedef TOutputImage                                  OutputImageType;
  typedef typename OutputImageType::IndexType           IndexType;
  typedef typename OutputImageType::PixelType           OutputPixelType;

  typedef typename Superclass::TimeStepType             TimeStepType;
  typedef typename Superclass::FiniteDifferenceFunctionType
                                                        FiniteDifferenceFunctionType;

  typedef typename Superclass::IdCellType               IdCellType;

  typedef TFunction                                     FunctionType;
  typedef typename FunctionType::Pointer                FunctionPtr;

  typedef TSharedData                                   SharedDataType;
  typedef typename SharedDataType::Pointer              SharedDataPointer;

  typedef RegionOfInterestImageFilter<
    FeatureImageType, FeatureImageType >                ROIFilterType;
  typedef typename ROIFilterType::Pointer               ROIFilterPointer;

#ifdef ITK_USE_CONCEPT_CHECKING

  itkConceptMacro(OutputHasNumericTraitsCheck,
    (Concept::HasNumericTraits<OutputPixelType>) );

#endif

  virtual void SetFeatureImage( const FeatureImagePointer f )
    {
    this->SetInput( f );
    }

protected:
  ScalarChanAndVeseSparseLevelSetImageFilter()
    {
    this->SetNumberOfLayers(2); // Narrow-band usage
    this->m_SharedData = SharedDataType::New();
    }

  ~ScalarChanAndVeseSparseLevelSetImageFilter() {}

  SharedDataPointer     m_SharedData;

  virtual void Initialize();
  virtual void InitializeIteration();
  virtual void UpdatePixel( unsigned int functionIndex,
    unsigned int idx, NeighborhoodIterator< InputImageType > &iterator,
    ValueType &newValue, bool &status );
};

} //end namespace itk

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

#endif