32 #include <mln/core/image/image3d.hh>
33 #include <mln/core/alias/neighb3d.hh>
34 #include <mln/make/box3d.hh>
36 #include <mln/core/image/vmorph/fun_image.hh>
37 #include <mln/fun/v2b/lnot.hh>
39 #include <mln/border/fill.hh>
41 #include <mln/labeling/blobs.hh>
43 #include <mln/topo/skeleton/breadth_first_thinning.hh>
45 #include <mln/topo/is_simple_point3d.hh>
46 #include <mln/topo/detach_point.hh>
47 #include <mln/topo/is_not_end_point.hh>
49 #include "apps/data.hh"
51 #include "image3d-skel.hh"
76 I input = make_triple_torus(20);
78 std::cerr << input.domain() << std::endl;
79 save_raw_3d(input,
"image3d-skel-with-end-points-input.raw");
80 save_vtk_polyhedrons(input,
"image3d-skel-with-end-points-input.vtk");
87 std::cerr <<
"n_fg_comps = " << n_fg_comps << std::endl;
88 std::cerr <<
"n_bg_comps = " << n_bg_comps << std::endl;
112 save_raw_3d(output,
"image3d-skel-with-end-points-skel.raw");
113 save_vtk_polyhedrons(output,
"image3d-skel-with-end-points-skel.vtk");