Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2009 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #ifndef MLN_MAKE_P_EDGES_WITH_MASS_CENTERS_HH 00027 # define MLN_MAKE_P_EDGES_WITH_MASS_CENTERS_HH 00028 00036 00037 00038 # include <mln/core/concept/image.hh> 00039 # include <mln/core/concept/graph.hh> 00040 # include <mln/core/site_set/p_edges.hh> 00041 00042 # include <mln/labeling/compute.hh> 00043 00044 # include <mln/accu/center.hh> 00045 00046 # include <mln/fun/i2v/array.hh> 00047 00048 # include <mln/util/site_pair.hh> 00049 00050 00051 namespace mln 00052 { 00053 00054 namespace make 00055 { 00056 00067 // 00068 template <typename W, typename G> 00069 inline 00070 p_edges<G, fun::i2v::array<util::site_pair<mln_site(W)> > > 00071 p_edges_with_mass_centers(const Image<W>& wst_, 00072 const Graph<G>& g_); 00073 00074 00075 # ifndef MLN_INCLUDE_ONLY 00076 00077 00078 template <typename W, typename G> 00079 inline 00080 p_edges<G, fun::i2v::array<util::site_pair<mln_site(W)> > > 00081 p_edges_with_mass_centers(const Image<W>& wst_, 00082 const Graph<G>& g_) 00083 { 00084 trace::entering("make::p_edges_with_mass_centers"); 00085 00086 const W& wst = exact(wst_); 00087 const G& g = exact(g_); 00088 mln_precondition(wst.is_valid()); 00089 mln_precondition(g.is_valid()); 00090 00091 typedef mln_site(W) P; 00092 00093 // nbasins does not include the background label. 00094 mln_value(W) nbasins = g.v_nmax() - 1; 00095 util::array<mln_vec(P)> 00096 mass_centers = labeling::compute(accu::center<P>(), wst, nbasins); 00097 00098 typedef fun::i2v::array< util::site_pair<P> > edge_sites_t; 00099 edge_sites_t edge_sites(g.e_nmax()); 00100 00101 mln_edge_iter(G) e(g); 00102 for_all(e) 00103 edge_sites(e.id()) = util::site_pair<P>(mass_centers[e.v1()], 00104 mass_centers[e.v2()]); 00105 00106 p_edges<G, edge_sites_t> pe(g, edge_sites); 00107 00108 trace::exiting("make::p_edges_with_mass_centers"); 00109 return pe; 00110 } 00111 00112 # endif // ! MLN_INCLUDE_ONLY 00113 00114 } // end of namespace mln::make 00115 00116 } // end of namespace mln 00117 00118 #endif // ! MLN_MAKE_P_EDGES_WITH_MASS_CENTERS_HH