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

on image boundaries, construct point away from the surface 

this avoids creating point clusters on surfaces
This commit is contained in:
Jane Tournois 2022-10-20 12:07:27 +02:00
parent 27cece3914
commit db6b051973
1 changed files with 16 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
Mesh_3/include/CGAL/Mesh_3/initialize_triangulation_from_labeled_image.h
View File

@ -33,6 +33,7 @@ struct Get_point
{ {
const double vx, vy, vz; const double vx, vy, vz;
const double tx, ty, tz; const double tx, ty, tz;
const std::size_t xdim, ydim, zdim;
Get_point(const CGAL::Image_3* image) Get_point(const CGAL::Image_3* image)
: vx(image->vx()) : vx(image->vx())
, vy(image->vy()) , vy(image->vy())
@ -40,15 +41,27 @@ struct Get_point
, tx(image->tx()) , tx(image->tx())
, ty(image->ty()) , ty(image->ty())
, tz(image->tz()) , tz(image->tz())
, xdim(image->xdim())
, ydim(image->ydim())
, zdim(image->zdim())
{} {}
Point operator()(const std::size_t i, Point operator()(const std::size_t i,
const std::size_t j, const std::size_t j,
const std::size_t k) const const std::size_t k) const
{ {
return Point(double(i) * vx + tx, double x = double(i) * vx + tx;
double(j) * vy + ty, double y = double(j) * vy + ty;
double(k) * vz + tz); double z = double(k) * vz + tz;
if (i == 0) x += 1. / 6. * vx;
else if (i == xdim - 1) x -= 1. / 6. * vx;
if (j == 0) y += 1. / 6. * vy;
else if (j == ydim - 1) y -= 1. / 6. * vy;
if (k == 0) z += 1. / 6. * vz;
else if (k == zdim - 1) z -= 1. / 6. * vz;
return Point(x, y, z);
} }
}; };
template<class C3T3, class MeshDomain, class MeshCriteria> template<class C3T3, class MeshDomain, class MeshCriteria>