handle non-manifold vertices when copying a mesh

This commit is contained in:
Sébastien Loriot 2017-07-26 11:59:02 +02:00
parent 594bed8350
commit bb175cdc66
1 changed files with 22 additions and 7 deletions

View File

@ -61,6 +61,7 @@ void copy_face_graph_impl(const SourceMesh& sm, TargetMesh& tm,
conv;
std::vector<tm_halfedge_descriptor> tm_border_halfedges;
std::vector<sm_halfedge_descriptor> sm_border_halfedges;
tm_face_descriptor tm_null_face = boost::graph_traits<TargetMesh>::null_face();
@ -78,12 +79,16 @@ void copy_face_graph_impl(const SourceMesh& sm, TargetMesh& tm,
if ( is_border(sm_h, sm) ){
tm_border_halfedges.push_back( tm_h );
sm_border_halfedges.push_back( sm_h );
set_face(tm_h, tm_null_face, tm);
CGAL_assertion(next(tm_h, tm) == boost::graph_traits<TargetMesh>::null_halfedge() );
}
if( is_border(sm_h_opp, sm) ){
tm_border_halfedges.push_back( tm_h_opp );
sm_border_halfedges.push_back( sm_h_opp );
set_face(tm_h_opp, tm_null_face, tm);
CGAL_assertion(next(tm_h_opp, tm) == boost::graph_traits<TargetMesh>::null_halfedge() );
}
//create a copy of interior vertices only once
@ -126,15 +131,25 @@ void copy_face_graph_impl(const SourceMesh& sm, TargetMesh& tm,
}
}
// update next/prev of tm border halfedges by visiting faces of the mesh
BOOST_FOREACH(tm_halfedge_descriptor h, tm_border_halfedges)
// update next/prev of tm border halfedges
std::size_t nb_border_hedges = tm_border_halfedges.size();
for (std::size_t i=0; i< nb_border_hedges; ++i)
{
tm_halfedge_descriptor candidate = opposite(h, tm);
do{
candidate = opposite( prev(candidate, tm), tm);
tm_halfedge_descriptor tm_h = tm_border_halfedges[i];
if ( next(tm_h, tm) != boost::graph_traits<TargetMesh>::null_halfedge() )
continue; //already set
tm_halfedge_descriptor tm_h_prev = tm_h;
CGAL_precondition(*halfedges_around_face(sm_border_halfedges[i], sm).first == sm_border_halfedges[i]);
BOOST_FOREACH(sm_halfedge_descriptor sm_h,
halfedges_around_face(next(sm_border_halfedges[i], sm), sm))
{
CGAL_assertion(next(tm_h_prev, tm) == boost::graph_traits<TargetMesh>::null_halfedge());
tm_h = get(hmap, sm_h);
set_next(tm_h_prev, tm_h, tm);
tm_h_prev=tm_h;
}
while( !is_border(candidate, tm) );
set_next(h, candidate, tm);
}
// update halfedge vertex of all but the vertex halfedge