Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
make_edge_image.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_INTER_PIXEL_DIM2_MAKE_EDGE_IMAGE_HH
27 # define MLN_WORLD_INTER_PIXEL_DIM2_MAKE_EDGE_IMAGE_HH
28 
32 
33 # include <mln/core/image/dmorph/extension_ima.hh>
34 # include <mln/core/image/dmorph/image_if.hh>
35 
36 # include <mln/core/routine/extend.hh>
37 
38 # include <mln/world/inter_pixel/neighb2d.hh>
39 # include <mln/world/inter_pixel/dim2/is_edge.hh>
40 # include <mln/data/fill.hh>
41 
42 namespace mln
43 {
44 
45  namespace world
46  {
47 
48  namespace inter_pixel
49  {
50 
51  namespace dim2
52  {
53 
54  template <typename I, typename F>
55  mln_concrete(I)
56  make_edge_image(const Image<I>& input_, const Function_vv2v<F>& f_)
57  {
58  trace::entering("world::inter_pixel::dim2::make_edge_image");
59 
60  const I& input = exact(input_);
61  const F& f = exact(f_);
62  mln_precondition(input.is_valid());
63 
64  typedef image_if<const image2d<value::int_u<8u> >,
65  world::inter_pixel::dim2::is_edge> edges_t;
66  edges_t edges = input | is_edge();
67  typedef extension_ima<const edges_t,const I> edges_ext_t;
68  edges_ext_t edges_ext = extend(edges, input);
69 
70  mln_piter(edges_ext_t) p(edges_ext.domain());
71  dbl_neighb2d nbh = e2v();
72  mln_niter(dbl_neighb2d) n(nbh, p);
73  mln_value(I) vs[2];
74 
75  mln_concrete(I) output;
76  initialize(output, input);
77  data::fill(output, literal::zero);
78 
79  for_all(p)
80  {
81  unsigned i = 0;
82  for_all(n)
83  {
84  mln_assertion(i < 2);
85  vs[i++] = input(n);
86  }
87  output(p) = f(vs[0], vs[1]);
88  }
89 
90  trace::exiting("world::inter_pixel::dim2::make_edge_image");
91  return output;
92  }
93 
94  } // end of namespace mln::world::inter_pixel::dim2
95 
96  } // end of namespace mln::world::inter_pixel
97 
98  } // end of namespace mln::world
99 
100 } // end of namespace mln
101 
102 #endif // ! MLN_WORLD_INTER_PIXEL_DIM2_MAKE_EDGE_IMAGE_HH