[Olena-patches] milena r2222: Clean morpho::erosion

URL: https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena ChangeLog: 2008-09-11 Matthieu Garrigues &lt;garrigues(a)lrde.epita.fr&gt; Clean morpho::erosion. * mln/morpho/erosion.hh, * mln/morpho/erosion.spe.hh: Clean the specialized versions of erosion. The new versions have the same signature than the facade. Add a new dispatch. * mln/win/rectangle2d.hh: Update an old include --- morpho/erosion.hh | 45 ++---- morpho/erosion.spe.hh | 325 ++++++++++++++------------------------------------ win/rectangle2d.hh | 2 3 files changed, 114 insertions(+), 258 deletions(-) Index: branches/cleanup-2008/milena/mln/win/rectangle2d.hh =================================================================== --- branches/cleanup-2008/milena/mln/win/rectangle2d.hh (revision 2221) +++ branches/cleanup-2008/milena/mln/win/rectangle2d.hh (revision 2222) @@ -35,7 +35,7 @@ * \todo Reactivate includes at EOF. */ -# include <mln/core/window.hh> +# include <mln/core/alias/window2d.hh> namespace mln Index: branches/cleanup-2008/milena/mln/morpho/erosion.spe.hh =================================================================== --- branches/cleanup-2008/milena/mln/morpho/erosion.spe.hh (revision 2221) +++ branches/cleanup-2008/milena/mln/morpho/erosion.spe.hh (revision 2222) @@ -40,31 +40,6 @@ * * \brief Specialization for mln::morpho::erosion. * - * Dispatch of erosion(input, win, output): - * - * If win = rectangle - * -> particular case - * Else If win = octagon - * -> particular case - * Else - * (If win is quite agglomerated (disjonction) - * -> using incremental algorithm - * Else - * -> using conventional algorithm) - * If input value kind = logical - * -> If image is Fast - * -> particular case on sets - * Else - * -> general case on sets - * Else - * If input value are high quantified - * -> using conventional accumulator - * Else - * -> using histogram based accumulator - * If image is Fast - * -> particular case on functions - * Else - * -> general case on functions * */ @@ -82,6 +57,27 @@ erosion(const Image<I>& input, const Window<W>& win); + namespace internal + { + template <typename I, typename W> + mln_concrete(I) + erosion_dispatch(const Image<I>& input, const Window<W>& win) + { + if (mlc_equal(mln_trait_image_kind(I)(), + trait::image::kind::logic)::value == true) + if (mlc_equal(mln_trait_image_speed(I)(), + trait::image::speed::fastest)::value == true) + impl::erosion_on_set_fastest(input, win); + else + impl::generic::erosion_on_set(input, win); + else + if (mlc_equal(mln_trait_image_speed(I)(), + trait::image::speed::fastest)::value == true) + impl::erosion_on_function_fastest(input, win); + else + impl::generic::erosion_on_function(input, win); + } + namespace impl { @@ -91,43 +87,27 @@ // Implementation is in mln/morpho/erosion.hh. template <typename I, typename W> mln_concrete(I) - erosion_on_function_(const I& input, const W& win); + erosion_on_function(const Image<I>& input_, const Window<W>& win_); // Fwd decl. // Implementation is in mln/morpho/erosion.hh. template <typename I, typename W> mln_concrete(I) - erosion_on_set_(const I& input, const W& win); + erosion_on_set(const Image<I>& input_, const Window<W>& win_); } - template <typename I, typename W> - inline mln_concrete(I) - erosion_iterative_(trait::image::kind::any, - trait::image::speed::any, - const I& input, const W& win) + erosion_on_function_fastest(const Image<I>& input_, const Window<W>& win_) { - return generic::erosion_on_function_(input, win); - } - - - template <typename I, typename W> - inline - mln_concrete(I) - erosion_iterative_(trait::image::kind::any, - trait::image::speed::fastest, - const I& input, const W& win) - { - trace::entering("morpho::impl::erosion_iterative_(kind::any, speed::fastest)"); - - border::adjust(input, win.delta()); + trace::entering("morpho::impl::erosion_on_function_fastest"); typedef mln_concrete(I) O; + const I& input = exact(input_); + const W& win = exact(win_); O output; - initialize(output, input); - mln_assertion(border::get(output) == border::get(input)); + initialize(output, input); border::fill(input, mln_max(mln_value(I))); mln_pixter(const I) p(input); @@ -141,230 +121,119 @@ min.take(q.val()); o.val() = min.to_result(); } - trace::exiting("morpho::impl::erosion_iterative_(kind::any, speed::fastest)"); - return output; - } + trace::exiting("morpho::impl::erosion_on_function_fastest"); - template <typename I, typename W> - inline - mln_concrete(I) - erosion_iterative_(trait::image::kind::logic, - trait::image::speed::any, - const I& input, const W& win) - { - return generic::erosion_on_set_(input, win); + return output; } - template <typename I, typename W> - inline mln_concrete(I) - erosion_iterative_(trait::image::kind::logic, - trait::image::speed::fastest, - const I& input, const W& win) + erosion_on_set_fastest(const Image<I>& input_, const Window<W>& win_) { - trace::entering("morpho::impl::erosion_iterative_(kind::logic, speed::fastest)"); + trace::entering("morpho::impl::erosion_on_set_fastest"); - border::adjust(input, win.delta()); - border::fill(input, true); typedef mln_concrete(I) O; + const I& input = exact(input_); + const W& win = exact(win_); O output; - mln_pixter(const I) p(input); - mln_qixter(const I, W) q(p, win); - - if (win.is_centered()) - { - output = clone(input); - - mln_pixter(O) p_out(output); - for_all_2(p, p_out) - // Consider P only if its value is `true' (the erosion - // is anti-extensive). - if (p.val()) - for_all(q) - if (! q.val()) - { - // Assign `false' to the value of P_OUT as soon - // as a `false' value is found in its vicinity. - p_out.val() = false; - break; - } - } - else - { initialize(output, input); - level::fill(output, input); + border::fill(input, true); + mln_pixter(const I) p(input); mln_pixter(O) p_out(output); for_all_2(p, p_out) - for_all(q) { - // Stop as soon as a `false' value is found the - // window. - if (! q.val()) + for_all(q) + if (q.val() == false) break; - // If a `false' value was not found, then Q is now - // invalid (past the end of the window), and vice - // versa. Use this information to assign P_OUT. p_out.val() = ! q.is_valid(); } - } - trace::exiting("morpho::impl::erosion_iterative_(kind::logic, speed::fastest)"); + trace::exiting("morpho::impl::erosion_on_set_fastest"); + return output; } - // Facade. template <typename I, typename W> - inline mln_concrete(I) - erosion_iterative_(const I& input, const W& win) + erosion_on_set_centered(const Image<I>& input_, const Window<W>& win_) { - return erosion_iterative_(mln_trait_image_kind(I)(), - mln_trait_image_speed(I)(), - input, win); - } - + trace::entering("morpho::impl::erosion_on_set_centered"); + typedef mln_concrete(I) O; + const I& input = exact(input_); + const W& win = exact(win_); + O output; + border::fill(input, true); - // // Old code. - - - - // // Stage 5: dispatch w.r.t. fast property - // // | - // // V - - // template <typename I, typename W, typename O> - // void erosion_set_wrt_fast(trait::image::speed::any, const I& input, - // const W& win, O& output) - // { - // generic::erosion_on_set(input, win, output); - // } - - // template <typename I, typename W, typename O> - // void erosion_set_wrt_fast(trait::image::speed::fastest, const I& input, - // const W& win, O& output) - // { - // impl::erosion_on_set_fast(input, win, output); - // } - - // template <typename I, typename W, typename O, typename A> - // void erosion_fun_wrt_fast(trait::image::speed::any, const I& input, - // const W& win, O& output, A& min) - // { - // generic::erosion_on_function(input, win, output, min); - // } - - // template <typename I, typename W, typename O, typename A> - // void erosion_fun_wrt_fast(trait::image::speed::fastest, const I& input, - // const W& win, O& output, A& min) - // { - // impl::erosion_on_function_fast(input, win, output, min); - // } - - // // ^ - // // | - // // end of stage 5 (dispatch w.r.t. fast property) - - // // Stage 4: dispatch w.r.t. data quantification - // // | - // // V - - // template <typename I, typename W, typename O> - // void erosion_wrt_data(trait::image::quant::high, const I& input, - // const W& win, O& output) - // { - // accu::min_<mln_value(I)> min; - // impl::erosion_fun_wrt_fast(mln_trait_image_speed(I)(), input, - // win, output, min); - // } - - // template <typename I, typename W, typename O> - // void erosion_wrt_data(trait::image::quant::low, const I& input, - // const W& win, O& output) - // { - // accu::min_h<mln_vset(I)> min; - // impl::erosion_fun_wrt_fast(mln_trait_image_speed(I)(), input, - // win, output, min); - // } - - // // ^ - // // | - // // end of stage 4 (dispatch w.r.t. the data quantification) - - // // Stage 3: dispatch w.r.t. the value type - // // | - // // V - - // template <typename I, typename W, typename O> - // void erosion_wrt_value(const I& input, const W& win, O& output) - // { - // if (mlc_is(mln_trait_image_kind(I), trait::image::kind::logic)::value) - // impl::erosion_set_wrt_fast(mln_trait_image_speed(I)(), input, - // win, output); - // // | - // // `--> call stage 5: dispatch w.r.t. fast property - // else - // impl::erosion_wrt_data(mln_trait_image_quant(I)(), input, - // win, output); - // // | - // // `--> call stage 4: dispatch w.r.t. the data quantification - // } - - // // ^ - // // | - // // end of stage 3 (dispatch w.r.t. the value type) - - - - // // Stage 2: dispatch w.r.t. the window morphology - // // | - // // V - - // template <typename I, typename W, typename O> - // void erosion_wrt_mor(const I& input, const W& win, O& output) - // { - // // FIXME : Choose the right algorithm between : - // impl::erosion_wrt_value(input, win, output); - // // and : - // // impl::erosion_incr_wrt_value(input, win, output); - // } - - // // ^ - // // | - // // end of stage 2 (dispatch w.r.t. the window morphology) + output = clone(input); + mln_piter(I) p(input.domain()); + mln_qiter(W) q(win, p); + for_all(p) + if (input(p) == true) + for_all(q) if (input.has(q)) + if (input(q) == false) + { + output(p) = false; + break; + } + trace::exiting("morpho::impl::erosion_on_set_centered"); + return output; + } + template <typename I, typename W> + mln_concrete(I) + erosion_on_set_centered_fastest(const Image<I>& input_, const Window<W>& win_) + { + trace::entering("morpho::impl::erosion_on_set_centered_fastest"); + typedef mln_concrete(I) O; + const I& input = exact(input_); + const W& win = exact(win_); + O output; + border::fill(input, true); - template <typename I, typename W> - inline - mln_concrete(I) - erosion_(const I& input, const W& win) + output = clone(input); + mln_pixter(const I) p(input); + mln_qixter(const I, W) q(p, win); + mln_pixter(O) p_out(output); + for_all_2(p, p_out) + if (p.val() == true) + for_all(q) + if (q.val() == false) { - // FIXME - return impl::erosion_iterative_(input, win); + p_out.val() = false; + break; + } + trace::exiting("morpho::impl::erosion_on_set_centered_fastest"); + + return output; } template <typename I> inline mln_concrete(I) - erosion_(const I& input, const win::rectangle2d& win) + erosion_rectangle2d(const Image<I>& input_, const Window<W>& win_) { - trace::entering("morpho::impl::erosion_(win::rectangle2d)"); + mlc_equal(W, win::rectangle2d)::check(); + + trace::entering("morpho::impl::erosion_rectangle2d"); + + const I& input = exact(input_); + const W& win = exact(win_); mln_concrete(I) temp, output; temp = morpho::erosion(input, win::hline2d(win.width())); output = morpho::erosion(temp, win::vline2d(win.height())); - trace::exiting("morpho::impl::erosion_(win::rectangle2d)"); + trace::exiting("morpho::impl::erosion_rectangle2d"); return output; } @@ -372,9 +241,11 @@ template <typename I> inline mln_concrete(I) - erosion_(const I& input, const win::octagon2d& win) + erosion_octagon2d(const I& input, const win::octagon2d& win) { - trace::entering("morpho::impl::erosion_(win::octagon2d)"); + mlc_equal(W, win::octagon2d)::check(); + + trace::entering("morpho::impl::erosion_octagon2d"); const unsigned len = win.length() / 3 + 1; mln_concrete(I) temp_1, temp_2, output; @@ -387,8 +258,6 @@ return output; } - - } // end of namespace mln::morpho::impl } // end of namespace mln::morpho Index: branches/cleanup-2008/milena/mln/morpho/erosion.hh =================================================================== --- branches/cleanup-2008/milena/mln/morpho/erosion.hh (revision 2221) +++ branches/cleanup-2008/milena/mln/morpho/erosion.hh (revision 2222) @@ -63,9 +63,12 @@ template <typename I, typename W> inline mln_concrete(I) - erosion_on_function_(const I& input, const W& win) + erosion_on_function(const Image<I>& input_, const Window<W>& win_); { - trace::entering("morpho::impl::generic::erosion_on_function_"); + trace::entering("morpho::impl::generic::erosion_on_function"); + + const I& input = exact(input_); + const W& win = exact(win_); mln_concrete(I) output; initialize(output, input); @@ -76,12 +79,12 @@ for_all(p) { min.init(); - for_all(q) if (input.owns_(q)) + for_all(q) if (input.has(q)) min.take(input(q)); output(p) = min; } - trace::exiting("morpho::impl::generic::erosion_on_function_"); + trace::exiting("morpho::impl::generic::erosion_on_function"); return output; } @@ -90,41 +93,25 @@ template <typename I, typename W> inline mln_concrete(I) - erosion_on_set_(const I& input, const W& win) + erosion_on_set(const Image<I>& input_, const Window<W>& win_); { - trace::entering("morpho::impl::generic::erosion_on_set_"); + trace::entering("morpho::impl::generic::erosion_on_set"); + + const I& input = exact(input_); + const W& win = exact(win_); mln_piter(I) p(input.domain()); mln_qiter(W) q(win, p); mln_concrete(I) output; - - if (win.is_centered()) - { - output = clone(input); - for_all(p) - if (input(p)) - for_all(q) if (input.owns_(q)) - if (! input(q)) - { - output(p) = false; - break; - } - } - else - { - initialize(output, input); for_all(p) { - for_all(q) if (input.owns_(q)) - if (! input(q)) + for_all(q) if (input.has(q)) + if (input(q) == false) break; - // If there was no break (so q is not valid) then - // output(p) <- true; otherwise, output(p) <- false. output(p) = ! q.is_valid(); } - } - trace::exiting("morpho::impl::generic::erosion_on_set_"); + trace::exiting("morpho::impl::generic::erosion_on_set"); return output; } @@ -145,7 +132,7 @@ mln_precondition(exact(input).has_data()); mln_precondition(! exact(win).is_empty()); - mln_concrete(I) output = impl::erosion_(exact(input), exact(win)); + mln_concrete(I) output = internal::erosion_dispatch(input, win); if (exact(win).is_centered()) mln_postcondition(output <= input);