Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2010 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_HISTO_EQUALIZATION_HH 00027 # define MLN_HISTO_EQUALIZATION_HH 00028 00029 # include <mln/core/concept/image.hh> 00030 # include <mln/histo/all.hh> 00031 00035 00036 namespace mln 00037 { 00038 00039 namespace histo 00040 { 00041 00046 template <typename I> 00047 mln_concrete(I) 00048 equalize(const Image<I>& input); 00049 00050 00051 # ifndef MLN_INCLUDE_ONLY 00052 00053 template <typename I> 00054 mln_concrete(I) 00055 equalize(const Image<I>& input_) 00056 { 00057 trace::entering("histo::equalize"); 00058 00059 const I& input = exact(input_); 00060 mln_concrete(I) out_image; 00061 00062 array<mln_value(I)> histogram = compute(input); 00063 array<mln_value(I)> histogram_correction; 00064 00065 unsigned cumulation = 0; 00066 int number_of_pixels = input.nsites(); 00067 //int number_of_colors=histogram.nvalues(); 00068 //int number_of_colors=mln_card(mln_value(I)); 00069 int max_color = mln_max(mln_value(I)); 00070 mln_piter(I) p(input.domain()); 00071 00072 mln_viter(mln::value::set<mln_value(I)>) v(histogram.vset()); 00073 for_all(v) 00074 if (histogram(v) != 0) 00075 { 00076 cumulation += histogram(v); 00077 histogram_correction(v) = (/*number_of_colors-1*/max_color) 00078 * cumulation / number_of_pixels; 00079 } 00080 00081 initialize(out_image, input); 00082 for_all(p) 00083 out_image(p) = histogram_correction[input(p)]; 00084 00085 trace::exiting("histo::equalize"); 00086 return out_image; 00087 } 00088 00089 # endif // ! MLN_INCLUDE_ONLY 00090 00091 } // end of namespace mln::histo 00092 00093 } // end of namespace mln 00094 00095 00096 00097 #endif // ! MLN_HISTO_EQUALIZATION_HH