itkShapedFloodFilledFunctionConditionalConstIterator.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright NumFOCUS
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         http://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
#ifndef itkShapedFloodFilledFunctionConditionalConstIterator_h
#define itkShapedFloodFilledFunctionConditionalConstIterator_h
 
#include <queue>
#include <vector>
 
#include "itkSize.h"
#include "itkConditionalConstIterator.h"
#include "itkConnectedComponentAlgorithm.h"
 
namespace itk
{
template <typename TImage, typename TFunction>
class ITK_TEMPLATE_EXPORT ShapedFloodFilledFunctionConditionalConstIterator : public ConditionalConstIterator<TImage>
{
public:
 
  using Self = ShapedFloodFilledFunctionConditionalConstIterator;
 
  using FunctionType = TFunction;
 
  using FunctionInputType = typename TFunction::InputType;
 
  using IndexType = typename TImage::IndexType;
 
  using SeedsContainerType = typename std::vector<IndexType>;
 
  using OffsetType = typename TImage::OffsetType;
 
  using SizeType = typename TImage::SizeType;
 
  using RegionType = typename TImage::RegionType;
 
  using ImageType = TImage;
 
  using InternalPixelType = typename TImage::InternalPixelType;
 
  using PixelType = typename TImage::PixelType;
 
  using NeighborhoodIteratorType = typename itk::ShapedNeighborhoodIterator<ImageType>;
 
  static constexpr unsigned int NDimensions = TImage::ImageDimension;
 
  ShapedFloodFilledFunctionConditionalConstIterator(const ImageType * imagePtr,
                                                    FunctionType *    fnPtr,
                                                    IndexType         startIndex);
 
  ShapedFloodFilledFunctionConditionalConstIterator(const ImageType *        imagePtr,
                                                    FunctionType *           fnPtr,
                                                    std::vector<IndexType> & startIndex);
 
  ShapedFloodFilledFunctionConditionalConstIterator(const ImageType * imagePtr, FunctionType * fnPtr);
 
  void
  FindSeedPixel();
 
  void
  FindSeedPixels();
 
  void
  InitializeIterator();
 
  ~ShapedFloodFilledFunctionConditionalConstIterator() override = default;
 
  bool
  IsPixelIncluded(const IndexType & index) const override = 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() override
  {
    return m_IndexStack.front();
  }
 
  const PixelType
  Get() const override
  {
    return this->m_Image->GetPixel(m_IndexStack.front());
  }
 
  bool
  IsAtEnd() const override
  {
    return this->m_IsAtEnd;
  }
 
  void
  AddSeed(const IndexType seed)
  {
    m_Seeds.push_back(seed);
  }
 
  void
  ClearSeeds()
  {
    m_Seeds.clear();
  }
 
  void
  GoToBegin()
  {
    // Clear the queue
    while (!m_IndexStack.empty())
    {
      m_IndexStack.pop();
    }
 
    this->m_IsAtEnd = true;
    // Initialize the temporary image
    m_TempPtr->FillBuffer(NumericTraits<typename TTempImage::PixelType>::ZeroValue());
 
    for (unsigned int i = 0; i < m_Seeds.size(); i++)
    {
      if (this->m_Image->GetBufferedRegion().IsInside(m_Seeds[i]) && this->IsPixelIncluded(m_Seeds[i]))
      {
        // Push the seed onto the queue
        m_IndexStack.push(m_Seeds[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_TempPtr->SetPixel(m_Seeds[i], 2);
      }
    }
  }
 
  void
  operator++() override
  {
    this->DoFloodStep();
  }
 
  void
  DoFloodStep();
 
  virtual SmartPointer<FunctionType>
  GetFunction() const
  {
    return m_Function;
  }
 
  void
  SetFullyConnected(const bool _arg);
 
  bool
  GetFullyConnected() const;
 
  itkBooleanMacro(FullyConnected);
 
  virtual const SeedsContainerType &
  GetSeeds() const
  {
    return m_Seeds;
  }
 
protected: // made protected so other iterators can access
  SmartPointer<FunctionType> m_Function;
 
  using TTempImage = Image<unsigned char, Self::NDimensions>;
 
  typename TTempImage::Pointer m_TempPtr;
 
  SeedsContainerType m_Seeds;
 
  typename ImageType::PointType m_ImageOrigin;
 
  typename ImageType::SpacingType m_ImageSpacing;
 
  NeighborhoodIteratorType m_NeighborhoodIterator;
 
  RegionType m_ImageRegion;
 
  std::queue<IndexType> m_IndexStack;
 
  FunctionInputType m_LocationVector;
 
  bool m_FoundUncheckedNeighbor;
 
  bool m_IsValidIndex;
 
  bool m_FullyConnected;
};
} // end namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkShapedFloodFilledFunctionConditionalConstIterator.hxx"
#endif
 
#endif