Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 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_LINEAR_GAUSSIAN_1D_HH 00027 # define MLN_LINEAR_GAUSSIAN_1D_HH 00028 00034 00035 #include <mln/core/image/image1d.hh> 00036 #include <mln/linear/gaussian_directional_2d.hh> 00037 00038 00039 00040 namespace mln 00041 { 00042 00043 namespace linear 00044 { 00045 00046 00047 template <typename I> 00048 mln_concrete(I) 00049 gaussian_1d(const Image<I>& input, 00050 double sigma, 00051 const mln_value(I)& bdr); 00052 00053 00054 00055 # ifndef MLN_INCLUDE_ONLY 00056 00057 template <typename I> 00058 inline 00059 mln_concrete(I) 00060 gaussian_1d(const Image<I>& input_, 00061 double sigma, 00062 const mln_value(I)& bdr) 00063 { 00064 trace::entering("linear::gaussian_1d"); 00065 00066 typedef mln_site(I) P; 00067 mlc_bool(P::dim == 1)::check(); 00068 00069 const I& input = exact(input_); 00070 mln_precondition(input.is_valid()); 00071 00072 my::recursivefilter_coef_ coef(1.68f, 3.735f, 00073 1.783f, 1.723f, 00074 -0.6803f, -0.2598f, 00075 0.6318f, 1.997f, 00076 sigma, 00077 my::recursivefilter_coef_::DericheGaussian); 00078 00079 extension::adjust_fill(input, 5 * int(sigma + .50001) + 1, bdr); 00080 mln_concrete(I) output = duplicate(input); 00081 00082 if (sigma < 0.006) 00083 return output; 00084 00085 int 00086 ninds = geom::ninds(input), 00087 b = input.border(); 00088 00089 recursivefilter_directional_fastest(output, coef, 00090 point1d(- b), 00091 point1d(ninds - 1 + b), 00092 ninds + 2 * b, 00093 dpoint1d(1), 00094 bdr); 00095 00096 trace::exiting("linear::gaussian_1d"); 00097 return output; 00098 } 00099 00100 # endif // ! MLN_INCLUDE_ONLY 00101 00102 } // end of namespace mln::linear 00103 00104 } // end of namespace mln 00105 00106 00107 #endif // ! MLN_LINEAR_GAUSSIAN_1D_HH