--- milena/sandbox/ChangeLog | 4 ++ milena/sandbox/lazzara/skeleton_crest/main.cc | 52 +++++++++++++++++++++++++ 2 files changed, 56 insertions(+), 0 deletions(-) create mode 100644 milena/sandbox/lazzara/skeleton_crest/main.cc diff --git a/milena/sandbox/ChangeLog b/milena/sandbox/ChangeLog index bed29e7..0fc24be 100644 --- a/milena/sandbox/ChangeLog +++ b/milena/sandbox/ChangeLog @@ -1,3 +1,7 @@ +2010-02-17 Guillaume Lazzara <z@lrde.epita.fr> + + * lazzara/skeleton_crest/main.cc: New tool. + 2010-02-16 Guillaume Lazzara <z@lrde.epita.fr> Add new tools in bin sandbox. diff --git a/milena/sandbox/lazzara/skeleton_crest/main.cc b/milena/sandbox/lazzara/skeleton_crest/main.cc new file mode 100644 index 0000000..e2bf808 --- /dev/null +++ b/milena/sandbox/lazzara/skeleton_crest/main.cc @@ -0,0 +1,52 @@ +# include <mln/core/image/image2d.hh> +# include <mln/core/alias/neighb2d.hh> +# include <mln/topo/skeleton/crest.hh> +# include <mln/topo/skeleton/is_simple_point.hh> +# include <mln/transform/distance_front.hh> +# include <mln/logical/not.hh> +# include <mln/make/w_window2d_int.hh> +# include <mln/value/int_u8.hh> +# include <mln/arith/revert.hh> +# include <mln/io/pgm/save.hh> +# include <mln/io/pgm/load.hh> +# include <mln/io/pbm/save.hh> +# include <mln/io/pbm/load.hh> +# include <mln/morpho/skeleton_constrained.hh> +# include <mln/core/routine/extend.hh> + +int main(int argc, char *argv[]) +{ + using namespace mln; + + image2d<bool> input; + io::pbm::load(input, argv[1]); + + int vals[] = { 0, 9, 0, 9, 0, + 9, 6, 4, 6, 9, + 0, 4, 0, 4, 0, // Values of distances. + 9, 6, 4, 6, 9, + 0, 9, 0, 9, 0 }; + + image2d<value::int_u8> + dist_map_n = transform::distance_front(logical::not_(input), c8(), + make::w_window2d_int(vals), + mln_max(value::int_u8)); + + image2d<value::int_u8> dist_map = arith::revert(dist_map_n); + + io::pgm::save(dist_map, "distance.pgm"); + io::pgm::save(dist_map_n, "distance_n.pgm"); + + // Constraint + image2d<bool> K = topo::skeleton::crest(input, dist_map_n, c8(), atoi(argv[2])); + io::pbm::save(K, "constraint.pbm"); + + typedef image2d<bool> I; + + I skel = + morpho::skeleton_constrained(input, c8(), + topo::skeleton::is_simple_point<I,neighb2d>, + extend(K, false), dist_map); + + io::pbm::save(skel, "skeleton.pbm"); +} -- 1.5.6.5