$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
flooding.cc
1 // Copyright (C) 2008, 2009, 2010 EPITA Research and Development
2 // Laboratory (LRDE)
3 //
4 // This file is part of Olena.
5 //
6 // Olena is free software: you can redistribute it and/or modify it under
7 // the terms of the GNU General Public License as published by the Free
8 // Software Foundation, version 2 of the License.
9 //
10 // Olena is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 // General Public License for more details.
14 //
15 // You should have received a copy of the GNU General Public License
16 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
17 //
18 // As a special exception, you may use this file as part of a free
19 // software project without restriction. Specifically, if other files
20 // instantiate templates or use macros or inline functions from this
21 // file, or you compile this file and link it with other files to produce
22 // an executable, this file does not by itself cause the resulting
23 // executable to be covered by the GNU General Public License. This
24 // exception does not however invalidate any other reasons why the
25 // executable file might be covered by the GNU General Public License.
26 
27 #include <iostream>
28 
29 #include <mln/core/image/image2d.hh>
30 #include <mln/core/alias/window2d.hh>
31 #include <mln/core/alias/neighb2d.hh>
32 
33 #include <mln/value/int_u8.hh>
34 #include <mln/value/int_u16.hh>
35 
36 #include <mln/morpho/watershed/flooding.hh>
37 #include <mln/data/transform.hh>
38 
39 #include <mln/io/pgm/load.hh>
40 #include <mln/io/pgm/save.hh>
41 
42 #include <mln/util/timer.hh>
43 
44 #include "tests/data.hh"
45 
46 struct f_16_to_8 : mln::Function_v2v< f_16_to_8 >
47 {
48  typedef mln::value::int_u8 result;
49  result operator()(const mln::value::int_u16& v) const
50  {
51  if (v == 0)
52  return 0;
53  return 1 + (v - 1) % 255;
54  }
55 };
56 
57 
58 int main()
59 {
60  using namespace mln;
61  using value::int_u8;
62  using value::int_u16;
63 
64  image2d<int_u8> input;
65  io::pgm::load(input, MLN_IMG_DIR "/squares.pgm");
66 
67  typedef int_u16 L;
68  L n_basins;
69  {
70  util::timer t;
71  t.start();
72  image2d<L> output =
74  std::cout << "gen: " << t << std::endl;
75  io::pgm::save(data::transform(output, f_16_to_8()),
76  "flooding-tmp_ref.pgm");
77  }
78  {
79  util::timer t;
80  t.start();
81  image2d<L> output =
82  morpho::watershed::impl::flooding_fastest(input, c4(), n_basins);
83  std::cout << "fast: " << t << std::endl;
84  io::pgm::save(data::transform(output, f_16_to_8()),
85  "flooding-tmp_out.pgm");
86  }
87 }