Milena (Olena)
User documentation 2.0a Id
|
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 }