$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
median_fast.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/win/rectangle2d.hh>
29 
30 #include <mln/io/pgm/load.hh>
31 #include <mln/io/pgm/save.hh>
32 
33 #include <mln/value/int_u8.hh>
34 #include <mln/debug/iota.hh>
35 #include <mln/debug/println.hh>
36 #include <mln/data/fast_median.hh>
37 
38 #include <mln/core/dpoints_pixter.hh>
39 #include <mln/core/pixel.hh>
40 
41 #include "tests/data.hh"
42 
43 
44 namespace mln
45 {
46 
47  template <typename I, typename W>
48  void test(I& input, const W& win)
49  {
50  mln_site(I) p;
51  p.row() = p.col() = 1;
52 
53  {
54  mln_qixter(I, W) qix(input, win, p);
55  for_all(qix)
56  std::cout << qix.val() << ' ';
57  std::cout << " : " << qix.center_val() << std::endl;
58  }
59 
60  {
61  pixel<I> pix(input, p);
62  mln_qixter(I, W) qix(input, win, pix);
63  for_all(qix)
64  std::cout << qix.val() << ' ';
65  std::cout << " : " << qix.center_val() << std::endl;
66  }
67  }
68 
69 }
70 
71 
72 
73 int main()
74 {
75  using namespace mln;
76  using value::int_u8;
77 
78  // FIXME: Dead code.
79 // {
80 // win::rectangle2d rect(3, 3);
81 // border::thickness = 4;
82 // image2d<int_u8> ima(3, 3);
83 // debug::iota(ima);
84 // debug::println(ima);
85 // test(ima, rect);
86 // }
87 
88 
89  {
90  win::rectangle2d rect(51, 51);
91  border::thickness = 52;
92 
93  image2d<int_u8> lena;
94  io::pgm::load(lena, MLN_IMG_DIR "/lena.pgm");
95  image2d<int_u8> out(lena.domain());
96 
97  data::fast_median(lena, rect, out);
98  io::pgm::save(out, "median_fast-out.pgm");
99  }
100 
101 }