+
+ // Create branches
+ TSkeleton* sk = this->GetOutput( );
+ typename std::vector< TIndex >::const_iterator eIt = end_points.begin( );
+ std::map< TIndex, unsigned long, typename TIndex::LexicographicCompare > tags;
+ unsigned long t = 1;
+ for( ; eIt != end_points.end( ); ++eIt )
+ {
+ // Tag path
+ TIndex it = *eIt;
+ TIndex p = mst->GetParent( it );
+ while( it != p )
+ {
+ tags[ it ] = t;
+ it = p;
+ p = mst->GetParent( it );
+
+ } // elihw
+ // TODO: More on it
+
+ } // rof