$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
flat_zones.cc
1 // Copyright (C) 2007, 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 <mln/core/image/image2d.hh>
28 #include <mln/value/int_u8.hh>
29 #include <mln/io/pgm/load.hh>
30 #include <mln/core/alias/neighb2d.hh>
31 #include <mln/labeling/flat_zones.hh>
32 
33 #include <mln/labeling/blobs.hh>
34 #include <mln/pw/all.hh>
35 #include <mln/data/compare.hh>
36 
37 #include "tests/data.hh"
38 
39 
40 int main()
41 {
42  using namespace mln;
43  using value::int_u8;
44 
45  image2d<int_u8> lena = io::pgm::load<int_u8>(MLN_IMG_DIR "/tiny.pgm");
46 
47 
48  // Bench code:
49 
50 // for (unsigned i = 0; i < 10; ++i)
51 // {
52 // typedef image2d<int_u8> I;
53 // labeling::impl::flat_zones_functor<I> f(lena);
54 // unsigned n;
55 // canvas::impl::labeling_video_fastest(lena,
56 // c4(),
57 // n,
58 // f);
59 // }
60 
61 
62  unsigned n;
63  image2d<unsigned> labels = labeling::flat_zones(lena, c4(), n);
64  mln_assertion(n == 247);
65 
66  {
67  typedef image2d<int_u8> I;
68  labeling::impl::flat_zones_functor<I,unsigned> f(lena);
69 
70  unsigned nlabels_generic, nlabels_fastest;
72  c4(),
73  nlabels_generic,
74  lena.domain(),
75  f)
76  ==
77  canvas::labeling::impl::video_fastest(lena,
78  c4(),
79  nlabels_fastest,
80  f));
81  mln_assertion(nlabels_generic == nlabels_fastest);
82  }
83 
84  {
85  unsigned n_ = 0;
86  for (unsigned i = 0; i <= 255; ++i)
87  {
88  unsigned n_i;
89  labeling::blobs((pw::value(lena) == pw::cst(i)) | lena.domain(),
90  c4(), n_i);
91  n_ += n_i;
92  }
93  mln_assertion(n_ == n);
94  }
95 
96 }