Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
world/inter_pixel/compute.cc
1 // Copyright (C) 2009 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 
26 #include <cmath>
27 
28 #include <mln/core/var.hh>
29 #include <mln/core/image/image2d.hh>
30 #include <mln/data/compare.hh>
31 #include <mln/world/inter_pixel/immerse.hh>
32 #include <mln/world/inter_pixel/compute.hh>
33 
34 #include <mln/world/inter_pixel/neighb2d.hh>
35 #include <mln/morpho/watershed/flooding.hh>
36 #include <mln/debug/println.hh>
37 
38 
39 
40 struct d_t : mln::Function_vv2v<d_t>
41 {
42  typedef int result;
43 
44  int operator()(int i1, int i2) const
45  {
46  return std::abs(i2 - i1);
47  }
48 }
49  d;
50 
51 
52 
53 int main()
54 {
55  using namespace mln;
56 
57  int vals[] = { 3, 4, 5,
58  1, 3, 6 ,
59  8, 7, 3 } ;
60 
61  typedef image2d<int> I;
62  I ima = make::image2d(vals);
63 
64  using namespace world::inter_pixel;
65  typedef image_if<I, is_pixel> Ix;
66  Ix imax = immerse(ima);
67 
68 
69  int refs[] = { 0, 1, 0, 1, 0,
70  2, 0, 1, 0, 1,
71  0, 2, 0, 3, 0,
72  7, 0, 4, 0, 3,
73  0, 1, 0, 4, 0 };
74 
75  mln_assertion(compute(imax, d) == (make::image2d(refs) | is_separator()));
76 
77  mln_VAR(g, compute(imax, d));
78 
79  unsigned n_basins;
80  debug::println(morpho::watershed::flooding(g, e2e(), n_basins));
81 }