mitkItkNonUniformBSpline.h

Go to the documentation of this file.

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

The Medical Imaging Interaction Toolkit (MITK)

Copyright (c) German Cancer Research Center,
Division of Medical and Biological Informatics.
All rights reserved.

This software is distributed WITHOUT ANY WARRANTY; without
even the implied warranty of MERCHANTABILITY or FITNESS FOR
A PARTICULAR PURPOSE.

See LICENSE.txt or http://www.mitk.org for details.

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

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

This file is based heavily on a corresponding ITK filter.

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

/*********************************
  This file was taken from ITK, CVS version 1.8 to circumvent a bug in ITK release 3.18 (see http://public.kitware.com/Bug/view.php?id=10633
 *********************************/

#ifndef __itkNonUniformBSpline_h
#define __itkNonUniformBSpline_h

#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif

#include <vector>

#include "itkPoint.h"
#include "itkObject.h"
#include "itkObjectFactory.h"
#include "itkArray.h"

namespace itk {

template < unsigned int TDimension = 3 >
class NonUniformBSpline
  : public Object
{
public:
  typedef NonUniformBSpline                            Self;
  typedef Object                                       Superclass;
  typedef SmartPointer < Self >                        Pointer;
  typedef SmartPointer < const Self >                  ConstPointer;
  typedef double                                       ScalarType;
  typedef itk::Point< ScalarType, TDimension >         PointType;
  typedef std::vector < PointType >                    PointListType;
  typedef PointListType *                              PointListPointer;
  typedef std::vector < double >                       KnotListType;
  typedef std::vector<double>                          CoordinateListType;
  typedef itk::Point<double, TDimension >              ControlPointType;
  typedef std::vector< ControlPointType >              ControlPointListType;
  typedef ControlPointListType *                       ControlPointListPointer;
  typedef std::vector<double>                          ChordLengthListType;

  itkFactorylessNewMacro(Self)
  itkCloneMacro(Self)

  itkTypeMacro( NonUniformBSpline, Object );

  void SetPoints( PointListType & newPoints );

  const PointListType & GetPoints() const;

  void SetKnots( KnotListType & newKnots);

  const KnotListType & GetKnots() const;

  void ComputeChordLengths();

  PointType  EvaluateSpline(const Array<double> & p) const;
  PointType  EvaluateSpline( double t ) const;

  void ComputeControlPoints();

  void SetControlPoints( ControlPointListType& ctrlpts );

  const ControlPointListType & GetControlPoints() const;

  double NonUniformBSplineFunctionRecursive(unsigned int order, unsigned int i, double t) const;

  itkSetMacro( SplineOrder, unsigned int );
  itkGetConstReferenceMacro( SplineOrder, unsigned int );


protected:

  NonUniformBSpline();

  virtual ~NonUniformBSpline();

  virtual void PrintSelf( std::ostream& os, Indent indent ) const override;

  PointListType         m_Points;

  KnotListType          m_Knots;

  ControlPointListType  m_ControlPoints;

  ChordLengthListType   m_ChordLength;

  ChordLengthListType   m_CumulativeChordLength;

  unsigned int          m_SplineOrder;

  unsigned int          m_SpatialDimension;

};

} // end namespace itk

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


#endif // __itkNonUniformBSpline_h