00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkNeighborhoodOperator_h
00018 #define __itkNeighborhoodOperator_h
00019
00020 #include "itkNeighborhood.h"
00021 #include "itkExceptionObject.h"
00022 #include "itkNumericTraits.h"
00023 #include <vector>
00024
00025 namespace itk {
00062 template< class TPixel, unsigned int VDimension,
00063 class TAllocator = NeighborhoodAllocator<TPixel> >
00064 class ITK_EXPORT NeighborhoodOperator
00065 : public Neighborhood<TPixel, VDimension, TAllocator>
00066 {
00067 public:
00069 typedef NeighborhoodOperator Self;
00070 typedef Neighborhood<TPixel, VDimension, TAllocator> Superclass;
00071
00073 typedef typename Superclass::SizeType SizeType;
00074
00076 typedef TPixel PixelType;
00077
00079 typedef SliceIterator<TPixel, Self> SliceIteratorType;
00080
00082 NeighborhoodOperator()
00083 { m_Direction = 0; }
00084
00086 NeighborhoodOperator(const Self &orig)
00087 : Neighborhood<TPixel, VDimension, TAllocator>(orig)
00088 { m_Direction = orig.m_Direction; }
00089
00091 Self &operator=( const Self &orig )
00092 {
00093 Superclass::operator=(orig);
00094 m_Direction = orig.m_Direction;
00095 return *this;
00096 }
00098
00100 void SetDirection(const unsigned long &direction)
00101 { m_Direction = direction; }
00102
00104 unsigned long GetDirection() const
00105 { return m_Direction; }
00106
00111 virtual void CreateDirectional();
00112
00117 virtual void CreateToRadius(const SizeType &);
00118
00123 virtual void CreateToRadius(const unsigned long);
00124
00127 virtual void FlipAxes();
00128
00130 virtual void PrintSelf(std::ostream& os, Indent i) const
00131 {
00132 os << i << "NeighborhoodOperator { this=" << this
00133 << " Direction = " << m_Direction << " }" << std::endl;
00134 Superclass::PrintSelf( os, i.GetNextIndent() );
00135 }
00137
00138 typedef typename NumericTraits< TPixel >::RealType PixelRealType;
00139
00141 void ScaleCoefficients( PixelRealType );
00142
00143 protected:
00146 typedef std::vector<double> CoefficientVector;
00147
00150 virtual CoefficientVector GenerateCoefficients() = 0;
00151
00154 virtual void Fill(const CoefficientVector &) = 0;
00155
00161 virtual void FillCenteredDirectional(const CoefficientVector &);
00162
00164 void InitializeToZero()
00165 {
00166 for (unsigned int i = 0; i< this->Size(); ++i)
00167 {
00168 this->operator[](i) = NumericTraits<PixelType>::Zero;
00169 }
00170 }
00172
00173 private:
00175 unsigned long m_Direction;
00176 };
00177
00178 }
00179
00180
00181 #define ITK_TEMPLATE_NeighborhoodOperator(_, EXPORT, x, y) namespace itk { \
00182 _(2(class EXPORT NeighborhoodOperator< ITK_TEMPLATE_2 x >)) \
00183 namespace Templates { typedef NeighborhoodOperator< ITK_TEMPLATE_2 x > \
00184 NeighborhoodOperator##y; } \
00185 }
00186
00187 #if ITK_TEMPLATE_EXPLICIT
00188 # include "Templates/itkNeighborhoodOperator+-.h"
00189 #endif
00190
00191 #if ITK_TEMPLATE_TXX
00192 # include "itkNeighborhoodOperator.txx"
00193 #endif
00194
00195
00196
00197
00198
00199
00200 #endif
00201