00001 // Copyright (C) 2008, 2009 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #include <mln/core/image/image2d.hh> 00027 #include <mln/core/alias/neighb2d.hh> 00028 #include <mln/pw/image.hh> 00029 00030 #include <mln/debug/println.hh> 00031 #include <mln/debug/iota.hh> 00032 #include <mln/morpho/elementary/dilation.hh> 00033 00034 #include <mln/morpho/tree/utils.hh> 00035 #include <mln/morpho/tree/max.hh> 00036 00037 00038 int main() 00039 { 00040 using namespace mln; 00041 00042 image2d<unsigned> ima(3, 3); 00043 debug::iota(ima); 00044 00045 ima = morpho::elementary::dilation(ima, c8()); 00046 debug::println(ima); 00047 00048 image2d<point2d> par = morpho::tree::max(ima, c4()); 00049 debug::println(par); 00050 00051 00052 // using morpho::tree::is_root; 00053 // using morpho::tree::is_a_node; 00054 00055 // mln_piter_(box2d) p(ima.domain()); 00056 // for_all(p) 00057 // if (is_root(par, p, ima)) 00058 // std::cout << "R "; 00059 // else if (is_a_node(par, p, ima)) 00060 // std::cout << "n "; 00061 // else 00062 // std::cout << ". "; 00063 // std::cout << std::endl; 00064 00065 p_array<point2d> s = data::sort_psites_increasing(ima); 00066 std::cout << morpho::tree::nodes(par, ima, s) << std::endl; 00067 }