itkHigherOrderAccurateDerivativeOperator.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright Insight Software Consortium
 *
 *  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 itkHigherOrderAccurateDerivativeOperator_h
#define itkHigherOrderAccurateDerivativeOperator_h

#include "itkNeighborhoodOperator.h"

namespace itk
{

template< class TPixel, unsigned int VDimension = 2,
          class TAllocator = NeighborhoodAllocator< TPixel > >
class HigherOrderAccurateDerivativeOperator:
  public NeighborhoodOperator< TPixel, VDimension, TAllocator >
{
public:

  typedef HigherOrderAccurateDerivativeOperator Self;
  typedef NeighborhoodOperator<
    TPixel, VDimension, TAllocator >           Superclass;

  typedef typename Superclass::PixelType     PixelType;
  typedef typename Superclass::PixelRealType PixelRealType;

  itkTypeMacro(HigherOrderAccurateDerivativeOperator, NeighborhoodOperator);

  HigherOrderAccurateDerivativeOperator():
    m_Order(1),
    m_OrderOfAccuracy(2)
  {}

  HigherOrderAccurateDerivativeOperator(const Self & other):
    NeighborhoodOperator< TPixel, VDimension, TAllocator >(other)
  { m_Order = other.m_Order;  }

  Self & operator=(const Self & other)
  {
    Superclass::operator=(other);
    m_Order = other.m_Order;
    return *this;
  }

  void SetOrder(const unsigned int & order)
  {
    this->m_Order = order;
  }

  unsigned int GetOrder() const { return m_Order; }

  void SetOrderOfAccuracy(const unsigned int & order)
  {
    this->m_OrderOfAccuracy = order;
  }

  unsigned int GetOrderOfAccuracy() const { return m_OrderOfAccuracy; }

  virtual void PrintSelf(std::ostream & os, Indent i) const
  {
    os << i << "HigherOrderAccurateDerivativeOperator { this=" << this
       << ", m_Order = " << m_Order
       << ", m_OrderOfAccuracy = " << m_OrderOfAccuracy << "}" << std::endl;
    Superclass::PrintSelf( os, i.GetNextIndent() );
  }

protected:
  typedef typename Superclass::CoefficientVector CoefficientVector;

  CoefficientVector GenerateCoefficients();

  void Fill(const CoefficientVector & coeff)
  {   Superclass::FillCenteredDirectional(coeff);  }
private:
  CoefficientVector GenerateFirstOrderCoefficients();

  unsigned int m_Order;

  unsigned int m_OrderOfAccuracy;
};

} // namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkHigherOrderAccurateDerivativeOperator.hxx"
#endif

#endif