Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2007, 2008, 2009, 2011 EPITA Research and Development 00002 // Laboratory (LRDE) 00003 // 00004 // This file is part of Olena. 00005 // 00006 // Olena is free software: you can redistribute it and/or modify it under 00007 // the terms of the GNU General Public License as published by the Free 00008 // Software Foundation, version 2 of the License. 00009 // 00010 // Olena is distributed in the hope that it will be useful, 00011 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00012 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00013 // General Public License for more details. 00014 // 00015 // You should have received a copy of the GNU General Public License 00016 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00017 // 00018 // As a special exception, you may use this file as part of a free 00019 // software project without restriction. Specifically, if other files 00020 // instantiate templates or use macros or inline functions from this 00021 // file, or you compile this file and link it with other files to produce 00022 // an executable, this file does not by itself cause the resulting 00023 // executable to be covered by the GNU General Public License. This 00024 // exception does not however invalidate any other reasons why the 00025 // executable file might be covered by the GNU General Public License. 00026 00027 #ifndef MLN_LABELING_FILL_HOLES_HH 00028 # define MLN_LABELING_FILL_HOLES_HH 00029 00033 00034 # include <mln/labeling/background.hh> 00035 # include <mln/labeling/compute.hh> 00036 00037 # include <mln/data/transform.hh> 00038 00039 # include <mln/core/image/dmorph/image_if.hh> 00040 # include <mln/accu/math/count.hh> 00041 00042 00043 namespace mln 00044 { 00045 00046 namespace labeling 00047 { 00048 00062 template <typename I, typename N, typename L> 00063 mln_concrete(I) 00064 fill_holes(const Image<I>& input, const Neighborhood<N>& nbh, 00065 L& nlabels); 00066 00067 00068 # ifndef MLN_INCLUDE_ONLY 00069 00070 template <typename I, typename N, typename L> 00071 inline 00072 mln_concrete(I) 00073 fill_holes(const Image<I>& input, const Neighborhood<N>& nbh, 00074 L& nlabels) 00075 { 00076 trace::entering("labeling::fill_holes"); 00077 00078 mlc_equal(mln_trait_image_kind(I), 00079 mln::trait::image::kind::binary)::check(); 00080 mln_precondition(exact(input).is_valid()); 00081 mln_precondition(exact(nbh).is_valid()); 00082 00083 mln_ch_value(I, L) lbls = labeling::background(input, nbh, nlabels); 00084 00085 accu::math::count<mln_value(I)> a_; 00086 util::array<unsigned> arr = labeling::compute(a_, input, lbls, nlabels); 00087 00088 unsigned bg_count = 0; 00089 unsigned bg_lbl = 0; 00090 00091 // We start at 1 to ignore the object. 00092 for (unsigned i = 1; i < arr.nelements(); ++i) 00093 { 00094 if (arr[i] > bg_count) 00095 { 00096 bg_count = arr[i]; 00097 bg_lbl = i; 00098 } 00099 } 00100 00101 util::array<bool> bg_relbl(arr.nelements(), true); 00102 bg_relbl(bg_lbl) = false; 00103 mln_ch_value(I, bool) output = data::transform(lbls, bg_relbl); 00104 00105 trace::exiting("labeling::fill_holes"); 00106 return output; 00107 } 00108 00109 # endif // ! MLN_INCLUDE_ONLY 00110 00111 } // end of namespace mln::labeling 00112 00113 } // end of namespace mln 00114 00115 00116 #endif // ! MLN_LABELING_FILL_HOLES_HH