Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2007, 2008, 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_DATA_NAIVE_MEDIAN_HH 00027 # define MLN_DATA_NAIVE_MEDIAN_HH 00028 00032 00033 # include <mln/core/concept/image.hh> 00034 # include <mln/core/alias/window2d.hh> 00035 # include <mln/data/median.hh> 00036 00037 00038 namespace mln 00039 { 00040 00041 namespace data 00042 { 00043 00044 namespace naive 00045 { 00046 00049 00061 template <typename I, typename W, typename O> 00062 void median(const Image<I>& input, const Window<W>& win, 00063 Image<O>& output); 00064 00065 00066 # ifndef MLN_INCLUDE_ONLY 00067 00068 00069 namespace impl 00070 { 00071 00072 template <typename I, typename W, typename O> 00073 inline 00074 void median(const I& input, const W& win, 00075 O& output) 00076 { 00077 mln_piter(I) p(input.domain()); 00078 mln_qiter(W) q(win, p); 00079 accu::stat::median_h<mln_value(I)> med; 00080 00081 for_all(p) 00082 { 00083 med.init(); 00084 for_all(q) if (input.has(q)) 00085 med.take(input(q)); 00086 output(p) = med.to_result(); 00087 } 00088 } 00089 00090 00091 } // end of namespace mln::data::naive::impl 00092 00093 00094 00095 // facade 00096 00097 template <typename I, typename W, typename O> 00098 inline 00099 void median(const Image<I>& input, const Window<W>& win, 00100 Image<O>& output) 00101 { 00102 impl::median(exact(input), exact(win), exact(output)); 00103 } 00104 00105 # endif // ! MLN_INCLUDE_ONLY 00106 00107 } // end of namespace mln::data::naive 00108 00109 } // end of namespace mln::data 00110 00111 } // end of namespace mln 00112 00113 00114 #endif // ! MLN_DATA_NAIVE_MEDIAN_HH