itkNonUniformBSpline.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkNonUniformBSpline.h,v $
  Language:  C++
  Date:      $Date: 2009-02-06 20:53:13 $
  Version:   $Revision: 1.5 $

  Copyright (c) Insight Software 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 __itkNonUniformBSpline_h
#define __itkNonUniformBSpline_h

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

#include <vector>

#include "itkPoint.h"
#include "itkObject.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;

  itkNewMacro( Self );

  itkTypeMacro( NonUniformBSpline, Object );

  void SetPoints( PointListType & newPoints );

  PointListType& GetPoints()
    {
    return m_Points;
    }

  void SetKnots( KnotListType & newKnots);

  KnotListType& GetKnots();

  void ComputeChordLengths();

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

  void ComputeControlPoints();

  void SetControlPoints( ControlPointListType& ctrlpts );

  ControlPointListType& GetControlPoints()
    { 
    return m_ControlPoints;
    }

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

  void SetSplineOrder(unsigned int order)
    {
    m_SplineOrder = order;
    this->Modified();
    }

  unsigned int GetSplineOrder()
    {
    return m_SplineOrder;
    }

protected:

  NonUniformBSpline();

  virtual ~NonUniformBSpline();

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

  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 "itkNonUniformBSpline.txx" 
#endif 


#endif // __itkNonUniformBSpline_h