00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
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 }
00087
00088 }
00089
00090 }
00091
00092 #endif // ! MLN_WORLD_BINARY_2D_PROJECTED_HISTO_HH