$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
fill_holes.hh
1 // Copyright (C) 2007, 2008, 2009, 2011 EPITA Research and Development
2 // Laboratory (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 MLN_LABELING_FILL_HOLES_HH
28 # define MLN_LABELING_FILL_HOLES_HH
29 
33 
34 # include <mln/labeling/background.hh>
35 # include <mln/labeling/compute.hh>
36 
37 # include <mln/data/transform.hh>
38 
39 # include <mln/core/image/dmorph/image_if.hh>
40 # include <mln/accu/math/count.hh>
41 
42 
43 namespace mln
44 {
45 
46  namespace labeling
47  {
48 
62  template <typename I, typename N, typename L>
63  mln_concrete(I)
64  fill_holes(const Image<I>& input, const Neighborhood<N>& nbh,
65  L& nlabels);
66 
67 
68 # ifndef MLN_INCLUDE_ONLY
69 
70  template <typename I, typename N, typename L>
71  inline
72  mln_concrete(I)
73  fill_holes(const Image<I>& input, const Neighborhood<N>& nbh,
74  L& nlabels)
75  {
76  mln_trace("labeling::fill_holes");
77 
78  mlc_equal(mln_trait_image_kind(I),
79  mln::trait::image::kind::binary)::check();
80  mln_precondition(exact(input).is_valid());
81  mln_precondition(exact(nbh).is_valid());
82 
83  mln_ch_value(I, L) lbls = labeling::background(input, nbh, nlabels);
84 
85  accu::math::count<mln_value(I)> a_;
86  util::array<unsigned> arr = labeling::compute(a_, input, lbls, nlabels);
87 
88  unsigned bg_count = 0;
89  unsigned bg_lbl = 0;
90 
91  // We start at 1 to ignore the object.
92  for (unsigned i = 1; i < arr.nelements(); ++i)
93  {
94  if (arr[i] > bg_count)
95  {
96  bg_count = arr[i];
97  bg_lbl = i;
98  }
99  }
100 
101  util::array<bool> bg_relbl(arr.nelements(), true);
102  bg_relbl(bg_lbl) = false;
103  mln_ch_value(I, bool) output = data::transform(lbls, bg_relbl);
104 
105  return output;
106  }
107 
108 # endif // ! MLN_INCLUDE_ONLY
109 
110  } // end of namespace mln::labeling
111 
112 } // end of namespace mln
113 
114 
115 #endif // ! MLN_LABELING_FILL_HOLES_HH