mitkItkNonUniformBSpline.h

Go to the documentation of this file.

/*============================================================================
 
The Medical Imaging Interaction Toolkit (MITK)
 
Copyright (c) German Cancer Research Center (DKFZ)
All rights reserved.
 
Use of this source code is governed by a 3-clause BSD license that can be
found in the LICENSE file.
 
============================================================================*/
 
#ifndef mitkItkNonUniformBSpline_h
#define mitkItkNonUniformBSpline_h
 
/*===================================================================
 
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
 *********************************/
 
#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();
 
  ~NonUniformBSpline() override;
 
  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