Milena (Olena): skeleton_constrained.cc Source File

00001 // Copyright (C) 2008, 2009, 2010, 2011 EPITA Research and Development
00002 // Laboratory (LRDE)
00003 //
00004 // This file is part of Olena.
00005 //
00006 // Olena is free software: you can redistribute it and/or modify it under
00007 // the terms of the GNU General Public License as published by the Free
00008 // Software Foundation, version 2 of the License.
00009 //
00010 // Olena is distributed in the hope that it will be useful,
00011 // but WITHOUT ANY WARRANTY; without even the implied warranty of
00012 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00013 // General Public License for more details.
00014 //
00015 // You should have received a copy of the GNU General Public License
00016 // along with Olena.  If not, see <http://www.gnu.org/licenses/>.
00017 //
00018 // As a special exception, you may use this file as part of a free
00019 // software project without restriction.  Specifically, if other files
00020 // instantiate templates or use macros or inline functions from this
00021 // file, or you compile this file and link it with other files to produce
00022 // an executable, this file does not by itself cause the resulting
00023 // executable to be covered by the GNU General Public License.  This
00024 // exception does not however invalidate any other reasons why the
00025 // executable file might be covered by the GNU General Public License.
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 }