itkMultiphaseDenseFiniteDifferenceImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkMultiphaseDenseFiniteDifferenceImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2009-07-29 15:13:55 $
  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 __itkMultiphaseDenseFiniteDifferenceImageFilter_h
#define __itkMultiphaseDenseFiniteDifferenceImageFilter_h

#include "itkMultiphaseFiniteDifferenceImageFilter.h"
#include "itkBinaryThresholdImageFilter.h"
#include "itkSignedMaurerDistanceMapImageFilter.h"
#include "itkImageRegionConstIterator.h"
#include "itkImageRegionIterator.h"
#include "itkNumericTraits.h"
#include "itkNeighborhoodAlgorithm.h"

#include <list>


namespace itk {

template < class TInputImage, class TFeatureImage, class TOutputImage, class TFunction,
  typename TIdCell = unsigned int >
class ITK_EXPORT MultiphaseDenseFiniteDifferenceImageFilter :
  public MultiphaseFiniteDifferenceImageFilter< TInputImage, TFeatureImage,
  TOutputImage, TFunction, TIdCell >
{
public:

  typedef MultiphaseDenseFiniteDifferenceImageFilter        Self;
  typedef MultiphaseFiniteDifferenceImageFilter< TInputImage,
    TFeatureImage, TOutputImage, TFunction >                Superclass;
  typedef SmartPointer<Self>                                Pointer;
  typedef SmartPointer<const Self>                          ConstPointer;

  itkTypeMacro( MultiphaseDenseFiniteDifferenceImageFilter, ImageToImageFilter );

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

  typedef typename Superclass::InputImageType           InputImageType;
  typedef typename Superclass::InputImagePointer        InputImagePointer;
  typedef typename Superclass::InputRegionType          InputRegionType;
  typedef typename Superclass::InputSizeType            InputSizeType;
  typedef typename Superclass::InputSizeValueType       InputSizeValueType;
  typedef typename Superclass::InputIndexType           InputIndexType;
  typedef typename Superclass::InputIndexValueType      InputIndexValueType;
  typedef typename Superclass::InputPixelType           InputPixelType;
  typedef typename Superclass::InputPointType           InputPointType;
  typedef typename Superclass::InputSpacingType         InputSpacingType;
  typedef typename InputImageType::ValueType            ValueType;

  typedef typename Superclass::FeatureImageType         FeatureImageType;
  typedef typename Superclass::FeatureSizeType          FeatureSizeType;
  typedef typename Superclass::FeatureImagePointer      FeatureImagePointer;
  typedef typename Superclass::FeatureRegionType        FeatureRegionType;
  typedef typename Superclass::FeatureSpacingType       FeatureSpacingType;
  typedef typename Superclass::FeaturePointType         FeaturePointType;
  typedef typename Superclass::FeaturePixelType         FeaturePixelType;

  typedef typename Superclass::OutputImageType          OutputImageType;
  typedef typename Superclass::OutputImagePointer       OutputImagePointer;
  typedef typename Superclass::OutputRegionType         OutputRegionType;
  typedef typename Superclass::OutputSizeType           OutputSizeType;
  typedef typename Superclass::OutputSizeValueType      SizeValueType;
  typedef typename Superclass::OutputIndexType          OutputIndexType;
  typedef typename Superclass::OutputIndexValueType     OutputIndexValueType;
  typedef typename OutputImageType::PixelType           OutputPixelType;

  typedef typename Superclass::IdCellType               IdCellType;

  typedef BinaryThresholdImageFilter< InputImageType, InputImageType >
                                                        ThresholdFilterType;
  typedef typename ThresholdFilterType::Pointer
                                                        ThresholdFilterPointer;

  typedef SignedMaurerDistanceMapImageFilter< InputImageType, InputImageType >
                                                        MaurerType;
  typedef typename MaurerType::Pointer                  MaurerPointer;

  typedef typename Superclass::FiniteDifferenceFunctionType
                                                        FiniteDifferenceFunctionType;
  typedef typename Superclass::FiniteDifferenceFunctionPointer
                                                        FiniteDifferenceFunctionPointer;
  typedef typename FiniteDifferenceFunctionType::NeighborhoodType
                                                        NeighborhoodIteratorType;

  typedef typename Superclass::TimeStepType TimeStepType;


  typedef NeighborhoodAlgorithm::ImageBoundaryFacesCalculator< InputImageType >
                                                        FaceCalculatorType;
  typedef typename FaceCalculatorType::FaceListType     FaceListType;

  void SetFunctionCount( const IdCellType& n );

#ifdef ITK_USE_CONCEPT_CHECKING

  itkConceptMacro(OutputTimesDoubleCheck,
    (Concept::MultiplyOperator<OutputPixelType, double>));
  itkConceptMacro(OutputAdditiveOperatorsCheck,
    (Concept::AdditiveOperators<OutputPixelType>));
  itkConceptMacro(InputConvertibleToOutputCheck,
    (Concept::Convertible<FeaturePixelType, OutputPixelType>));

#endif

protected:
  MultiphaseDenseFiniteDifferenceImageFilter()
    {
    this->m_ReinitializeCounter = 1;
    this->m_UpdateCounter = 0;        // FIXME: Should this be a bool ?
    }

  ~MultiphaseDenseFiniteDifferenceImageFilter() {}

  virtual void CopyInputToOutput();
  virtual void PostProcessOutput();

  virtual void ApplyUpdate(TimeStepType dt);

  unsigned int m_ReinitializeCounter;  // FIXME: Should this be a boolean ?
  unsigned int m_UpdateCounter;        // FIXME: Should this be a boolean ?

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

  virtual void AllocateUpdateBuffer();

  virtual TimeStepType CalculateChange();

  std::vector< InputImagePointer > m_UpdateBuffers;
};


}// end namespace itk

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

#endif