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 <cstdlib> 00028 00029 #include <mln/core/image/image3d.hh> 00030 #include <mln/core/alias/neighb3d.hh> 00031 #include <mln/data/fill.hh> 00032 #include <mln/opt/at.hh> 00033 #include <mln/transform/distance_and_closest_point_geodesic.hh> 00034 #include <mln/value/int_u8.hh> 00035 #include <mln/util/couple.hh> 00036 00037 int main() 00038 { 00039 using namespace mln; 00040 using value::int_u8; 00041 00042 const unsigned 00043 nslices = 100, 00044 nrows = 250, 00045 ncols = 250; 00046 image3d<bool> input(nslices, nrows, ncols); 00047 data::fill(input, false); 00048 for (unsigned i = 0; i < 100; ++i) 00049 opt::at(input, 00050 std::rand() % nslices, 00051 std::rand() % nrows, 00052 std::rand() % ncols) = true; 00053 00054 util::couple<image3d<unsigned>, image3d<point3d> > output = 00055 transform::distance_and_closest_point_geodesic(input, 00056 c6(), 00057 mln_max(unsigned)); 00058 }