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