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_WORLD_BINARY_2D_PROJECTED_HISTO_HH 00027 # define MLN_WORLD_BINARY_2D_PROJECTED_HISTO_HH 00028 00032 00033 # include <mln/core/concept/image.hh> 00034 # include <mln/core/image/image1d.hh> 00035 # include <mln/core/image/image2d.hh> 00036 # include <mln/core/alias/dpoint2d.hh> 00037 # include <mln/geom/all.hh> 00038 # include <mln/opt/at.hh> 00039 # include <mln/value/int_u12.hh> 00040 # include <mln/util/couple.hh> 00041 00042 namespace mln 00043 { 00044 00045 namespace world 00046 { 00047 00048 namespace binary_2d 00049 { 00050 00051 template <typename I> 00052 util::couple<image1d<float>, image1d<float> > 00053 projected_histo(const Image<I>& input_, bool value = true) 00054 { 00055 mln_ch_value(I, bool) input = exact(input_); 00056 00057 const unsigned nrows = geom::nrows(input); 00058 image1d<float> row_histo(nrows); 00059 data::fill(row_histo, 0); 00060 00061 const unsigned ncols = geom::ncols(input); 00062 image1d<float> col_histo(ncols); 00063 data::fill(col_histo, 0); 00064 00065 for (unsigned row = 0; row < nrows; ++row) 00066 { 00067 for (unsigned col = 0; col < ncols; ++col) 00068 { 00069 if (opt::at(input, geom::bbox(input).min_row() + row, geom::bbox(input).min_col() + col) == value) 00070 { 00071 ++row_histo.at_(row); 00072 ++col_histo.at_(col); 00073 } 00074 } 00075 } 00076 00077 for (unsigned i = 0; i < row_histo.ninds(); ++i) 00078 row_histo.at_(i) = row_histo.at_(i) * 100 / ncols; 00079 00080 for (unsigned i = 0; i < col_histo.ninds(); ++i) 00081 col_histo.at_(i) = col_histo.at_(i) * 100 / nrows; 00082 00083 return util::couple<image1d<float>, image1d<float> > (row_histo, col_histo); 00084 } 00085 00086 } // end of namespace mln::world::binary_2d 00087 00088 } // end of namspace mln::world 00089 00090 } // end of namespace mln 00091 00092 #endif // ! MLN_WORLD_BINARY_2D_PROJECTED_HISTO_HH