songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

replace PCA of points by PCA of triangles 

to make projection more precise and avoid moving a point inside a
protecting ball
This commit is contained in:
Jane Tournois 2020-11-26 16:25:55 +01:00
parent a8f143d6a8
commit 9666b3cb73
1 changed files with 15 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
Mesh_3/include/CGAL/Mesh_3/C3T3_helpers.h
View File

@ -3441,6 +3441,8 @@ get_least_square_surface_plane(const Vertex_handle& v,
// Get adjacent surface points
std::vector<Bare_point> surface_point_vector;
std::vector<typename C3T3::Triangulation::Triangle> triangles;
typename Facet_vector::iterator fit = facets.begin(), fend = facets.end();
for ( ; fit != fend; ++fit )
{
@ -3460,6 +3462,8 @@ get_least_square_surface_plane(const Vertex_handle& v,
const Bare_point& bp = tr_.get_closest_point(cp(position),
cell->get_facet_surface_center(i));
surface_point_vector.push_back(bp);
triangles.push_back(c3t3_.triangulation().triangle(*fit));
}
}
@ -3470,11 +3474,19 @@ get_least_square_surface_plane(const Vertex_handle& v,
// Compute least square fitting plane
Plane_3 plane;
Bare_point point;
CGAL::linear_least_squares_fitting_3(surface_point_vector.begin(),
surface_point_vector.end(),
// CGAL::linear_least_squares_fitting_3(surface_point_vector.begin(),
// surface_point_vector.end(),
// plane,
// point,
// Dimension_tag<0>(),
// tr_.geom_traits(),
// Default_diagonalize_traits<FT, 3>());
CGAL::linear_least_squares_fitting_3(triangles.begin(),
triangles.end(),
plane,
point,
Dimension_tag<0>(),
Dimension_tag<2>(),
tr_.geom_traits(),
Default_diagonalize_traits<FT, 3>());