55unsigned char findMainAxis (
const T & curve,
const long double & t )
58 value[0] = std::abs ( curve.xp ( t )[0] );
59 value[1] = std::abs ( curve.xp ( t )[1] );
60 value[2] = std::abs ( curve.xp ( t )[2] );
62 if ( value[0] >= value[1] && value[0] >= value[2] )
64 else if ( value[1] >= value[0] && value[1] >= value[2] )
71int main(
int argc,
char** argv )
74 typedef EllipticHelix < Space > MyHelix;
75 typedef NaiveParametricCurveDigitizer3D < MyHelix > DigitizerHelix;
79 trace.info() <<
"exampleParamCurve3dDigitization" << endl;
84 MyDigitalCurve digitalCurve;
86 MyHelix helix( 15, 10, 1 );
88 digitize.init ( M_PI / 2., ( MyHelix::getPeriod() * 10. ) + M_PI / 2., 0.0001 );
93 digitize.digitize( back_insert_iterator < MyDigitalCurve> ( digitalCurve ), back_insert_iterator < MyMetaData > ( metaData ) );
96 trace.info() <<
"Number of points: " << digitalCurve.size () <<
" number of metadata: " << metaData.size () << endl;
99 for (
unsigned int i = 0; i < digitalCurve.size ( ); i++ )
101 if (
findMainAxis ( helix, metaData.at ( i ).first ) == 0 )
102 viewer <<
Color ( 255, 0, 0, 128 );
103 if (
findMainAxis ( helix, metaData.at ( i ).first ) == 1 )
104 viewer <<
Color ( 0, 255, 0, 128 );
105 if (
findMainAxis ( helix, metaData.at ( i ).first ) == 2 )
106 viewer <<
Color ( 0, 0, 255, 128 );
107 viewer << digitalCurve.at ( i );
void digitize(Shape &shape, std::vector< SCell > &sCells0, std::vector< SCell > &sCells1, KSpace &kspace, const double h)