[Olena-patches] proto-1.0 173: New implementation of geodesic_dilation.hh

2005-04-28 Nicolas Widynski &lt;widyns_n(a)ptaouchnok.lrde.epita.fr&gt; * oln/morpho/geodesic_dilation.hh: New type of implementation, like geodesic_erosion.hh +++ oln/morpho/geodesic_dilation.hh (revision 173) @@ -0,0 +1,111 @@ +// Copyright (C) 2001, 2002, 2003, 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_MORPHO_GEODESIC_DILATION_HH +# define OLENA_MORPHO_GEODESIC_DILATION_HH + +# include <mlc/cmp.hh> +# include <mlc/contract.hh> + +# include <oln/utils/record.hh> +# include <oln/core/pw/all.hh> +# include <oln/core/gen/image_with_nbh.hh> // FIXME: should be in core/abstract/ (?) +# include <oln/morpho/elementary_dilation.hh> + + +// FIXME: such routines should be written somewhere else... + +template <typename L, typename R> +bool operator <= (const oln::abstract::image<L>& lhs, + const oln::abstract::image<R>& rhs) +{ + return oln::pw::check(oln::p_value(lhs) <= oln::p_value(rhs)); +} + + + +namespace oln { + + namespace morpho { + + + + /// Generic geodesic dilation. + + template<typename K, typename I, typename II> + oln_type_of(I, concrete) geodesic_dilation(const tag::kind<K>& kind, + const abstract::image_with_nbh<I>& marker, + const abstract::image<II>& mask) + { + // FIXME: later, test the short code: + // FIXME: return max(elementary_dilation(kind, marker), mask) + + // FIXME: and even later: + // FIXME: return elementary_dilation(kind, marker) _max_ mask + + mlc::eq<oln_type_of(I, grid), oln_type_of(II, grid)>::ensure(); + + entering("morpho::geodesic_dilation"); + + registering(marker, "marker"); + registering(mask, "mask"); + precondition(marker.size() == mask.size() + and marker <= mask); + + oln_type_of(I, concrete) + dil("dil"), + output(marker.size(), "output"); + + dil = elementary_dilation(kind, marker); + + // FIXME: replace code below by a pw::max... + oln_type_of(I, fwd_piter) p(marker.size()); + for_all_p (p) + output[p] = std::min(dil[p].value(), mask[p].value()); + + exiting("morpho::geodesic_dilation"); + return output; + } + + + /// Generic classical geodesic dilation. + + template<typename I, typename II> + oln_type_of(I, concrete) geodesic_dilation(const abstract::image_with_nbh<I>& marker, + const abstract::image<II>& mask) + { + return geodesic_dilation(tag::classical, marker, mask); + } + + + + } // end of namespace oln::morpho + +} // end of namespace oln + + +#endif // OLENA_MORPHO_GEODESIC_DILATION_HH