itkFloodFilledFunctionConditionalConstIterator.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkFloodFilledFunctionConditionalConstIterator.h,v $
  Language:  C++
  Date:      $Date: 2009-02-05 19:04:56 $
  Version:   $Revision: 1.20 $

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

#include <queue>
#include <vector>

#include "itkIndex.h"
#include "itkSize.h"
#include "itkConditionalConstIterator.h"
#include "itkImage.h"

namespace itk
{

template<class TImage, class TFunction>
class ITK_EXPORT FloodFilledFunctionConditionalConstIterator:
    public ConditionalConstIterator<TImage>
{
public:

  typedef FloodFilledFunctionConditionalConstIterator Self;

  typedef TFunction FunctionType;

  typedef typename TFunction::InputType FunctionInputType;

  typedef typename TImage::IndexType  IndexType;

  typedef typename TImage::SizeType    SizeType;

  typedef typename TImage::RegionType    RegionType;

  typedef TImage   ImageType;

  typedef typename TImage::InternalPixelType   InternalPixelType;

  typedef typename TImage::PixelType   PixelType;

  itkStaticConstMacro(NDimensions, unsigned int, TImage::ImageDimension);

  FloodFilledFunctionConditionalConstIterator(const ImageType *imagePtr,
                                     FunctionType *fnPtr,
                                     IndexType startIndex);

  FloodFilledFunctionConditionalConstIterator(const ImageType *imagePtr,
                                     FunctionType *fnPtr,
                                     std::vector<IndexType> & startIndices);

  FloodFilledFunctionConditionalConstIterator(const ImageType *imagePtr,
                                              FunctionType *fnPtr);

  void FindSeedPixel();

  void FindSeedPixels();

  void InitializeIterator();

  virtual ~FloodFilledFunctionConditionalConstIterator() {};

  virtual bool IsPixelIncluded(const IndexType & index) const = 0;

  Self &operator=(const Self& it)
    {
    this->m_Image = it.m_Image;     // copy the smart pointer
    this->m_Region = it.m_Region;
    return *this;
    } 

  static unsigned int GetIteratorDimension() 
    {return TImage::ImageDimension;}

  const IndexType GetIndex()
    { return m_IndexStack.front();}

  const PixelType & Get(void) const
    { return this->m_Image->GetPixel(m_IndexStack.front() ); }

  bool IsAtEnd()
    { return this->m_IsAtEnd; }

  void AddSeed ( const IndexType seed )
    {
    m_StartIndices.push_back ( seed );
    }

  void ClearSeeds ()
    {
    m_StartIndices.clear();
    }

  void GoToBegin()
    {
    // Clear the queue
    while (!m_IndexStack.empty())
      {
      m_IndexStack.pop();
      }

    this->m_IsAtEnd = true;
    // Initialize the temporary image
    m_TemporaryPointer->FillBuffer(
      NumericTraits<ITK_TYPENAME TTempImage::PixelType>::Zero
      );
    
    for ( unsigned int i = 0; i < m_StartIndices.size(); i++ )
      {
      if( this->m_Image->GetBufferedRegion().IsInside ( m_StartIndices[i] ) &&
          this->IsPixelIncluded(m_StartIndices[i]) )
        {
        // Push the seed onto the queue
        m_IndexStack.push(m_StartIndices[i]);
        
        // Obviously, we're at the beginning
        this->m_IsAtEnd = false;
        
        // Mark the start index in the temp image as inside the
        // function, neighbor check incomplete
        m_TemporaryPointer->SetPixel(m_StartIndices[i], 2);
        }
      }
    }

  void operator++()
    { this->DoFloodStep(); }

  void DoFloodStep();
  
  virtual SmartPointer<FunctionType> GetFunction() const
    {
    return m_Function;
    }


protected: //made protected so other iterators can access 
  SmartPointer<FunctionType> m_Function;

  typedef Image<unsigned char, itkGetStaticConstMacro(NDimensions)> TTempImage;
  typename TTempImage::Pointer m_TemporaryPointer;

  std::vector<IndexType> m_StartIndices;

  typename ImageType::PointType m_ImageOrigin;

  typename ImageType::SpacingType m_ImageSpacing;

  RegionType   m_ImageRegion;

  std::queue<IndexType> m_IndexStack;

  FunctionInputType m_LocationVector;

  bool m_FoundUncheckedNeighbor;

  bool m_IsValidIndex;
};

} // end namespace itk

// Define instantiation macro for this template.
#define ITK_TEMPLATE_FloodFilledFunctionConditionalConstIterator(_, EXPORT, x, y) namespace itk { \
  _(2(class EXPORT FloodFilledFunctionConditionalConstIterator< ITK_TEMPLATE_2 x >)) \
  namespace Templates { typedef FloodFilledFunctionConditionalConstIterator< ITK_TEMPLATE_2 x > \
                        FloodFilledFunctionConditionalConstIterator##y; } \
  }

#if ITK_TEMPLATE_EXPLICIT
# include "Templates/itkFloodFilledFunctionConditionalConstIterator+-.h"
#endif

#if ITK_TEMPLATE_TXX
# include "itkFloodFilledFunctionConditionalConstIterator.txx"
#endif

#endif