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_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
00068
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) = (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 }
00092
00093 }
00094
00095
00096
00097 #endif // ! MLN_HISTO_EQUALIZATION_HH