Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2007, 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 #ifndef MLN_MORPHO_LAPLACIAN_HH 00027 # define MLN_MORPHO_LAPLACIAN_HH 00028 00036 # include <mln/morpho/includes.hh> 00037 00038 00039 namespace mln 00040 { 00041 00042 namespace morpho 00043 { 00044 00049 template <typename I, typename W, typename O> 00050 void laplacian(const Image<I>& input, const Window<W>& win, 00051 Image<O>& output); 00052 00053 00054 template <typename I, typename W> 00055 mln_trait_op_minus_twice(mln_concrete(I)) 00056 laplacian(const Image<I>& input, const Window<W>& win); 00057 00058 00059 # ifndef MLN_INCLUDE_ONLY 00060 00061 template <typename I, typename W, typename O> 00062 inline 00063 void laplacian(const Image<I>& input, const Window<W>& win, Image<O>& output) 00064 { 00065 trace::entering("morpho::laplacian"); 00066 mln_precondition(exact(output).domain() == exact(input).domain()); 00067 mln_precondition(! exact(win).is_empty()); 00068 00069 mln_concrete(I) 00070 d_I = morpho::minus(dilation(input, win), input), 00071 e_I = morpho::minus(input, erosion(input, win)); 00072 data::fill(output, d_I - e_I); 00073 00074 trace::exiting("morpho::laplacian"); 00075 } 00076 00077 template <typename I, typename W> 00078 inline 00079 mln_trait_op_minus_twice(mln_concrete(I)) 00080 laplacian(const Image<I>& input, const Window<W>& win) 00081 { 00082 trace::entering("morpho::laplacian"); 00083 mln_precondition(exact(input).is_valid()); 00084 mln_precondition(! exact(win).is_empty()); 00085 00086 mln_trait_op_minus_twice(mln_concrete(I)) output; 00087 initialize(output, input); 00088 laplacian(input, win, output); // Calls previous version. 00089 00090 trace::exiting("morpho::laplacian"); 00091 return output; 00092 } 00093 00094 # endif // ! MLN_INCLUDE_ONLY 00095 00096 } // end of namespace mln::morpho 00097 00098 } // end of namespace mln 00099 00100 00101 #endif // ! MLN_MORPHO_LAPLACIAN_HH