itkNeighborhood.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkNeighborhood.h,v $
  Language:  C++
  Date:      $Date: 2002/09/11 19:57:13 $
  Version:   $Revision: 1.33 $

  Copyright (c) 2002 Insight 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 __itkNeighborhood_h
#define __itkNeighborhood_h

#include <iostream>
#include "itkNeighborhoodAllocator.h"
#include "itkSize.h"
#include "itkIndent.h"
#include "itkSliceIterator.h"
#include "vnl/vnl_vector.h"

namespace itk {
  
template<class TPixel, unsigned int VDimension = 2,
         class TAllocator = NeighborhoodAllocator<TPixel> >
class ITK_EXPORT Neighborhood 
{
public:
  typedef Neighborhood Self;

  typedef TAllocator AllocatorType;

  itkStaticConstMacro(NeighborhoodDimension, unsigned int, VDimension);
  
  typedef TPixel PixelType;
  
  typedef typename AllocatorType::iterator Iterator;
  typedef typename AllocatorType::const_iterator ConstIterator;
  
  typedef Size<VDimension> SizeType;
  typedef typename SizeType::SizeValueType SizeValueType;
  
  typedef Size<VDimension> RadiusType;

  typedef SliceIterator<TPixel, Self> SliceIteratorType;
  
  Neighborhood() { m_Radius.Fill(0); m_Size.Fill(0);  }

  virtual ~Neighborhood() {}
    
  Neighborhood(const Self& other)
  {
    m_Radius     = other.m_Radius;
    m_Size       = other.m_Size;
    m_DataBuffer = other.m_DataBuffer;
  }

  Self &operator=(const Self& other)
  {
    m_Radius     = other.m_Radius;
    m_Size       = other.m_Size;
    m_DataBuffer = other.m_DataBuffer;
    return *this;
  }

  bool
  operator==(const Self& other) const
  {
    return (m_Radius == other.m_Radius &&
            m_Size   == other.m_Size &&
            m_DataBuffer == other.m_DataBuffer);
  }

  bool operator!=(const Self& other) const
  {
    return (m_Radius != other.m_Radius ||
            m_Size   != other.m_Size ||
            m_DataBuffer != other.m_DataBuffer);
  }

  const SizeType GetRadius() const
    { return m_Radius; }

  unsigned long GetRadius(const unsigned long n) const
    { return m_Radius[n]; }

  unsigned long GetSize(const unsigned long n) const
    { return m_Size[n]; }

  SizeType GetSize() const
    { return m_Size; }
  
  unsigned long GetStride(const unsigned long) const;
  
  Iterator End()
    { return m_DataBuffer.end(); }
  Iterator Begin()
    { return m_DataBuffer.begin(); }
  ConstIterator End() const
    { return m_DataBuffer.end(); }
  ConstIterator Begin() const
    { return m_DataBuffer.begin(); }
  
  unsigned int Size() const
    { return m_DataBuffer.size(); }
  
  TPixel &operator[](unsigned int i)
    { return m_DataBuffer[i]; }
  const TPixel &operator[](unsigned int i) const
  { return m_DataBuffer[i]; }

  TPixel GetCenterValue() const
  {    return (this->operator[]((this->Size())>>1)); }

  void SetRadius(const SizeType &);

  void SetRadius(const unsigned long *rad)
  {
    SizeType s;
    memcpy(s.m_Size, rad, sizeof(unsigned long) * VDimension);
    this->SetRadius(s);
  }
  
  void SetRadius(const unsigned long);

  void Print(std::ostream& os) const
    { this->PrintSelf(os, Indent(0));  }

  AllocatorType &GetBufferReference()
    { return m_DataBuffer; }
  const AllocatorType &GetBufferReference() const
    { return m_DataBuffer; }
  
protected:
  void SetSize()
  {
    for (unsigned int i=0; i<VDimension; ++i)
      { m_Size[i] = m_Radius[i]*2+1; }
  }

  virtual void Allocate(unsigned int i)
    { m_DataBuffer.resize(i); }

  virtual void PrintSelf(std::ostream&, Indent) const;
  
private:
  SizeType m_Radius;

  SizeType m_Size;

  AllocatorType m_DataBuffer;
};

template <class TPixel, unsigned int VDimension, class TContainer>
std::ostream & operator<<(std::ostream &os, const Neighborhood<TPixel,VDimension,TContainer> &neighborhood)
{
  os << "Neighborhood:" << std::endl;
  os << "    Radius:" << neighborhood.GetRadius() << std::endl;
  os << "    Size:" << neighborhood.GetSize() << std::endl;
  os << "    DataBuffer:" << neighborhood.GetBufferReference() << std::endl;

  return os;
}

} // namespace itk

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

#endif

---