$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
w_window3d.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_MAKE_W_WINDOW3D_HH
27 # define MLN_MAKE_W_WINDOW3D_HH
28 
34 # include <cmath>
35 
36 # include <mln/core/w_window.hh>
37 # include <mln/core/alias/dpoint3d.hh>
38 
39 #include <iostream>
40 
41 namespace mln
42 {
43 
44  namespace make
45  {
46 
55  template <typename W, unsigned M>
57 
58 
59 # ifndef MLN_INCLUDE_ONLY
60 
61  template <typename W, unsigned M>
62  inline
64  w_window3d(W (&weights)[M])
65  {
66  int h = unsigned(std::pow(float(M), float(1. / 3))) / 2;
67  mln_precondition((2 * h + 1) * (2 * h + 1) * (2 * h + 1) == M);
69  unsigned i = 0;
70  for (int sli = - h; sli <= h; ++sli)
71  for (int row = - h; row <= h; ++row)
72  for (int col = - h; col <= h; ++col)
73  {
74  const W& cur = weights[i];
75 
76  if (cur != 0)
77  tmp.insert(cur, dpoint3d(sli, row, col));
78  i++;
79  }
80  return tmp;
81  }
82 
83 # endif // ! MLN_INCLUDE_ONLY
84 
85  } // end of namespace mln::make
86 
87 } // end of namespace mln
88 
89 
90 #endif // ! MLN_MAKE_W_WINDOW3D_HH