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
00027 #include <iostream>
00028
00029 #include <mln/core/image/dmorph/image_if.hh>
00030 #include <mln/core/image/image2d.hh>
00031 #include <mln/core/alias/window2d.hh>
00032 #include <mln/core/alias/neighb2d.hh>
00033
00034 #include <mln/value/int_u8.hh>
00035 #include <mln/value/label_16.hh>
00036
00037 #include <mln/morpho/gradient.hh>
00038 #include <mln/morpho/closing/area.hh>
00039 #include <mln/morpho/meyer_wst.hh>
00040
00041 #include <mln/convert/to_window.hh>
00042
00043 #include <mln/pw/cst.hh>
00044 #include <mln/pw/value.hh>
00045
00046 #include <mln/io/pgm/load.hh>
00047 #include <mln/io/ppm/save.hh>
00048 #include <mln/labeling/colorize.hh>
00049
00050 #include "tests/data.hh"
00051
00052
00053 int main()
00054 {
00055 using namespace mln;
00056 using value::int_u8;
00057
00058 image2d<int_u8> input;
00059 io::pgm::load(input, MLN_IMG_DIR "/lena.pgm");
00060
00061 image2d<int_u8>
00062 grad = morpho::gradient(input, c4().win()),
00063 clo = morpho::closing::area(grad, c4(), 200);
00064
00065
00066 typedef value::label_16 L;
00067 L nbasins;
00068 image2d<L> ws = morpho::meyer_wst(clo, c4(), nbasins);
00069
00070
00071 io::ppm::save(labeling::colorize(value::rgb8(), ws, nbasins),
00072 "meyer_wst_long-out.ppm");
00073 }