Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
laplacian.hh
1 // Copyright (C) 2007, 2008, 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_MORPHO_LAPLACIAN_HH
27 # define MLN_MORPHO_LAPLACIAN_HH
28 
36 # include <mln/morpho/includes.hh>
37 
38 
39 namespace mln
40 {
41 
42  namespace morpho
43  {
44 
49  template <typename I, typename W, typename O>
50  void laplacian(const Image<I>& input, const Window<W>& win,
51  Image<O>& output);
52 
53 
54  template <typename I, typename W>
55  mln_trait_op_minus_twice(mln_concrete(I))
56  laplacian(const Image<I>& input, const Window<W>& win);
57 
58 
59 # ifndef MLN_INCLUDE_ONLY
60 
61  template <typename I, typename W, typename O>
62  inline
63  void laplacian(const Image<I>& input, const Window<W>& win, Image<O>& output)
64  {
65  trace::entering("morpho::laplacian");
66  mln_precondition(exact(output).domain() == exact(input).domain());
67  mln_precondition(! exact(win).is_empty());
68 
69  mln_concrete(I)
70  d_I = morpho::minus(dilation(input, win), input),
71  e_I = morpho::minus(input, erosion(input, win));
72  data::fill(output, d_I - e_I);
73 
74  trace::exiting("morpho::laplacian");
75  }
76 
77  template <typename I, typename W>
78  inline
79  mln_trait_op_minus_twice(mln_concrete(I))
80  laplacian(const Image<I>& input, const Window<W>& win)
81  {
82  trace::entering("morpho::laplacian");
83  mln_precondition(exact(input).is_valid());
84  mln_precondition(! exact(win).is_empty());
85 
86  mln_trait_op_minus_twice(mln_concrete(I)) output;
87  initialize(output, input);
88  laplacian(input, win, output); // Calls previous version.
89 
90  trace::exiting("morpho::laplacian");
91  return output;
92  }
93 
94 # endif // ! MLN_INCLUDE_ONLY
95 
96  } // end of namespace mln::morpho
97 
98 } // end of namespace mln
99 
100 
101 #endif // ! MLN_MORPHO_LAPLACIAN_HH