$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
distance_and_closest_point_geodesic.cc
1 // Copyright (C) 2008, 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 <mln/core/image/image2d.hh>
27 #include <mln/core/alias/neighb2d.hh>
28 #include <mln/value/int_u8.hh>
29 #include <mln/data/fill.hh>
30 #include <mln/debug/println.hh>
31 #include <mln/opt/at.hh>
32 #include <mln/data/compare.hh>
33 
34 #include <mln/transform/distance_and_closest_point_geodesic.hh>
35 
36 unsigned dmap_ref[] = { 4, 4, 4, 3, 2, 3, 4, 4, 4,
37  4, 4, 3, 2, 1, 2, 3, 4, 4,
38  4, 3, 2, 1, 0, 1, 2, 3, 4,
39  3, 2, 1, 2, 1, 2, 3, 4, 4,
40  2, 1, 0, 1, 2, 3, 4, 4, 4,
41  3, 2, 1, 2, 3, 4, 4, 4, 4,
42  4, 3, 2, 3, 4, 4, 4, 4, 4,
43  4, 4, 3, 4, 4, 4, 4, 4, 4,
44  4, 4, 4, 4, 4, 4, 4, 4, 4 };
45 
46 
47 unsigned cp_idx_ref[] = { 2, 2, 0, 0, 0, 0, 0, 2, 2,
48  2, 0, 0, 0, 0, 0, 0, 0, 2,
49  1, 0, 0, 0, 0, 0, 0, 0, 0,
50  1, 1, 1, 0, 0, 0, 0, 0, 2,
51  1, 1, 1, 1, 0, 0, 0, 2, 2,
52  1, 1, 1, 1, 0, 0, 2, 2, 2,
53  1, 1, 1, 1, 1, 2, 2, 2, 2,
54  2, 1, 1, 1, 2, 2, 2, 2, 2,
55  2, 2, 1, 2, 2, 2, 2, 2, 2 };
56 
57 
58 int main()
59 {
60  using namespace mln;
61  using value::int_u8;
62 
63  image2d<bool> input(9, 9);
64  data::fill(input, false);
65  opt::at(input, 2, 4) = true;
66  opt::at(input, 4, 2) = true;
67 
68  {
70  output_t output = transform::distance_and_closest_point_geodesic(input,
71  c4(),
72  int_u8(4));
73  mln_assertion(output.first() == make::image2d(dmap_ref));
75  }
76 
77  {
78  p_array<point2d> arr;
79  arr.insert(point2d(2, 4));
80  arr.insert(point2d(4, 2));
81  typedef util::couple<image2d<int_u8>, image2d<unsigned> > output_t;
83  input.domain(),
84  c4(),
85  int_u8(4));
86  mln_assertion(output.first() == make::image2d(dmap_ref));
87  mln_assertion(output.second() == make::image2d(cp_idx_ref));
88  }
89 }