October 2008 - Olena-patches

cleanup-2008 2743: Fix wrong include.
by Ugo Jardonnet 31 Oct '08

31 Oct '08

https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena Index: ChangeLog from Ugo Jardonnet <jardonnet(a)lrde.epita.fr> Fix wrong include. * mln/core/site_set/p_graph.hh, * mln/core/image/graph_image.hh, * mln/display/remove.hh, * mln/display/save.hh, * mln/display/show.hh, * mln/util/greater_point.hh: Include image/image_if.hh instead of image_if_value.hh. core/image/graph_image.hh | 2 +- core/site_set/p_graph.hh | 1 - display/remove.hh | 2 +- display/save.hh | 2 +- display/show.hh | 2 +- util/greater_point.hh | 2 +- 6 files changed, 5 insertions(+), 6 deletions(-) Index: mln/core/site_set/p_graph.hh --- mln/core/site_set/p_graph.hh (revision 2742) +++ mln/core/site_set/p_graph.hh (working copy) @@ -35,7 +35,6 @@ # include <mln/util/graph.hh> # include <mln/util/tracked_ptr.hh> -# include <mln/core/image/graph_psite.hh> # include <mln/core/site_set/p_graph_piter.hh> /* FIXME: This class shares a lot with p_line_graph. Factor as much Index: mln/core/image/graph_image.hh --- mln/core/image/graph_image.hh (revision 2742) +++ mln/core/image/graph_image.hh (working copy) @@ -38,7 +38,7 @@ # include <mln/core/internal/image_primary.hh> # include <mln/algebra/vec.hh> # include <mln/core/site_set/p_graph.hh> -# include <mln/core/image/graph_psite.hh> +# include <mln/core/image/bgraph_psite.hh> # include <mln/value/set.hh> /* FIXME: This class shares a lot with line_graph_image. Factor as Index: mln/display/remove.hh --- mln/display/remove.hh (revision 2742) +++ mln/display/remove.hh (working copy) @@ -35,7 +35,7 @@ */ # include <mln/trait/image_from_grid.hh> -# include <mln/core/image_if_value.hh> +# include <mln/core/image/image_if.hh> # include <mln/core/image/image2d.hh> # include <mln/display/save.hh> Index: mln/display/save.hh --- mln/display/save.hh (revision 2742) +++ mln/display/save.hh (working copy) @@ -35,7 +35,7 @@ */ # include <mln/trait/image_from_grid.hh> -# include <mln/core/image_if_value.hh> +# include <mln/core/image/image_if.hh> # include <mln/core/image/image2d.hh> # include <mln/value/rgb8.hh> # include <mln/level/fill.hh> Index: mln/display/show.hh --- mln/display/show.hh (revision 2742) +++ mln/display/show.hh (working copy) @@ -35,7 +35,7 @@ */ # include <mln/trait/image_from_grid.hh> -# include <mln/core/image_if_value.hh> +# include <mln/core/image/image_if.hh> # include <mln/core/image/image2d.hh> # include <mln/display/save.hh> Index: mln/util/greater_point.hh --- mln/util/greater_point.hh (revision 2742) +++ mln/util/greater_point.hh (working copy) @@ -1,4 +1,4 @@ -:/ Copyright (C) 2005, 2008 EPITA Research and Development Laboratory (LRDE) +// Copyright (C) 2005, 2008 EPITA Research and Development Laboratory (LRDE) // // This file is part of the Olena Library. This library is free // software; you can redistribute it and/or modify it under the terms

1 0

cleanup-2008 2743: Fix wrong include.
by Ugo Jardonnet 30 Oct '08

30 Oct '08

https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena Index: ChangeLog from Ugo Jardonnet <jardonnet(a)lrde.epita.fr> Fix wrong include. * mln/core/site_set/p_graph.hh, * mln/core/image/graph_image.hh, * mln/display/remove.hh, * mln/display/save.hh, * mln/display/show.hh, * mln/util/greater_point.hh: Include image/image_if.hh instead of image_if_value.hh. core/image/graph_image.hh | 2 +- core/site_set/p_graph.hh | 1 - display/remove.hh | 2 +- display/save.hh | 2 +- display/show.hh | 2 +- util/greater_point.hh | 2 +- 6 files changed, 5 insertions(+), 6 deletions(-) Index: mln/core/site_set/p_graph.hh --- mln/core/site_set/p_graph.hh (revision 2742) +++ mln/core/site_set/p_graph.hh (working copy) @@ -35,7 +35,6 @@ # include <mln/util/graph.hh> # include <mln/util/tracked_ptr.hh> -# include <mln/core/image/graph_psite.hh> # include <mln/core/site_set/p_graph_piter.hh> /* FIXME: This class shares a lot with p_line_graph. Factor as much Index: mln/core/image/graph_image.hh --- mln/core/image/graph_image.hh (revision 2742) +++ mln/core/image/graph_image.hh (working copy) @@ -38,7 +38,7 @@ # include <mln/core/internal/image_primary.hh> # include <mln/algebra/vec.hh> # include <mln/core/site_set/p_graph.hh> -# include <mln/core/image/graph_psite.hh> +# include <mln/core/image/bgraph_psite.hh> # include <mln/value/set.hh> /* FIXME: This class shares a lot with line_graph_image. Factor as Index: mln/display/remove.hh --- mln/display/remove.hh (revision 2742) +++ mln/display/remove.hh (working copy) @@ -35,7 +35,7 @@ */ # include <mln/trait/image_from_grid.hh> -# include <mln/core/image_if_value.hh> +# include <mln/core/image/image_if.hh> # include <mln/core/image/image2d.hh> # include <mln/display/save.hh> Index: mln/display/save.hh --- mln/display/save.hh (revision 2742) +++ mln/display/save.hh (working copy) @@ -35,7 +35,7 @@ */ # include <mln/trait/image_from_grid.hh> -# include <mln/core/image_if_value.hh> +# include <mln/core/image/image_if.hh> # include <mln/core/image/image2d.hh> # include <mln/value/rgb8.hh> # include <mln/level/fill.hh> Index: mln/display/show.hh --- mln/display/show.hh (revision 2742) +++ mln/display/show.hh (working copy) @@ -35,7 +35,7 @@ */ # include <mln/trait/image_from_grid.hh> -# include <mln/core/image_if_value.hh> +# include <mln/core/image/image_if.hh> # include <mln/core/image/image2d.hh> # include <mln/display/save.hh> Index: mln/util/greater_point.hh --- mln/util/greater_point.hh (revision 2742) +++ mln/util/greater_point.hh (working copy) @@ -1,4 +1,4 @@ -:/ Copyright (C) 2005, 2008 EPITA Research and Development Laboratory (LRDE) +// Copyright (C) 2005, 2008 EPITA Research and Development Laboratory (LRDE) // // This file is part of the Olena Library. This library is free // software; you can redistribute it and/or modify it under the terms

1 0

cleanup-2008 2741: Add majority accu.
by Ugo Jardonnet 30 Oct '08

30 Oct '08

https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena Index: ChangeLog from Ugo Jardonnet <ugo.jardonnet(a)lrde.epita.fr> Add majority accu. * mln/accu/maj_h.hh: Add majority accumulator. * mln/accu/histo.hh: Remove bad use of size_t. * mln/accu/mean.hh: Temporary fix. * mln/fun/x2x/rotation.hh: Fix matrix parameters. accu/histo.hh | 24 ++--- accu/maj_h.hh | 213 ++++++++++++++++++++++++++++++++++++++++++++++++++++ accu/mean.hh | 9 ++ fun/x2x/rotation.hh | 2 4 files changed, 235 insertions(+), 13 deletions(-) Index: mln/accu/histo.hh --- mln/accu/histo.hh (revision 2740) +++ mln/accu/histo.hh (working copy) @@ -56,7 +56,7 @@ */ template <typename V> struct histo : - public mln::accu::internal::base<const std::vector<std::size_t>& , + public mln::accu::internal::base<const std::vector<unsigned>& , histo<V> > { histo(); @@ -70,16 +70,16 @@ void untake(const argument& t); void init(); - std::size_t operator()(const argument& t) const; - std::size_t operator[](unsigned i) const; + unsigned operator()(const argument& t) const; + unsigned operator[](unsigned i) const; unsigned nvalues() const; - std::size_t sum() const; + unsigned sum() const; /// \} /// Get the value of the accumulator. /// \{ - const std::vector<std::size_t>& vect() const; - const std::vector<std::size_t>& to_result() const; + const std::vector<unsigned>& vect() const; + const std::vector<unsigned>& to_result() const; /// \} const value::set<V>& vset() const; @@ -91,7 +91,7 @@ protected: mln::histo::data<V> h_; - std::size_t sum_; + unsigned sum_; }; template <typename V> @@ -165,7 +165,7 @@ template <typename V> inline - std::size_t + unsigned histo<V>::operator()(const argument& t) const { return h_[h_.vset().index_of(t)]; @@ -173,7 +173,7 @@ template <typename V> inline - std::size_t + unsigned histo<V>::operator[](unsigned i) const { mln_precondition(i < h_.vset().nvalues()); @@ -190,7 +190,7 @@ template <typename V> inline - std::size_t + unsigned histo<V>::sum() const { return sum_; @@ -198,7 +198,7 @@ template <typename V> inline - const std::vector<std::size_t>& + const std::vector<unsigned>& histo<V>::vect() const { return h_.vect(); @@ -206,7 +206,7 @@ template <typename V> inline - const std::vector<std::size_t>& + const std::vector<unsigned>& histo<V>::to_result() const { return this->vect(); Index: mln/accu/maj_h.hh --- mln/accu/maj_h.hh (revision 0) +++ mln/accu/maj_h.hh (revision 0) @@ -0,0 +1,213 @@ +// Copyright (C) 2007, 2008 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, 51 Franklin Street, Fifth Floor, +// 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 MLN_ACCU_MAJ_H_HH +# define MLN_ACCU_MAJ_H_HH + +/*! \file mln/accu/maj_h.hh + * + * \brief Define an accumulator that computes a maj_h. + */ + +# include <mln/core/concept/meta_accumulator.hh> +# include <mln/accu/internal/base.hh> +# include <mln/trait/value_.hh> +# include <mln/util/pix.hh> +# include <mln/accu/histo.hh> +# include <vector> + + +namespace mln +{ + + namespace accu + { + + + /*! \brief Generic maj_h accumulator class. + * + * The parameter \c T is the type of values. + */ + template <typename T> + struct maj_h : public mln::accu::internal::base< const T& , maj_h<T> > + { + typedef T argument; + + maj_h(); + + /// Manipulators. + /// \{ + void init(); + void take(const argument& t); + void untake(const argument& t); + void take(const maj_h<T>& other); + /// \} + + /// Get the value of the accumulator. + const T& to_result() const; + + //operator T () const; + + /// Check whether this accu is able to return a result. + /// Always true here. + bool is_valid() const; + + void update_() const; + + mutable bool valid_; + + const accu::histo<T>& histo() const; + + protected: + + mutable argument t_; + mutable accu::histo<T> h_; + }; + + + template <typename I> struct maj_h< util::pix >; + + + namespace meta + { + + /// Meta accumulator for maj_h. + + struct maj_h : public Meta_Accumulator< maj_h > + { + template <typename T> + struct with + { + typedef accu::maj_h<T> ret; + }; + }; + + } // end of namespace mln::accu::meta + + + +# ifndef MLN_INCLUDE_ONLY + + template <typename T> + inline + maj_h<T>::maj_h() + { + init(); + valid_ = true; + } + + template <typename T> + inline + void + maj_h<T>::init() + { + h_.init(); + } + + template <typename T> + inline + void + maj_h<T>::take(const argument& t) + { + h_.take(t); + + //update return value + if (h_(t) > h_(t_)) + t_ = t; + } + + template <typename T> + inline + void + maj_h<T>::untake(const argument& t) + { + h_.untake(t); + + if (valid_) + valid_ = false; + } + + template <typename T> + inline + void + maj_h<T>::take(const maj_h<T>& other) + { + h_.take(other.h_); + + //FIXME: t_ is wrong then + } + + template <typename T> + inline + void + maj_h<T>::update_() const + { + const std::vector<unsigned>& v = h_.vect(); + + for(unsigned i = 0; i != v.size(); i++) + { + // if nb referents of occurrence i > nb referents of t_ + if (v[i] > h_(t_)) + t_ = h_.vset()[i]; // t_ <- current + } + valid_ = true; + } + + template <typename T> + inline + const T& + maj_h<T>::to_result() const + { + if (not valid_) + update_(); + return t_; + } + + template <typename T> + inline + bool + maj_h<T>::is_valid() const + { + return true; + } + + template <typename V> + inline + const accu::histo<V>& + maj_h<V>::histo() const + { + return h_; + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::accu + +} // end of namespace mln + + +#endif // ! MLN_ACCU_MAJ_H _HH Index: mln/accu/mean.hh --- mln/accu/mean.hh (revision 2740) +++ mln/accu/mean.hh (working copy) @@ -61,6 +61,7 @@ struct mean : public mln::accu::internal::base< M , mean<T,S,M> > { typedef T argument; + typedef M result; mean(); @@ -73,6 +74,7 @@ /// Get the value of the accumulator. M to_result() const; + operator M () const; /// Check whether this accu is able to return a result. /// Always true here. @@ -153,6 +155,13 @@ template <typename T, typename S, typename M> inline + mean<T,S,M>::operator M() const + { + return M(sum_.to_result() / count_.to_result()); + } + + template <typename T, typename S, typename M> + inline bool mean<T,S,M>::is_valid() const { Index: mln/fun/x2x/rotation.hh --- mln/fun/x2x/rotation.hh (revision 2740) +++ mln/fun/x2x/rotation.hh (working copy) @@ -189,7 +189,7 @@ xy + zw, 1.f - x2 - z2, yz - xw, 0, xz - yw, yz + xw, 1.f - x2 - y2, 0, 0, 0, 0, 1}; - this->m_(make::mat<4,4,16,float>(t)); + this->m_(make::mat<4,4>(t)); } template <unsigned n, typename C>

1 0

cleanup-2008 2740: INIM: Classif: Make use of maj_h. Fix Makefile.
by Ugo Jardonnet 30 Oct '08

30 Oct '08

https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena/sandbox Index: ChangeLog from Ugo Jardonnet <ugo.jardonnet(a)lrde.epita.fr> INIM: Classif: Make use of maj_h. Fix Makefile. * classif/proj.hh: Make use of majority. * classif/Makefile: Proper dependencies. Makefile | 21 ++++++++++++++++++--- proj.hh | 10 ++++++---- 2 files changed, 24 insertions(+), 7 deletions(-) Index: classif/proj.hh --- classif/proj.hh (revision 2739) +++ classif/proj.hh (working copy) @@ -32,7 +32,7 @@ #include <mln/io/pgm/save.hh> #include <mln/level/paste.hh> #include <mln/accu/mean.hh> -#include <mln/accu/max.hh> +#include <mln/accu/maj_h.hh> #include <mln/literal/white.hh> namespace mln @@ -59,11 +59,13 @@ { //accu::mean<unsigned, unsigned long, value::int_u8> mean_accu; - accu::max<value::int_u8> max_1;//FIXME: use majoritaire accu + accu::maj_h<value::int_u8> max_1; image2d<value::int_u8> hproj = proj(histo, max_1); - accu::max max_2;//FIXME: use majoritaire accu - image2d proj_class = proj(ws, max_2); + accu::maj_h<value::int_u8> max_2; + image2d<value::int_u8> proj_class = proj(ws, max_2); + + //std::cout << histo; image2d<value::rgb8> out(proj_class.domain()); Index: classif/Makefile --- classif/Makefile (revision 2739) +++ classif/Makefile (working copy) @@ -12,7 +12,7 @@ # V2 -V2_INCLUDES=-I../.. -I../garrigues/fllt +V2_INCLUDES=-I../.. V2_SRC=v2.cc max_tree.hh V2=v2 V2_DBG=v2_dbg @@ -32,11 +32,11 @@ g++ $(V2_INCLUDES) -O1 -DNDEBUG $(V2_SRC) -o $(V2) -$(V2_DBG): $(V2_SRC) +$(V2_DBG): $(V2_SRC) .dep_v2 g++ $(V2_INCLUDES) -O0 -g3 -ggdb3 $(V2_SRC) -o $(V2_DBG) -.PHONY:clean check check-debug valgrind +.PHONY:clean check check-debug valgrind .dep_iccvg .dep_v2 clean: rm -f $(ICCVG) @@ -44,6 +44,7 @@ rm -f $(V2) rm -f $(V2_DBG) rm -f *.log + rm -f .dep_iccvg .dep_v2 check-debug: $(ICCVG_DBG) ./iccvg_dbg $(IMG) $(DIV) $(LAMBDA) $(LOG) @@ -70,3 +71,17 @@ v2-gdb: $(V2_DBG) echo "run ../../img/lena.ppm $(DIV) $(LAMBDA) $(LOG)" > gdb.cmd gdb $(V2_DBG) -x gdb.cmd + + +.dep_iccvg: + g++ $(ICCVG_INCLUDES) -MM iccvg04.cc > $@ + @sed -ir s/iccvg04.cc// .dep_iccvg + @sed -ir s/iccvg04.o/iccvg/ .dep_iccvg + +.dep_v2: + g++ $(V2_INCLUDES) -MM v2.cc > $@ + @sed -ir s/v2.cc// .dep_v2 + @sed -ir s/v2.o/v2/ .dep_v2 + +-include .dep_iccvg +-include .dep_v2 \ No newline at end of file

1 0

cleanup-2008 2739: Add a generic discrete distance transform.
by Thierry Geraud 30 Oct '08

30 Oct '08

https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008 Index: ChangeLog from Thierry Geraud <thierry.geraud(a)lrde.epita.fr> Add a generic discrete distance transform. * milena/mln/transform: New. * milena/mln/transform/distance.hh: New. * milena/tests/transform: New directory. * milena/tests/transform/Makefile.am: New. * milena/tests/transform/distance.cc: New. * milena/sandbox/geraud/skel.cc: Layout. * milena/sandbox/geraud/dmap.cc: Update. Slight change in w_window. * milena/mln/core/w_window.hh (from_to): Prepare new overload. For that, add some new meta-code: * milena/mln/metal/math/root.hh: New. * milena/mln/metal/math/all.hh: Update. * milena/mln/metal/math/sqrt.hh (todo): New. * milena/tests/metal/math/root.cc: New. * milena/tests/metal/math/pow.cc: Remove echo. * milena/tests/metal/math/Makefile.am: Update. * milena/tests/morpho/Makefile.am (SUBDIRS): Fix missing. mln/core/w_window.hh | 17 +++ mln/metal/math/all.hh | 3 mln/metal/math/root.hh | 71 +++++---------- mln/metal/math/sqrt.hh | 2 mln/transform/distance.hh | 197 +++++++++++++++++++++++++++++++++++++++++++ sandbox/geraud/dmap.cc | 33 +++---- sandbox/geraud/skel.cc | 4 tests/metal/math/Makefile.am | 4 tests/metal/math/pow.cc | 5 - tests/metal/math/root.cc | 15 +-- tests/morpho/Makefile.am | 2 tests/transform/Makefile.am | 10 ++ tests/transform/distance.cc | 61 +++++++++++++ 13 files changed, 349 insertions(+), 75 deletions(-) Index: milena/tests/metal/math/root.cc --- milena/tests/metal/math/root.cc (revision 2732) +++ milena/tests/metal/math/root.cc (working copy) @@ -1,4 +1,4 @@ -// Copyright (C) 2007 EPITA Research and Development Laboratory +// Copyright (C) 2007, 2008 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 @@ -25,23 +25,20 @@ // reasons why the executable file might be covered by the GNU General // Public License. -/*! \file tests/metal/math/pow.cc +/*! \file tests/metal/math/root.cc * - * \brief Test on mln::metal::math::pow. + * \brief Test on mln::metal::math::root. */ #include <iostream> #include <mln/core/contract.hh> -#include <mln/metal/math/pow.hh> +#include <mln/metal/math/root.hh> int main() { using namespace mln; - using namespace mln::metal; - int res = metal::math::pow_int<2,3>::value; - mln_assertion(res == 8); - - std::cout << metal::math::pow< int_<2>, int_<3> >::ret().name() << std::endl; + int res = metal::math::root<3,8>::value; + mln_assertion(res == 2); } Index: milena/tests/metal/math/pow.cc --- milena/tests/metal/math/pow.cc (revision 2738) +++ milena/tests/metal/math/pow.cc (working copy) @@ -1,4 +1,4 @@ -// Copyright (C) 2007 EPITA Research and Development Laboratory +// Copyright (C) 2007, 2008 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 @@ -43,5 +43,6 @@ int res = metal::math::pow_int<2,3>::value; mln_assertion(res == 8); - std::cout << metal::math::pow< int_<2>, int_<3> >::ret().name() << std::endl; + std::string s = metal::math::pow< int_<2>, int_<3> >::ret().name(); + mln_assertion(s == "metal::int_<8>"); } Index: milena/tests/metal/math/Makefile.am --- milena/tests/metal/math/Makefile.am (revision 2738) +++ milena/tests/metal/math/Makefile.am (working copy) @@ -3,8 +3,10 @@ include $(top_srcdir)/milena/tests/tests.mk check_PROGRAMS = \ - pow + pow \ + root pow_SOURCES = pow.cc +root_SOURCES = root.cc TESTS = $(check_PROGRAMS) Index: milena/tests/transform/Makefile.am --- milena/tests/transform/Makefile.am (revision 0) +++ milena/tests/transform/Makefile.am (revision 0) @@ -0,0 +1,10 @@ +## Process this file through Automake to create Makefile.in -*- Makefile -*- + +include $(top_srcdir)/milena/tests/tests.mk + +check_PROGRAMS = \ + distance + +distance_SOURCES = distance.cc + +TESTS = $(check_PROGRAMS) Index: milena/tests/transform/distance.cc --- milena/tests/transform/distance.cc (revision 0) +++ milena/tests/transform/distance.cc (revision 0) @@ -0,0 +1,61 @@ +// Copyright (C) 2008 EPITA Research and Development Laboratory (LRDE) +// +// 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, 51 Franklin Street, Fifth Floor, +// 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. + +/*! \file tests/transform/distance.cc + * + * \brief Test on mln::transform::distance. + */ + +#include <mln/core/image/image2d.hh> +#include <mln/core/alias/neighb2d.hh> +#include <mln/make/w_window2d_int.hh> +#include <mln/value/int_u8.hh> +#include <mln/level/fill.hh> +#include <mln/debug/println.hh> + +#include <mln/transform/distance.hh> + + +int main() +{ + using namespace mln; + using value::int_u8; + + image2d<bool> input(9, 9); + level::fill(input, false); + input.at(4, 4) = true; + + int vals[] = { 0, 9, 0, 9, 0, + 9, 6, 4, 6, 9, + 0, 4, 0, 4, 0, + 9, 6, 4, 6, 9, + 0, 9, 0, 9, 0 }; + + image2d<int_u8> output = transform::distance(int_u8(), input, + c4(), make::w_window2d_int(vals)); + debug::println(output); +} Index: milena/tests/morpho/Makefile.am --- milena/tests/morpho/Makefile.am (revision 2738) +++ milena/tests/morpho/Makefile.am (working copy) @@ -2,6 +2,8 @@ include $(top_srcdir)/milena/tests/tests.mk +SUBDIRS = elementary + check_PROGRAMS = \ artificial_line_graph_image_wst \ closing_area \ Index: milena/mln/core/w_window.hh --- milena/mln/core/w_window.hh (revision 2738) +++ milena/mln/core/w_window.hh (working copy) @@ -148,6 +148,10 @@ void from_to(const Image& from, w_window<D,W>& to); + template <typename V, unsigned S, typename D, typename W> + void + from_to(const V (&values)[S], w_window<D,W>& to); + } // end of namespace mln::convert @@ -339,6 +343,19 @@ to.insert(ima(p), convert::to<D>(p)); } +// template <typename V, unsigned S, typename D, typename W> +// void +// from_to(const V (&values)[S], w_window<D,W>& to) +// { +// enum { d = D::dim, +// s = mlc_root(d,S)::value / 2 }; +// metal::bool_<(mlc_pow_int(2 * s + 1, d) == S)>::check(); +// to.clear(); +// D dp; +// dp.set_all(-s); +// FIXME +// } + } // end of namespace mln::convert Index: milena/mln/metal/math/root.hh --- milena/mln/metal/math/root.hh (revision 2732) +++ milena/mln/metal/math/root.hh (working copy) @@ -25,19 +25,20 @@ // reasons why the executable file might be covered by the GNU General // Public License. -#ifndef MLN_METAL_MATH_SQRT_HH -# define MLN_METAL_MATH_SQRT_HH +#ifndef MLN_METAL_MATH_ROOT_HH +# define MLN_METAL_MATH_ROOT_HH -/*! \file mln/metal/math/sqrt.hh +/*! \file mln/metal/math/root.hh * - * \brief Definition of the 'sqrt' static function. + * \brief Definition of the 'nth-root' static function. + * + * \todo Have consistent writing of math meta-routines. */ -# include <mln/metal/bool.hh> -# include <mln/metal/int.hh> +# include <mln/metal/math/pow.hh> -# define mlc_sqrt_int(N) mln::metal::math::sqrt_int<( N )>::value +# define mlc_root(N,X) mln::metal::math::root<( N ),( X )> namespace mln @@ -49,63 +50,43 @@ namespace math { - // sqrt_int<x> - namespace impl { - template <int n, int lo = 1, int hi = n> - struct sqrt_int_ + template <unsigned n, unsigned x, + unsigned lo = 1, unsigned hi = x> + struct root { enum { mid = (lo + hi + 1) / 2, - val_lo = sqrt_int_<n, lo, mid-1>::value, - val_hi = sqrt_int_<n, mid, hi>::value + val_lo = root<n, x, lo, mid-1>::value, + val_hi = root<n, x, mid, hi>::value }; - enum { value = n < mid * mid ? val_lo : val_hi }; + enum { value = x < mlc_pow_int(mid, n) ? val_lo : val_hi }; }; - template<int n, int m> - struct sqrt_int_<n, m, m> - { - enum { value = m }; - }; - - // Entry. - - template <int n, bool b> - struct sqrt_int_if_ : sqrt_int_<n> - { - enum { value_ = sqrt_int_<n>::value, - reminder_ = n - value_ * value_ }; - // FIXME: Check that reminder_ == 0. - }; - - template <int n> - struct sqrt_int_if_< n, false > + template<unsigned n, unsigned x, unsigned m> + struct root<n, x, m, m> { + enum { value = m }; // Found. }; } // end of namespace mln::metal::math::impl - template <int n> - struct sqrt_int : impl::sqrt_int_if_< n, (n >= 0) > + template <unsigned n, unsigned x> + struct root : bool_<(n != 0)>::check_t { + enum { value = impl::root<n,x>::value, + reminder = x - mlc_pow_int(value, n) }; }; - - // sqrt<N> - - template <typename N> - struct sqrt; - - template <int n> - struct sqrt< int_<n> > + template <unsigned n> + struct root<n, 0> : bool_<(n != 0)>::check_t { - typedef int_< sqrt_int<n>::value > ret; + enum { value = 0, + reminder = 0}; }; - } // end of namespace mln::metal::math } // end of namespace mln::metal @@ -113,4 +94,4 @@ } // end of namespace mln -#endif // ! MLN_METAL_MATH_SQRT_HH +#endif // ! MLN_METAL_MATH_ROOT_HH Index: milena/mln/metal/math/all.hh --- milena/mln/metal/math/all.hh (revision 2738) +++ milena/mln/metal/math/all.hh (working copy) @@ -1,4 +1,4 @@ -// Copyright (C) 2007 EPITA Research and Development Laboratory +// Copyright (C) 2007, 2008 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 @@ -54,6 +54,7 @@ # include <mln/metal/math/pow.hh> +# include <mln/metal/math/root.hh> # include <mln/metal/math/sqrt.hh> # include <mln/metal/math/max.hh> // ... Index: milena/mln/metal/math/sqrt.hh --- milena/mln/metal/math/sqrt.hh (revision 2738) +++ milena/mln/metal/math/sqrt.hh (working copy) @@ -31,6 +31,8 @@ /*! \file mln/metal/math/sqrt.hh * * \brief Definition of the 'sqrt' static function. + * + * \todo Rely on nth-root. */ # include <mln/metal/bool.hh> Index: milena/mln/transform/distance.hh --- milena/mln/transform/distance.hh (revision 0) +++ milena/mln/transform/distance.hh (revision 0) @@ -0,0 +1,197 @@ +// Copyright (C) 2008 EPITA Research and Development Laboratory (LRDE) +// +// 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, 51 Franklin Street, Fifth Floor, +// 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 MLN_TRANSFORM_DISTANCE_HH +# define MLN_TRANSFORM_DISTANCE_HH + +/*! \file mln/transform/distance.hh + * + * \brief Discrete distance transform. + */ + +# include <vector> +# include <mln/core/concept/image.hh> +# include <mln/core/concept/neighborhood.hh> +# include <mln/core/concept/weighted_window.hh> +# include <mln/level/fill.hh> +# include <mln/accu/max.hh> + + + +namespace mln +{ + + namespace transform + { + + /*! Discrete distance transform. + * + * FIXME: doc. + */ + template <typename D, typename I, typename N, typename W> + mln_ch_value(I, D) + distance(D, const Image& input, + const Neighborhood<N>& nbh, const Weighted_Window<W>& w_win); + + +# ifndef MLN_INCLUDE_ONLY + + namespace impl + { + + namespace generic + { + + template <typename D, typename I, typename N, typename W> + mln_ch_value(I, D) + distance(D, const Image& input_, + const Neighborhood<N>& nbh_, const Weighted_Window<W>& w_win_) + { + trace::entering("transform::impl::generic::distance"); + + const I& input = exact(input_); + const N& nbh = exact(nbh_); + const W& w_win = exact(w_win_); + + mln_precondition(input.has_data()); + + const D Max = mln_max(D); + + typedef mln_site(I) P; + typedef std::vector bucket_t; + + mln_ch_value(I, D) output; + initialize(output, input); + level::fill(output, Max); + + // Modulus determination. + unsigned mod; + { + mln::accu::max<unsigned> accu; + P p; + mln_qiter(W) q(w_win, p); + for_all(q) + accu.take(q.w()); + mod = accu.to_result() + 1; + } + + std::vector<bucket_t> bucket; + bucket.resize(mod); + unsigned bucket_size = 0; + + // Initialization. + { + mln_piter(I) p(input.domain()); + mln_niter(N) n(nbh, p); + for_all(p) + if (input(p) == true) + { + output(p) = literal::zero; + for_all(n) + if (input.domain().has(n) && input(n) == false) + { + bucket[0].push_back(p); + ++bucket_size; + break; + } + } + } + + // Propagation. + { + P p; + mln_qiter(W) q(w_win, p); + + for (unsigned d = 0; bucket_size != 0; ++d) + { + bucket_t& bucket_d = bucket[d % mod]; + for (unsigned i = 0; i < bucket_d.size(); ++i) + { + p = bucket_d[i]; + if (output(p) < d) + continue; + + for_all(q) + if (output.domain().has(q) && output(q) > d) + { + if (unsigned(Max - q.w()) < d) // Saturation => Stop! + { + // trace::warning... + trace::exiting("transform::impl::generic::distance"); + return output; + } + unsigned d_ = d + q.w(); + + if (d_ < output(q)) + { + output(q) = d_; + bucket[d_ % mod].push_back(q); + ++bucket_size; + } + } + } + bucket_size -= bucket_d.size(); + bucket_d.clear(); + } + } + + trace::exiting("transform::impl::generic::distance"); + return output; + + } + + } // end of namespace mln::transform::impl::generic + + } // end of namespace mln::transform::impl + + + // Facade. + + template <typename D, typename I, typename N, typename W> + inline + mln_ch_value(I, D) + distance(D, const Image& input, + const Neighborhood<N>& nbh, const Weighted_Window<W>& w_win) + { + trace::entering("transform::distance"); + + // FIXME: tests. + mln_ch_value(I, D) output = impl::generic::distance(D(), input, + nbh, w_win); + + trace::exiting("transform::distance"); + return output; + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::transform + +} // end of namespace mln + + +#endif // ! MLN_TRANSFORM_DISTANCE_HH Index: milena/sandbox/geraud/skel.cc --- milena/sandbox/geraud/skel.cc (revision 2738) +++ milena/sandbox/geraud/skel.cc (working copy) @@ -12,7 +12,6 @@ # include <mln/debug/println.hh> - namespace mln { @@ -75,6 +74,9 @@ return output; } +*/ + + } // mln Index: milena/sandbox/geraud/dmap.cc --- milena/sandbox/geraud/dmap.cc (revision 2738) +++ milena/sandbox/geraud/dmap.cc (working copy) @@ -21,14 +21,14 @@ const W& win = exact(win_); mln_precondition(input.has_data()); - typedef std::vector<mln_point(I)> bucket_t; + typedef std::vector<mln_site(I)> bucket_t; // Initialization of distance. mln_ch_value(I, unsigned) distance; initialize(distance, input); // Mod determination. - mln_accu_with_(accu::max, unsigned) max; + accu::max<unsigned> max; for (unsigned i = 0; i < win.size(); ++i) max.take(win.w(i)); unsigned mod = unsigned(max) + 1; @@ -59,7 +59,7 @@ // Propagation. { - mln_point(I) p; + mln_site(I) p; mln_qiter(W) q(win, p); for (unsigned d = 0; bucket_size != 0; ++d) @@ -104,26 +104,27 @@ using namespace mln; using value::int_u8; - const unsigned n = 256; + const unsigned n = 9; image2d<bool> ima(n, n); level::fill(ima, false); ima.at(n / 2, n / 2) = true; - int ws[] = { 3, 2, 3, - 2, 0, 2, - 3, 2, 3 }; - -// int ws[] = { 0, 9, 0, 9, 0, -// 9, 6, 4, 6, 9, -// 0, 4, 0, 4, 0, -// 9, 6, 4, 6, 9, -// 0, 9, 0, 9, 0 }; // coef = 4.1203 +// int ws[] = { 3, 2, 3, +// 2, 0, 2, +// 3, 2, 3 }; + + int ws[] = { 0, 9, 0, 9, 0, + 9, 6, 4, 6, 9, + 0, 4, 0, 4, 0, + 9, 6, 4, 6, 9, + 0, 9, 0, 9, 0 }; // coef = 4.1203 w_window2d_int win = make::w_window2d(ws); image2d<unsigned> d = dmap(ima, win); + debug::println(d); - image2d<int_u8> out(d.domain()); - level::stretch(d, out); - io::pgm::save(out, "out.pgm"); +// image2d<int_u8> out; +// level::stretch(d, out); +// io::pgm::save(out, "out.pgm"); }

1 0

scool r123: Fix the way verbatim C++ is handled
by Maxime van Noppen 30 Oct '08

30 Oct '08

From: Maxime van Noppen <yabo(a)lrde.epita.fr> To: transformers-patches(a)lrde.epita.fr, olena-patches(a)lrde.epita.fr Subject: scool r123: Fix the way verbatim C++ is handled URL: https://svn.lrde.epita.fr/svn/scool/branches/scool-ng ChangeLog: 2008-10-30 Maxime van Noppen <yabo(a)lrde.epita.fr> Fix the way verbatim C++ is handled. * pp-cxx/CxxDecl.str: Handle the CxxProgram constructor. * scl-syn/Lexical.sdf: Fix an ambiguity on layout. * scl-syn/Statement.sdf: Add verbatim C++ as statements. * scoolt/Expression.str: Don't un-double-quote strings. * scoolt/Statement.str: Handle the CxxProgram constructor as a declaration. --- pp-cxx/CxxDecl.str | 3 +++ scl-syn/Lexical.sdf | 1 + scl-syn/Statement.sdf | 4 ++++ scoolt/Expression.str | 3 ++- scoolt/Statement.str | 7 ++++++- 5 files changed, 16 insertions(+), 2 deletions(-) Index: branches/scool-ng/src/pp-cxx/CxxDecl.str =================================================================== --- branches/scool-ng/src/pp-cxx/CxxDecl.str (revision 122) +++ branches/scool-ng/src/pp-cxx/CxxDecl.str (revision 123) @@ -57,6 +57,9 @@ CxxComment(txt) -> H hs=1 [ "//" ~S(txt) ] + CxxTopLevelToAbox: + CxxProgram(c) -> S(c) + CxxBodyToAbox: CxxPublic(b) -> V is=2 [ "public:" V [ ~*<map(CxxTopLevelToAbox)> b ] ] Index: branches/scool-ng/src/scoolt/Expression.str =================================================================== --- branches/scool-ng/src/scoolt/Expression.str (revision 122) +++ branches/scool-ng/src/scoolt/Expression.str (revision 123) @@ -25,7 +25,8 @@ <BaseExpressionToCxx> bexp => cxx_bexp ExpressionToCxx: - String(str) -> <un-double-quote> str + String(str) -> str + //String(str) -> <un-double-quote> str ExpressionToCxx: Sum(e1, e2) -> CxxSum(cxx_e1, cxx_e2) Index: branches/scool-ng/src/scoolt/Statement.str =================================================================== --- branches/scool-ng/src/scoolt/Statement.str (revision 122) +++ branches/scool-ng/src/scoolt/Statement.str (revision 123) @@ -2,6 +2,11 @@ imports libstratego-lib Cxx Scool Declaration +rules + + RawCxx: + Cxx(cxx_code) -> CxxProgram(cxx_code) + strategies - StatementToCxx = DeclarationToCxx + StatementToCxx = DeclarationToCxx <+ RawCxx Index: branches/scool-ng/src/scl-syn/Statement.sdf =================================================================== --- branches/scool-ng/src/scl-syn/Statement.sdf (revision 122) +++ branches/scool-ng/src/scl-syn/Statement.sdf (revision 123) @@ -9,6 +9,10 @@ Declaration -> Statement + Cxx -> CxxStatement + CxxStatement -> Statement + CxxStatement -> FunctionStatement + Expression ";" -> FunctionStatement {cons("ExpressionStatement")} SimpleDeclaration -> FunctionStatement "->" Expression ";" -> FunctionStatement {cons("ReturnStatement")} Index: branches/scool-ng/src/scl-syn/Lexical.sdf =================================================================== --- branches/scool-ng/src/scl-syn/Lexical.sdf (revision 122) +++ branches/scool-ng/src/scl-syn/Lexical.sdf (revision 123) @@ -63,6 +63,7 @@ Id -/- [A-Za-z0-9\_] Asterisk -/- [\/] Slash -/- [\*] + RawCxx -/- [\ \t\n] context-free restrictions LAYOUT? -/- [\ \t\n]

1 0

scool r122: Add reference files for tests on functions
by Maxime van Noppen 30 Oct '08

30 Oct '08

From: Maxime van Noppen <yabo(a)lrde.epita.fr> To: transformers-patches(a)lrde.epita.fr, olena-patches(a)lrde.epita.fr Subject: scool r122: Add reference files for tests on functions URL: https://svn.lrde.epita.fr/svn/scool/branches/scool-ng ChangeLog: 2008-10-30 Maxime van Noppen <yabo(a)lrde.epita.fr> Add reference files for tests on functions. * function/c++_build/call_001.g++: New. * function/c++_build/call_002.g++: New. * function/c++_build/call_003.g++: New. * function/c++_build/call_004.g++: New. * function/c++_build/call_005.g++: New. * function/c++_build/call_006.g++: New. * function/c++_build/declaration_003.g++: New. * function/c++_build/declaration_004.g++: New. * function/c++_build/definition_002.g++: New. * function/c++_ast/call_001.aterm, * function/c++_ast/call_002.aterm, * function/c++_ast/call_003.aterm, * function/c++_ast/call_004.aterm, * function/c++_ast/call_005.aterm, * function/c++_ast/call_006.aterm, * function/c++_ast/declaration_003.aterm, * function/c++_ast/declaration_004.aterm, * function/c++_ast/definition_002.aterm, * function/c++_src/call_001.cc, * function/c++_src/call_002.cc, * function/c++_src/call_003.cc, * function/c++_src/call_004.cc, * function/c++_src/call_005.cc, * function/c++_src/call_006.cc, * function/c++_src/declaration_003.cc, * function/c++_src/declaration_004.cc, * function/c++_src/definition_002.cc, * function/scool_ast/call_001.aterm, * function/scool_ast/call_002.aterm, * function/scool_ast/call_003.aterm, * function/scool_ast/call_004.aterm, * function/scool_ast/call_005.aterm, * function/scool_ast/call_006.aterm, * function/scool_ast/declaration_003.aterm, * function/scool_ast/declaration_004.aterm, * function/scool_ast/definition_002.aterm, * function/scool_src/call_001.scl, * function/scool_src/call_002.scl, * function/scool_src/call_003.scl, * function/scool_src/call_004.scl, * function/scool_src/call_005.scl, * function/scool_src/call_006.scl, * function/scool_src/declaration_003.scl, * function/scool_src/declaration_004.scl, * function/scool_src/definition_002.scl: Here. --- c++_ast/call_001.aterm | 2 +- c++_ast/call_002.aterm | 2 +- c++_ast/call_003.aterm | 2 +- c++_ast/call_004.aterm | 2 +- c++_ast/call_005.aterm | 2 +- c++_ast/call_006.aterm | 2 +- c++_ast/declaration_003.aterm | 2 +- c++_ast/declaration_004.aterm | 2 +- c++_ast/definition_002.aterm | 2 +- c++_src/call_001.cc | 2 +- c++_src/call_002.cc | 4 ++-- c++_src/call_003.cc | 8 ++++++-- c++_src/call_004.cc | 2 ++ c++_src/call_005.cc | 2 ++ c++_src/call_006.cc | 2 ++ c++_src/declaration_003.cc | 1 + c++_src/declaration_004.cc | 1 + c++_src/definition_002.cc | 2 +- scool_ast/call_001.aterm | 2 +- scool_ast/call_002.aterm | 2 +- scool_ast/call_003.aterm | 2 +- scool_ast/call_004.aterm | 2 +- scool_ast/call_005.aterm | 2 +- scool_ast/call_006.aterm | 2 +- scool_ast/declaration_003.aterm | 2 +- scool_ast/declaration_004.aterm | 2 +- scool_ast/definition_002.aterm | 2 +- scool_src/call_001.scl | 2 +- scool_src/call_002.scl | 4 ++-- scool_src/call_003.scl | 10 ++++++++-- scool_src/call_004.scl | 2 ++ scool_src/call_005.scl | 2 ++ scool_src/call_006.scl | 2 ++ scool_src/declaration_003.scl | 2 ++ scool_src/declaration_004.scl | 2 ++ scool_src/definition_002.scl | 2 +- 36 files changed, 58 insertions(+), 30 deletions(-) Index: branches/scool-ng/tests/function/c++_ast/declaration_004.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/declaration_004.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/declaration_004.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("foo"),None,[]),CxxFun([],CxxType("float"),CxxId("foo"),None,[(CxxType("int"),CxxId("arg1"))]),CxxFun([],CxxType("char"),CxxId("foo"),None,[(CxxConstType(CxxRefType(CxxType("toto"))),CxxId("arg1")),(CxxRefType(CxxType("double")),CxxId("arg2"))])]) +CxxProgram([CxxClassDecl([],CxxId("toto"),[]),CxxFun([],CxxType("void"),CxxId("foo"),None,[]),CxxFun([],CxxType("float"),CxxId("foo"),None,[(CxxType("int"),CxxId("arg1"))]),CxxFun([],CxxType("char"),CxxId("foo"),None,[(CxxConstType(CxxRefType(CxxType("toto"))),CxxId("arg1")),(CxxRefType(CxxType("double")),CxxId("arg2"))])]) Index: branches/scool-ng/tests/function/c++_ast/definition_002.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/definition_002.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/definition_002.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("answer"),[],None,[CxxDecl(CxxConstType(CxxType("int")),CxxId("i"))])]) +CxxProgram([CxxFun([],CxxType("void"),CxxId("answer"),[],None,[CxxDecl(CxxConstType(CxxType("int")),CxxId("i"),CxxInt("42"))])]) Index: branches/scool-ng/tests/function/c++_ast/call_001.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/call_001.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/call_001.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([],CxxId("foo"),[]))])]) +CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([],CxxId("doit"),[]))])]) Index: branches/scool-ng/tests/function/c++_ast/call_002.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/call_002.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/call_002.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([],CxxId("foo"),[CxxInt("42")]))])]) +CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[(CxxType("int"),CxxId("i"))],None,[CxxExpStm(CxxFunCall([],CxxId("doit"),[CxxInt("42")]))])]) Index: branches/scool-ng/tests/function/c++_ast/call_003.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/call_003.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/call_003.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([],CxxId("foo"),["bar",CxxInt("42")]))])]) +CxxProgram([CxxProgram("#include <string>\n\nusing std::string;\n"),CxxFun([],CxxType("void"),CxxId("doit"),[(CxxConstType(CxxRefType(CxxType("string"))),CxxId("s")),(CxxType("int"),CxxId("i"))],None,[CxxExpStm(CxxFunCall([],CxxId("doit"),["\"bar\"",CxxInt("42")]))])]) Index: branches/scool-ng/tests/function/c++_ast/call_004.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/call_004.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/call_004.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([CxxType("int")],CxxId("foo"),[]))])]) +CxxProgram([CxxFun([(CxxType("typename"),CxxId("T"))],CxxType("void"),CxxId("foo"),None,[]),CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([CxxType("int")],CxxId("foo"),[]))])]) Index: branches/scool-ng/tests/function/c++_ast/call_005.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/call_005.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/call_005.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([CxxInt("42")],CxxId("foo"),[]))])]) +CxxProgram([CxxFun([(CxxType("int"),CxxId("i"))],CxxType("void"),CxxId("foo"),None,[]),CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([CxxInt("42")],CxxId("foo"),[]))])]) Index: branches/scool-ng/tests/function/c++_ast/call_006.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/call_006.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/call_006.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([CxxType("int"),CxxInt("42")],CxxId("foo"),[CxxInt("42")]))])]) +CxxProgram([CxxFun([(CxxType("typename"),CxxId("T")),(CxxType("int"),CxxId("i"))],CxxType("void"),CxxId("foo"),None,[(CxxType("int"),CxxId("j"))]),CxxFun([],CxxType("void"),CxxId("doit"),[],None,[CxxExpStm(CxxFunCall([CxxType("int"),CxxInt("42")],CxxId("foo"),[CxxInt("42")]))])]) Index: branches/scool-ng/tests/function/c++_ast/declaration_003.aterm =================================================================== --- branches/scool-ng/tests/function/c++_ast/declaration_003.aterm (revision 121) +++ branches/scool-ng/tests/function/c++_ast/declaration_003.aterm (revision 122) @@ -1 +1 @@ -CxxProgram([CxxFun([],CxxType("void"),CxxId("foo"),None,[(CxxType("int"),CxxId("a")),(CxxConstType(CxxRefType(CxxType("bar"))),CxxId("b"))])]) +CxxProgram([CxxClassDecl([],CxxId("bar"),[]),CxxFun([],CxxType("void"),CxxId("foo"),None,[(CxxType("int"),CxxId("a")),(CxxConstType(CxxRefType(CxxType("bar"))),CxxId("b"))])]) Index: branches/scool-ng/tests/function/c++_src/call_001.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/call_001.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/call_001.cc (revision 122) @@ -1,4 +1,4 @@ void doit() { - foo(); + doit(); } Index: branches/scool-ng/tests/function/c++_src/call_002.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/call_002.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/call_002.cc (revision 122) @@ -1,4 +1,4 @@ -void doit() +void doit(int i) { - foo(42); + doit(42); } Index: branches/scool-ng/tests/function/c++_src/call_003.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/call_003.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/call_003.cc (revision 122) @@ -1,4 +1,8 @@ -void doit() +#include <string> + +using std::string; + +void doit(const string& s, int i) { - foo(bar, 42); + doit("bar", 42); } Index: branches/scool-ng/tests/function/c++_src/call_004.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/call_004.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/call_004.cc (revision 122) @@ -1,3 +1,5 @@ +template < typename T > +void foo(); void doit() { foo< int >(); Index: branches/scool-ng/tests/function/c++_src/call_005.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/call_005.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/call_005.cc (revision 122) @@ -1,3 +1,5 @@ +template < int i > +void foo(); void doit() { foo< 42 >(); Index: branches/scool-ng/tests/function/c++_src/call_006.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/call_006.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/call_006.cc (revision 122) @@ -1,3 +1,5 @@ +template < typename T, int i > +void foo(int j); void doit() { foo< int, 42 >(42); Index: branches/scool-ng/tests/function/c++_src/declaration_003.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/declaration_003.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/declaration_003.cc (revision 122) @@ -1 +1,2 @@ +class bar; void foo(int a, const bar& b); Index: branches/scool-ng/tests/function/c++_src/declaration_004.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/declaration_004.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/declaration_004.cc (revision 122) @@ -1,3 +1,4 @@ +class toto; void foo(); float foo(int arg1); char foo(const toto& arg1, double& arg2); Index: branches/scool-ng/tests/function/c++_src/definition_002.cc =================================================================== --- branches/scool-ng/tests/function/c++_src/definition_002.cc (revision 121) +++ branches/scool-ng/tests/function/c++_src/definition_002.cc (revision 122) @@ -1,4 +1,4 @@ void answer() { - const int i; + const int i = 42; } Index: branches/scool-ng/tests/function/scool_ast/declaration_004.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/declaration_004.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/declaration_004.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void"))),FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("arg1"),SimpleType("int"))]),SimpleType("float"))),FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("arg1"),SimpleType("toto")),TypedId(Identifier("arg2"),ReferenceType(SimpleType("double")))]),SimpleType("char")))]) +Program([ClassDeclaration(Some("final"),Identifier("toto"),None),FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void"))),FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("arg1"),SimpleType("int"))]),SimpleType("float"))),FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("arg1"),SimpleType("toto")),TypedId(Identifier("arg2"),ReferenceType(SimpleType("double")))]),SimpleType("char")))]) Index: branches/scool-ng/tests/function/scool_ast/definition_002.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/definition_002.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/definition_002.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("answer"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([SimpleDeclaration(None,Identifier("i"),SimpleType("int"),None)]))]) +Program([FunctionDefinition(None,Identifier("answer"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([SimpleDeclaration(None,Identifier("i"),SimpleType("int"),Some(Initialiser(Integer("42"))))]))]) Index: branches/scool-ng/tests/function/scool_ast/call_001.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/call_001.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/call_001.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),None,Arguments([])))]))]) +Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("doit"),None,Arguments([])))]))]) Index: branches/scool-ng/tests/function/scool_ast/call_002.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/call_002.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/call_002.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),None,Arguments([Integer("42")])))]))]) +Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("i"),SimpleType("int"))]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("doit"),None,Arguments([Integer("42")])))]))]) Index: branches/scool-ng/tests/function/scool_ast/call_003.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/call_003.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/call_003.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),None,Arguments([String("\"bar\""),Integer("42")])))]))]) +Program([Cxx("#include <string>\n\nusing std::string;\n"),FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("s"),SimpleType("string")),TypedId(Identifier("i"),SimpleType("int"))]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("doit"),None,Arguments([String("\"bar\""),Integer("42")])))]))]) Index: branches/scool-ng/tests/function/scool_ast/call_004.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/call_004.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/call_004.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),Some(Parameters([SimpleType("int")])),Arguments([])))]))]) +Program([FunctionDeclaration(None,Identifier("foo"),FunctionType(Some(ParametersDeclaration([TypedId(Identifier("T"),"type")])),ArgumentsDeclaration([]),SimpleType("void"))),FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),Some(Parameters([SimpleType("int")])),Arguments([])))]))]) Index: branches/scool-ng/tests/function/scool_ast/call_005.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/call_005.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/call_005.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),Some(Parameters([Integer("42")])),Arguments([])))]))]) +Program([FunctionDeclaration(None,Identifier("foo"),FunctionType(Some(ParametersDeclaration([TypedId(Identifier("i"),SimpleType("int"))])),ArgumentsDeclaration([]),SimpleType("void"))),FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),Some(Parameters([Integer("42")])),Arguments([])))]))]) Index: branches/scool-ng/tests/function/scool_ast/call_006.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/call_006.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/call_006.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),Some(Parameters([SimpleType("int"),Integer("42")])),Arguments([Integer("42")])))]))]) +Program([FunctionDeclaration(None,Identifier("foo"),FunctionType(Some(ParametersDeclaration([TypedId(Identifier("T"),"type"),TypedId(Identifier("i"),SimpleType("int"))])),ArgumentsDeclaration([TypedId(Identifier("j"),SimpleType("int"))]),SimpleType("void"))),FunctionDefinition(None,Identifier("doit"),FunctionType(None,ArgumentsDeclaration([]),SimpleType("void")),FunctionBlock([ExpressionStatement(FunctionCall(Identifier("foo"),Some(Parameters([SimpleType("int"),Integer("42")])),Arguments([Integer("42")])))]))]) Index: branches/scool-ng/tests/function/scool_ast/declaration_003.aterm =================================================================== --- branches/scool-ng/tests/function/scool_ast/declaration_003.aterm (revision 121) +++ branches/scool-ng/tests/function/scool_ast/declaration_003.aterm (revision 122) @@ -1 +1 @@ -Program([FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("a"),SimpleType("int")),TypedId(Identifier("b"),SimpleType("bar"))]),SimpleType("void")))]) +Program([ClassDeclaration(Some("final"),Identifier("bar"),None),FunctionDeclaration(None,Identifier("foo"),FunctionType(None,ArgumentsDeclaration([TypedId(Identifier("a"),SimpleType("int")),TypedId(Identifier("b"),SimpleType("bar"))]),SimpleType("void")))]) Index: branches/scool-ng/tests/function/scool_src/call_005.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/call_005.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/call_005.scl (revision 122) @@ -1,3 +1,5 @@ +foo : [i : int]() -> void; + doit : () -> void = { foo[42](); Index: branches/scool-ng/tests/function/scool_src/call_006.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/call_006.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/call_006.scl (revision 122) @@ -1,3 +1,5 @@ +foo : [T : type, i : int](j : int) -> void; + doit : () -> void = { foo[int, 42](42); Index: branches/scool-ng/tests/function/scool_src/declaration_003.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/declaration_003.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/declaration_003.scl (revision 122) @@ -1 +1,3 @@ +final bar : class; + foo : (a : int, b : bar) -> void; Index: branches/scool-ng/tests/function/scool_src/declaration_004.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/declaration_004.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/declaration_004.scl (revision 122) @@ -1,3 +1,5 @@ +final toto : class; + foo : () -> void; foo : (arg1 : int) -> float; foo : (arg1 : toto, arg2 : ref double) -> char; Index: branches/scool-ng/tests/function/scool_src/definition_002.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/definition_002.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/definition_002.scl (revision 122) @@ -1,4 +1,4 @@ answer : () -> void = { - i : int; + i : int = 42; } Index: branches/scool-ng/tests/function/scool_src/call_001.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/call_001.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/call_001.scl (revision 122) @@ -1,4 +1,4 @@ doit : () -> void = { - foo(); + doit(); } Index: branches/scool-ng/tests/function/scool_src/call_002.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/call_002.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/call_002.scl (revision 122) @@ -1,4 +1,4 @@ -doit : () -> void = +doit : (i : int) -> void = { - foo(42); + doit(42); } Index: branches/scool-ng/tests/function/scool_src/call_003.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/call_003.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/call_003.scl (revision 122) @@ -1,4 +1,10 @@ -doit : () -> void = +|[ +#include <string> + +using std::string; +]| + +doit : (s : string, i : int) -> void = { - foo("bar", 42); + doit("bar", 42); } Index: branches/scool-ng/tests/function/scool_src/call_004.scl =================================================================== --- branches/scool-ng/tests/function/scool_src/call_004.scl (revision 121) +++ branches/scool-ng/tests/function/scool_src/call_004.scl (revision 122) @@ -1,3 +1,5 @@ +foo : [T : type] () -> void; + doit : () -> void = { foo[int](); Index: branches/scool-ng/tests/function/c++_build/declaration_003.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/declaration_004.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/definition_002.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/call_001.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/call_002.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/call_003.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/call_004.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/call_005.g++ =================================================================== Index: branches/scool-ng/tests/function/c++_build/call_006.g++ ===================================================================

1 0

cleanup-2008 2738: Add co-occurence matrix computing.
by Alexandre Abraham 30 Oct '08

30 Oct '08

https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena/sandbox Index: ChangeLog from Alexandre Abraham <abraham(a)lrde.epita.fr> Add co-occurence matrix computing. * nature/co_occurence.hh: New (mln::co_occurence) : compute co-occurence matrix. * nature/mco.cc: New Dummy test. co_occurence.hh | 75 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ mco.cc | 45 +++++++++++++++++++++++++++++++++ 2 files changed, 120 insertions(+) Index: nature/co_occurence.hh --- nature/co_occurence.hh (revision 0) +++ nature/co_occurence.hh (revision 0) @@ -0,0 +1,75 @@ +// Copyright (C) 2007, 2008 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, 51 Franklin Street, Fifth Floor, +// 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 MLN_CO_OCCURENCE_HH +# define MLN_CO_OCCURENCE_HH + +/*! \file ?? + * + * \brief FIXME + */ + +# include <mln/core/concept/image.hh> +# include <mln/core/dpoint.hh> +# include <mln/core/image/image2d.hh> +# include <mln/level/fill.hh> + + +namespace mln +{ + + + template <typename I, typename D> + // image2d<unsigned> co_occurence (const Image&, const Dpoint<D>&); + image2d<unsigned> co_occurence (const Image&, const Gdpoint<D>&); + +# ifndef MLN_INCLUDE_ONLY + + template <typename I, typename D> + image2d<unsigned> co_occurence (const Image &ima_, const Gdpoint<D> &dp_) + { + mln::metal::equal<mln_psite(I), mln_psite(D)>::check(); + + const I &ima = exact(ima_); + const D &dp = exact(dp_); + image2d<unsigned> mco(mln_card(mln_value(I)), mln_card(mln_value(I)), 0); + level::fill(mco, 0); + + + mln_piter(I) p(ima.domain()); + for_all(p) + if (ima.domain().has(p + dp)) + mco(point2d(ima(p), ima(p + dp)))++; + + return mco; + } + +#endif // MLN_INCLUDE_ONLY + +} + +#endif // MLN_CO_OCCURENCE_HH Index: nature/mco.cc --- nature/mco.cc (revision 0) +++ nature/mco.cc (revision 0) @@ -0,0 +1,45 @@ +#include <iostream> + +#include <mln/value/int_u.hh> +#include <mln/core/image/image2d.hh> +#include "co_occurence.hh" +#include <mln/core/alias/dpoint2d.hh> + +int main () +{ + using namespace mln; + using namespace value; + + typedef image2d< int_u<3> > I; + + int_u<3> vs[6][5] = { + + { 3, 3, 4, 4, 4 }, + { 2, 1, 1, 1, 1 }, + { 1, 4, 4, 4, 6 }, + { 1, 4, 3, 4, 1 }, + { 7, 4, 5, 3, 1 }, + { 7, 7, 1, 1, 0 } + + }; + + I ima(make::image2d(vs)); + + dpoint2d d(0, 1); + + // std::cout << co_occurence(ima, d) << std::endl; + + image2d<unsigned> co(co_occurence(ima, d)); + + unsigned cpt = 0; + + mln_piter_(image2d<unsigned>) p(co.domain()); + for_all(p) + { + std::cout << p << " = " << co(p) << std::endl; + cpt += co(p); + } + + std::cout << cpt << std::endl; + +}

1 0

r2737: Add skeleton algorithm
by Jimmy Ma 30 Oct '08

30 Oct '08

From: Jimmy Ma <jimmy.ma(a)lrde.epita.fr> To: olena-patches(a)lrde.epita.fr Subject: r2737: Add skeleton algorithm URL: https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena/sandbox ChangeLog: 2008-10-30 Jimmy Ma <jimmy.ma(a)lrde.epita.fr> Add skeleton algorithm. * garrigues/ocr/skel.cc: New. --- skel.cc | 155 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 155 insertions(+) Index: branches/cleanup-2008/milena/sandbox/garrigues/ocr/skel.cc =================================================================== --- branches/cleanup-2008/milena/sandbox/garrigues/ocr/skel.cc (revision 0) +++ branches/cleanup-2008/milena/sandbox/garrigues/ocr/skel.cc (revision 2737) @@ -0,0 +1,155 @@ +# include <mln/core/var.hh> + +# include <mln/core/image/image2d.hh> +# include <mln/core/alias/neighb2d.hh> +# include <mln/core/site_set/p_queue_fast.hh> +# include <mln/core/site_set/p_priority.hh> + +# include <mln/io/pbm/load.hh> +# include <mln/io/pbm/save.hh> +# include <mln/value/int_u8.hh> + +# include <mln/level/fill.hh> + +# include <mln/debug/println.hh> + +# include "simple_point.hh" + + +namespace mln +{ + + template <typename D, typename I, typename N> + mln_ch_value(I, D) + distance(D, const I& input, const N& nbh) + { + const D M = mln_max(D); + mln_ch_value(I, D) output; + initialize(output, input); + + mln_ch_value(I, bool) deja_vu; + initialize(deja_vu, input); + level::fill(deja_vu, input); + + typedef mln_site(I) P; + p_queue_fast q; + + // Initialization. + { + level::fill(output, M); + mln_piter(I) p(input.domain()); + mln_niter(N) n(nbh, p); + for_all(p) + if (input(p) == true) // p in object + { + output(p) = 0; + for_all(n) + if (input.domain().has(n) && input(n) == false) // n in background + { + q.push(p); + break; + } + } + } + + // Propagation. + { + P p; + mln_niter(N) n(nbh, p); + while (! q.is_empty()) + { + p = q.pop_front(); + for_all(n) + if (input.domain().has(n) && output(n) == M) + { + output(n) = output(p) + 1; + if (output(n) == M) + { + // Saturation so stop. + q.clear(); + break; + } + q.push(n); + deja_vu(n) = true; + } + } + } + + return output; + } + + image2d<bool> + skel(const image2d<bool>& input, unsigned nbh_i) + { + neighb2d nbh = int_to_neighb(nbh_i); + image2d<bool> output; + initialize(output, input); + + image2d<value::int_u16> dist_map = distance(value::int_u16(), input, nbh); + + typedef mln_site_(image2d<bool>) P; + p_priority<value::int_u16, p_queue_fast > q; + + // Initialization. + { + level::fill(output, input); + mln_piter_(image2d<bool>) p(output.domain()); + for_all(p) + if (output(p) && + is_simple_point2d(output, nbh_i, p) && // p is a not simple point of object + !is_curve_extremum(output, nbh_i, p)) + { + output(p) = false; // Remove p from object + q.push(dist_map(p), p); + } + } + + // Propagation. + { + P p; + mln_niter_(neighb2d) n(nbh, p); + while (! q.is_empty()) + { + p = q.pop_front(); + + for_all(n) + if (output.domain().has(n) && + output(n) && + is_simple_point2d(output, nbh_i, n) && // n is not simple + !is_curve_extremum(output, nbh_i, n)) + { + output(n) = false; // Remove p from object + q.push(dist_map(n), n); + } + } + } + + return output; + } + +} // mln + + + +void usage(char* argv[]) +{ + std::cerr << "usage: " << argv[0] << " input.pbm output.pbm" << std::endl; + abort(); +} + + + +int main(int argc, char* argv[]) +{ + using namespace mln; + using value::int_u8; + + if (argc != 3) + usage(argv); + + image2d<bool> input; + io::pbm::load(input, argv[1]); + + image2d<bool> output = skel(input, 4); + io::pbm::save(output, argv[2]); +}

1 0

r2736: Add is_curve_extremum to detect curve extrema
by Matthieu Garrigues 30 Oct '08

30 Oct '08

URL: https://svn.lrde.epita.fr/svn/oln/branches/cleanup-2008/milena/sandbox ChangeLog: 2008-10-30 Matthieu Garrigues <garrigues(a)lrde.epita.fr> Add is_curve_extremum to detect curve extrema. * garrigues/ocr/simple_point.hh: (is_curve_extremum) Tell if a point is a curve extremum. --- simple_point.hh | 32 ++++++++++++++++++++++++-------- 1 file changed, 24 insertions(+), 8 deletions(-) Index: branches/cleanup-2008/milena/sandbox/garrigues/ocr/simple_point.hh =================================================================== --- branches/cleanup-2008/milena/sandbox/garrigues/ocr/simple_point.hh (revision 2735) +++ branches/cleanup-2008/milena/sandbox/garrigues/ocr/simple_point.hh (revision 2736) @@ -30,7 +30,7 @@ /*! \file simple_point.hh * - * \brief simple_point tell if a point is simple or not (Cf + * \brief is_simple_point tell if a point is simple or not (Cf * bertrand.07.chap). * */ @@ -57,10 +57,12 @@ * */ - bool simple_point(const image2d<bool>& ima, const neighb2d& nbh, const point2d& p); + bool is_simple_point(const image2d<bool>& ima, const neighb2d& nbh, const point2d& p); unsigned nb_connexity2d(const image2d<bool>& ima, const neighb2d& nbh, const point2d& p); + bool is_curve_extremum(const image2d<bool>& ima, unsigned nbh, const point2d& p); + # ifndef MLN_INCLUDE_ONLY static const unsigned char nb_connexity_c8[256] = @@ -129,15 +131,15 @@ return 0; } - unsigned nb_connexity2d(const image2d<bool>& ima, unsigned nbh, const point2d& p, bool complemented) + unsigned nb_connexity2d(const image2d<bool>& ima, unsigned nbh, const point2d& p, bool object) { unsigned res = 0; - mln_fwd_niter_(neighb2d) n(c8() , p); + mln_bkd_niter_(neighb2d) n(c8() , p); for_all(n) { res = (res << 1); - if (ima(n) == complemented) + if (ima(n) == object) res = res | 1; } @@ -150,12 +152,26 @@ } } - bool simple_point2d(const image2d<bool>& ima, unsigned nbh, const point2d& p) + bool is_curve_extremum(const image2d<bool>& ima, unsigned nbh, const point2d& p) + { + unsigned cpt = 0; + mln_bkd_niter_(neighb2d) n(c8() , p); + + for_all(n) + { + if (ima(n) == true) + cpt++; + } + + return cpt == 1; + } + + bool is_simple_point2d(const image2d<bool>& ima, unsigned nbh, const point2d& p) { mln_assertion(nbh == 4 || nbh == 8); - return nb_connexity2d(ima, nbh, p, true) == 1 && - nb_connexity2d(ima, complement_neighb(nbh), p, false) == 1; + return (nb_connexity2d(ima, nbh, p, true) == 1) && + (nb_connexity2d(ima, complement_neighb(nbh), p, false) == 1); } # endif // MLN_INCLUDE_ONLY

1 0