$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
singh_functor.hh
1 // Copyright (C) 2012, 2013 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_BINARIZATION_INTERNAL_SINGH_FUNCTOR_HH
28 # define SCRIBO_BINARIZATION_INTERNAL_SINGH_FUNCTOR_HH
29 
33 
34 # include <mln/core/image/image2d.hh>
35 # include <mln/core/alias/neighb2d.hh>
36 # include <mln/extension/fill.hh>
37 
38 # include <scribo/binarization/internal/singh_formula.hh>
39 
40 # ifdef SCRIBO_LOCAL_THRESHOLD_DEBUG
41 # include <scribo/binarization/internal/local_threshold_debug.hh>
42 # endif // ! SCRIBO_LOCAL_THRESHOLD_DEBUG
43 
44 
45 namespace scribo
46 {
47 
48  namespace binarization
49  {
50 
51  namespace internal
52  {
53 
54  using namespace mln;
55 
56 
57  template <typename I>
59  {
60  typedef mln_value(I) V;
61 
62  // Moves in input and output images are made using "step"
63  // pixels. It corresponds to the scale ratio between the input
64  // image and the integral image used to give the statistics
65  // values.
66  enum { step = 3 };
67 
68  singh_functor(const Image<I>& input, double K);
69 
70  void init();
71 
72  // Run every 4 pixels.
73  void exec(double mean, double stddev);
74 
75  void end_of_row(int);
76 
77  void finalize();
78 
79 
80  const I input;
81  mln_ch_value(I,bool) output;
82 
83  const V* pi;
84  bool* po;
85 
86  double K_;
87 
89 
90  unsigned next_line3;
91  unsigned offset1;
92  unsigned offset2;
93  };
94 
95 #ifndef MLN_INCLUDE_ONLY
96 
97  template <typename I>
98  singh_functor<I>::singh_functor(const Image<I>& input_, double K)
99  : input(exact(input_)),
100  pi(&input(input.domain().pmin())),
101  K_(K)
102  {
103  mln_precondition(exact(input).is_valid());
104  mln_precondition(K > 0.);
105  }
106 
107  template <typename I>
108  void
109  singh_functor<I>::init()
110  {
111  // This initialization MUST be done here since input image
112  // borders may have changed!
113 
114  // Since we iterate from a smaller image in the largest ones
115  // and image at scale 1 does not always have a size which can
116  // be divided by 3, some sites in the border may not be
117  // processed and we must skip them.
118  int more_offset = - (3 - input.ncols() % 3);
119  if (more_offset == - 3)
120  more_offset = 0; // No offset needed.
121 
122  next_line3 = input.delta_offset(dpoint2d(+2,0))
123  + 2 * input.border() + more_offset;
124 
125  offset1 = input.delta_offset(dpoint2d(+1,0));
126  offset2 = input.delta_offset(dpoint2d(+2,0));
127 
128  initialize(output, input);
129  po = &output(output.domain().pmin());
130  }
131 
132  template <typename I>
133  void
134  singh_functor<I>::exec(double mean, double stddev)
135  {
136  mln_precondition(input.border() == output.border());
137  (void) stddev;
138 
139 
140  for (int i = 0; i < step; ++i, ++po, ++pi)
141  {
142  double th = formula_(*pi, mean, K_);
143  *po = (*pi <= th);
144  th = formula_(*(pi + offset1), mean, K_);
145  *(po + offset1) = (*(pi + offset1) <= th);
146  th = formula_(*(pi + offset2), mean, K_);
147  *(po + offset2) = (*(pi + offset2) <= th);
148  }
149 
150 # ifdef SCRIBO_LOCAL_THRESHOLD_DEBUG
151  // Store local mean
152  unsigned index = pi - input.buffer();
153 
154  debug_mean.element(index) = mean * mean_debug_factor;
155  debug_threshold.element(index) = th;
156 # endif // ! SCRIBO_LOCAL_THRESHOLD_DEBUG
157 
158  }
159 
160  template <typename I>
161  void
162  singh_functor<I>::end_of_row(int)
163  {
164  po += next_line3;
165  pi += next_line3;
166  }
167 
168  template <typename I>
169  void
170  singh_functor<I>::finalize()
171  {
172  }
173 
174 #endif // ! MLN_INCLUDE_ONLY
175 
176  } // end of namespace scribo::binarization::internal
177 
178  } // end of namespace scribo::binarization
179 
180 } // end of namespace scribo
181 
182 #endif // SCRIBO_BINARIZATION_INTERNAL_SINGH_FUNCTOR_HH