$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
init_integral_image.hh
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 #ifndef SCRIBO_CORE_INIT_INTEGRAL_IMAGE_HH
28 # define SCRIBO_CORE_INIT_INTEGRAL_IMAGE_HH
29 
33 
34 # include <mln/core/image/image2d.hh>
35 
36 namespace scribo
37 {
38  using namespace mln;
39 
40 
41  template <typename I, typename F>
42  mln_ch_value(I,double)
43  init_integral_image(const Image<I>& input_, F& func);
44 
45 
46 # ifndef MLN_INCLUDE_ONLY
47 
48  namespace internal
49  {
50 
51  inline
52  double square_(const double& val)
53  {
54  double v = static_cast<double>(val);
55  return v * v;
56  }
57 
58  inline
59  double identity_(const double& val)
60  {
61  return static_cast<double>(val);
62  }
63 
64  } // end of namespace scribo::internal
65 
66 
67  // Facade
68 
69  template <typename I, typename F>
70  mln_ch_value(I,double)
71  init_integral_image(const Image<I>& input_, F& func)
72  {
73  mln_trace("scribo::init_integral_image");
74 
75  const I& input = exact(input_);
76  mln_precondition(input.is_valid());
77  mln_precondition(input.domain().pmin() == literal::origin);
78 
79  mln_ch_value(I,double) output;
80  initialize(output, input);
81 
82  unsigned
83  nrows_ = input.nrows(),
84  ncols_ = input.ncols();
85 
86  output.at_(0,0) = func(input.at_(0, 0));
87 
88  for (unsigned row = 1; row < nrows_; ++row)
89  output.at_(row,0) = output.at_(row - 1, 0) + func(input.at_(row, 0));
90 
91  for (unsigned col = 1; col < ncols_; ++col)
92  output.at_(0,col) = output.at_(0, col - 1)
93  + func(input.at_(0, col));
94 
95  for (unsigned row = 1; row < nrows_; ++row)
96  for (unsigned col = 1; col < ncols_; ++col)
97  output.at_(row, col) = output.at_(row - 1, col)
98  + output.at_(row, col - 1)
99  - output.at_(row - 1, col - 1)
100  + func(input.at_(row, col));
101 
102  return output;
103  }
104 
105 #endif // ! MLN_INCLUDE_ONLY
106 
107 } // end of namespace scribo
108 
109 #endif // ! SCRIBO_CORE_INIT_INTEGRAL_IMAGE_HH