generic_salembier_int_u8.cc

00001 // Copyright (C) 2006  EPITA Research and Development Laboratory
00002 //
00003 // This file is part of the Olena Library.  This library is free
00004 // software; you can redistribute it and/or modify it under the terms
00005 // of the GNU General Public License version 2 as published by the
00006 // Free Software Foundation.
00007 //
00008 // This library is distributed in the hope that it will be useful,
00009 // but WITHOUT ANY WARRANTY; without even the implied warranty of
00010 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00011 // General Public License for more details.
00012 //
00013 // You should have received a copy of the GNU General Public License
00014 // along with this library; see the file COPYING.  If not, write to
00015 // the Free Software Foundation, 51 Franklin Street, Fifth Floor,
00016 // Boston, Boston, MA 02110-1301, USA.
00017 //
00018 // As a special exception, you may use this file as part of a free
00019 // software library 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
00022 // produce an executable, this file does not by itself cause the
00023 // resulting executable to be covered by the GNU General Public
00024 // License.  This exception does not however invalidate any other
00025 // reasons why the executable file might be covered by the GNU General
00026 // Public License.
00027 
00028 #include <cstdlib>
00029 
00030 #include <ntg/int.hh>
00031 #include <oln/basics2d.hh>
00032 #include <oln/utils/timer.hh>
00033 
00034 #include <oln/lrde/ufmt/generic_salembier.hh>
00035 
00036 
00037 void usage(char* argv[])
00038 {
00039   std::cerr << "generic max-tree computation with salembier's flooding"
00040             << " (int_u8 version)." << std::endl
00041             << "usage: " << argv[0] << " input.pgm c" << std::endl;
00042   exit(1);
00043 }
00044 
00045 
00046 int main(int argc, char* argv[])
00047 {
00048   if (argc != 3)
00049     usage(argv);
00050 
00051   using namespace oln;
00052   typedef image2d<ntg::int_u8> image_t;
00053 
00054   image_t input = io::load(argv[1]);
00055   assert (input.has_impl());
00056 
00057   int c = atoi(argv[2]);
00058   if (not (c == 2 or c == 4 or c == 8))
00059     usage(argv);
00060 
00061   neighborhood2d nbh;
00062   if (c == 2)
00063     nbh.add(dpoint2d(0, 1));
00064   else
00065     nbh = c == 4 ? neighb_c4() : neighb_c8();
00066 
00067   typedef lrde::ufmt::generic_salembier<image_t> algorithm_t;
00068 
00069   algorithm_t run(input, nbh);
00070   utils::timer t;
00071   t.start();
00072   run.go();
00073   t.stop();
00074   std::cout << "n level roots = " << run.n_level_roots() << std::endl
00075             << "elapsed time  = " << t.last_time() << std::endl;
00076 }

Generated on Tue Feb 20 20:19:43 2007 for Olena by  doxygen 1.5.1