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_INTER_PIXEL_COMPUTE_HH
00027 # define MLN_WORLD_INTER_PIXEL_COMPUTE_HH
00028
00032
00033 # include <mln/core/concept/image.hh>
00034 # include <mln/core/concept/function.hh>
00035 # include <mln/world/inter_pixel/is_separator.hh>
00036 # include <mln/world/inter_pixel/separator_to_pixels.hh>
00037
00038
00039 namespace mln
00040 {
00041
00042 namespace world
00043 {
00044
00045 namespace inter_pixel
00046 {
00047
00048
00049 template <typename I, typename F>
00050 image_if<mln_ch_value(mln_unmorph(I), mln_result(F)), is_separator>
00051 compute(const Image<I>& input, const Function_vv2v<F>& f);
00052
00053
00054
00055 # ifndef MLN_INCLUDE_ONLY
00056
00057 template <typename I, typename F>
00058 inline
00059 image_if<mln_ch_value(mln_unmorph(I), mln_result(F)), is_separator>
00060 compute(const Image<I>& input_, const Function_vv2v<F>& f_)
00061 {
00062 trace::entering("world::inter_pixel::compute");
00063
00064 const I& input = exact(input_);
00065 const F& f = exact(f_);
00066
00067 mln_precondition(input.is_valid());
00068
00069 typedef mln_unmorph(I) I_;
00070 typedef mln_ch_value(I_, mln_result(F)) O_;
00071 O_ output_;
00072 initialize(output_, input.unmorph_());
00073
00074 typedef image_if<O_, is_separator> O;
00075 O output(output_, is_separator());
00076
00077 mln_piter(O) e(output.domain());
00078 for_all(e)
00079 {
00080 mln_site(I) p1, p2;
00081 separator_to_pixels(e, p1, p2);
00082 output(e) = f(input(p1), input(p2));
00083 }
00084
00085 trace::exiting("world::inter_pixel::compute");
00086 return output;
00087 }
00088
00089 # endif // ! MLN_INCLUDE_ONLY
00090
00091 }
00092
00093 }
00094
00095 }
00096
00097 #endif // ! MLN_WORLD_INTER_PIXEL_COMPUTE_HH