CGAL users discussion list

["Laurent Rineau (CGAL/GeometryFactory)" <>]

Le Friday 13 June 2014 09:59:24 lednakashim a écrit :
> Hi Folks,
> 
> I was wondering if somebody could kindly point me in the right direction.
> 
> I want to generate a volume from an 3D image stack (`
> Mesh_3/mesh_3D_image.cpp
> <http://doc.cgal.org/latest/Mesh_3/Mesh_3_2mesh_3D_image_8cpp-example.html>
> `), modify the mesh, and get back binary mask to extract the values inside
> my mesh. I want to use the mesh as a kind of `3D` threshold.
> 
> How can I go from a mesh back to an 3D image stack?

Given a 3D point, the function c3t3.triangulation().locate(point) can locale 
the triangulation cell that contains the point. Then, your can read the 
domain 
index of the cell with:
  int domain_index = c3t3.subdomain_index(cell_handle);

You can reconstruct a new image using that trick on a regular 3D grid of 
points.

-- 
Laurent Rineau, PhD
R&D Engineer at GeometryFactory           http://www.geometryfactory.com/
Release Manager of the CGAL Project       http://www.cgal.org/