MeshNavigationT.cc 3.38 KB
Newer Older
Jan Möbius's avatar
 
Jan Möbius committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109
//=============================================================================
//
//                               OpenFlipper
//        Copyright (C) 2008 by Computer Graphics Group, RWTH Aachen
//                           www.openflipper.org
//
//-----------------------------------------------------------------------------
//
//                                License
//
//  OpenFlipper is free software: you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or
//  (at your option) any later version.
// 
//  OpenFlipper is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU Lesser General Public License for more details.
// 
//  You should have received a copy of the GNU Lesser General Public License
//  along with OpenFlipper.  If not, see <http://www.gnu.org/licenses/>.
//
//-----------------------------------------------------------------------------
//
//   $Revision$
//   $Author$
//   $Date$
//
//=============================================================================




//=============================================================================
//
//  IMPLEMENTATION
//
//=============================================================================

#define MESHNAVIGATIONT_C

//== INCLUDES =================================================================

#include "MeshNavigationT.hh"

#include <OpenMesh/Core/Mesh/TriMesh_ArrayKernelT.hh>
#include <OpenMesh/Core/Geometry/MathDefs.hh>

#include <queue>
//== NAMESPACES ===============================================================

namespace MeshNavigation {

//== IMPLEMENTATION ==========================================================

template < typename MeshT >
inline
typename MeshT::HalfedgeHandle
opposite_halfedge(MeshT& _mesh , typename MeshT::HalfedgeHandle _he) {
   typename MeshT::VertexOHalfedgeIter opp_he( _mesh , _he);
   ++opp_he;
   ++opp_he;
   ++opp_he;
   return (opp_he.handle());
}

template < typename MeshT >
inline
typename MeshT::VertexHandle
findClosestBoundary(MeshT* _mesh , typename MeshT::VertexHandle _vh){

  //add visited property
  OpenMesh::VPropHandleT< bool > visited;
  _mesh->add_property(visited,"Visited Property" );
  
  //init visited property
  typename MeshT::VertexIter v_it, v_end = _mesh->vertices_end();
  for( v_it = _mesh->vertices_begin(); v_it != v_end; ++v_it )
    _mesh->property( visited, v_it ) = false;

  std::queue< typename MeshT::VertexHandle > queue;
  queue.push( _vh );

  while(!queue.empty()){
    typename MeshT::VertexHandle vh = queue.front();
    queue.pop();
    if (_mesh->property(visited, vh)) continue;

    for (typename MeshT::VertexOHalfedgeIter voh_it(*_mesh,vh); voh_it; ++voh_it){

      if ( _mesh->is_boundary(voh_it) ){
        _mesh->remove_property(visited);
        return _mesh->to_vertex_handle(voh_it);
      }else{
        queue.push( _mesh->to_vertex_handle(voh_it) );
      }
    }
    _mesh->property(visited, vh) = true;
  }

  _mesh->remove_property(visited);
  return typename MeshT::VertexHandle(-1);
}


//=============================================================================
} // MeshNavigation Namespace 
//=============================================================================