itkScalarChanAndVeseSparseLevelSetImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkScalarChanAndVeseSparseLevelSetImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2009-05-16 12:35:12 $
  Version:   $Revision: 1.4 $

  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"

namespace itk
{
template < class TInputImage, class TFeatureImage, class TOutputImage, class TFunction,
  class TSharedData, typename TIdCell = unsigned int >
class ITK_EXPORT ScalarChanAndVeseSparseLevelSetImageFilter :
public MultiphaseSparseFiniteDifferenceImageFilter< TInputImage, TOutputImage, TFunction, TIdCell >
{
public:
  typedef ScalarChanAndVeseSparseLevelSetImageFilter      Self;
  typedef MultiphaseSparseFiniteDifferenceImageFilter<
    TInputImage, 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 TFeatureImage                                 FeatureImageType;
  typedef typename FeatureImageType::Pointer            FeatureImagePtr;
  typedef typename FeatureImageType::PixelType          FeaturePixelType;
  typedef typename FeatureImageType::IndexType          FeatureIndexType;
  typedef typename FeatureIndexType::IndexValueType     FeatureIndexValueType;
  typedef typename FeatureImageType::RegionType         FeatureRegionType;

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

  typedef typename Superclass::OutputImageType          OutputImageType;
  typedef typename OutputImageType::ValueType           ValueType;
  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 FeatureImageType *f)
    {
    this->ProcessObject::SetNthInput( 0, const_cast< FeatureImageType *>(f) );
    }

  virtual const FeatureImageType * GetFeatureImage() const
    {
    return (static_cast< const FeatureImageType*>(this->ProcessObject::GetInput(0)));
    }

protected:
  ScalarChanAndVeseSparseLevelSetImageFilter()
    {
    this->SetNumberOfLayers(5); // 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< OutputImageType > &iterator,
    ValueType &newValue, bool &status );
};

} //end namespace itk

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

#endif

---