1 /*# ---------------------------------------------------------------------
3 # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
5 # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
6 # Previous Authors : Laurent Guigues, Jean-Pierre Roux
7 # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
9 # This software is governed by the CeCILL-B license under French law and
10 # abiding by the rules of distribution of free software. You can use,
11 # modify and/ or redistribute the software under the terms of the CeCILL-B
12 # license as circulated by CEA, CNRS and INRIA at the following URL
13 # http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
14 # or in the file LICENSE.txt.
16 # As a counterpart to the access to the source code and rights to copy,
17 # modify and redistribute granted by the license, users are provided only
18 # with a limited warranty and the software's author, the holder of the
19 # economic rights, and the successive licensors have only limited
22 # The fact that you are presently reading this means that you have had
23 # knowledge of the CeCILL-B license and that you accept its terms.
24 # ------------------------------------------------------------------------ */
26 /*=========================================================================
29 Module: $RCSfile: curve.hxx,v $
31 Date: $Date: 2012/11/15 14:15:31 $
32 Version: $Revision: 1.3 $
34 Copyright: (c) 2002, 2003
37 This software is distributed WITHOUT ANY WARRANTY; without even
38 the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
39 PURPOSE. See the above copyright notice for more information.
41 =========================================================================*/
43 #ifndef KGFO__CURVE__HXX
44 #define KGFO__CURVE__HXX
49 /*! \brief Parametrized n-dimensional curve.
51 * This class realizes a parametric curve model, given some control points.
53 class MAR_KERNEL_EXPORT kCurve
57 kCurve( size_t dim = 3, size_t sdim = 0 );
60 kCurve& operator=( const kCurve& r );
62 uint32_t getDimension( ) const
67 uint32_t getStateDimension( ) const
72 uint32_t getClosestControlPoint( double* p );
74 void getPoint( double* p, double s );
75 void getState( double* st, double s );
76 void getTangent( double* tg, double s );
77 void getNormal( double* n, double s );
78 void getBinormal( double* b, double s );
80 void addControlPoint( double* cp, double* sp );
81 void getControlPoint( int i, double* cp, double* sp );
82 void setControlPoint( int i, double* cp, double* sp );
83 uint32_t getNumberOfControlPoints( )
85 return( _controlPoints.size( ) );
88 double length( double step = 0.01 );
90 double projectOverControlPoints( double* pr, double* pt );
91 double projectOverCurve( double* pr, double* pt );
93 void evaluate( double* p, int i, double t );
94 void evalState( double* s, int i, double t );
96 void derivative1( double* d, int i, double t );
97 void derivative2( double* d, int i, double t );
103 static int MAX_DIMENSION;
104 static int MAX_STATE_DIMENSION;
108 std::vector< double* > _controlPoints;
109 std::vector< double* > _controlState;
115 void loadCatmullControlMatrix( int i );
116 void loadCatmullStateMatrix( int i );
117 void calculeSplineArguments( double s, int& i, double& t );
119 std::vector< double > _controlT;
120 std::vector< double > _controlL;
127 #endif // KGFO__CURVE__HXX