$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
crest.cc
1 // Copyright (C) 2009, 2010 EPITA Research and Development Laboratory
2 // (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 <mln/core/alias/neighb2d.hh>
28 # include <mln/core/image/image2d.hh>
29 # include <mln/data/compare.hh>
30 # include <mln/logical/not.hh>
31 # include <mln/make/image.hh>
32 # include <mln/make/w_window2d_int.hh>
33 # include <mln/topo/skeleton/crest.hh>
34 # include <mln/transform/distance_front.hh>
35 # include <mln/value/int_u8.hh>
36 
37 int main()
38 {
39  using namespace mln;
40 
41  bool ref_dat[][9] = { { 0, 0, 0, 0, 0, 0, 0, 0, 0 },
42  { 0, 0, 0, 0, 0, 1, 0, 0, 0 },
43  { 0, 0, 0, 0, 1, 1, 0, 0, 0 },
44  { 0, 0, 0, 1, 1, 1, 0, 0, 0 },
45  { 0, 0, 0, 0, 1, 0, 0, 0, 0 },
46  { 0, 0, 0, 0, 0, 0, 0, 0, 0 },
47  { 0, 0, 0, 0, 0, 0, 0, 0, 0 } };
48 
49  bool in_dat[][9] = { { 0, 0, 0, 0, 1, 1, 1, 0, 0 },
50  { 0, 0, 0, 1, 1, 1, 1, 0, 0 },
51  { 0, 0, 1, 1, 1, 1, 1, 1, 1 },
52  { 0, 1, 1, 1, 1, 1, 1, 1, 1 },
53  { 0, 0, 1, 1, 1, 1, 1, 0, 0 },
54  { 0, 0, 0, 1, 1, 1, 1, 0, 0 },
55  { 0, 0, 0, 0, 1, 0, 0, 0, 0 } };
56 
57  int vals[] = { 0, 9, 0, 9, 0,
58  9, 6, 4, 6, 9,
59  0, 4, 0, 4, 0, // Values of distances.
60  9, 6, 4, 6, 9,
61  0, 9, 0, 9, 0 };
62 
63  image2d<bool> input = make::image(in_dat);
64  image2d<bool> crest_ref = make::image(ref_dat);
65 
67  dist_map = mln::transform::distance_front(logical::not_(input), c8(),
69  mln_max(value::int_u8));
70  image2d<bool> crest_ima = topo::skeleton::crest(input, dist_map, c8());
71 
72  mln_assertion(crest_ima == crest_ref);
73 }