CGAL users discussion list

[DemoKioussis <>]

Hello, I am trying to use mesh simplification to eventually simplify large
meshes, however I am experimenting on a relatively small mesh (about 80mb)
for now. Eventually I will probably move to make my own version of the cost
/ placement functions already defined, but I am having trouble with the
Lindstrom-Turk placement function. I have no problem using the midpoint
placement function. The error doesn't seem to be in the cost function.

The error I am having is that when I simplify with Lindstrom-Turk placement
( with default parameters) I get many spikes all over the mesh, but
simplifying with midpoint placement results in a good simplification ( with
either edge-length cost or Lindstrom-Turk cost). 

Here are images describing the situation in better detail

Original mesh:
  <http://cgal-discuss.949826.n4.nabble.com/file/t375925/Original.png> 

Lindstrom-Turk placement with Lindstrom-Turk cost and both with default
parameters:
<http://cgal-discuss.949826.n4.nabble.com/file/t375925/LindstromTurk.png> 


Midpoint placement with Edge Length cost:
<http://cgal-discuss.949826.n4.nabble.com/file/t375925/Midpoint.png> 


Midpoint placement with Lindstrom-Turk cost: 
<http://cgal-discuss.949826.n4.nabble.com/file/t375925/Midpoint_LindstromCost.png>
 

The borders of the midpoint placements are preserved which is why they were
not simplified.

Here is how i am calling the simplification


///////////////////////////////////////////////////////
Lindstrom-Turk placement with Lindstrom-Turk cost and both with default
parameters:

int iRes = SMS::edge_collapse( surfaceMesh, stop,
                                        CGAL::parameters::vertex_index_map( 
get( CGAL::vertex_external_index,
surfaceMesh) )
                                        .halfedge_index_map( get( 
CGAL::halfedge_external_index, surfaceMesh) )
                                        .get_placement( 
SMS::LindstromTurk_placement<SurfaceMesh>() )
                                        .get_cost( SMS::LindstromTurk_cost 
<SurfaceMesh>() ) );

I believe this is equivalent to:
int iRes = SMS::edge_collapse( surfaceMesh, stop,
                                        CGAL::parameters::vertex_index_map( 
get( CGAL::vertex_external_index,
surfaceMesh) )
                                        .halfedge_index_map( get( 
CGAL::halfedge_external_index, surfaceMesh) )

which results in identical errors

///////////////////////////////////////////////////////
Midpoint placement with Edge Length cost:

                                iRes = SMS::edge_collapse( surfaceMesh, stop,
                                           
CGAL::parameters::vertex_index_map( get(
CGAL::vertex_external_index, surfaceMesh) )
                                           .halfedge_index_map( get( 
CGAL::halfedge_external_index,
surfaceMesh) )
                                           .edge_is_constrained_map( 
constraintMap )
                                           .get_placement(
SMS::Constrained_placement<SMS::Midpoint_placement&lt;SurfaceMesh>,
ContstrainedEdgeMap>( constraintMap ))
                                           .get_cost( 
SMS::Edge_length_cost<SurfaceMesh>() ) );

///////////////////////////////////////////////////////
ContstrainedEdgeMap is defined very similarly to how it was defined here:
https://doc.cgal.org/latest/Surface_mesh_simplification/Surface_mesh_simplification_2edge_collapse_OpenMesh_8cpp-example.html

I would assume that the surface mesh is valid because further up I call
CGAL::is_valid( surfaceMesh) && CGAL::is_triangle_mesh( surfaceMesh) before
simplification, and also the same mesh simplifies correctly with midpoint
placement.

I don't know what other details you may need, but I am willing give more



--
Sent from: http://cgal-discuss.949826.n4.nabble.com/