Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2008, 2009, 2010 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 <iostream> 00028 00029 #include <mln/core/image/image2d.hh> 00030 #include <mln/core/alias/window2d.hh> 00031 #include <mln/core/alias/neighb2d.hh> 00032 00033 #include <mln/value/int_u8.hh> 00034 #include <mln/value/int_u16.hh> 00035 00036 #include <mln/morpho/watershed/flooding.hh> 00037 #include <mln/data/transform.hh> 00038 00039 #include <mln/io/pgm/load.hh> 00040 #include <mln/io/pgm/save.hh> 00041 00042 #include <mln/util/timer.hh> 00043 00044 #include "tests/data.hh" 00045 00046 struct f_16_to_8 : mln::Function_v2v< f_16_to_8 > 00047 { 00048 typedef mln::value::int_u8 result; 00049 result operator()(const mln::value::int_u16& v) const 00050 { 00051 if (v == 0) 00052 return 0; 00053 return 1 + (v - 1) % 255; 00054 } 00055 }; 00056 00057 00058 int main() 00059 { 00060 using namespace mln; 00061 using value::int_u8; 00062 using value::int_u16; 00063 00064 image2d<int_u8> input; 00065 io::pgm::load(input, MLN_IMG_DIR "/squares.pgm"); 00066 00067 typedef int_u16 L; 00068 L n_basins; 00069 { 00070 util::timer t; 00071 t.start(); 00072 image2d<L> output = 00073 morpho::watershed::impl::generic::flooding(input, c4(), n_basins); 00074 std::cout << "gen: " << t << std::endl; 00075 io::pgm::save(data::transform(output, f_16_to_8()), 00076 "flooding-tmp_ref.pgm"); 00077 } 00078 { 00079 util::timer t; 00080 t.start(); 00081 image2d<L> output = 00082 morpho::watershed::impl::flooding_fastest(input, c4(), n_basins); 00083 std::cout << "fast: " << t << std::endl; 00084 io::pgm::save(data::transform(output, f_16_to_8()), 00085 "flooding-tmp_out.pgm"); 00086 } 00087 }