$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
image2d-skel-unconstrained.cc
1 // Copyright (C) 2010-2012 EPITA Research and Development Laboratory (LRDE)
2 //
3 // This file is part of Olena.
4 //
5 // Olena is free software: you can redistribute it and/or modify it under
6 // the terms of the GNU General Public License as published by the Free
7 // Software Foundation, version 2 of the License.
8 //
9 // Olena is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // General Public License for more details.
13 //
14 // You should have received a copy of the GNU General Public License
15 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
16 //
17 // As a special exception, you may use this file as part of a free
18 // software project without restriction. Specifically, if other files
19 // instantiate templates or use macros or inline functions from this
20 // file, or you compile this file and link it with other files to produce
21 // an executable, this file does not by itself cause the resulting
22 // executable to be covered by the GNU General Public License. This
23 // exception does not however invalidate any other reasons why the
24 // executable file might be covered by the GNU General Public License.
25 
28 
29 #include <mln/core/image/image2d.hh>
30 #include <mln/core/alias/neighb2d.hh>
31 
32 #include <mln/topo/skeleton/breadth_first_thinning.hh>
33 
34 #include <mln/topo/is_simple_point2d.hh>
35 #include <mln/topo/detach_point.hh>
36 
37 #include <mln/io/pbm/all.hh>
38 
39 #include <mln/util/timer.hh>
40 
41 
42 int
43 main(int argc, char* argv[])
44 {
45  using namespace mln;
46 
47  if (argc != 3)
48  {
49  std::cerr << "usage: " << argv[0] << " input.pbm output.pbm"
50  << std::endl;
51  std::exit(1);
52  }
53 
54  std::string input_filename = argv[1];
55  std::string output_filename = argv[2];
56 
57  typedef image2d<bool> I;
58  typedef neighb2d N;
59 
60  // Add a border of (at least) 1 pixel, to a guarantee a meaningful
61  // result of the computation of connectivity numbers (called within
62  // is_simple_2d); indeed, this computation always expects each pixel
63  // to have 8 neighboring sites.
65 
66  I input = io::pbm::load(input_filename);
67  std::cout << input.nrows() << " x " << input.ncols() << " = "
68  << input.nrows() * input.ncols() << " pixels"
69  << std::endl;
70 
71  // FIXME: Use a dual neighborhood instead?
72 
73  // Foreground neighborhood.
74  neighb2d nbh_fg = c4();
75  // Background neighborhood.
76  neighb2d nbh_bg = c8();
77 
78  // Simplicity criterion functor.
79  topo::is_simple_point2d<I, N> is_simple(nbh_fg, nbh_bg);
80  // Simple point detach procedure.
81  topo::detach_point<I> detach;
82 
83  util::timer t;
84  t.start();
85  I output = topo::skeleton::breadth_first_thinning(input, nbh_fg,
86  is_simple,
87  detach);
88  t.stop();
89  std::cout << t.read() << " s" << std::endl;
90 
91  io::pbm::save(output, output_filename);
92 }