Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2009, 2010 EPITA Research and Development Laboratory 00002 // (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/alias/neighb2d.hh> 00028 # include <mln/core/image/image2d.hh> 00029 # include <mln/data/compare.hh> 00030 # include <mln/logical/not.hh> 00031 # include <mln/make/image.hh> 00032 # include <mln/make/w_window2d_int.hh> 00033 # include <mln/topo/skeleton/crest.hh> 00034 # include <mln/transform/distance_front.hh> 00035 # include <mln/value/int_u8.hh> 00036 00037 int main() 00038 { 00039 using namespace mln; 00040 00041 bool ref_dat[][9] = { { 0, 0, 0, 0, 0, 0, 0, 0, 0 }, 00042 { 0, 0, 0, 0, 0, 1, 0, 0, 0 }, 00043 { 0, 0, 0, 0, 1, 1, 0, 0, 0 }, 00044 { 0, 0, 0, 1, 1, 1, 0, 0, 0 }, 00045 { 0, 0, 0, 0, 1, 0, 0, 0, 0 }, 00046 { 0, 0, 0, 0, 0, 0, 0, 0, 0 }, 00047 { 0, 0, 0, 0, 0, 0, 0, 0, 0 } }; 00048 00049 bool in_dat[][9] = { { 0, 0, 0, 0, 1, 1, 1, 0, 0 }, 00050 { 0, 0, 0, 1, 1, 1, 1, 0, 0 }, 00051 { 0, 0, 1, 1, 1, 1, 1, 1, 1 }, 00052 { 0, 1, 1, 1, 1, 1, 1, 1, 1 }, 00053 { 0, 0, 1, 1, 1, 1, 1, 0, 0 }, 00054 { 0, 0, 0, 1, 1, 1, 1, 0, 0 }, 00055 { 0, 0, 0, 0, 1, 0, 0, 0, 0 } }; 00056 00057 int vals[] = { 0, 9, 0, 9, 0, 00058 9, 6, 4, 6, 9, 00059 0, 4, 0, 4, 0, // Values of distances. 00060 9, 6, 4, 6, 9, 00061 0, 9, 0, 9, 0 }; 00062 00063 image2d<bool> input = make::image(in_dat); 00064 image2d<bool> crest_ref = make::image(ref_dat); 00065 00066 image2d<value::int_u8> 00067 dist_map = mln::transform::distance_front(logical::not_(input), c8(), 00068 make::w_window2d_int(vals), 00069 mln_max(value::int_u8)); 00070 image2d<bool> crest_ima = topo::skeleton::crest(input, dist_map, c8()); 00071 00072 mln_assertion(crest_ima == crest_ref); 00073 }