Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
subsampling.hh
1 // Copyright (C) 2008, 2009 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 #ifndef MLN_SUBSAMPLING_SUBSAMPLING_HH
28 # define MLN_SUBSAMPLING_SUBSAMPLING_HH
29 
33 
34 # include <mln/geom/ncols.hh>
35 # include <mln/geom/nrows.hh>
36 
37 
38 namespace mln
39 {
40 
41  namespace subsampling
42  {
43 
45  template <typename I>
46  inline
47  mln_concrete(I)
48  subsampling(const Image<I>& input,
49  const mln_deduce(I, site, delta)& first_p,
50  const mln_deduce(I, site, coord)& gap);
51 
52 # ifndef MLN_INCLUDE_ONLY
53 
54  namespace impl
55  {
56 
57  //FIXME : add version for every image types.
58 
59  template <typename T>
60  inline
61  mln_concrete(image2d<T>)
62  subsampling_(const image2d<T>& input,
63  const mln_deduce(image2d<T>, site, delta)& first_p,
64  const mln_deduce(image2d<T>, site, coord)& gap)
65  {
66  trace::entering("subsampling::impl::subsampling_");
67  mln_concrete(image2d<T>) output(geom::nrows(input) / gap,
68  geom::ncols(input) / gap);
69 
70  for (def::coord j = geom::min_col(output);
71  j <= geom::max_col(output); ++j)
72  for (def::coord i = geom::min_row(output);
73  i <= geom::max_row(output); ++i)
74  {
75  point2d p1(i, j);
76  point2d p2(first_p[0] + i * gap, first_p[1] + j * gap);
77 
78  output(p1) = input(p2);
79  }
80 
81  trace::exiting("subsampling::impl::subsampling_");
82  return output;
83  }
84 
85  } // end of namespace mln::subsampling::impl
86 
87 
88  template <typename I>
89  inline
90  mln_concrete(I)
91  subsampling(const Image<I>& input,
92  const mln_deduce(I, site, delta)& first_p,
93  const mln_deduce(I, site, coord)& gap)
94  {
95  trace::entering("subsampling::subsampling");
96  mln_precondition(exact(input).is_valid());
97 
98  mln_concrete(I) output(geom::nrows(input) / gap,
99  geom::ncols(input) / gap); // FIXME: only 2d here
100 
101  output = impl::subsampling_(exact(input), first_p, gap);
102 
103  trace::exiting("subsampling::subsampling");
104  return output;
105  }
106 
107 # endif // ! MLN_INCLUDE_ONLY
108 
109  } // end of namespace mln::subsampling
110 
111 } // end of namespace mln
112 
113 
114 #endif // ! MLN_SUBSAMPLING_SUBSAMPLING_HH