CGAL users discussion list

[Yann Cobigo <>]

Hi everyone,

I try to produce a tetrahedral mesh (M3D) based on a file file.inr (I3D). In this purpose, I use the example Mesh_3/mesh_3D_image.cpp.

My image, I3D, is divided into n domains with indexes: {i_{d}}_{d = 0 .. n-1}.

If my domains' indexes are consecutive, i.e. i_{0} = 0, i_{1} = 1, i_{2} = 2, ..., i_{n-1} = n-1, the tetrahedrization produces a mesh M3D with the same index nomenclature for the domains: m_{1} = 1, m_{2} = 2, ..., m_{n-1} = n-1.

On the other hand, if ths image I3D domains' indexes are not consecutive, i.e. i_{0} = 0, i_{1} = 1, i_{2} = 13, ..., i_{n-1} = N, the indexes of the mesh no longer respect the same I3D nomenclature: m_{1} = 1, m_{2} = 2, ..., m_{n-1} = n-1, instead of m_{1} = 1, m_{2} = 13, ..., m_{n-1} = N.

Is there a way to ensure that the mesh domains' indexes m_{d} be the same as the image domains' indexes i_ {d}?

Best regards,
Yann