Milena (Olena)
User documentation 2.0a Id
|
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/value/int_u8.hh> 00029 #include <mln/data/fill.hh> 00030 #include <mln/debug/println.hh> 00031 #include <mln/opt/at.hh> 00032 #include <mln/data/compare.hh> 00033 00034 #include <mln/transform/distance_and_closest_point_geodesic.hh> 00035 00036 unsigned dmap_ref[] = { 4, 4, 4, 3, 2, 3, 4, 4, 4, 00037 4, 4, 3, 2, 1, 2, 3, 4, 4, 00038 4, 3, 2, 1, 0, 1, 2, 3, 4, 00039 3, 2, 1, 2, 1, 2, 3, 4, 4, 00040 2, 1, 0, 1, 2, 3, 4, 4, 4, 00041 3, 2, 1, 2, 3, 4, 4, 4, 4, 00042 4, 3, 2, 3, 4, 4, 4, 4, 4, 00043 4, 4, 3, 4, 4, 4, 4, 4, 4, 00044 4, 4, 4, 4, 4, 4, 4, 4, 4 }; 00045 00046 00047 unsigned cp_idx_ref[] = { 2, 2, 0, 0, 0, 0, 0, 2, 2, 00048 2, 0, 0, 0, 0, 0, 0, 0, 2, 00049 1, 0, 0, 0, 0, 0, 0, 0, 0, 00050 1, 1, 1, 0, 0, 0, 0, 0, 2, 00051 1, 1, 1, 1, 0, 0, 0, 2, 2, 00052 1, 1, 1, 1, 0, 0, 2, 2, 2, 00053 1, 1, 1, 1, 1, 2, 2, 2, 2, 00054 2, 1, 1, 1, 2, 2, 2, 2, 2, 00055 2, 2, 1, 2, 2, 2, 2, 2, 2 }; 00056 00057 00058 int main() 00059 { 00060 using namespace mln; 00061 using value::int_u8; 00062 00063 image2d<bool> input(9, 9); 00064 data::fill(input, false); 00065 opt::at(input, 2, 4) = true; 00066 opt::at(input, 4, 2) = true; 00067 00068 { 00069 typedef util::couple<image2d<int_u8>, image2d<point2d> > output_t; 00070 output_t output = transform::distance_and_closest_point_geodesic(input, 00071 c4(), 00072 int_u8(4)); 00073 mln_assertion(output.first() == make::image2d(dmap_ref)); 00075 } 00076 00077 { 00078 p_array<point2d> arr; 00079 arr.insert(point2d(2, 4)); 00080 arr.insert(point2d(4, 2)); 00081 typedef util::couple<image2d<int_u8>, image2d<unsigned> > output_t; 00082 output_t output = transform::distance_and_closest_point_geodesic(arr, 00083 input.domain(), 00084 c4(), 00085 int_u8(4)); 00086 mln_assertion(output.first() == make::image2d(dmap_ref)); 00087 mln_assertion(output.second() == make::image2d(cp_idx_ref)); 00088 } 00089 }