CGAL users discussion list

["Laurent Rineau (GeometryFactory)" <>]

On Wednesday 11 March 2009 22:23:58 Qianqian Fang wrote:
> The attached zip file contains a 3D binary image stored in INR format.
> I found that CGAL surface mesher produced strange results with this
> dataset.
>
> By using a threadhold value=0.5, bounding sphere center=(122.0, 82.0,
> 93.0),
> radius^2=32258, and CGAL::Surface_mesh_default_criteria_3<Tr>
> criteria(30,4,4),
> I run the mesh_a_3d_gray_image example with this data, and unfortunately,
> the program had ran forever and exhausted the system memory.

Hi,

You should give us the usual needed informations to get help:
  http://www.cgal.org/bug_report.html
At the very least: which CGAL version, which operating system, which 
compiler, 
which compiler version?

As for your problem, my crystal ball tells me that you should try to change 
the manifold tag, in the following line:
CGAL::make_surface_mesh(c2t3, surface, criteria, CGAL::Manifold_tag());
Change "Manifold_tag" into "Non_manifold_tag". The levelset 0.5 of your image 
is probably a non manifold surface.

Please also post a working .cpp file that shows your problem. You tell us 
that 
it is a modified version of examples/Surface_mesher/mesh_a_3d_gray_image.cpp, 
but I prefere to see all your modifications at once.

-- 
Laurent Rineau, PhD
Engineer at GeometryFactory
http://www.geometryfactory.com/