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 #include <mln/core/image/image2d.hh>
00027 #include <mln/core/alias/neighb2d.hh>
00028 #include <mln/core/var.hh>
00029 #include <mln/pw/all.hh>
00030
00031 #include <mln/value/int_u8.hh>
00032 #include <mln/io/pgm/load.hh>
00033 #include <mln/io/pgm/save.hh>
00034
00035 #include <mln/morpho/tree/component_tree.hh>
00036 #include <mln/morpho/tree/compute_attribute_image.hh>
00037 #include <mln/morpho/attribute/card.hh>
00038
00039 #include <mln/morpho/tree/filter/all.hh>
00040 #include <mln/morpho/tree/propagate_representative.hh>
00041
00042 #include <mln/morpho/closing/area.hh>
00043
00044 #include "tests/data.hh"
00045
00046 #include <iostream>
00047
00048 int main()
00049 {
00050 using namespace mln;
00051 using value::int_u8;
00052
00053
00054 typedef image2d<int_u8> I;
00055 I input;
00056 io::pgm::load(input, MLN_IMG_DIR "/tiny.pgm");
00057
00058 typedef p_array< mln_site_(I) > S;
00059 typedef morpho::tree::data<I,S> tree_t;
00060
00061 tree_t tree = morpho::tree::min_tree(input, c4());
00062
00063
00064 {
00065 typedef morpho::attribute::card<I> attribute_t;
00066 typedef mln_ch_value_(I, unsigned) A;
00067
00068 A a = morpho::tree::compute_attribute_image(attribute_t (), tree);
00069
00070 unsigned lambda = 5;
00071 mln_VAR(predicate, pw::value(a) >= pw::cst(lambda));
00072 I ref = morpho::closing::area(input, c4(), lambda);
00073
00074 {
00075 I out = duplicate(input);
00076 morpho::tree::filter::min(tree, out, predicate);
00077 morpho::tree::propagate_representative(tree, out);
00078 mln_assertion(out == ref);
00079 }
00080
00081 {
00082 I out = duplicate(input);
00083 morpho::tree::filter::max(tree, out, predicate);
00084 morpho::tree::propagate_representative(tree, out);
00085 mln_assertion(out == ref);
00086 }
00087
00088 {
00089 I out = duplicate(input);
00090 morpho::tree::filter::direct(tree, out, predicate);
00091 morpho::tree::propagate_representative(tree, out);
00092 mln_assertion(out == ref);
00093 }
00094
00095
00096 }
00097 }