$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
projected_histo.hh
1 // Copyright (C) 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 #ifndef MLN_WORLD_BINARY_2D_PROJECTED_HISTO_HH
27 # define MLN_WORLD_BINARY_2D_PROJECTED_HISTO_HH
28 
32 
33 # include <mln/core/concept/image.hh>
34 # include <mln/core/image/image1d.hh>
35 # include <mln/core/image/image2d.hh>
36 # include <mln/core/alias/dpoint2d.hh>
37 # include <mln/geom/all.hh>
38 # include <mln/opt/at.hh>
39 # include <mln/value/int_u12.hh>
40 # include <mln/util/couple.hh>
41 
42 namespace mln
43 {
44 
45  namespace world
46  {
47 
48  namespace binary_2d
49  {
50 
51  template <typename I>
52  util::couple<image1d<float>, image1d<float> >
53  projected_histo(const Image<I>& input_, bool value = true)
54  {
55  mln_ch_value(I, bool) input = exact(input_);
56 
57  const unsigned nrows = geom::nrows(input);
58  image1d<float> row_histo(nrows);
59  data::fill(row_histo, 0);
60 
61  const unsigned ncols = geom::ncols(input);
62  image1d<float> col_histo(ncols);
63  data::fill(col_histo, 0);
64 
65  for (unsigned row = 0; row < nrows; ++row)
66  {
67  for (unsigned col = 0; col < ncols; ++col)
68  {
69  if (opt::at(input, geom::bbox(input).min_row() + row, geom::bbox(input).min_col() + col) == value)
70  {
71  ++row_histo.at_(row);
72  ++col_histo.at_(col);
73  }
74  }
75  }
76 
77  for (unsigned i = 0; i < row_histo.ninds(); ++i)
78  row_histo.at_(i) = row_histo.at_(i) * 100 / ncols;
79 
80  for (unsigned i = 0; i < col_histo.ninds(); ++i)
81  col_histo.at_(i) = col_histo.at_(i) * 100 / nrows;
82 
83  return util::couple<image1d<float>, image1d<float> > (row_histo, col_histo);
84  }
85 
86  } // end of namespace mln::world::binary_2d
87 
88  } // end of namspace mln::world
89 
90 } // end of namespace mln
91 
92 #endif // ! MLN_WORLD_BINARY_2D_PROJECTED_HISTO_HH