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 <mln/core/image/image2d.hh>
00028 #include <mln/core/alias/neighb2d.hh>
00029 #include <mln/core/var.hh>
00030 #include <mln/value/int_u8.hh>
00031
00032 #include <mln/make/dual_neighb.hh>
00033 #include <mln/topo/is_simple_2d.hh>
00034 #include <mln/morpho/skeleton_constrained.hh>
00035
00036 #include <mln/io/pbm/load.hh>
00037 #include <mln/data/fill.hh>
00038
00039 #include <mln/pw/value.hh>
00040 #include <mln/core/image/dmorph/image_if.hh>
00041 #include <mln/extension/adjust_duplicate.hh>
00042
00043 #include <mln/logical/not.hh>
00044 #include <mln/arith/revert.hh>
00045 #include <mln/transform/distance_geodesic.hh>
00046
00047 #include <mln/data/compare.hh>
00048
00049 #include "tests/data.hh"
00050
00051
00052 bool skl_ref_[22][18] =
00053 { { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00054 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00055 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00056 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00057 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00058 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00059 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00060 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00061 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00062 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00063 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00064 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00065 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00066 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00067 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00068 { 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 },
00069 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0 },
00070 { 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0 },
00071 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0 },
00072 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0 },
00073 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0, 0 },
00074 { 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0 } };
00075
00076
00077 int main()
00078 {
00079 using namespace mln;
00080 using value::int_u8;
00081
00082 border::thickness = 1;
00083
00084 image2d<bool> pic;
00085 io::pbm::load(pic, MLN_IMG_DIR "/tiny.pbm");
00086
00087 typedef
00088 neighb<win::multiple_size<2u, window<dpoint<grid::square, short int> >,
00089 pw::value_<image2d<bool> > > >
00090 nbh_t;
00091 nbh_t nbh = make::dual_neighb(pic, c4(), c8());
00092
00093
00094 image2d<bool> K(pic.domain());
00095 data::fill(K, false);
00096
00097
00098 image2d<int_u8> dmap = transform::distance_geodesic(logical::not_(pic),
00099 nbh.foreground(),
00100 mln_max(int_u8));
00101 dmap = arith::revert(dmap);
00102
00103 image2d<bool> skl = morpho::impl::generic::skeleton_constrained(pic,
00104 nbh,
00105 mln::topo::is_simple_2d_t<nbh_t>(nbh),
00106 K, dmap);
00107 image2d<bool> skl_ref = make::image(skl_ref_);
00108 mln_assertion(skl == skl_ref);
00109 }