- Olena-patches

1029: Fix ChangeLogs.
by Roland Levillain 23 Jul '07

23 Jul '07

Index: ChangeLog =================================================================== --- ChangeLog (revision 1028) +++ ChangeLog (revision 1029) @@ -1,15 +1,3 @@ -2007-07-23 Roland Levillain <roland(a)lrde.epita.fr> - - Clean up stc/scoop.hxx a bit. - - * stc/scoop.hxx: Add a copyright header. - Typos. - Re-indent some parts. - (helper_check): Fix documentation. - Actually use check_final_stm and check_delegatee_inherited as - described in the documentation. - (helper_find): Add some documentation. - 2007-06-29 Roland Levillain <roland(a)lrde.epita.fr> Add support for Doxygen in Milena. Index: static/ChangeLog =================================================================== --- static/ChangeLog (revision 1028) +++ static/ChangeLog (revision 1029) @@ -1,3 +1,15 @@ +2007-07-23 Roland Levillain <roland(a)lrde.epita.fr> + + Clean up stc/scoop.hxx a bit. + + * stc/scoop.hxx: Add a copyright header. + Typos. + Re-indent some parts. + (helper_check): Fix documentation. + Actually use check_final_stm and check_delegatee_inherited as + described in the documentation. + (helper_find): Add some documentation. + 2007-07-17 Roland Levillain <roland(a)lrde.epita.fr> Update stc/scoop.{hh,hxx} with the latest version from Olena.

1 0

1027: Renaming some concepts with more explicit and longer names.
by Thierry Geraud 23 Jul '07

23 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Thierry Geraud <thierry.geraud(a)lrde.epita.fr> Renaming some concepts with more explicit and longer names. * mln/core/concept/piter.hh, * mln/core/concept/genpoint.hh, * mln/core/concept/psite.hh, * mln/core/concept/viter.hh, * mln/core/concept/doc/genpoint.hh, * mln/core/concept/doc/piter.hh, * mln/core/concept/doc/viter.hh, * mln/core/concept/doc/genpixel.hh, * mln/core/concept/genpixel.hh: Rename respectively as... * mln/core/concept/point_iterator.hh, * mln/core/concept/generalized_point.hh, * mln/core/concept/point_site.hh, * mln/core/concept/value_iterator.hh, * mln/core/concept/doc/generalized_point.hh, * mln/core/concept/doc/point_iterator.hh, * mln/core/concept/doc/value_iterator.hh, * mln/core/concept/doc/generalized_pixel.hh, * mln/core/concept/generalized_pixel.hh: ...these. * mln/convert/to_dpoint.hh, * mln/geom/pmin_pmax.hh, * mln/core/dpoints_piter.hh, * mln/core/psubset.hh, * mln/core/box.hh, * mln/core/psubset_piter.hh, * mln/core/rectangle2d.hh, * mln/core/box_piter.hh, * mln/core/window.hh, * mln/core/dpoints_pixter.hh, * mln/core/internal/image_adaptor.hh, * mln/core/internal/image_base.hh, * mln/core/internal/piter_adaptor.hh, * mln/core/pixel.hh, * mln/core/hline2d.hh, * mln/core/neighb.hh, * mln/core/concept/value_set.hh, * mln/core/concept/point.hh, * mln/core/concept/pixel_iterator.hh, * mln/core/concept/point_set.hh, * mln/core/concept/doc/image.hh, * mln/core/concept/doc/value_set.hh, * mln/core/concept/doc/pixel_iterator.hh, * mln/core/concept/doc/point_set.hh, * mln/core/concept/doc/window.hh, * mln/core/concept/doc/neighborhood.hh, * mln/core/vline2d.hh, * mln/level/approx/median.hh, * mln/value/lut_vec.hh, * mln/value/viter.hh, * mln/value/set.hh: Update. * TODO: Update. TODO | 17 +++++++-- mln/convert/to_dpoint.hh | 6 +-- mln/core/box.hh | 6 +-- mln/core/box_piter.hh | 10 ++--- mln/core/concept/doc/generalized_pixel.hh | 12 +++--- mln/core/concept/doc/generalized_point.hh | 10 ++--- mln/core/concept/doc/image.hh | 8 ++-- mln/core/concept/doc/neighborhood.hh | 6 +-- mln/core/concept/doc/pixel_iterator.hh | 2 - mln/core/concept/doc/point_iterator.hh | 16 ++++---- mln/core/concept/doc/point_set.hh | 6 +-- mln/core/concept/doc/value_iterator.hh | 10 ++--- mln/core/concept/doc/value_set.hh | 4 +- mln/core/concept/doc/window.hh | 6 +-- mln/core/concept/generalized_pixel.hh | 14 +++---- mln/core/concept/generalized_point.hh | 54 +++++++++++++++--------------- mln/core/concept/pixel_iterator.hh | 4 +- mln/core/concept/point.hh | 4 +- mln/core/concept/point_iterator.hh | 16 ++++---- mln/core/concept/point_set.hh | 2 - mln/core/concept/point_site.hh | 16 ++++---- mln/core/concept/value_iterator.hh | 20 +++++------ mln/core/concept/value_set.hh | 2 - mln/core/dpoints_piter.hh | 12 +++--- mln/core/dpoints_pixter.hh | 10 ++--- mln/core/hline2d.hh | 4 +- mln/core/internal/image_adaptor.hh | 2 - mln/core/internal/image_base.hh | 6 +-- mln/core/internal/piter_adaptor.hh | 6 +-- mln/core/neighb.hh | 4 +- mln/core/pixel.hh | 4 +- mln/core/psubset.hh | 6 +-- mln/core/psubset_piter.hh | 2 - mln/core/rectangle2d.hh | 4 +- mln/core/vline2d.hh | 4 +- mln/core/window.hh | 6 +-- mln/geom/pmin_pmax.hh | 4 +- mln/level/approx/median.hh | 2 - mln/value/lut_vec.hh | 4 +- mln/value/set.hh | 4 +- mln/value/viter.hh | 6 +-- 41 files changed, 175 insertions(+), 166 deletions(-) Index: TODO --- TODO (revision 1026) +++ TODO (working copy) @@ -25,17 +25,13 @@ in Fast_Image: memcpy and memset-like methods (?) .offset() into GenPixel -notion of low-quantization -vset and .values() in Image built-in op objects -> reverse lhs/rhs * renaming -GenPixel, GenPoint, Psite, Piter: longer name (?) mlc into metal point-wise material: pw::value_of, pw::constant, and pw::variable -kind-related: mln_value_kind(I) and mln_kind(V) * clean-up @@ -52,3 +48,16 @@ border::* and no-op versions if not fast image arith::inplace_plus et al. linear:: for convolutions + + +svn mv genpoint.hh generalized_point.hh +svn mv genpixel.hh generalized_pixel.hh +svn mv psite.hh point_site.hh +svn mv viter.hh value_iterator.hh +svn mv piter.hh point_iterator.hh + + +svn mv doc/genpoint.hh doc/generalized_point.hh +svn mv doc/genpixel.hh doc/generalized_pixel.hh +svn mv doc/viter.hh doc/value_iterator.hh +svn mv doc/piter.hh doc/point_iterator.hh Index: mln/convert/to_dpoint.hh --- mln/convert/to_dpoint.hh (revision 1026) +++ mln/convert/to_dpoint.hh (working copy) @@ -33,7 +33,7 @@ * \brief Convertions to mln::Dpoint. */ -# include <mln/core/concept/genpoint.hh> +# include <mln/core/concept/generalized_point.hh> namespace mln @@ -44,13 +44,13 @@ /// Convert a generalized point \p p into a delta-point. template <typename P> - mln_dpoint(P) to_dpoint(const GenPoint& p); + mln_dpoint(P) to_dpoint(const Generalized_Point& p); # ifndef MLN_INCLUDE_ONLY template <typename P> - mln_dpoint(P) to_dpoint(const GenPoint& p_) + mln_dpoint(P) to_dpoint(const Generalized_Point& p_) { const P& p = internal::force_exact(p_); mln_dpoint(P) tmp; Index: mln/geom/pmin_pmax.hh --- mln/geom/pmin_pmax.hh (revision 1026) +++ mln/geom/pmin_pmax.hh (working copy) @@ -48,7 +48,7 @@ /// iterator \p p. template <typename I> std::pair<mln_point(I), mln_point(I)> - pmin_pmax(const Piter& p); + pmin_pmax(const Point_Iterator& p); @@ -56,7 +56,7 @@ template <typename I> std::pair<mln_point(I), mln_point(I)> - pmin_pmax(Piter& p_) + pmin_pmax(Point_Iterator& p_) { I& p = exact(p_); Index: mln/core/dpoints_piter.hh --- mln/core/dpoints_piter.hh (revision 1026) +++ mln/core/dpoints_piter.hh (working copy) @@ -33,8 +33,8 @@ * \brief Definition of mln::dpoints_fwd_piter and mln::dpoints_bkd_piter. */ -# include <mln/core/concept/piter.hh> -# include <mln/core/concept/genpoint.hh> +# include <mln/core/concept/point_iterator.hh> +# include <mln/core/concept/generalized_point.hh> namespace mln @@ -46,7 +46,7 @@ * The parameter \c D is the type of delta-points. */ template <typename D> - class dpoints_fwd_piter : public Piter< dpoints_fwd_piter<D> > + class dpoints_fwd_piter : public Point_Iterator< dpoints_fwd_piter<D> > { public: @@ -59,7 +59,7 @@ /// Point associated type. typedef mln_point(D) point; - /// Psite associated type. + /// Point_Site associated type. typedef point psite; /// Coordinate associated type. @@ -72,7 +72,7 @@ */ template <typename Dps, typename Pref> dpoints_fwd_piter(const Dps& dps, // FIXME: explicitly set_of_<D>? - const GenPoint<Pref>& p_ref); + const Generalized_Point<Pref>& p_ref); /// Convertion to point. operator mln_point(D) () const; @@ -117,7 +117,7 @@ template <typename D> template <typename Dps, typename Pref> dpoints_fwd_piter<D>::dpoints_fwd_piter(const Dps& dps, - const GenPoint<Pref>& p_ref) + const Generalized_Point<Pref>& p_ref) : dps_(exact(dps).vec()), p_ref_(* internal::force_exact<Pref>(p_ref).pointer_()) { Index: mln/core/psubset.hh --- mln/core/psubset.hh (revision 1026) +++ mln/core/psubset.hh (working copy) @@ -67,17 +67,17 @@ template <typename S, typename F> struct psubset : public Point_Set< psubset<S,F> > { - /// Psite associated type. + /// Point_Site associated type. typedef mln_psite(S) psite; /// Point associated type. typedef mln_point(S) point; - /// Forward Piter associated type. + /// Forward Point_Iterator associated type. typedef psubset_fwd_piter_<S,F> fwd_piter; - /// Backward Piter associated type. + /// Backward Point_Iterator associated type. typedef internal::fixme bkd_piter; Index: mln/core/box.hh --- mln/core/box.hh (revision 1026) +++ mln/core/box.hh (working copy) @@ -54,7 +54,7 @@ struct box_ : public Box< box_ >, public internal::box_impl_< P::dim, mln_coord(P), box_ > { - /*! \brief Psite associated type. + /*! \brief Point_Site associated type. */ typedef P psite; @@ -62,11 +62,11 @@ */ typedef P point; - /*! \brief Forward Piter associated type. + /*! \brief Forward Point_Iterator associated type. */ typedef box_fwd_piter_ fwd_piter; - /*! \brief Backward Piter associated type. + /*! \brief Backward Point_Iterator associated type. */ typedef box_bkd_piter_ bkd_piter; Index: mln/core/psubset_piter.hh --- mln/core/psubset_piter.hh (revision 1026) +++ mln/core/psubset_piter.hh (working copy) @@ -33,7 +33,7 @@ * \brief Definition of iterators on points of psubsetes. */ -# include <mln/core/concept/piter.hh> +# include <mln/core/concept/point_iterator.hh> # include <mln/core/internal/piter_adaptor.hh> # include <mln/core/psubset.hh> Index: mln/core/rectangle2d.hh --- mln/core/rectangle2d.hh (revision 1026) +++ mln/core/rectangle2d.hh (working copy) @@ -62,12 +62,12 @@ /// Dpoint associated type. typedef dpoint2d dpoint; - /*! \brief Piter type to browse a rectangle such as: "for each row + /*! \brief Point_Iterator type to browse a rectangle such as: "for each row * (increasing), for each column (increasing)." */ typedef dpoints_fwd_piter<dpoint2d> fwd_qiter; - /*! \brief Piter type to browse a rectangle such as: "for each row + /*! \brief Point_Iterator type to browse a rectangle such as: "for each row * (decreasing), for each column (decreasing)." */ typedef internal::fixme bkd_qiter; Index: mln/core/box_piter.hh --- mln/core/box_piter.hh (revision 1026) +++ mln/core/box_piter.hh (working copy) @@ -33,7 +33,7 @@ * \brief Definition of iterators on points of boxes. */ -# include <mln/core/concept/piter.hh> +# include <mln/core/concept/point_iterator.hh> # include <mln/core/concept/box.hh> @@ -47,14 +47,14 @@ * \see mln::box_ */ template <typename P> - class box_fwd_piter_ : public Piter< box_fwd_piter_ > + class box_fwd_piter_ : public Point_Iterator< box_fwd_piter_ > { public: /// Space dimension. enum { dim = P::dim }; - /// Psite associated type. + /// Point_Site associated type. typedef P psite; /// Point associated type. @@ -107,14 +107,14 @@ * \see mln::box_ */ template <typename P> - class box_bkd_piter_ : public Piter< box_bkd_piter_ > + class box_bkd_piter_ : public Point_Iterator< box_bkd_piter_ > { public: /// Space dimension. enum { dim = P::dim }; - /// Psite associated type. + /// Point_Site associated type. typedef P psite; /// Point associated type. Index: mln/core/window.hh --- mln/core/window.hh (revision 1026) +++ mln/core/window.hh (working copy) @@ -34,7 +34,7 @@ */ # include <mln/core/concept/window.hh> -# include <mln/core/concept/genpoint.hh> +# include <mln/core/concept/generalized_point.hh> # include <mln/core/internal/set_of.hh> # include <mln/core/dpoint.hh> # include <mln/core/box.hh> @@ -67,12 +67,12 @@ /// Dpoint associated type. typedef D dpoint; - /*! \brief Piter type to browse the points of a generic window + /*! \brief Point_Iterator type to browse the points of a generic window * w.r.t. the ordering of delta-points. */ typedef dpoints_fwd_piter<D> fwd_qiter; - /*! \brief Piter type to browse the points of a generic window + /*! \brief Point_Iterator type to browse the points of a generic window * w.r.t. the reverse ordering of delta-points. */ typedef internal::fixme bkd_qiter; Index: mln/core/dpoints_pixter.hh --- mln/core/dpoints_pixter.hh (revision 1026) +++ mln/core/dpoints_pixter.hh (working copy) @@ -38,7 +38,7 @@ # include <vector> # include <mln/core/concept/pixel_iterator.hh> -# include <mln/core/concept/genpixel.hh> +# include <mln/core/concept/generalized_pixel.hh> # include <mln/core/internal/pixel_impl.hh> @@ -70,7 +70,7 @@ template <typename Dps, typename Pref> dpoints_fwd_pixter(I& image, const Dps& dps, - const GenPoint<Pref>& p_ref); + const Generalized_Point<Pref>& p_ref); /*! \brief Constructor. * @@ -81,7 +81,7 @@ template <typename Dps, typename Pref> dpoints_fwd_pixter(I& image, const Dps& dps, - const GenPixel<Pref>& p_ref); + const Generalized_Pixel<Pref>& p_ref); /// Start an iteration. void start(); @@ -134,7 +134,7 @@ template <typename Dps, typename Pref> dpoints_fwd_pixter::dpoints_fwd_pixter(I& image, const Dps& dps, - const GenPoint<Pref>& p_ref) + const Generalized_Point<Pref>& p_ref) : super(image) { mln_precondition(image.has_data()); @@ -148,7 +148,7 @@ template <typename Dps, typename Pref> dpoints_fwd_pixter::dpoints_fwd_pixter(I& image, const Dps& dps, - const GenPixel<Pref>& p_ref) + const Generalized_Pixel<Pref>& p_ref) : super(image) { mln_precondition(image.has_data()); Index: mln/core/internal/image_adaptor.hh --- mln/core/internal/image_adaptor.hh (revision 1026) +++ mln/core/internal/image_adaptor.hh (working copy) @@ -51,7 +51,7 @@ typename S = mln_pset(I) > struct image_adaptor_ : public internal::image_base_< S, E > { - /// Psite associated type. + /// Point_Site associated type. typedef mln_psite(S) psite; /// Point_Set associated type. Index: mln/core/internal/image_base.hh --- mln/core/internal/image_base.hh (revision 1026) +++ mln/core/internal/image_base.hh (working copy) @@ -70,7 +70,7 @@ /// Point_Set associated type. typedef S pset; - /// Psite associated type. + /// Point_Site associated type. typedef mln_psite(S) psite; /// Point associated type. @@ -84,10 +84,10 @@ typedef mln_coord(point) coord; - /// Forward Piter associated type. + /// Forward Point_Iterator associated type. typedef mln_fwd_piter(S) fwd_piter; - /// Backward Piter associated type. + /// Backward Point_Iterator associated type. typedef mln_bkd_piter(S) bkd_piter; Index: mln/core/internal/piter_adaptor.hh --- mln/core/internal/piter_adaptor.hh (revision 1026) +++ mln/core/internal/piter_adaptor.hh (working copy) @@ -33,7 +33,7 @@ * \brief Definition of iterators on points of boxes. */ -# include <mln/core/concept/piter.hh> +# include <mln/core/concept/point_iterator.hh> # include <mln/core/concept/box.hh> @@ -49,14 +49,14 @@ * parameter E is the exact type. */ template <typename Pi, typename E> - class piter_adaptor_ : public Piter<E> + class piter_adaptor_ : public Point_Iterator<E> { public: /// Space dimension. enum { dim = Pi::dim }; - /// Psite associated type. + /// Point_Site associated type. typedef mln_psite(Pi) psite; /// Point associated type. Index: mln/core/pixel.hh --- mln/core/pixel.hh (revision 1026) +++ mln/core/pixel.hh (working copy) @@ -33,7 +33,7 @@ * \brief Definition of the generic pixel class mln::pixel. */ -# include <mln/core/concept/genpixel.hh> +# include <mln/core/concept/generalized_pixel.hh> # include <mln/core/internal/pixel_impl.hh> @@ -47,7 +47,7 @@ */ template <typename I> struct pixel : public Object< pixel >, - public GenPixel< pixel >, + public Generalized_Pixel< pixel >, internal::pixel_impl_< I, pixel > { /// Constructor. Index: mln/core/hline2d.hh --- mln/core/hline2d.hh (revision 1026) +++ mln/core/hline2d.hh (working copy) @@ -60,12 +60,12 @@ /// Dpoint associated type. typedef dpoint2d dpoint; - /*! \brief Piter type to browse a hline such as: "for each row + /*! \brief Point_Iterator type to browse a hline such as: "for each row * (increasing), for each column (increasing)." */ typedef dpoints_fwd_piter<dpoint2d> fwd_qiter; - /*! \brief Piter type to browse a hline such as: "for each row + /*! \brief Point_Iterator type to browse a hline such as: "for each row * (decreasing), for each column (decreasing)." */ typedef internal::fixme bkd_qiter; Index: mln/core/neighb.hh --- mln/core/neighb.hh (revision 1026) +++ mln/core/neighb.hh (working copy) @@ -58,12 +58,12 @@ /// Dpoint associated type. typedef D dpoint; - /*! \brief Piter type to browse the points of a generic + /*! \brief Point_Iterator type to browse the points of a generic * neighborhood w.r.t. the ordering of delta-points. */ typedef dpoints_fwd_piter<D> fwd_niter; - /*! \brief Piter type to browse the points of a generic + /*! \brief Point_Iterator type to browse the points of a generic * neighborhood w.r.t. the reverse ordering of delta-points. */ typedef dpoints_bkd_piter<D> bkd_niter; Index: mln/core/concept/value_iterator.hh --- mln/core/concept/value_iterator.hh (revision 1026) +++ mln/core/concept/value_iterator.hh (working copy) @@ -28,13 +28,13 @@ #ifndef MLN_CORE_CONCEPT_VITER_HH # define MLN_CORE_CONCEPT_VITER_HH -/*! \file mln/core/concept/viter.hh +/*! \file mln/core/concept/value_iterator.hh * - * \brief Definition of the concept of mln::Viter. + * \brief Definition of the concept of mln::Value_Iterator. */ # include <mln/core/concept/iterator.hh> -# include <mln/core/concept/genpoint.hh> +# include <mln/core/concept/generalized_point.hh> namespace mln @@ -46,11 +46,11 @@ * An iterator on values is an iterator that browse over a set of * values. * - * \see mln::doc::Viter for a complete documentation of this class + * \see mln::doc::Value_Iterator for a complete documentation of this class * contents. */ template <typename E> - struct Viter : public Iterator<E> + struct Value_Iterator : public Iterator<E> { /* typedef value; @@ -58,7 +58,7 @@ */ protected: - Viter(); + Value_Iterator(); }; @@ -72,17 +72,17 @@ * * \return The modified output stream \p ostr. * - * \relates mln::Viter + * \relates mln::Value_Iterator */ template <typename E> - std::ostream& operator<<(std::ostream& ostr, const Viter<E>& v); + std::ostream& operator<<(std::ostream& ostr, const Value_Iterator<E>& v); # ifndef MLN_INCLUDE_ONLY template <typename E> - Viter<E>::Viter() + Value_Iterator<E>::Value_Iterator() { typedef mln_value(E) value; value (E::*m)() const = & E::operator value; @@ -90,7 +90,7 @@ } template <typename E> - std::ostream& operator<<(std::ostream& ostr, const Viter<E>& v) + std::ostream& operator<<(std::ostream& ostr, const Value_Iterator<E>& v) { mln_precondition(exact(v).is_valid()); return ostr << mln_value(E)(exact(v)); Index: mln/core/concept/value_set.hh --- mln/core/concept/value_set.hh (revision 1026) +++ mln/core/concept/value_set.hh (working copy) @@ -32,7 +32,7 @@ * \brief Definition of the concept of mln::Value_Set. */ -# include <mln/core/concept/viter.hh> +# include <mln/core/concept/value_iterator.hh> namespace mln Index: mln/core/concept/point.hh --- mln/core/concept/point.hh (revision 1026) +++ mln/core/concept/point.hh (working copy) @@ -33,7 +33,7 @@ * \brief Definition of the concept of mln::Point. */ -# include <mln/core/concept/psite.hh> +# include <mln/core/concept/point_site.hh> namespace mln @@ -47,7 +47,7 @@ * discrete square grid of the 2D plane. */ template <typename P> - struct Point : public Psite + struct Point : public Point_Site { /*! \brief The associated point type is itself. */ Index: mln/core/concept/pixel_iterator.hh --- mln/core/concept/pixel_iterator.hh (revision 1026) +++ mln/core/concept/pixel_iterator.hh (working copy) @@ -34,7 +34,7 @@ */ # include <mln/core/concept/iterator.hh> -# include <mln/core/concept/genpixel.hh> +# include <mln/core/concept/generalized_pixel.hh> namespace mln @@ -51,7 +51,7 @@ */ template <typename E> struct Pixel_Iterator : public Iterator<E>, - public GenPixel<E> + public Generalized_Pixel<E> { /* typedef lvalue; Index: mln/core/concept/generalized_pixel.hh --- mln/core/concept/generalized_pixel.hh (revision 1026) +++ mln/core/concept/generalized_pixel.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_CORE_CONCEPT_GENPIXEL_HH # define MLN_CORE_CONCEPT_GENPIXEL_HH -/*! \file mln/core/concept/genpixel.hh +/*! \file mln/core/concept/generalized_pixel.hh * * \brief Definition of the concept of mln::Genpixel. */ @@ -40,21 +40,21 @@ namespace mln { - // FIXME: \class GenPixel GenPixel.hh "mln/core/concept/doc/GenPixel.hh" + // FIXME: \class Generalized_Pixel Generalized_Pixel.hh "mln/core/concept/doc/Generalized_Pixel.hh" /*! \brief Base class for implementation classes that are pixels or that * have the behavior of pixels. * - * "GenPixel" is "Generalized Pixel" for short. + * "Generalized_Pixel" is "Generalized Pixel" for short. * * \warning This class does \em not derive from mln::Object; it is * for use as a parallel hierarchy. * - * \see mln::doc::GenPixel for a complete documentation of this + * \see mln::doc::Generalized_Pixel for a complete documentation of this * class contents. */ template <typename E> - struct GenPixel + struct Generalized_Pixel { /* typedef value; @@ -64,14 +64,14 @@ value** address_() const; */ protected: - GenPixel(); + Generalized_Pixel(); }; # ifndef MLN_INCLUDE_ONLY template <typename E> - GenPixel<E>::GenPixel() + Generalized_Pixel<E>::Generalized_Pixel() { typedef mln_value(E) value; typedef mln_rvalue(E) rvalue; Index: mln/core/concept/generalized_point.hh --- mln/core/concept/generalized_point.hh (revision 1026) +++ mln/core/concept/generalized_point.hh (working copy) @@ -28,8 +28,8 @@ #ifndef MLN_CORE_CONCEPT_GENPOINT_HH # define MLN_CORE_CONCEPT_GENPOINT_HH -/*! \file mln/core/concept/genpoint.hh - * \brief Definition of the concept of mln::GenPoint. +/*! \file mln/core/concept/generalized_point.hh + * \brief Definition of the concept of mln::Generalized_Point. */ # include <mlc/equal.hh> # include <mlc/same_point.hh> @@ -42,21 +42,21 @@ namespace mln { - // FIXME: \class GenPoint genpoint.hh "mln/core/concept/doc/genpoint.hh" + // FIXME: \class Generalized_Point generalized_point.hh "mln/core/concept/doc/generalized_point.hh" /*! \brief Base class for implementation classes that are points or that * have the behavior of points. * - * "GenPoint" is "Generalized Point" for short. + * "Generalized_Point" is "Generalized Point" for short. * * \warning This class does \em not derive from mln::Object; it is * for use as a parallel hierarchy. * - * \see mln::doc::GenPoint for a complete documentation of this + * \see mln::doc::Generalized_Point for a complete documentation of this * class contents. */ template <typename E> - struct GenPoint + struct Generalized_Point { /* @@ -75,7 +75,7 @@ */ protected: - GenPoint(); + Generalized_Point(); }; @@ -93,10 +93,10 @@ * \return True if both generalized points have the same * coordinates, otherwise false. * - * \relates mln::GenPoint + * \relates mln::Generalized_Point */ template <typename Pl, typename Pr> - bool operator=(const GenPoint<Pl>& lhs, const GenPoint<Pr>& rhs); + bool operator=(const Generalized_Point<Pl>& lhs, const Generalized_Point<Pr>& rhs); /*! \brief Ordering "less than" comparison between a couple of generalized @@ -117,10 +117,10 @@ * \return True if \p lhs is before \p rhs in the sense of the * coordinates lexicographic comparison, otherwise false. * - * \relates mln::GenPoint + * \relates mln::Generalized_Point */ template <typename Pl, typename Pr> - bool operator<(const GenPoint<Pl>& lhs, const GenPoint<Pr>& rhs); + bool operator<(const Generalized_Point<Pl>& lhs, const Generalized_Point<Pr>& rhs); /*! \brief Difference between a couple of generalized point \p lhs and \p @@ -142,11 +142,11 @@ * \return A delta-point (temporary object). * * \see mln::Dpoint - * \relates mln::GenPoint + * \relates mln::Generalized_Point */ template <typename Pl, typename Pr> mln_dpoint(Pl) - operator-(const GenPoint<Pl>& lhs, const GenPoint<Pr>& rhs); + operator-(const Generalized_Point<Pl>& lhs, const Generalized_Point<Pr>& rhs); /*! \brief Add a delta-point \p rhs to a generalized point \p lhs. @@ -160,14 +160,14 @@ * \return A point (temporary object). * * \see mln::Dpoint - * \relates mln::GenPoint + * \relates mln::Generalized_Point * * \todo Introduce the notion of "generalized dpoint" and - * add the more general extra operator+(GenPoint, GenDpoint). + * add the more general extra operator+(Generalized_Point, GenDpoint). */ template <typename P> mln_point(P) - operator+(const GenPoint& lhs, const mln_dpoint(P)& rhs); + operator+(const Generalized_Point& lhs, const mln_dpoint(P)& rhs); /*! \brief Substract a delta-point \p rhs to a generalized point \p lhs. @@ -181,11 +181,11 @@ * \return A point (temporary object). * * \see mln::Dpoint - * \relates mln::GenPoint + * \relates mln::Generalized_Point */ template <typename P> mln_point(P) - operator-(const GenPoint& lhs, const mln_dpoint(P)& rhs); + operator-(const Generalized_Point& lhs, const mln_dpoint(P)& rhs); /*! \brief Print a generalized point \p p into the output stream \p @@ -196,10 +196,10 @@ * * \return The modified output stream \p ostr. * - * \relates mln::GenPoint + * \relates mln::Generalized_Point */ template <typename P> - std::ostream& operator<<(std::ostream& ostr, const GenPoint& p); + std::ostream& operator<<(std::ostream& ostr, const Generalized_Point& p); @@ -207,7 +207,7 @@ template <typename E> - GenPoint<E>::GenPoint() + Generalized_Point<E>::Generalized_Point() { int dim = E::dim; mln_invariant(dim > 0); @@ -223,7 +223,7 @@ template <typename Pl, typename Pr> - bool operator=(const GenPoint<Pl>& lhs, const GenPoint<Pr>& rhs) + bool operator=(const Generalized_Point<Pl>& lhs, const Generalized_Point<Pr>& rhs) { // FIXME: mlc::same_grid<Pl, Pr>::check(); const Pl& lhs_ = internal::force_exact<Pl>(lhs); @@ -236,7 +236,7 @@ } template <typename Pl, typename Pr> - bool operator<(const GenPoint<Pl>& lhs, const GenPoint<Pr>& rhs) + bool operator<(const Generalized_Point<Pl>& lhs, const Generalized_Point<Pr>& rhs) { // FIXME: mlc::same_grid<Pl, Pr>::check(); const Pl& lhs_ = internal::force_exact<Pl>(lhs); @@ -252,7 +252,7 @@ template <typename Pl, typename Pr> mln_dpoint(Pl) - operator-(const GenPoint<Pl>& lhs, const GenPoint<Pr>& rhs) + operator-(const Generalized_Point<Pl>& lhs, const Generalized_Point<Pr>& rhs) { mlc::equal<mln_dpoint(Pl), mln_dpoint(Pr)>::check(); // FIXME: mlc::same_grid<Pl, Pr>::check(); @@ -268,7 +268,7 @@ template <typename P> mln_point(P) - operator+(const GenPoint& lhs, const mln_dpoint(P)& rhs) + operator+(const Generalized_Point& lhs, const mln_dpoint(P)& rhs) { const P& lhs_ = internal::force_exact(lhs); mln_point(P) tmp; @@ -279,7 +279,7 @@ template <typename P> mln_point(P) - operator-(const GenPoint& lhs, const mln_dpoint(P)& rhs) + operator-(const Generalized_Point& lhs, const mln_dpoint(P)& rhs) { const P& lhs_ = internal::force_exact(lhs); mln_point(P) tmp; @@ -289,7 +289,7 @@ } template <typename P> - std::ostream& operator<<(std::ostream& ostr, const GenPoint& p) + std::ostream& operator<<(std::ostream& ostr, const Generalized_Point& p) { const P& p_ = internal::force_exact(p); Index: mln/core/concept/point_iterator.hh --- mln/core/concept/point_iterator.hh (revision 1026) +++ mln/core/concept/point_iterator.hh (working copy) @@ -28,13 +28,13 @@ #ifndef MLN_CORE_CONCEPT_PITER_HH # define MLN_CORE_CONCEPT_PITER_HH -/*! \file mln/core/concept/piter.hh +/*! \file mln/core/concept/point_iterator.hh * - * \brief Definition of the concept of mln::Piter. + * \brief Definition of the concept of mln::Point_Iterator. */ # include <mln/core/concept/iterator.hh> -# include <mln/core/concept/genpoint.hh> +# include <mln/core/concept/generalized_point.hh> namespace mln @@ -46,12 +46,12 @@ * An iterator on points is an iterator that browse over a set of * points. * - * \see mln::doc::Piter for a complete documentation of this class + * \see mln::doc::Point_Iterator for a complete documentation of this class * contents. */ template <typename E> - struct Piter : public Iterator<E>, - public GenPoint<E> + struct Point_Iterator : public Iterator<E>, + public Generalized_Point<E> { /* typedef psite; @@ -59,14 +59,14 @@ */ protected: - Piter(); + Point_Iterator(); }; # ifndef MLN_INCLUDE_ONLY template <typename E> - Piter<E>::Piter() + Point_Iterator<E>::Point_Iterator() { typedef mln_psite(E) psite; psite (E::*m)() const = & E::operator psite; Index: mln/core/concept/point_set.hh --- mln/core/concept/point_set.hh (revision 1026) +++ mln/core/concept/point_set.hh (working copy) @@ -34,7 +34,7 @@ */ # include <mln/core/concept/point.hh> -# include <mln/core/concept/piter.hh> +# include <mln/core/concept/point_iterator.hh> namespace mln Index: mln/core/concept/doc/image.hh --- mln/core/concept/doc/image.hh (revision 1026) +++ mln/core/concept/doc/image.hh (working copy) @@ -127,8 +127,8 @@ */ typedef void point; - /*! \brief Psite associated type. - * \invariant This type has to derive from mln::Psite. + /*! \brief Point_Site associated type. + * \invariant This type has to derive from mln::Point_Site. */ typedef void psite; @@ -144,12 +144,12 @@ /*! \brief Forward point iterator associated type. - * \invariant This type has to derive from mln::Piter. + * \invariant This type has to derive from mln::Point_Iterator. */ typedef void fwd_piter; /*! \brief Backward point iterator associated type. - * \invariant This type has to derive from mln::Piter. + * \invariant This type has to derive from mln::Point_Iterator. */ typedef void bkd_piter; Index: mln/core/concept/doc/value_iterator.hh --- mln/core/concept/doc/value_iterator.hh (revision 1026) +++ mln/core/concept/doc/value_iterator.hh (working copy) @@ -25,8 +25,8 @@ // reasons why the executable file might be covered by the GNU General // Public License. -/*! \file mln/core/concept/doc/viter.hh - * \brief This file documents the concept of mln::Viter. +/*! \file mln/core/concept/doc/value_iterator.hh + * \brief This file documents the concept of mln::Value_Iterator. */ namespace mln @@ -35,12 +35,12 @@ namespace doc { - /*! \brief Documentation class for mln::Viter. + /*! \brief Documentation class for mln::Value_Iterator. * - * \see mln::Viter + * \see mln::Value_Iterator */ template <typename E> - struct Viter : public Iterator<E> + struct Value_Iterator : public Iterator<E> { /*! \brief Value associated type. */ Index: mln/core/concept/doc/value_set.hh --- mln/core/concept/doc/value_set.hh (revision 1026) +++ mln/core/concept/doc/value_set.hh (working copy) @@ -46,11 +46,11 @@ */ typedef void value; - /*! \brief Forward Viter associated type. + /*! \brief Forward Value_Iterator associated type. */ typedef void fwd_viter; - /*! \brief Backward Viter associated type. + /*! \brief Backward Value_Iterator associated type. */ typedef void bkd_viter; Index: mln/core/concept/doc/generalized_pixel.hh --- mln/core/concept/doc/generalized_pixel.hh (revision 1026) +++ mln/core/concept/doc/generalized_pixel.hh (working copy) @@ -28,8 +28,8 @@ #ifndef MLN_CORE_CONCEPT_DOC_GENPIXEL_HH # define MLN_CORE_CONCEPT_DOC_GENPIXEL_HH -/*! \file mln/core/concept/doc/genpixel.hh - * \brief This file documents the concept of mln::GenPixel. +/*! \file mln/core/concept/doc/generalized_pixel.hh + * \brief This file documents the concept of mln::Generalized_Pixel. */ @@ -39,12 +39,12 @@ namespace doc { - /*! \brief Documentation class for mln::GenPixel. + /*! \brief Documentation class for mln::Generalized_Pixel. * - * \see mln::GenPixel + * \see mln::Generalized_Pixel */ template <typename E> - struct GenPixel + struct Generalized_Pixel { /// Value associated type. @@ -70,7 +70,7 @@ value** address_() const; protected: - GenPixel(); + Generalized_Pixel(); }; } // end of namespace mln::doc Index: mln/core/concept/doc/pixel_iterator.hh --- mln/core/concept/doc/pixel_iterator.hh (revision 1026) +++ mln/core/concept/doc/pixel_iterator.hh (working copy) @@ -41,7 +41,7 @@ */ template <typename E> struct Pixel_Iterator : public Iterator<E>, - public GenPixel<E> + public Generalized_Pixel<E> { /*! \brief Type returned by the read-write dereference operator. */ Index: mln/core/concept/doc/point_iterator.hh --- mln/core/concept/doc/point_iterator.hh (revision 1026) +++ mln/core/concept/doc/point_iterator.hh (working copy) @@ -25,8 +25,8 @@ // reasons why the executable file might be covered by the GNU General // Public License. -/*! \file mln/core/concept/doc/piter.hh - * \brief This file documents the concept of mln::Piter. +/*! \file mln/core/concept/doc/point_iterator.hh + * \brief This file documents the concept of mln::Point_Iterator. */ namespace mln @@ -35,16 +35,16 @@ namespace doc { - /*! \brief Documentation class for mln::Piter. + /*! \brief Documentation class for mln::Point_Iterator. * - * \see mln::Piter + * \see mln::Point_Iterator */ template <typename E> - struct Piter : public Iterator<E>, - public GenPoint<E> + struct Point_Iterator : public Iterator<E>, + public Generalized_Point<E> { - /*! \brief Psite associated type. - * \invariant This type has to derive from mln::Psite. + /*! \brief Point_Site associated type. + * \invariant This type has to derive from mln::Point_Site. */ typedef void psite; Index: mln/core/concept/doc/generalized_point.hh --- mln/core/concept/doc/generalized_point.hh (revision 1026) +++ mln/core/concept/doc/generalized_point.hh (working copy) @@ -25,8 +25,8 @@ // reasons why the executable file might be covered by the GNU General // Public License. -/*! \file mln/core/concept/doc/genpoint.hh - * \brief This file documents the concept of mln::GenPoint. +/*! \file mln/core/concept/doc/generalized_point.hh + * \brief This file documents the concept of mln::Generalized_Point. */ namespace mln @@ -35,12 +35,12 @@ namespace doc { - /*! \brief Documentation class for mln::GenPoint. + /*! \brief Documentation class for mln::Generalized_Point. * - * \see mln::GenPoint + * \see mln::Generalized_Point */ template <typename E> - struct GenPoint + struct Generalized_Point { /*! \var dim * \brief Dimension of the space. Index: mln/core/concept/doc/point_set.hh --- mln/core/concept/doc/point_set.hh (revision 1026) +++ mln/core/concept/doc/point_set.hh (working copy) @@ -46,15 +46,15 @@ */ typedef void point; - /*! \brief Psite associated type. + /*! \brief Point_Site associated type. */ typedef void psite; - /*! \brief Forward Piter associated type. + /*! \brief Forward Point_Iterator associated type. */ typedef void fwd_piter; - /*! \brief Backward Piter associated type. + /*! \brief Backward Point_Iterator associated type. */ typedef void bkd_piter; Index: mln/core/concept/doc/window.hh --- mln/core/concept/doc/window.hh (revision 1026) +++ mln/core/concept/doc/window.hh (working copy) @@ -45,17 +45,17 @@ template <typename E> struct Window : public Object<E> { - /*! \brief Piter type associated to this window to browse its + /*! \brief Point_Iterator type associated to this window to browse its * points. */ typedef void qiter; - /*! \brief Piter type associated to this window to browse its + /*! \brief Point_Iterator type associated to this window to browse its * points in a forward way. */ typedef void fwd_qiter; - /*! \brief Piter type associated to this window to browse its + /*! \brief Point_Iterator type associated to this window to browse its * points in a backward way. */ typedef void bkd_qiter; Index: mln/core/concept/doc/neighborhood.hh --- mln/core/concept/doc/neighborhood.hh (revision 1026) +++ mln/core/concept/doc/neighborhood.hh (working copy) @@ -42,17 +42,17 @@ template <typename E> struct Neighborhood : public Object<E> { - /*! \brief Piter type associated to this neighborhood to browse + /*! \brief Point_Iterator type associated to this neighborhood to browse * neighbors. */ typedef void niter; - /*! \brief Piter type associated to this neighborhood to browse + /*! \brief Point_Iterator type associated to this neighborhood to browse * neighbors in a forward way. */ typedef void fwd_niter; - /*! \brief Piter type associated to this neighborhood to browse + /*! \brief Point_Iterator type associated to this neighborhood to browse * neighbors in a backward way. */ typedef void bkd_niter; Index: mln/core/concept/point_site.hh --- mln/core/concept/point_site.hh (revision 1026) +++ mln/core/concept/point_site.hh (working copy) @@ -28,11 +28,11 @@ #ifndef MLN_CORE_CONCEPT_PSITE_HH # define MLN_CORE_CONCEPT_PSITE_HH -/*! \file mln/core/concept/psite.hh - * \brief Definition of the concept of mln::Psite. +/*! \file mln/core/concept/point_site.hh + * \brief Definition of the concept of mln::Point_Site. */ -# include <mln/core/concept/genpoint.hh> +# include <mln/core/concept/generalized_point.hh> namespace mln @@ -46,7 +46,7 @@ * is either a point or designates a point: regular points, deriving * from mln::Point, are point sites, yet some point sites are not * points.) A point site has the behavior expected from every - * point; see mln::GenPoint. + * point; see mln::Generalized_Point. * * When a point site is not a point, it is automatically convertible * to the point it designates. @@ -61,8 +61,8 @@ * proper point site is a pair (index of run, index within the run). */ template <typename E> - struct Psite : public Object<E>, - public GenPoint<E> + struct Point_Site : public Object<E>, + public Generalized_Point<E> { /* const point* pointer_() const @@ -72,14 +72,14 @@ */ protected: - Psite(); + Point_Site(); }; # ifndef MLN_INCLUDE_ONLY template <typename E> - Psite<E>::Psite() + Point_Site<E>::Point_Site() { } Index: mln/core/vline2d.hh --- mln/core/vline2d.hh (revision 1026) +++ mln/core/vline2d.hh (working copy) @@ -62,12 +62,12 @@ /// Dpoint associated type. typedef dpoint2d dpoint; - /*! \brief Piter type to browse a vline such as: "for each row + /*! \brief Point_Iterator type to browse a vline such as: "for each row * (increasing), for each column (increasing)." */ typedef dpoints_fwd_piter<dpoint2d> fwd_qiter; - /*! \brief Piter type to browse a vline such as: "for each row + /*! \brief Point_Iterator type to browse a vline such as: "for each row * (decreasing), for each column (decreasing)." */ typedef internal::fixme bkd_qiter; Index: mln/level/approx/median.hh --- mln/level/approx/median.hh (revision 1026) +++ mln/level/approx/median.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_LEVEL_APPROX_MEDIAN_HH # define MLN_LEVEL_APPROX_MEDIAN_HH -/*! \file mln/level/median.hh +/*! \file mln/level/approx/median.hh * * \brief Approximates of some median filters of an image. */ Index: mln/value/lut_vec.hh --- mln/value/lut_vec.hh (revision 1026) +++ mln/value/lut_vec.hh (working copy) @@ -63,10 +63,10 @@ /// Value associated type. typedef T value; - /// Forward Viter associated type. + /// Forward Value_Iterator associated type. typedef fwd_viter_< lut_vec<S,T> > fwd_viter; - /// Backward Viter associated type. + /// Backward Value_Iterator associated type. typedef bkd_viter_< lut_vec<S,T> > bkd_viter; /// Give the \p i-th value. Index: mln/value/viter.hh --- mln/value/viter.hh (revision 1026) +++ mln/value/viter.hh (working copy) @@ -33,7 +33,7 @@ * \brief Definition of iterators on value sets. */ -# include <mln/core/concept/viter.hh> +# include <mln/core/concept/value_iterator.hh> @@ -49,7 +49,7 @@ * The parameter \c S is the type of value set. */ template <typename S> - struct fwd_viter_ : public Viter< fwd_viter_<S> > + struct fwd_viter_ : public Value_Iterator< fwd_viter_<S> > { /// Value associated type. typedef mln_value(S) value; @@ -86,7 +86,7 @@ * The parameter \c S is the type of value set. */ template <typename S> - struct bkd_viter_ : public Viter< bkd_viter_<S> > + struct bkd_viter_ : public Value_Iterator< bkd_viter_<S> > { /// Value associated type. typedef mln_value(S) value; Index: mln/value/set.hh --- mln/value/set.hh (revision 1026) +++ mln/value/set.hh (working copy) @@ -58,10 +58,10 @@ /// Value associated type. typedef T value; - /// Forward Viter associated type. + /// Forward Value_Iterator associated type. typedef fwd_viter_< set_<T> > fwd_viter; - /// Backward Viter associated type. + /// Backward Value_Iterator associated type. typedef bkd_viter_< set_<T> > bkd_viter; /// Test if \p v belongs to this set: always true!

1 0

1026: Handle low quantization and fast properties.
by Thierry Geraud 20 Jul '07

20 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Thierry Geraud <thierry.geraud(a)lrde.epita.fr> Handle low quantization and fast properties. * tests/transform.cc: New. * tests/fill.cc: New. * mln/core/macros.hh (mln_value_kind, mln_vset): New. * mln/core/trait/is_fast.hh: New. * mln/core/trait/all.hh: New. * mln/core/trait/is_lowq.hh: New. * mln/core/internal/image_base.hh (select_image_concept_): New. * mln/core/concept/image.hh (vset, values): New. Update. * mln/core/concept/fast_image.hh (operator[]): New. * mln/metal/math.hh: New. * mln/metal/bool.hh: New. * mln/metal/bexpr.hh: New. * mln/value/lut_vec.hh: New. * mln/value/int_u.hh (nbits, card): New. * tests/README (g++-2.95): New hints. * tests/pixter_dpoint2d.cc: Fix warning. * tests/erosion.cc: Fix doc. * doc/Doxyfile.in: Update. * mln/debug/iota.hh: New overload. * mln/fun/pw_value.hh (function_): Rename as... (select_function_): ...this. * mln/core/trait/pixter.hh (pixter, qixter, nixter): Remove; obsolete. * mln/core/concept/value_set.hh (viter): Likewise. * mln/core/concept/doc/image.hh, * mln/core/concept/doc/value_set.hh, * mln/core/concept/doc/fast_image.hh: Update. * mln/core/image2d_b.hh: Update. * mln/metal/none.hh: Fix doc. * mln/morpho/erosion.hh: Update. * mln/level/fill.hh: Add todos. New overload. * mln/level/fast_median.hh: Fix doc. * mln/level/transform.hh: New overload. * mln/value/props.hh (min): Change to var. * mln/value/kind.hh: Fix doc. * mln/value/set.hh (viter): Remove; obsolete. doc/Doxyfile.in | 1 mln/core/concept/doc/fast_image.hh | 25 +++++ mln/core/concept/doc/image.hh | 17 +++ mln/core/concept/doc/value_set.hh | 4 mln/core/concept/fast_image.hh | 20 +++- mln/core/concept/image.hh | 6 + mln/core/concept/value_set.hh | 2 mln/core/image2d_b.hh | 65 +++++++++++++- mln/core/internal/image_base.hh | 20 ++++ mln/core/macros.hh | 12 ++ mln/core/tags.hh | 1 mln/core/trait/all.hh | 52 +++++++++++ mln/core/trait/is_fast.hh | 60 +++++++++++++ mln/core/trait/is_lowq.hh | 65 ++++++++++++++ mln/core/trait/pixter.hh | 15 --- mln/debug/iota.hh | 23 +++++ mln/fun/pw_value.hh | 6 - mln/level/fast_median.hh | 2 mln/level/fill.hh | 75 ++++++++++++++-- mln/level/transform.hh | 56 +++++++++++- mln/metal/bexpr.hh | 102 ++++++++++++++++++++++ mln/metal/bool.hh | 76 +++++++++++++++++ mln/metal/math.hh | 63 ++++++++++++++ mln/metal/none.hh | 2 mln/morpho/erosion.hh | 2 mln/value/int_u.hh | 3 mln/value/kind.hh | 2 mln/value/lut_vec.hh | 165 +++++++++++++++++++++++++++++++++++++ mln/value/props.hh | 85 ++++++++++--------- mln/value/set.hh | 11 -- tests/README | 4 tests/erosion.cc | 4 tests/fill.cc | 63 ++++++++++++++ tests/pixter_dpoint2d.cc | 2 tests/transform.cc | 66 ++++++++++++++ 35 files changed, 1066 insertions(+), 111 deletions(-) Index: tests/README --- tests/README (revision 1025) +++ tests/README (working copy) @@ -1 +1,3 @@ -g++-4.1 -I.. -ansi -pedantic -W -Wall -Wextra -Wconversion sample.cc +g++-4.1 -ansi -pedantic -W -Wall -Wextra -Wconversion -I.. sample.cc + +g++-2.95 -ansi -pedantic -W -Wall -Wconversion -ftemplate-depth-51 -I.. sample.cc \ No newline at end of file Index: tests/transform.cc --- tests/transform.cc (revision 0) +++ tests/transform.cc (revision 0) @@ -0,0 +1,66 @@ +// Copyright (C) 2007 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. + +/*! \file tests/transform.cc + * + * \brief Tests on mln::level::transform + */ + +#include <cmath> + +#include <mln/core/image2d_b.hh> +#include <mln/level/transform.hh> +#include <mln/debug/iota.hh> + + +struct mysqrt : mln::Function_v2v<mysqrt> +{ + typedef unsigned short result; + result operator()(unsigned short c) const + { + return result( std::sqrt(float(c)) ); + } +}; + + + +int main() +{ + using namespace mln; + + const unsigned size = 10000; + image2d_b<unsigned short> + ima(size, size); + + (std::cout << "iota... ").flush(); + debug::iota(ima); + std::cout << "done" << std::endl; + + (std::cout << "transform... ").flush(); + level::transform(ima, mysqrt(), ima); + std::cout << "done" << std::endl; +} Index: tests/pixter_dpoint2d.cc --- tests/pixter_dpoint2d.cc (revision 1025) +++ tests/pixter_dpoint2d.cc (working copy) @@ -67,7 +67,7 @@ qixter qix(ima, win, p); for_all(p) - if (p[0] > 0 && p[1] > 0 && p[0] < size - 1 && p[1] < size - 1) + if (p[0] > 0 && p[1] > 0 && p[0] < int(size - 1) && p[1] < int(size - 1)) for_all(qix) mln_assertion(*qix = value); } Index: tests/erosion.cc --- tests/erosion.cc (revision 1025) +++ tests/erosion.cc (working copy) @@ -25,9 +25,9 @@ // reasons why the executable file might be covered by the GNU General // Public License. -/*! \file tests/median.cc +/*! \file tests/erosion.cc * - * \brief Test on mln::level::median. + * \brief Test on mln::morpho::erosion. */ #include <mln/core/image2d_b.hh> Index: tests/fill.cc --- tests/fill.cc (revision 0) +++ tests/fill.cc (revision 0) @@ -0,0 +1,63 @@ +// Copyright (C) 2007 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. + +/*! \file tests/fill.cc + * + * \brief Tests on mln::level::fill + */ + +#include <mln/core/image2d_b.hh> +#include <mln/level/fill.hh> + +#include <mln/debug/println.hh> +#include <mln/value/props.hh> + + +int main() +{ + using namespace mln; + + + unsigned u = 300; + unsigned char uc = u; + mln_assertion(uc = 44); + +// { +// const unsigned size = 3; +// image2d_b<unsigned> ima(size, size); +// level::fill(ima, u); +// debug::println(ima); +// } + + { + const unsigned size = 10000; + image2d_b<unsigned char> ima(size, size); + for (unsigned i = 0; i < 5; ++i) + level::fill(ima, uc); + } + +} Index: doc/Doxyfile.in --- doc/Doxyfile.in (revision 1025) +++ doc/Doxyfile.in (working copy) @@ -1055,6 +1055,7 @@ "mln_fwd_viter(T)=typename T::fwd_viter" \ "mln_bkd_viter(T)=typename T::bkd_viter" \ "mln_value(T)=typename T::value" \ + "mln_vset(T)=typename T::vset" \ "mln_rvalue(T)=typename T::rvalue" \ "mln_lvalue(T)=typename T::lvalue" \ "mln_coord(T)=typename T::coord" \ Index: mln/debug/iota.hh --- mln/debug/iota.hh (revision 1025) +++ mln/debug/iota.hh (working copy) @@ -49,6 +49,9 @@ # ifndef MLN_INCLUDE_ONLY + namespace impl + { + template <typename I> void iota(Image& input_) { @@ -59,6 +62,26 @@ input(p) = ++i; } + template <typename I> + void iota(Fast_Image& input_) + { + unsigned i = 0; + I& input = exact(input_); + mln_pixter(I) p(input); + for_all(p) + *p = ++i; + } + + } // end of namespace mln::debug::impl + + + template <typename I> + void iota(Image& input) + { + mln_precondition(exact(input).has_data()); + impl::iota(exact(input)); + } + # endif // ! MLN_INCLUDE_ONLY } // end of namespace mln::debug Index: mln/fun/pw_value.hh --- mln/fun/pw_value.hh (revision 1025) +++ mln/fun/pw_value.hh (working copy) @@ -63,11 +63,11 @@ { template <typename K, typename E> - struct function_ : Function_p2v<E> + struct select_function_ : Function_p2v<E> {}; template <typename E> - struct function_< value::binary_kind, E > : Function_p2b<E> + struct select_function_< value::binary_kind, E > : Function_p2b<E> {}; } // end of namespace mln::fun::internal @@ -76,7 +76,7 @@ // FIXME: Doc! template <typename I> - struct pw_value : public internal::function_< mln_kind(I), pw_value > + struct pw_value : public internal::select_function_< mln_value_kind(I), pw_value > { typedef mln_value(I) result; pw_value(const I& ima); Index: mln/core/macros.hh --- mln/core/macros.hh (revision 1025) +++ mln/core/macros.hh (working copy) @@ -105,14 +105,14 @@ // p -/// Shortcut to access the piter type associated to T. +/// Shortcut to access the type of point iterator (piter) associated to T. # define mln_piter(T) typename T::fwd_piter # define mln_piter_(T) T::fwd_piter -/// Shortcut to access the pset type associated to T. +/// Shortcut to access the type of point set (pset) associated to T. # define mln_pset(T) typename T::pset -/// Shortcut to access the psite type associated to T. +/// Shortcut to access the type of point site (psite) associated to T. # define mln_psite(T) typename T::psite /// Shortcut to access the point type associated to T. @@ -145,6 +145,12 @@ /// Shortcut to access the value type associated to T. # define mln_value(T) typename T::value +/// Shortcut to the kind of values for an image with type \c I. +# define mln_value_kind(T) typename mln::value::props< mln_value(I) >::kind + +/// Shortcut to access the type of value set (vset) associated to T. +# define mln_vset(T) typename T::vset + /// Shortcut to access the viter type associated to T. # define mln_viter(T) typename T::fwd_viter Index: mln/core/trait/is_fast.hh --- mln/core/trait/is_fast.hh (revision 0) +++ mln/core/trait/is_fast.hh (revision 0) @@ -0,0 +1,60 @@ +// Copyright (C) 2007 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_CORE_TRAIT_IS_FAST_HH +# define MLN_CORE_TRAIT_IS_FAST_HH + +/*! \file mln/core/trait/is_fast.hh + * + * \brief Definition of the is_fast image trait. + */ + +# include <mln/metal/bool.hh> + + +namespace mln +{ + + namespace trait + { + + + // FIXME: Doc! + + template <typename I> + struct is_fast + { + typedef metal::false_ ret; + }; + + + } // end of namespace mln::trait + +} // end of namespace mln + + +#endif // ! MLN_CORE_TRAIT_IS_FAST_HH Index: mln/core/trait/all.hh --- mln/core/trait/all.hh (revision 0) +++ mln/core/trait/all.hh (revision 0) @@ -0,0 +1,52 @@ +// Copyright (C) 2007 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_CORE_TRAIT_ALL_HH +# define MLN_CORE_TRAIT_ALL_HH + +/*! \file mln/core/trait/all.hh + * + * \brief File that includes all traits. + */ + + +namespace mln +{ + + /*! Namespace for image traits. + */ + namespace trait {} + +} + + +# include <mln/core/trait/is_fast.hh> +# include <mln/core/trait/is_lowq.hh> +# include <mln/core/trait/pixter.hh> + + +#endif // ! MLN_CORE_TRAIT_ALL_HH Index: mln/core/trait/pixter.hh --- mln/core/trait/pixter.hh (revision 1025) +++ mln/core/trait/pixter.hh (working copy) @@ -57,11 +57,6 @@ typedef metal::none ret; }; - template <typename I> - struct pixter : fwd_pixter - { - }; - // qixter @@ -77,11 +72,6 @@ typedef metal::none ret; }; - template <typename I, typename W> - struct qixter : fwd_qixter<I, W> - { - }; - // nixter @@ -97,11 +87,6 @@ typedef metal::none ret; }; - template <typename I, typename N> - struct nixter : fwd_nixter<I, N> - { - }; - } // end of namespace mln::trait Index: mln/core/trait/is_lowq.hh --- mln/core/trait/is_lowq.hh (revision 0) +++ mln/core/trait/is_lowq.hh (revision 0) @@ -0,0 +1,65 @@ +// Copyright (C) 2007 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_CORE_TRAIT_IS_LOWQ_HH +# define MLN_CORE_TRAIT_IS_LOWQ_HH + +/*! \file mln/core/trait/is_lowq.hh + * + * \brief Definition of the is_lowq image trait. + */ + +# include <mln/metal/bool.hh> +# include <mln/value/props.hh> + + +# define mln_is_lowq(I) typename mln::trait::is_lowq::ret + + + +namespace mln +{ + + namespace trait + { + + + // FIXME: Doc! + + template <typename I> + struct is_lowq + { + typedef typename metal::bool_<( mln_card(mln_value(I)) != 0 )>::type ret; + }; + + + } // end of namespace mln::trait + +} // end of namespace mln + + +#endif // ! MLN_CORE_TRAIT_IS_LOWQ_HH Index: mln/core/internal/image_base.hh --- mln/core/internal/image_base.hh (revision 1025) +++ mln/core/internal/image_base.hh (working copy) @@ -42,12 +42,30 @@ namespace internal { + + + template <typename Is_fast, typename E> + struct select_image_concept_; + + template <typename E> + struct select_image_concept_< metal::true_, E > + : public Fast_Image<E> + {}; + + template <typename E> + struct select_image_concept_< metal::false_, E > + : public Image<E> + {}; + + + /*! \brief A base class for images. * * \internal */ template <typename S, typename E> - struct image_base_ : public Image<E> + struct image_base_ : public select_image_concept_< typename trait::is_fast<E>::ret, + E > { /// Point_Set associated type. typedef S pset; Index: mln/core/concept/image.hh --- mln/core/concept/image.hh (revision 1025) +++ mln/core/concept/image.hh (working copy) @@ -33,6 +33,7 @@ */ # include <mln/core/concept/point_set.hh> +# include <mln/core/trait/all.hh> namespace mln @@ -52,8 +53,10 @@ typedef value; typedef rvalue; typedef lvalue; + typedef vset; bool has_data() const; + const vset& values() const; bool owns_(const psite& p) const; const pset& domain() const; @@ -138,4 +141,7 @@ } // end of namespace mln +# include <mln/core/concept/fast_image.hh> + + #endif // ! MLN_CORE_CONCEPT_IMAGE_HH Index: mln/core/concept/value_set.hh --- mln/core/concept/value_set.hh (revision 1025) +++ mln/core/concept/value_set.hh (working copy) @@ -48,7 +48,6 @@ { /* typedef value; - typedef viter; typedef fwd_viter; typedef bkd_viter; @@ -71,7 +70,6 @@ Value_Set<E>::Value_Set() { typedef mln_value(E) value; - typedef mln_viter(E) viter; typedef mln_fwd_viter(E) fwd_viter; typedef mln_bkd_viter(E) bkd_viter; Index: mln/core/concept/fast_image.hh --- mln/core/concept/fast_image.hh (revision 1025) +++ mln/core/concept/fast_image.hh (working copy) @@ -33,7 +33,6 @@ */ # include <mln/core/concept/image.hh> -# include <mln/core/trait/pixter.hh> namespace mln @@ -55,6 +54,9 @@ point point_at_offset(unsigned o) const; const value* buffer() const; + + rvalue operator[](unsigned o) const; + lvalue operator[](unsigned o); */ protected: @@ -69,7 +71,6 @@ { typedef mln_point(E) point; typedef mln_dpoint(E) dpoint; - typedef mln_value(E) value; typedef mln_fwd_pixter(E) fwd_pixter; typedef mln_bkd_pixter(E) bkd_pixter; @@ -80,11 +81,22 @@ m2 = 0; point (E::*m3)(unsigned) const = & E::point_at_offset; m3 = 0; - const value* (E::*m4)() const = & E::buffer; + unsigned (E::*m4)() const = & E::border; m4 = 0; - unsigned (E::*m5)() const = & E::border; + + typedef mln_value(E) value; + + const value* (E::*m5)() const = & E::buffer; m5 = 0; + typedef mln_rvalue(E) rvalue; + typedef mln_lvalue(E) lvalue; + + rvalue (E::*m6)(unsigned) const = & E::operator[]; + m6 = 0; + lvalue (E::*m7)(unsigned) = & E::operator[]; + m7 = 0; + // FIXME: how to check that qixter are defined when W is unknown! } Index: mln/core/concept/doc/image.hh --- mln/core/concept/doc/image.hh (revision 1025) +++ mln/core/concept/doc/image.hh (working copy) @@ -56,6 +56,11 @@ */ typedef void lvalue; + /*! \brief Value set associated type. + * \invariant This type has to derive from mln::Value_Set. + */ + typedef void vset; + /*! \brief Test if the image have been initialized. */ bool has_data() const; @@ -73,13 +78,19 @@ */ const pset& domain() const; - /*! \brief Read access to the image value located at \p p. + /*! \brief Give the set of values of the image. + * + * \return A reference to the value set. + */ + const vset& values() const; + + /*! \brief Read-only access to the image value located at \p p. * * \param[in] p A point site. * * \pre The image has to own the site \p p. * - * \return The value at \p p. + * \return The value at \p p (not assignable). */ rvalue operator()(const psite& p) const; @@ -89,7 +100,7 @@ * * \pre The image has to own the site \p p. * - * \return The value at \p p. + * \return The value at \p p (assignable). */ lvalue operator()(const psite& p); Index: mln/core/concept/doc/value_set.hh --- mln/core/concept/doc/value_set.hh (revision 1025) +++ mln/core/concept/doc/value_set.hh (working copy) @@ -46,10 +46,6 @@ */ typedef void value; - /*! \brief Viter associated type. - */ - typedef void viter; - /*! \brief Forward Viter associated type. */ typedef void fwd_viter; Index: mln/core/concept/doc/fast_image.hh --- mln/core/concept/doc/fast_image.hh (revision 1025) +++ mln/core/concept/doc/fast_image.hh (working copy) @@ -51,7 +51,7 @@ /*! \brief Give the offset corresponding to the delta-point \p * dp. * - * \param[in] p A delta-point. + * \param[in] dp A delta-point. * * \pre The image has to be initialized. */ @@ -69,7 +69,7 @@ /*! \brief Give the point at offset \p o. * - * \param[in] p An offset. + * \param[in] o An offset. * * \pre The image has to be initialized. * \pre o < ncells() @@ -85,6 +85,27 @@ */ const value* buffer() const; + /*! \brief Read-only access to the image value at offset \p o. + * + * \param[in] o An offset. + * + * \pre o < ncells() + * + * \return The value at \p o (not assignable). + */ + rvalue operator[](unsigned o) const; + + + /*! \brief Read-write access to the image value at offset \p o. + * + * \param[in] o An offset. + * + * \pre o < ncells() + * + * \return The value at \p o (assignable). + */ + lvalue operator[](unsigned o); + }; } // end of namespace mln::doc Index: mln/core/tags.hh --- mln/core/tags.hh (revision 1025) +++ mln/core/tags.hh (working copy) @@ -42,6 +42,7 @@ // FIXME: ... + } // end of namespace mln::core } // end of namespace mln Index: mln/core/image2d_b.hh --- mln/core/image2d_b.hh (revision 1025) +++ mln/core/image2d_b.hh (working copy) @@ -38,6 +38,7 @@ # include <mln/core/box2d.hh> # include <mln/border/thickness.hh> +# include <mln/value/set.hh> # include <mln/fun/all.hh> @@ -50,6 +51,23 @@ namespace mln { + // Fwd decl. + template <typename T> struct image2d_b; + + + namespace trait + { + + template <typename T> + struct is_fast< image2d_b<T> > + { + typedef metal::true_ ret; + }; + + } // end of mln::trait + + + /*! \brief Basic 2D image class. * * The parameter \c T is the type of pixel values. This image class @@ -71,6 +89,7 @@ // end of warning + // FIXME: // /// Forward pixel iterator associated to image2d @@ -99,6 +118,9 @@ typedef image2d_b ret; }; + /// Value_Set associated type. + typedef mln::value::set_<T> vset; + /// Constructor without argument. image2d_b(); @@ -134,6 +156,9 @@ /// Test if this image has been initialized. bool has_data() const; + /// Give the set of values of the image. + const vset& values() const; + /// Give the definition domain. const box2d& domain() const; @@ -143,12 +168,18 @@ /// Give the number of cells (points including border ones). std::size_t ncells() const; - /// Read-only access to the image value located at \p p. + /// Read-only access to the image value located at point \p p. const T& operator()(const point2d& p) const; - /// Read-write access to the image value located at \p p. + /// Read-write access to the image value located at point \p p. T& operator()(const point2d& p); + /// Read-only access to the image value located at offset \p o. + const T& operator[](unsigned o) const; + + /// Read-write access to the image value located at offset \p o. + T& operator[](unsigned o); + /// Read-only access to the image value located at (\p row, \p col). const T& at(int row, int col) const; @@ -279,6 +310,13 @@ } template <typename T> + const typename image2d_b<T>::vset& + image2d_b<T>::values() const + { + return vset::the(); + } + + template <typename T> const box2d& image2d_b<T>::domain() const { @@ -328,6 +366,22 @@ template <typename T> const T& + image2d_b<T>::operator[](unsigned o) const + { + mln_precondition(o < ncells()); + return *(buffer_ + o); + } + + template <typename T> + T& + image2d_b<T>::operator[](unsigned o) + { + mln_precondition(o < ncells()); + return *(buffer_ + o); + } + + template <typename T> + const T& image2d_b<T>::at(int row, int col) const { mln_precondition(this->owns_(make::point2d(row, col))); @@ -380,7 +434,6 @@ image2d_b<T>::point_at_offset(unsigned o) const { mln_precondition(o < ncells()); - mln_precondition(this->has_data()); point2d p = make::point2d(o / vb_.len(1) + vb_.min_row(), o % vb_.len(1) + vb_.min_col()); mln_postcondition(& this->operator()(p) = this->buffer_ + o); @@ -462,6 +515,12 @@ }; template <typename T> + struct fwd_pixter< const image2d_b<T> > + { + typedef fwd_pixter2d_b< const image2d_b<T> > ret; + }; + + template <typename T> struct bkd_pixter< image2d_b<T> > { typedef internal::fixme ret; Index: mln/metal/math.hh --- mln/metal/math.hh (revision 0) +++ mln/metal/math.hh (revision 0) @@ -0,0 +1,63 @@ +// Copyright (C) 2007 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_METAL_MATH_HH +# define MLN_METAL_MATH_HH + +/*! \file mln/metal/math.hh + * + * \brief Definition of some mathematical static functions. + */ + + +namespace mln +{ + + namespace metal + { + + // pow<x, n> + + template <int x, unsigned n> + struct pow + { + enum { value = x * pow<x, n-1>::value }; + }; + + template <int x> + struct pow< x, 0 > + { + enum { value = 1 }; + }; + + + } // end of namespace mln::metal + +} // end of namespace mln + + +#endif // ! MLN_METAL_MATH_HH Index: mln/metal/bool.hh --- mln/metal/bool.hh (revision 0) +++ mln/metal/bool.hh (revision 0) @@ -0,0 +1,76 @@ +// Copyright (C) 2007 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_METAL_BOOL_HH +# define MLN_METAL_BOOL_HH + +/*! \file mln/metal/bool.hh + * + * \brief Definition of a Boolean value type. + */ + + +namespace mln +{ + + namespace metal + { + + // Fwd decls. + struct true_; + struct false_; + + + + // FIXME: Doc! + + template <bool b> struct bool_; + + template <> + struct bool_< true > + { + typedef true_ type; + static const bool value = true; + }; + + template <> + struct bool_< false > + { + typedef false_ type; + static const bool value = false; + }; + + + } // end of namespace mln::metal + +} // end of namespace mln + + +# include <mln/metal/bexpr.hh> + + +#endif // ! MLN_METAL_BOOL_HH Index: mln/metal/none.hh --- mln/metal/none.hh (revision 1025) +++ mln/metal/none.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_CORE_METAL_NONE_HH # define MLN_CORE_METAL_NONE_HH -/*! \file mln/core/metal/none.hh +/*! \file mln/metal/none.hh * * \brief Definition of a type that means "none". */ Index: mln/metal/bexpr.hh --- mln/metal/bexpr.hh (revision 0) +++ mln/metal/bexpr.hh (revision 0) @@ -0,0 +1,102 @@ +// Copyright (C) 2007 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_METAL_BEXPR_HH +# define MLN_METAL_BEXPR_HH + +/*! \file mln/metal/bexpr.hh + * + * \brief Definition of types for static Boolean expressions. + */ + +# include <mln/metal/bool.hh> + + + +namespace mln +{ + + namespace metal + { + + + /// "true" type. + struct true_ + { + static void check(); + typedef true_ eval; + enum { to_bool = true }; + }; + + + /// "false" type. + struct false_ + { + typedef false_ eval; + enum { to_bool = false }; + }; + + + /// Negate type. + template <typename B> + struct not_ : bool_<( ! B::to_bool )>::type + {}; + + + /// And type. + template <typename L, typename R> + struct and_ : bool_<( L::to_bool && R::to_bool )>::type + {}; + + + /// Or type. + template <typename L, typename R> + struct or_ : bool_<( L::to_bool || R::to_bool )>::type + {}; + + + /// Xor type. + template <typename L, typename R> + struct xor_ : bool_<( L::to_bool ^ R::to_bool )>::type + {}; + + + +# ifndef MLN_INCLUDE_ONLY + + void true_::check() + { + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::metal + +} // end of namespace mln + + +#endif // ! MLN_METAL_BOOL_HH Index: mln/morpho/erosion.hh --- mln/morpho/erosion.hh (revision 1025) +++ mln/morpho/erosion.hh (working copy) @@ -132,7 +132,7 @@ template <typename I, typename W, typename O> void erosion_wrt_win(const Image& input, const Window<W>& win, Image<O>& output) { - erosion_wrt_value(mln_kind(I)(), exact(input), exact(win), output); + erosion_wrt_value(mln_value_kind(I)(), exact(input), exact(win), output); // | // --> call stage 2: dispatch w.r.t. the value kind } Index: mln/level/fill.hh --- mln/level/fill.hh (revision 1025) +++ mln/level/fill.hh (working copy) @@ -33,6 +33,8 @@ * \brief Fill an image, that is, set pixel values. */ +# include <cstring> + # include <mln/core/concept/image.hh> # include <mln/core/concept/function.hh> @@ -49,6 +51,8 @@ * \param[in] v The value to assign to all pixels. * * \pre \p ima has to be initialized. + * + * \todo Optimize when \p ima is large and sizeof(mln_value(I)) > 1. */ template <typename I> void fill(Image& ima, const mln_value(I)& v); @@ -60,6 +64,8 @@ * \param[in] f The function. * * \pre \p ima has to be initialized. + * + * \todo Take benefit from quantization when possible. */ template <typename I, typename F> void fill(Image& ima, const Function_p2v<F>& f); @@ -74,6 +80,8 @@ * " value f(const point& p) " * * \pre \p ima has to be initialized. + * + * \todo Take benefit from quantization when possible. */ template <typename I> void fill(Image& ima, mln_value(I) (*f)(const mln_point(I)& p)); @@ -103,28 +111,70 @@ * the one of \p data. * * \pre \p ima.domain <= \p data.domain. + * + * \todo Use memcpy when possible. */ template <typename I, typename J> void fill(Image& ima, const Image<J>& data); + # ifndef MLN_INCLUDE_ONLY + namespace impl + { + + + // fill_with_value + template <typename I> - void fill(Image& ima_, - const mln_value(I)& value) + void fill_with_value(Image& ima_, const mln_value(I)& value, + bool force_call = false) { + force_call = false; // just to avoid warning ("unused param") I& ima = exact(ima_); - mln_precondition(ima.has_data()); mln_piter(I) p(ima.domain()); for_all(p) ima(p) = value; } + template <typename I> + void fill_with_value(Fast_Image& ima_, const mln_value(I)& value) + { + I& ima = exact(ima_); + if (sizeof(mln_value(I)) = 1) + { + std::memset((void*)(ima.buffer()), + value, + sizeof(mln_value(I)) * ima.ncells()); + } + else + { + // FIXME: Use memcpy on a chunck when image size is large! + fill_with_value(ima, value, true); + } + } + + + } // end of namespace mln::level::impl + + + + // with: value + + template <typename I> + void fill(Image& ima, const mln_value(I)& value) + { + mln_precondition(exact(ima).has_data()); + impl::fill_with_value(exact(ima), value); + } + + + // with: Function_p2v<F> + template <typename I, typename F> - void fill(Image& ima_, - const Function_p2v<F>& f_) + void fill(Image& ima_, const Function_p2v<F>& f_) { I& ima = exact(ima_); mln_precondition(ima.has_data()); @@ -134,6 +184,9 @@ ima(p) = f(p); } + + // with: value f(const point&) + template <typename I> void fill(Image& ima_, mln_value(I) (*f)(const mln_point(I)& p)) @@ -145,9 +198,11 @@ ima(p) = f(p); } + + // with: value arr[N] + template <typename I, unsigned N> - void fill(Image& ima_, - mln_value(I) (&arr)[N]) + void fill(Image& ima_, mln_value(I) (&arr)[N]) { I& ima = exact(ima_); mln_precondition(ima.has_data()); @@ -158,9 +213,11 @@ ima(p) = arr[i++]; } + + // with: Image<J> + template <typename I, typename J> - void fill(Image& ima_, - const Image<J>& data_) + void fill(Image& ima_, const Image<J>& data_) { I& ima = exact(ima_); const J& data = exact(data_); Index: mln/level/fast_median.hh --- mln/level/fast_median.hh (revision 1025) +++ mln/level/fast_median.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_LEVEL_FAST_MEDIAN_HH # define MLN_LEVEL_FAST_MEDIAN_HH -/*! \file mln/level/median.hh +/*! \file mln/level/fast_median.hh * * \brief Fast Median filtering of an image. */ Index: mln/level/transform.hh --- mln/level/transform.hh (revision 1025) +++ mln/level/transform.hh (working copy) @@ -36,6 +36,9 @@ # include <mln/core/concept/image.hh> # include <mln/core/concept/function.hh> +# include <mln/value/set.hh> +# include <mln/value/lut_vec.hh> + namespace mln { @@ -61,20 +64,67 @@ # ifndef MLN_INCLUDE_ONLY + namespace impl + { + template <typename I, typename F, typename O> - void transform(const Image& input_, const Function_v2v<F>& f_, Image<O>& output_) + void transform(metal::false_, // general case + const Image& input_, const Function_v2v<F>& f_, Image<O>& output_) { const I& input = exact(input_); const F& f = exact(f_); O& output = exact(output_); - mln_precondition(output.domain() >= input.domain()); - mln_piter(I) p(input.domain()); for_all(p) output(p) = f( input(p) ); } + template <typename I, typename F, typename O> + void transform(metal::true_, // low quantization + const Image& input_, const Function_v2v<F>& f_, Image<O>& output_) + { + const I& input = exact(input_); + const F& f = exact(f_); + O& output = exact(output_); + + value::lut_vec<mln_vset(I), mln_result(F)> lut(input.values(), f); + mln_piter(I) p(input.domain()); + for_all(p) + output(p) = lut(input(p)); + } + +// template <typename I, typename F, typename O> +// void transform(metal::true_, // low quantization +// const Fast_Image& input_, const Function_v2v<F>& f_, Image<O>& output_) +// { +// const I& input = exact(input_); +// const F& f = exact(f_); +// O& output = exact(output_); + +// value::lut_vec<mln_vset(I), mln_result(F)> lut(input.values(), f); +// mln_pixter(const I) pi(input); // FIXME +// mln_pixter(O) po(output); +// po.start(); +// for_all(pi) +// { +// *po = lut(*pi); +// po.next(); +// } +// } + + } // end of namespace mln::level::impl + + + template <typename I, typename F, typename O> + void transform(const Image& input, const Function_v2v<F>& f, Image<O>& output) + { + mln_precondition(exact(output).domain() >= exact(input).domain()); + impl::transform(mln_is_lowq(I)(), + exact(input), exact(f), exact(output)); + } + + # endif // ! MLN_INCLUDE_ONLY } // end of namespace mln::level Index: mln/value/lut_vec.hh --- mln/value/lut_vec.hh (revision 0) +++ mln/value/lut_vec.hh (revision 0) @@ -0,0 +1,165 @@ +// Copyright (C) 2007 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_VALUE_LUT_VEC_HH +# define MLN_VALUE_LUT_VEC_HH + +/*! \file mln/value/lut_vec.hh + * + * \brief Define some basic lut_vecs of values from value types. + */ + +# include <vector> + +# include <mln/core/concept/value_set.hh> +# include <mln/core/concept/function.hh> +# include <mln/value/props.hh> + + +namespace mln +{ + + namespace value + { + + // Fwd decls. + template <typename S> struct fwd_viter_; + template <typename S> struct bkd_viter_; + + + /*! Class that defines FIXME + * + * \warning This is a multi-set!!! + * FIXME + * + */ + template <typename S, typename T> + struct lut_vec : public Value_Set< lut_vec<S,T> > + { + /// Value associated type. + typedef T value; + + /// Forward Viter associated type. + typedef fwd_viter_< lut_vec<S,T> > fwd_viter; + + /// Backward Viter associated type. + typedef bkd_viter_< lut_vec<S,T> > bkd_viter; + + /// Give the \p i-th value. + T operator[](std::size_t i) const; + + /// Give the number of values. + std::size_t nvalues() const; + + // Apply the look-up-table. FIXME: Doc! + T operator()(const mln_value(S)& val) const; + + /// Test if \p v belongs to this set. + bool has(const value& v) const; + + /// Give the index of value \p v in this set. + std::size_t index_of(const value& v) const; + + /// Ctor. FIXME! + template <typename F> + lut_vec(const S& vset, const Function_v2v<F>& f); + + protected: + + const S& vset_; + std::vector<T> vec_; + unsigned n_; + }; + + + +# ifndef MLN_INCLUDE_ONLY + + template <typename S, typename T> + bool + lut_vec<S,T>::has(const T&) const + { + mln_invariant(0); // FIXME + return false; + } + + template <typename S, typename T> + std::size_t + lut_vec<S,T>::index_of(const T& v) const + { + mln_invariant(0); // FIXME + return 0; + } + + template <typename S, typename T> + template <typename F> + lut_vec<S,T>::lut_vec(const S& vset, const Function_v2v<F>& f) + : vset_(vset) + { + const F& f_ = exact(f); + n_ = vset.nvalues(); + vec_.reserve(n_); + for (unsigned i = 0; i < n_; ++i) + vec_[i] = f_(vset[i]); + } + + template <typename S, typename T> + T + lut_vec<S,T>::operator()(const mln_value(S)& val) const + { + unsigned i = vset_.index_of(val); + mln_precondition(i < n_); + return vec_[vset_.index_of(val)]; + } + + template <typename S, typename T> + T + lut_vec<S,T>::operator[](std::size_t i) const + { + mln_precondition(i < nvalues()); + return vec_[i]; + } + + template <typename S, typename T> + std::size_t + lut_vec<S,T>::nvalues() const + { + return vec_.size(); + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::value + + +} // end of namespace mln + + +# include <mln/value/viter.hh> + + +#endif // ! MLN_VALUE_LUT_VEC_HH Index: mln/value/props.hh --- mln/value/props.hh (revision 1025) +++ mln/value/props.hh (working copy) @@ -41,19 +41,19 @@ /// Get the minimum value of type \c T. -# define mln_min(T) mln::value::props<T>::min() +# define mln_min(T) mln::value::props< T >::min /// Get the maximum value of type \c T. -# define mln_max(T) mln::value::props<T>::max() +# define mln_max(T) mln::value::props< T >::max /// Get the number of values for value type \c T. -# define mln_card(T) mln::value::props<T>::card() +# define mln_card(T) mln::value::props< T >::card -/// Get the kind of the value type of image \c I. -# define mln_kind(I) typename mln::value::props< mln_value(I) >::kind +/// Get the kind of value type \c T. +# define mln_kind(T) typename mln::value::props< T >::kind @@ -71,119 +71,124 @@ struct props { /// Minimum value for type \c T. - static T min(); + static const T min; /// Maximum value for type \c T. - static T max(); + static const T max; /// Number of values for type \c T. - static std::size_t card(); + static const std::size_t card; }; -# ifndef MLN_INCLUDE_ONLY + + // bool template <> struct props<bool> { - static bool min() { return false; } - static bool max() { return true; } - static std::size_t card() { return 2; } + static const bool min = false; + static const bool max = true; + static const std::size_t card = 2; typedef binary_kind kind; }; + // integers template <> struct props<unsigned char> { - static unsigned char min() { return 0; } - static unsigned char max() { return 255; } - static std::size_t card() { return 256; } + static const unsigned char min = 0; + static const unsigned char max = 255; + static const std::size_t card = 256; typedef data_kind kind; }; template <> struct props<signed char> { - static signed char min() { return -128; } - static signed char max() { return 127; } - static std::size_t card() { return 256; } + static const signed char min = -128; + static const signed char max = 127; + static const std::size_t card = 256; typedef data_kind kind; }; template <> struct props<unsigned short> { - static unsigned short min() { return 0; } - static unsigned short max() { return 65535; } - static std::size_t card() { return 65536; } + static const unsigned short min = 0; + static const unsigned short max = 65535; + static const std::size_t card = 65536; typedef data_kind kind; }; template <> struct props<signed short> { - static signed short min() { return -32768; } - static signed short max() { return 32767; } - static std::size_t card() { return 655356; } + static const signed short min = -32768; + static const signed short max = 32767; + static const std::size_t card = 655356; typedef data_kind kind; }; template <> struct props<unsigned int> { - static unsigned int min() { return 0; } - static unsigned int max() { return UINT_MAX; } - static std::size_t card() { return std::size_t(UINT_MAX) + 1; } + static const unsigned int min = 0; + static const unsigned int max = UINT_MAX; typedef data_kind kind; + static const std::size_t card = 0; }; template <> struct props<signed int> { - static signed int min() { return INT_MIN; } - static signed int max() { return INT_MAX; } - static std::size_t card() { return std::size_t(UINT_MAX) + 1; } + static const signed int min = INT_MIN; + static const signed int max = INT_MAX; typedef data_kind kind; + static const std::size_t card = 0; }; template <> struct props<unsigned long int> { - static unsigned long int min() { return 0; } - static unsigned long int max() { return ULONG_MAX; } + static const unsigned long int min = 0; + static const unsigned long int max = ULONG_MAX; typedef data_kind kind; + static const std::size_t card = 0; }; template <> struct props<signed long int> { - static signed long int min() { return LONG_MIN; } - static signed long int max() { return LONG_MAX; } + static const signed long int min = LONG_MIN; + static const signed long int max = LONG_MAX; typedef data_kind kind; + static const std::size_t card = 0; }; + // floating template <> struct props<float> { - static float min() { return FLT_MIN; } - static float max() { return FLT_MAX; } + static const float min() { return FLT_MIN; } + static const float max() { return FLT_MAX; } typedef data_kind kind; + static const std::size_t card = 0; }; template <> struct props<double> { - static double min() { return DBL_MIN; } - static double max() { return DBL_MAX; } + static const double min() { return DBL_MIN; } + static const double max() { return DBL_MAX; } typedef data_kind kind; + static const std::size_t card = 0; }; -# endif // ! MLN_INCLUDE_ONLY - } // end of namespace mln::value } // end of namespace mln Index: mln/value/kind.hh --- mln/value/kind.hh (revision 1025) +++ mln/value/kind.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_VALUE_KIND_HH # define MLN_VALUE_KIND_HH -/*! \file mln/value/tags.hh +/*! \file mln/value/kind.hh * * \brief Define the set of kinds of value. */ Index: mln/value/set.hh --- mln/value/set.hh (revision 1025) +++ mln/value/set.hh (working copy) @@ -64,9 +64,6 @@ /// Backward Viter associated type. typedef bkd_viter_< set_<T> > bkd_viter; - /// Viter associated type. - typedef fwd_viter viter; - /// Test if \p v belongs to this set: always true! bool has(const T& v) const; @@ -129,14 +126,6 @@ } // end of namespace mln::value - - // FIXME: Turn definitions into objects! - typedef value::set_<bool> Booleans; - typedef value::set_<int> integers; - typedef value::set_<unsigned> naturals; - - - } // end of namespace mln Index: mln/value/int_u.hh --- mln/value/int_u.hh (revision 1025) +++ mln/value/int_u.hh (working copy) @@ -34,6 +34,7 @@ */ # include <mln/core/concept/value.hh> +# include <mln/metal/math.hh> # include <mln/value/internal/value_like.hh> # include <mln/value/props.hh> @@ -83,6 +84,8 @@ template <unsigned n> struct props< int_u_<n> > : public props< typename internal::encoding_<n>::ret > { + static const unsigned nbits = n; + static const std::size_t card = metal::pow<2, n>::value; };

1 0

1025: Clean-up and update pixel related material.
by Thierry Geraud 20 Jul '07

20 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Thierry Geraud <thierry.geraud(a)lrde.epita.fr> Clean-up and update pixel related material. * tests/README: New. * mln/core/trait: New. * mln/core/trait/pixter.hh: New. * mln/core/internal/pixel_impl.hh: New. * mln/core/internal/fixme.hh: New. * mln/core/concept/fast_image.hh: New. * mln/core/concept/doc/fast_iterator.hh: Rename as... * mln/core/concept/doc/pixel_iterator.hh: ...this new file. Update. * mln/core/concept/doc/fast_image.hh: New. * mln/core/concept/doc/genpixel.hh: New. * mln/core/tags.hh: New. * mln/metal: New. * mln/metal/none.hh: New. * mln/core/internal/lineary_pixel_iterator_base.hh: Remove. * mln/core/dpoints_piter.hh, * mln/core/box_piter.hh, * mln/core/internal/image_base.hh, * mln/core/internal/piter_adaptor.hh, * mln/core/concept/genpoint.hh, * mln/core/concept/psite.hh, * mln/core/concept/point.hh, * mln/core/concept/doc/genpoint.hh (pointer): Rename as... (pointer_): ...this. * tests/window2d.cc: Cosmetic change. * TODO: Update. * mln/core/macros.hh (mln_ima): Remove. (*ixter): Update. (mln_piter, mln_qiter, mln_niter): Update. * mln/core/psubset.hh (bkd_piter): Use fixme. (npoints): Update. * mln/core/box.hh (piter): Remove. * mln/core/rectangle2d.hh (bkd_piter): Use fixme. * mln/core/pixter2d_b.hh: Clean-up and update. * mln/core/window.hh (bkd_piter): Use fixme. * mln/core/dpoints_pixter.hh: Clean-up and update. * mln/core/internal/pixel_iterator_base.hh: Clean-up and update. * mln/core/internal/box_impl.hh: Touch doc. * mln/core/internal/image_base.hh (dpoint, coord): New. (pointer): Update. * mln/core/pixel.hh: Clean-up and update. * mln/core/hline2d.hh (bkd_piter): Use fixme. * mln/core/concept/image.hh (piter): Remove. (dpoint, coord): New. * mln/core/concept/pixel_iterator.hh: Update. * mln/core/concept/point_set.hh (piter): Remove. * mln/core/concept/object.hh: New include. * mln/core/concept/doc/image.hh: Update. * mln/core/concept/doc/point_set.hh: Update. * mln/core/concept/genpixel.hh: Clean-up and update. * mln/core/vline2d.hh (bkd_piter): Use fixme. * mln/core/image2d_b.hh: Clean-up and update. (buffer): New. (offset, point_at_offset): Fix. * mln/level/fast_median.hh: Update. * mln/level/median.hh: Prefer mln_coord to int. * mln/canvas/sbrowsing.hh: Update doc. * tests/fast_median.cc, * tests/pixter_point2d.cc, * tests/pixter_dpoint2d.cc, * tests/pixel.cc: Update. TODO | 55 +++++++-- mln/canvas/sbrowsing.hh | 18 ++- mln/core/box.hh | 4 mln/core/box_piter.hh | 8 - mln/core/concept/doc/fast_image.hh | 92 ++++++++++++++++ mln/core/concept/doc/genpixel.hh | 78 +++++++++++++ mln/core/concept/doc/genpoint.hh | 2 mln/core/concept/doc/image.hh | 17 +-- mln/core/concept/doc/pixel_iterator.hh | 24 ++-- mln/core/concept/doc/point_set.hh | 4 mln/core/concept/fast_image.hh | 96 +++++++++++++++++ mln/core/concept/genpixel.hh | 62 +---------- mln/core/concept/genpoint.hh | 7 - mln/core/concept/image.hh | 8 + mln/core/concept/object.hh | 1 mln/core/concept/pixel_iterator.hh | 32 ++--- mln/core/concept/point.hh | 4 mln/core/concept/point_set.hh | 4 mln/core/concept/psite.hh | 2 mln/core/dpoints_piter.hh | 8 - mln/core/dpoints_pixter.hh | 104 +++++++++++------- mln/core/hline2d.hh | 2 mln/core/image2d_b.hh | 128 +++++++++++++++++----- mln/core/internal/box_impl.hh | 2 mln/core/internal/fixme.hh | 55 +++++++++ mln/core/internal/image_base.hh | 10 + mln/core/internal/piter_adaptor.hh | 6 - mln/core/internal/pixel_impl.hh | 174 ++++++++++++++++++------------- mln/core/internal/pixel_iterator_base.hh | 135 ++++++------------------ mln/core/macros.hh | 50 +++++--- mln/core/pixel.hh | 104 +++++------------- mln/core/pixter2d_b.hh | 75 ++++++------- mln/core/psubset.hh | 7 - mln/core/rectangle2d.hh | 2 mln/core/tags.hh | 50 ++++++++ mln/core/trait/pixter.hh | 111 +++++++++++++++++++ mln/core/vline2d.hh | 2 mln/core/window.hh | 2 mln/level/fast_median.hh | 19 +-- mln/level/median.hh | 11 + mln/metal/none.hh | 57 ++++++++++ tests/README | 1 tests/fast_median.cc | 4 tests/pixel.cc | 7 - tests/pixter_dpoint2d.cc | 84 ++++++++++---- tests/pixter_point2d.cc | 28 ++-- tests/window2d.cc | 3 47 files changed, 1177 insertions(+), 582 deletions(-) Index: tests/window2d.cc --- tests/window2d.cc (revision 1024) +++ tests/window2d.cc (working copy) @@ -42,12 +42,15 @@ using namespace mln; window2d w; + mln_assertion(w.is_centered() = false); mln_assertion(w.is_symmetric() = true); w.insert(make::dpoint2d(-1,-1)); w.insert(make::dpoint2d( 1, 1)); + image2d_b<bool> ima = convert::to_image(w); debug::println(ima); + mln_assertion(w.delta() = 1); } Index: tests/README --- tests/README (revision 0) +++ tests/README (revision 0) @@ -0,0 +1 @@ +g++-4.1 -I.. -ansi -pedantic -W -Wall -Wextra -Wconversion sample.cc Index: tests/fast_median.cc --- tests/fast_median.cc (revision 1024) +++ tests/fast_median.cc (working copy) @@ -55,7 +55,7 @@ p.row() = p.col() = 1; { - mln_qixter(I) qix(win, p, input); + mln_qixter(I, W) qix(input, win, p); for_all(qix) std::cout << *qix << ' '; std::cout << " : " << qix.center_value() << std::endl; @@ -63,7 +63,7 @@ { pixel pix(input, p); - mln_qixter(I) qix(win, pix, input); + mln_qixter(I, W) qix(input, win, pix); for_all(qix) std::cout << *qix << ' '; std::cout << " : " << qix.center_value() << std::endl; Index: tests/pixter_point2d.cc --- tests/pixter_point2d.cc (revision 1024) +++ tests/pixter_point2d.cc (working copy) @@ -34,31 +34,29 @@ { using namespace mln; - image2d_b<int> ima(20, 20); + typedef image2d_b<int> I; + I ima(20, 20); - image2d_b<int>::piter p(ima.domain()); - image2d_b<int>::pixter f(ima); + mln_piter_(I) p(ima.domain()); + + mln_pixter_(I) pix(ima); int i = 0; for_all(p) - { ima(p) = i++; - } i = 0; - for_all(f) + for_all(pix) { - assert(*f = i ++); - *f = 5; + mln_assertion(*pix = i ++); + *pix = 5; } for_all(p) - { - assert(ima(p) = 5); - } + mln_assertion(ima(p) = 5); - f.start(); - assert(f.is_valid()); - f.invalidate(); - assert(!f.is_valid()); + pix.start(); + mln_assertion(pix.is_valid()); + pix.invalidate(); + mln_assertion(! pix.is_valid()); } Index: tests/pixter_dpoint2d.cc --- tests/pixter_dpoint2d.cc (revision 1024) +++ tests/pixter_dpoint2d.cc (working copy) @@ -1,37 +1,73 @@ +// Copyright (C) 2007 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. + +/*! \file tests/pixter_dpoint2d.cc + * + * \brief Test on mln::dpoints_fwd_pixter. + */ + #include <cassert> #include <iostream> + #include <mln/core/image2d_b.hh> #include <mln/core/window.hh> -#include <mln/level/fill.hh> #include <mln/core/dpoints_pixter.hh> -int -main() +#include <mln/level/fill.hh> + + +int main() { using namespace mln; - typedef image2d_b<int> I; + typedef image2d_b<int> I; + typedef I::dpoint D; + typedef window_<D> W; - I ima(20, 20); - level::fill(ima, 5); - I::piter f(ima.domain()); - window_<I::dpoint> set; - - set.insert(make::dpoint2d(0, -1)); - set.insert(make::dpoint2d(0, -1)); - set.insert(make::dpoint2d(1, 0)); - set.insert(make::dpoint2d(1, 0)); + typedef dpoints_fwd_pixter qixter; - for_all(f) - { - if (f[0] > 0 && f[1] > 0 && f[0] < 19 && f[1] < 19) - { - I::qixter dit(set, f, ima); + const unsigned size = 20; + I ima(size, size); - for_all(dit) - { - assert(*dit = 5); - } - } - } + const int value = 51; + level::fill(ima, value); + + W win; + win + .insert(make::dpoint2d(0, -1)) + .insert(make::dpoint2d(0, -1)) + .insert(make::dpoint2d(1, 0)) + .insert(make::dpoint2d(1, 0)); + + mln_piter_(I) p(ima.domain()); + qixter qix(ima, win, p); + + for_all(p) + if (p[0] > 0 && p[1] > 0 && p[0] < size - 1 && p[1] < size - 1) + for_all(qix) + mln_assertion(*qix = value); } Index: tests/pixel.cc --- tests/pixel.cc (revision 1024) +++ tests/pixel.cc (working copy) @@ -42,15 +42,13 @@ I ima(3, 3); { - pixel pxl(ima); - pxl.site() = make::point2d(1, 1); + pixel pxl(ima, make::point2d(1, 1)); *pxl = 51; mln_assertion(ima.at(1,1) = 51); } { - pixel<const I> pxl(ima); - pxl.site() = make::point2d(1, 1); + pixel<const I> pxl(ima, make::point2d(1, 1)); ima.at(1,1) = 51; mln_assertion(*pxl = 51); @@ -58,4 +56,5 @@ // *pxl = 0; // assignment of read-only location } + } Index: TODO --- TODO (revision 1024) +++ TODO (working copy) @@ -1,25 +1,54 @@ -*- outline -*- -* compilation +* large changes -g++-4.1 -I.. int_u8.cc -Wall -Wextra -Wconversion +introduce grids -* value +* adds-on -Value_Set -props { ... } -notion of lowq -+ values() : vset dans Image ? +** abstractions +Function_i2i +Weighted_Window -* function +** types -Function_p2v -Function_p2b -Function_v2w +value::proxy to dispatch read/write + the corresponding image type +a mean_value object { sum; count } and operator+ +value_cast<T>(image) +t_image to "transpose" the 0 and the i-th coordinates +image defined by f/pset -F apply(I, F) -O transform(I, F) +** extensions +in Fast_Image: memcpy and memset-like methods (?) +.offset() into GenPixel +notion of low-quantization +vset and .values() in Image +built-in op objects -> reverse lhs/rhs + + +* renaming + +GenPixel, GenPoint, Psite, Piter: longer name (?) +mlc into metal +point-wise material: pw::value_of, pw::constant, and pw::variable +kind-related: mln_value_kind(I) and mln_kind(V) + + +* clean-up + +accu::histo and median: remove inheritance + + +* processing routines + +reconstructions + their canvases +fast versions of level::fill and level::paste +histogram class +sorting points w.r.t. their value (array of offsets, psites, points) +border::* and no-op versions if not fast image +arith::inplace_plus et al. +linear:: for convolutions Index: mln/core/dpoints_piter.hh --- mln/core/dpoints_piter.hh (revision 1024) +++ mln/core/dpoints_piter.hh (working copy) @@ -78,7 +78,7 @@ operator mln_point(D) () const; /// Address of the point this iterator designates. - const point* pointer() const; + const point* pointer_() const; /// Test the iterator validity. bool is_valid() const; @@ -108,7 +108,7 @@ unsigned i_; point p_; // location of this iterator; p_ makes this iterator be - // itself a potential center point (Cf. the pointer() method). + // itself a potential center point (Cf. the pointer_() method). }; @@ -119,7 +119,7 @@ dpoints_fwd_piter<D>::dpoints_fwd_piter(const Dps& dps, const GenPoint<Pref>& p_ref) : dps_(exact(dps).vec()), - p_ref_(* internal::force_exact<Pref>(p_ref).pointer()) + p_ref_(* internal::force_exact<Pref>(p_ref).pointer_()) { invalidate(); } @@ -133,7 +133,7 @@ template <typename D> const mln_point(D)* - dpoints_fwd_piter<D>::pointer() const + dpoints_fwd_piter<D>::pointer_() const { return & p_; } Index: mln/core/macros.hh --- mln/core/macros.hh (revision 1024) +++ mln/core/macros.hh (working copy) @@ -81,29 +81,20 @@ /// Shortcut to access the fwd_piter type associated to T. # define mln_fwd_piter(T) typename T::fwd_piter -/// Shortcut to access the fwd_pixter type associated to T. -# define mln_fwd_pixter(T) typename T::fwd_pixter - /// Shortcut to access the bkd_piter type associated to T. # define mln_bkd_piter(T) typename T::bkd_piter /// Shortcut to access the fwd_qiter type associated to T. # define mln_fwd_qiter(T) typename T::fwd_qiter -/// Shortcut to access the fwd_qixter type associated to T. -# define mln_fwd_qixter(T) typename T::fwd_qixter - /// Shortcut to access the fwd_viter type associated to T. # define mln_fwd_viter(T) typename T::fwd_viter // i -/// Shortcut to access the fwd_viter type associated to T. +/// Shortcut to access the image type associated to T. # define mln_image(T) typename T::image -/// Shortcut to access the fwd_viter type associated to T. -# define mln_ima(T) typename T::ima - // l @@ -115,10 +106,8 @@ // p /// Shortcut to access the piter type associated to T. -# define mln_piter(T) typename T::piter - -/// Shortcut to access the pixter type associated to T. -# define mln_pixter(T) typename T::pixter +# define mln_piter(T) typename T::fwd_piter +# define mln_piter_(T) T::fwd_piter /// Shortcut to access the pset type associated to T. # define mln_pset(T) typename T::pset @@ -133,16 +122,13 @@ // q /// Shortcut to access the qiter type associated to T. -# define mln_qiter(T) typename T::qiter - -// Shortcut to access to the qixter type associated to T -# define mln_qixter(T) typename T::qixter +# define mln_qiter(T) typename T::fwd_qiter // n /// Shortcut to access the niter type associated to T. -# define mln_niter(T) typename T::niter +# define mln_niter(T) typename T::fwd_niter // r @@ -160,7 +146,9 @@ # define mln_value(T) typename T::value /// Shortcut to access the viter type associated to T. -# define mln_viter(T) typename T::viter +# define mln_viter(T) typename T::fwd_viter + + @@ -172,4 +160,26 @@ # define mln_ch_value(I, T) typename I::template change_value<T>::ret + +// FIXME: Doc! + +# define mln_fwd_pixter(I) typename mln::trait::fwd_pixter::ret +# define mln_fwd_pixter_(I) mln::trait::fwd_pixter::ret + +# define mln_bkd_pixter(I) typename mln::trait::bkd_pixter::ret +# define mln_pixter(I) mln_fwd_pixter(I) +# define mln_pixter_(I) mln_fwd_pixter_(I) + + +# define mln_fwd_qixter(I, W) typename mln::trait::fwd_qixter<I, W>::ret +# define mln_bkd_qixter(I, W) typename mln::trait::bkd_qixter<I, W>::ret +# define mln_qixter(I, W) mln_fwd_qixter(I, W) + + +# define mln_fwd_nixter(I, N) typename mln::trait::fwd_nixter<I, N>::ret +# define mln_bkd_nixter(I, N) typename mln::trait::bkd_nixter<I, N>::ret +# define mln_nixter(I, N) mln_fwd_nixter(I, N) + + + #endif // ! MLN_CORE_MACROS_HH Index: mln/core/trait/pixter.hh --- mln/core/trait/pixter.hh (revision 0) +++ mln/core/trait/pixter.hh (revision 0) @@ -0,0 +1,111 @@ +// Copyright (C) 2007 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_CORE_TRAIT_PIXTER_HH +# define MLN_CORE_TRAIT_PIXTER_HH + +/*! \file mln/core/trait/pixter.hh + * + * \brief Definition of some traits over image types. + */ + +# include <mln/metal/none.hh> + + +namespace mln +{ + + namespace trait + { + + + // pixter + + template <typename I> + struct fwd_pixter + { + typedef metal::none ret; + }; + + template <typename I> + struct bkd_pixter + { + typedef metal::none ret; + }; + + template <typename I> + struct pixter : fwd_pixter + { + }; + + + // qixter + + template <typename I, typename W> + struct fwd_qixter + { + typedef metal::none ret; + }; + + template <typename I, typename W> + struct bkd_qixter + { + typedef metal::none ret; + }; + + template <typename I, typename W> + struct qixter : fwd_qixter<I, W> + { + }; + + + // nixter + + template <typename I, typename N> + struct fwd_nixter + { + typedef metal::none ret; + }; + + template <typename I, typename N> + struct bkd_nixter + { + typedef metal::none ret; + }; + + template <typename I, typename N> + struct nixter : fwd_nixter<I, N> + { + }; + + + } // end of namespace mln::trait + +} // end of namespace mln + + +#endif // ! MLN_CORE_TRAIT_PIXTER_HH Index: mln/core/psubset.hh --- mln/core/psubset.hh (revision 1024) +++ mln/core/psubset.hh (working copy) @@ -78,10 +78,7 @@ typedef psubset_fwd_piter_<S,F> fwd_piter; /// Backward Piter associated type. - typedef psubset_fwd_piter_<S,F> bkd_piter; // FIXME: bkd!!! - - /// Piter associated type. - typedef fwd_piter piter; + typedef internal::fixme bkd_piter; /// Constructor without argument. @@ -180,7 +177,7 @@ psubset<S,F>::npoints() const { std::size_t n = 0; - piter p(*this); + fwd_piter p(*this); for_all(p) ++n; return n; Index: mln/core/box.hh --- mln/core/box.hh (revision 1024) +++ mln/core/box.hh (working copy) @@ -70,10 +70,6 @@ */ typedef box_bkd_piter_ bkd_piter; - /*! \brief Piter associated type. - */ - typedef fwd_piter piter; - /*! \brief Minimum point. */ P pmin() const; Index: mln/core/rectangle2d.hh --- mln/core/rectangle2d.hh (revision 1024) +++ mln/core/rectangle2d.hh (working copy) @@ -70,7 +70,7 @@ /*! \brief Piter type to browse a rectangle such as: "for each row * (decreasing), for each column (decreasing)." */ - typedef dpoints_fwd_piter<dpoint2d> bkd_qiter; // FIXME: bkd!!! + typedef internal::fixme bkd_qiter; /*! \brief Same as fwd_qiter. */ Index: mln/core/box_piter.hh --- mln/core/box_piter.hh (revision 1024) +++ mln/core/box_piter.hh (working copy) @@ -76,7 +76,7 @@ operator P() const; /// Address of the point. - const P* pointer() const; + const P* pointer_() const; /// Give the i-th coordinate. coord operator[](unsigned i) const; @@ -136,7 +136,7 @@ operator P() const; /// Address of the point. - const P* pointer() const; + const P* pointer_() const; /// Give the i-th coordinate. coord operator[](unsigned i) const; @@ -183,7 +183,7 @@ template <typename P> const P* - box_fwd_piter_::pointer() const + box_fwd_piter_::pointer_() const { return & p_; } @@ -254,7 +254,7 @@ template <typename P> const P* - box_bkd_piter_::pointer() const + box_bkd_piter_::pointer_() const { return & p_; } Index: mln/core/pixter2d_b.hh --- mln/core/pixter2d_b.hh (revision 1024) +++ mln/core/pixter2d_b.hh (working copy) @@ -28,80 +28,77 @@ #ifndef MLN_CORE_PIXTER2D_B_HH # define MLN_CORE_PIXTER2D_B_HH -# include <mln/core/internal/lineary_pixel_iterator_base.hh> -# include <mln/core/point2d.hh> -# include <iostream> - - -//FIXME comment /*! \file mln/core/pixter2d_b.hh * * \brief Pixel iterator class on a image 2d with border. */ +# include <mln/core/internal/pixel_iterator_base.hh> +# include <mln/core/point2d.hh> + namespace mln { - // Forward declaration - template <typename T> class image2d_b; - - template <typename T> - class fwd_pixter2d_b : - public internal::lineary_pixel_iterator_base_<image2d_b<T>, fwd_pixter2d_b<T> > + template <typename I> + class fwd_pixter2d_b : public internal::pixel_iterator_base_< I, fwd_pixter2d_b > { - typedef internal::lineary_pixel_iterator_base_<image2d_b<T>, fwd_pixter2d_b<T> > super; + typedef internal::pixel_iterator_base_< I, fwd_pixter2d_b > super_; + typedef mln_value(super_) value_; public: - /// Image pixel value type. - typedef mln_value(image2d_b<T>) value; /*! \brief Constructor. * - * \param[in] ima Image to iterate. + * \param[in] image Image to iterate over its pixels. */ - fwd_pixter2d_b(image2d_b<T>& image); - /// Move the iterator on the next element. + fwd_pixter2d_b(I& image); + + /// Go to the next pixel. void next_(); private: - /// Size of the image border. - unsigned border_size_; + + /// Twice the size of the image border. + unsigned border_x2_; + /// Row offset. unsigned row_offset_; - /// End of a row. - value *eor_; + + /// End of the current row. + value_* eor_; }; + + // FIXME: bkd_pixter2d_b + + #ifndef MLN_INCLUDE_ONLY - template <typename T> - fwd_pixter2d_b<T>::fwd_pixter2d_b(image2d_b<T>& image) : - super(image), - border_size_(image.border()), - row_offset_((image.domain().pmax()[1] - image.domain().pmin()[1]) - + 2 * border_size_ + 1) + template <typename I> + fwd_pixter2d_b::fwd_pixter2d_b(I& image) : + super_(image) { - this->start_ = &image(image.domain().pmin()); - this->eor_ = &image(make::point2d(image.domain().pmin()[0], image.domain().pmax()[1])) + 1; - this->eoi_ = &image(image.domain().pmax()) + 1; + mln_precondition(image.has_data()); + border_x2_ = 2 * image.border(); + row_offset_ = image.max_col() - image.min_col() + 1 + border_x2_; + eor_ = & image.at(image.min_row(), image.max_col()) + 1; } - template <typename T> - void fwd_pixter2d_b<T>::next_() + template <typename I> + void + fwd_pixter2d_b::next_() { - ++(this->value_ptr_); - - if (this->value_ptr_ = this->eor_ && this->value_ptr_ != this->eoi_) + ++this->value_ptr_; + if (this->value_ptr_ = eor_ && this->value_ptr_ != this->eoi_) { - this->value_ptr_ += 2 * this->border_size_; - this->eor_ += this->row_offset_; + this->value_ptr_ += border_x2_; + eor_ += row_offset_; } } #endif // ! MLN_INCLUDE_ONLY - } // end of namespace mln #endif // ! MLN_CORE_PIXTER2D_B_HH Index: mln/core/window.hh --- mln/core/window.hh (revision 1024) +++ mln/core/window.hh (working copy) @@ -75,7 +75,7 @@ /*! \brief Piter type to browse the points of a generic window * w.r.t. the reverse ordering of delta-points. */ - typedef dpoints_bkd_piter<D> bkd_qiter; + typedef internal::fixme bkd_qiter; /*! \brief Same as fwd_qiter. */ Index: mln/core/dpoints_pixter.hh --- mln/core/dpoints_pixter.hh (revision 1024) +++ mln/core/dpoints_pixter.hh (working copy) @@ -28,47 +28,60 @@ #ifndef MLN_CORE_DPOINTS_PIXTER_HH # define MLN_CORE_DPOINTS_PIXTER_HH +/*! \file mln/core/dpoints_pixter.hh + * + * \brief Definition of mln::dpoints_fwd_pixter and + * mln::dpoints_bkd_pixter. + */ + # include <cassert> # include <vector> -# include <mln/core/internal/pixel_iterator_base.hh> + +# include <mln/core/concept/pixel_iterator.hh> +# include <mln/core/concept/genpixel.hh> +# include <mln/core/internal/pixel_impl.hh> + namespace mln { - /*! \brief A generic forward iterator on image pixel of a windows and of + /*! \brief A generic forward iterator on pixels of windows and of * neighborhoods. * * The parameter \c I is the image type. */ template <typename I> - class dpoints_pixter : public internal::pixel_iterator_base_< I, dpoints_pixter > + class dpoints_fwd_pixter : public Pixel_Iterator< dpoints_fwd_pixter >, + public internal::pixel_impl_< I, dpoints_fwd_pixter > { - typedef typename internal::pixel_iterator_base_< I, dpoints_pixter > super; + typedef typename internal::pixel_impl_< I, dpoints_fwd_pixter > super; public: - /// Image pixel value - typedef mln_value(I) value; - /// Image pixel rvalue - typedef mln_rvalue(I) rvalue; - /// Image pixel lvalue - typedef mln_lvalue(I) lvalue; + /// Using super value type. + typedef mln_value(super) value; /*! \brief Constructor. * - * \param[in] dps Object that can provide an array of delta-points. - * \param[in] p_ref Center point to iterate around. - * \param[in] ima Image to iterate. + * \param[in] image Image subject to iteration. + * \param[in] dps Object that can provide a set of delta-points. + * \param[in] p_ref Center (generalized) point to iterate around. */ template <typename Dps, typename Pref> - dpoints_pixter(const Dps& dps, - const GenPoint<Pref>& p_ref, - I& image); + dpoints_fwd_pixter(I& image, + const Dps& dps, + const GenPoint<Pref>& p_ref); + /*! \brief Constructor. + * + * \param[in] image Image subject to iteration. + * \param[in] dps Object that can provide a set of delta-points. + * \param[in] p_ref Center (generalized) pixel to iterate around. + */ template <typename Dps, typename Pref> - dpoints_pixter(const Dps& dps, - const GenPixel<Pref>& p_ref, - I& image); + dpoints_fwd_pixter(I& image, + const Dps& dps, + const GenPixel<Pref>& p_ref); /// Start an iteration. void start(); @@ -87,7 +100,7 @@ void update(); /// The value around which this iterator moves. - const value& center_value() const; + const mln_value(I)& center_value() const; private: @@ -95,7 +108,7 @@ std::vector<int> offset_; /// current offset - int i_; + unsigned i_; /// reference pixel / point in the image value** value_ref_; @@ -108,16 +121,24 @@ }; + // FIXME: dpoints_bkd_pixter + + + #ifndef MLN_INCLUDE_ONLY + + // dpoints_fwd_pixter + template <typename I> template <typename Dps, typename Pref> - dpoints_pixter::dpoints_pixter(const Dps& dps, - const GenPoint<Pref>& p_ref, - I& image) + dpoints_fwd_pixter::dpoints_fwd_pixter(I& image, + const Dps& dps, + const GenPoint<Pref>& p_ref) : super(image) { - p_ref_ = internal::force_exact<Pref>(p_ref).pointer(); + mln_precondition(image.has_data()); + p_ref_ = internal::force_exact<Pref>(p_ref).pointer_(); value_ref_ = 0; init_(dps); } @@ -125,19 +146,20 @@ template <typename I> template <typename Dps, typename Pref> - dpoints_pixter::dpoints_pixter(const Dps& dps, - const GenPixel<Pref>& p_ref, - I& image) + dpoints_fwd_pixter::dpoints_fwd_pixter(I& image, + const Dps& dps, + const GenPixel<Pref>& p_ref) : super(image) { + mln_precondition(image.has_data()); p_ref_ = 0; - value_ref_ = internal::force_exact<Pref>(p_ref).address(); + value_ref_ = internal::force_exact<Pref>(p_ref).address_(); init_(dps); } template <typename I> const mln_value(I)& - dpoints_pixter::center_value() const + dpoints_fwd_pixter::center_value() const { mln_invariant(value_ref_ != 0 || p_ref_ != 0); if (p_ref_) @@ -149,21 +171,21 @@ template <typename I> template <typename Dps> void - dpoints_pixter::init_(const Dps& dps) + dpoints_fwd_pixter::init_(const Dps& dps) { for (unsigned i = 0; i < dps.nelements(); ++i) offset_.push_back(this->image_.offset(dps.element(i))); - + // offset_[0] is absolute + // other offsets are relative: if (dps.nelements() > 1) for (unsigned i = dps.nelements() - 1; i > 0; --i) offset_[i] -= offset_[i - 1]; - - // offset_[0] is absolute; other offsets are relative. invalidate(); } template <typename I> - void dpoints_pixter::update() + void + dpoints_fwd_pixter::update() { if (is_valid()) { @@ -175,27 +197,31 @@ } template <typename I> - void dpoints_pixter::start() + void + dpoints_fwd_pixter::start() { i_ = 0; update(); } template <typename I> - void dpoints_pixter::next_() + void + dpoints_fwd_pixter::next_() { ++i_; this->value_ptr_ += offset_[i_]; } template <typename I> - bool dpoints_pixter::is_valid() const + bool + dpoints_fwd_pixter::is_valid() const { return i_ != offset_.size(); } template <typename I> - void dpoints_pixter::invalidate() + void + dpoints_fwd_pixter::invalidate() { i_ = offset_.size(); } Index: mln/core/internal/pixel_impl.hh --- mln/core/internal/pixel_impl.hh (revision 1024) +++ mln/core/internal/pixel_impl.hh (working copy) @@ -25,15 +25,16 @@ // reasons why the executable file might be covered by the GNU General // Public License. -#ifndef MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH -# define MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH +#ifndef MLN_CORE_INTERNAL_PIXEL_IMPL_HH +# define MLN_CORE_INTERNAL_PIXEL_IMPL_HH -/*! \file mln/core/internal/pixel_iterator_base.hh +/*! \file mln/core/internal/pixel_impl.hh * - * \brief Base class for Pixel_Iterator concept implementation classes. + * \brief Define a couple of implementation classes to provide methods + * to classes of generalized pixels. */ -# include <mln/core/concept/pixel_iterator.hh> +# include <mln/core/internal/force_exact.hh> namespace mln @@ -42,11 +43,13 @@ namespace internal { - /*! \brief pixel_iterator_base_ class + /*! \brief Implementation class to equip generalized pixel + * classes based on mutable images. * + * \internal */ template <typename I, typename E> - class pixel_iterator_base_ : public Pixel_Iterator<E> + class pixel_impl_ { public: @@ -68,111 +71,138 @@ /// Address of the current iterator value/pixel. - value** address() const; + value** address_() const; + protected: - // FIXME: Inactivated: + /// Image associated to the iterator + I& image_; -// /// I type. -// typedef I ima; + /// Current pixel / value + value* value_ptr_; -// /// Image psite type. -// typedef mln_psite(I) psite; + /// Constructor + pixel_impl_(I& image); -// /// Get the image associated to the current pixel iterator. -// const ima& image() const; + private: + bool is_valid_() const; + }; -// /// psite associated to the iterator. -// const psite& site() const; -// /// psite associated to the iterator. -// psite& site(); + /*! \brief Implementation class to equip generalized pixel + * classes based on constant images. + * + * \internal + */ + template <typename I, typename E> + class pixel_impl_< const I, E > + { + public: -// /// Address of the current iterator value/pixel. -// value* address(); + /// Image value type. + typedef mln_value(I) value; + + /// Image rvalue type. + typedef mln_rvalue(I) rvalue; + + + /// Get the pixel iterator value. + rvalue operator*() const; + + + /// Address of the current iterator value/pixel. + value** address_() const; protected: /// Image associated to the iterator - I& image_; + const I& image_; /// Current pixel / value - value* value_ptr_; - - // FIXME: Inactivated: + const value* value_ptr_; -// /// Psite of the pixel -// psite p_; + /// Constructor + pixel_impl_(const I& image); - pixel_iterator_base_(I& image); + private: + bool is_valid_() const; }; + #ifndef MLN_INCLUDE_ONLY + // pixel_impl_<I, E> + template <typename I, typename E> - pixel_iterator_base_<I, E>::pixel_iterator_base_(I& image) : + bool + pixel_impl_<I, E>::is_valid_() const + { + return value_ptr_ != 0 && internal::force_exact<E>(*this).is_valid(); + } + + template <typename I, typename E> + pixel_impl_<I, E>::pixel_impl_(I& image) : image_(image), value_ptr_(0) { } -// template <typename I, typename E> -// const typename pixel_iterator_base_<I, E>::ima& -// pixel_iterator_base_<I, E>::image() const -// { -// return ima_; -// } - template <typename I, typename E> mln_lvalue(I) - pixel_iterator_base_<I, E>::operator* () + pixel_impl_<I, E>::operator*() + { + mln_precondition(is_valid_()); + return *value_ptr_; + } + + template <typename I, typename E> + mln_rvalue(I) + pixel_impl_<I, E>::operator*() const { - mln_precondition(exact(this)->is_valid()); - mln_precondition(value_ptr_ != 0); + mln_precondition(is_valid_()); return *value_ptr_; } template <typename I, typename E> + mln_value(I) ** + pixel_impl_<I, E>::address_() const + { + mln_precondition(is_valid_()); + return (value**)(& this->value_ptr_); + } + + + // pixel_impl_<const I, E> + + template <typename I, typename E> + bool + pixel_impl_<const I, E>::is_valid_() const + { + return value_ptr_ != 0 && internal::force_exact<E>(*this).is_valid(); + } + + template <typename I, typename E> + pixel_impl_<const I, E>::pixel_impl_(const I& image) : + image_(image), + value_ptr_(0) + { + } + + template <typename I, typename E> mln_rvalue(I) - pixel_iterator_base_<I, E>::operator* () const + pixel_impl_<const I, E>::operator*() const { - mln_precondition(exact(this)->is_valid()); - mln_precondition(value_ptr_ != 0); + mln_precondition(is_valid_()); return *value_ptr_; } template <typename I, typename E> mln_value(I) ** - pixel_iterator_base_<I, E>::address() const + pixel_impl_<const I, E>::address_() const { - mln_precondition(exact(this)->is_valid()); - mln_precondition(value_ptr_ != 0); - return & (mln_value(I)*)(value_ptr_); - } - -// template <typename I, typename E> -// mln_value(I)* -// pixel_iterator_base_<I, E>::address() -// { -// mln_precondition(exact(this)->is_valid()); -// return value_ptr_; -// } - -// template <typename I, typename E> -// const typename pixel_iterator_base_<I, E>::psite& -// pixel_iterator_base_<I, E>::site() const -// { -// //FIXME: update psite -// return p_; -// } - -// template <typename I, typename E> -// typename pixel_iterator_base_<I, E>::psite& -// pixel_iterator_base_<I, E>::site() -// { -// //FIXME: update psite -// return p_; -// } + mln_precondition(is_valid_()); + return (value**)(& this->value_ptr_); + } #endif // ! MLN_INCLUDE_ONLY @@ -181,4 +211,4 @@ } // end of namespace mln -#endif // ! MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH +#endif // ! MLN_CORE_INTERNAL_PIXEL_IMPL_HH Index: mln/core/internal/pixel_iterator_base.hh --- mln/core/internal/pixel_iterator_base.hh (revision 1024) +++ mln/core/internal/pixel_iterator_base.hh (working copy) @@ -30,10 +30,11 @@ /*! \file mln/core/internal/pixel_iterator_base.hh * - * \brief Base class for Pixel_Iterator concept implementation classes. + * \brief Base class to factor code for pixel iterator classes. */ # include <mln/core/concept/pixel_iterator.hh> +# include <mln/core/internal/pixel_impl.hh> namespace mln @@ -42,138 +43,74 @@ namespace internal { - /*! \brief pixel_iterator_base_ class + /*! \brief A base class for pixel iterators. * */ template <typename I, typename E> - class pixel_iterator_base_ : public Pixel_Iterator<E> + class pixel_iterator_base_ : public Pixel_Iterator<E>, + public internal::pixel_impl_<I, E> { - public: - - /// Image value type. - typedef mln_value(I) value; - - /// Image lvalue type. - typedef mln_lvalue(I) lvalue; - - /// Image rvalue type. - typedef mln_rvalue(I) rvalue; - - - /// pixel iterator value. - lvalue operator* (); - - /// Get the pixel iterator value. - rvalue operator* () const; - - - /// Address of the current iterator value/pixel. - value** address() const; - + typedef internal::pixel_impl_<I, E> super_; + typedef mln_value(super_) value_; - // FIXME: Inactivated: - -// /// I type. -// typedef I ima; - -// /// Image psite type. -// typedef mln_psite(I) psite; - -// /// Get the image associated to the current pixel iterator. -// const ima& image() const; + public: -// /// psite associated to the iterator. -// const psite& site() const; + /// Start an iteration. + void start(); -// /// psite associated to the iterator. -// psite& site(); + /// Invalidate the iterator. + void invalidate(); -// /// Address of the current iterator value/pixel. -// value* address(); + /// Test if the iterator is valid. + bool is_valid() const; protected: - /// Image associated to the iterator - I& image_; - - /// Current pixel / value - value* value_ptr_; + /// Beginning of the image. + value_* boi_; - // FIXME: Inactivated: - -// /// Psite of the pixel -// psite p_; + /// End of the image (past-the-end). + value_* eoi_; + /// Constructor. pixel_iterator_base_(I& image); }; + #ifndef MLN_INCLUDE_ONLY template <typename I, typename E> - pixel_iterator_base_<I, E>::pixel_iterator_base_(I& image) : - image_(image), - value_ptr_(0) + pixel_iterator_base_<I, E>::pixel_iterator_base_(I& image) + : super_(image) { + mln_precondition(image.has_data()); + I& ima = this->image_; + boi_ = & ima( ima.domain().pmin() ); + eoi_ = & ima( ima.domain().pmax() ) + 1; + invalidate(); } -// template <typename I, typename E> -// const typename pixel_iterator_base_<I, E>::ima& -// pixel_iterator_base_<I, E>::image() const -// { -// return ima_; -// } - template <typename I, typename E> - mln_lvalue(I) - pixel_iterator_base_<I, E>::operator* () + void + pixel_iterator_base_<I, E>::start() { - mln_precondition(exact(this)->is_valid()); - mln_precondition(value_ptr_ != 0); - return *value_ptr_; + this->value_ptr_ = boi_; } template <typename I, typename E> - mln_rvalue(I) - pixel_iterator_base_<I, E>::operator* () const + void + pixel_iterator_base_<I, E>::invalidate() { - mln_precondition(exact(this)->is_valid()); - mln_precondition(value_ptr_ != 0); - return *value_ptr_; + this->value_ptr_ = eoi_; } template <typename I, typename E> - mln_value(I) ** - pixel_iterator_base_<I, E>::address() const + bool + pixel_iterator_base_<I, E>::is_valid() const { - mln_precondition(exact(this)->is_valid()); - mln_precondition(value_ptr_ != 0); - return & (mln_value(I)*)(value_ptr_); + return this->value_ptr_ != eoi_; } -// template <typename I, typename E> -// mln_value(I)* -// pixel_iterator_base_<I, E>::address() -// { -// mln_precondition(exact(this)->is_valid()); -// return value_ptr_; -// } - -// template <typename I, typename E> -// const typename pixel_iterator_base_<I, E>::psite& -// pixel_iterator_base_<I, E>::site() const -// { -// //FIXME: update psite -// return p_; -// } - -// template <typename I, typename E> -// typename pixel_iterator_base_<I, E>::psite& -// pixel_iterator_base_<I, E>::site() -// { -// //FIXME: update psite -// return p_; -// } - #endif // ! MLN_INCLUDE_ONLY } // end of namespace internal Index: mln/core/internal/fixme.hh --- mln/core/internal/fixme.hh (revision 0) +++ mln/core/internal/fixme.hh (revision 0) @@ -0,0 +1,55 @@ +// Copyright (C) 2007 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_CORE_INTERNAL_FIXME_HH +# define MLN_CORE_INTERNAL_FIXME_HH + +/*! \file mln/core/internal/fixme.hh + * + * \brief Definition of mln::fixme for internal use only. + */ + + +namespace mln +{ + + namespace internal + { + + /*! \brief A FIXME class to make explicit in code that a type is + * not yet implemented. + * + * \internal + */ + struct fixme; + + } // end of namespace mln::internal + +} // end of namespace mln + + +#endif // ! MLN_CORE_INTERNAL_FIXME_HH Index: mln/core/internal/box_impl.hh --- mln/core/internal/box_impl.hh (revision 1024) +++ mln/core/internal/box_impl.hh (working copy) @@ -48,7 +48,7 @@ // box_impl /*! \brief Implementation class to equip objects having a bounding - * box with methods. + * box. * * \internal */ Index: mln/core/internal/image_base.hh --- mln/core/internal/image_base.hh (revision 1024) +++ mln/core/internal/image_base.hh (working copy) @@ -58,14 +58,20 @@ /// Point associated type. typedef mln_point(S) point; + + /// Dpoint associated type. + typedef mln_dpoint(point) dpoint; + + /// Coordinate associated type. + typedef mln_coord(point) coord; + + /// Forward Piter associated type. typedef mln_fwd_piter(S) fwd_piter; /// Backward Piter associated type. typedef mln_bkd_piter(S) bkd_piter; - /// Piter associated type. - typedef fwd_piter piter; /// Test if \p p belongs to the image domain. bool has(const psite& p) const; Index: mln/core/internal/piter_adaptor.hh --- mln/core/internal/piter_adaptor.hh (revision 1024) +++ mln/core/internal/piter_adaptor.hh (working copy) @@ -77,7 +77,7 @@ operator mln_point(Pi) () const; /// Address of the point. - const point* pointer() const; + const point* pointer_() const; /// Give the i-th coordinate. coord operator[](unsigned i) const; @@ -120,9 +120,9 @@ template <typename Pi, typename E> const mln_point(Pi)* - piter_adaptor_<Pi,E>::pointer() const + piter_adaptor_<Pi,E>::pointer_() const { - return piter_.pointer(); + return piter_.pointer_(); } template <typename Pi, typename E> Index: mln/core/pixel.hh --- mln/core/pixel.hh (revision 1024) +++ mln/core/pixel.hh (working copy) @@ -34,57 +34,36 @@ */ # include <mln/core/concept/genpixel.hh> +# include <mln/core/internal/pixel_impl.hh> namespace mln { - template <typename I> - struct pixel_lvalue - { - typedef mln_lvalue(I) ret; - }; - - template <typename I> - struct pixel_lvalue< const I > - { - typedef mln_rvalue(I) ret; - }; - - - /*! \brief Generic pixel class. * * The parameter is \c I the type of the image it belongs to. */ template <typename I> struct pixel : public Object< pixel >, - public GenPixel< pixel > + public GenPixel< pixel >, + internal::pixel_impl_< I, pixel > { - typedef mln_psite(I) psite; - typedef mln_value(I) value; - typedef mln_rvalue(I) rvalue; - typedef I ima; - - pixel(I& ima); - pixel(I& ima, const psite& p); - - const I& image() const; + /// Constructor. + pixel(I& image); - const psite& site() const; - psite& site(); + /// Constructor. + pixel(I& image, const mln_point(I)& p); - mln_rvalue(I) operator*() const; - mln_lvalue(I) operator*(); + /// Change the pixel to the one at point \p p. + void change_to(const mln_point(I)& p); - value** address() const; + /// Test if this pixel is valid. + bool is_valid() const; - protected: - - I& ima_; - psite p_; - value* value_ptr_; + private: + typedef internal::pixel_impl_<I, pixel > super; }; @@ -92,58 +71,37 @@ template <typename I> pixel::pixel(I& image) - : ima_(image) - { - } - - template <typename I> - pixel::pixel(I& image, const psite& p) - : ima_(image), - p_(p) - { - value_ptr_ = & ima_(p); - } - - template <typename I> - const I& - pixel::image() const - { - return ima_; - } - - template <typename I> - const mln_psite(I)& - pixel::site() const - { - return p_; - } - - template <typename I> - mln_psite(I)& - pixel::site() + : super(image) { - return p_; + this->value_ptr_ = 0; } template <typename I> - mln_rvalue(I) - pixel::operator*() const + pixel::pixel(I& image, const mln_point(I)& p) + : super(image) { - return ima_(p_); + mln_precondition(this->image_.owns_(p)); + change_to(p); } template <typename I> - mln_lvalue(I) - pixel::operator*() + void + pixel::change_to(const mln_point(I)& p) { - return ima_(p_); + mln_precondition(this->image_.owns_(p)); + this->value_ptr_ = & this->image_(p); } template <typename I> - mln_value(I)** - pixel::address() const + bool + pixel::is_valid() const { - return (mln_value(I)**)(& this->value_ptr_); + I& ima = this->image_; + if (this->value_ptr_ = 0 || ! ima.has_data()) + return false; + int o = this->value_ptr_ - ima.buffer(); + mln_point(I) p = ima.point_at_offset(o); + return ima.owns_(p); } # endif // ! MLN_INCLUDE_ONLY Index: mln/core/hline2d.hh --- mln/core/hline2d.hh (revision 1024) +++ mln/core/hline2d.hh (working copy) @@ -68,7 +68,7 @@ /*! \brief Piter type to browse a hline such as: "for each row * (decreasing), for each column (decreasing)." */ - typedef dpoints_fwd_piter<dpoint2d> bkd_qiter; // FIXME: bkd!!! + typedef internal::fixme bkd_qiter; /*! \brief Same as fwd_qiter. */ Index: mln/core/concept/image.hh --- mln/core/concept/image.hh (revision 1024) +++ mln/core/concept/image.hh (working copy) @@ -74,7 +74,9 @@ typedef point; typedef psite; - typedef piter; + typedef coord; + typedef dpoint; + typedef fwd_piter; typedef bkd_piter; @@ -99,7 +101,9 @@ typedef mln_point(E) point; typedef mln_psite(E) psite; - typedef mln_piter(E) piter; + typedef mln_coord(E) coord; + typedef mln_dpoint(E) dpoint; + typedef mln_fwd_piter(E) fwd_piter; typedef mln_bkd_piter(E) bkd_piter; Index: mln/core/concept/genpoint.hh --- mln/core/concept/genpoint.hh (revision 1024) +++ mln/core/concept/genpoint.hh (working copy) @@ -36,9 +36,6 @@ # include <mlc/same_coord.hh> # include <mln/core/concept/object.hh> - -# include <mln/core/macros.hh> -# include <mln/core/contract.hh> # include <mln/core/internal/force_exact.hh> @@ -73,7 +70,7 @@ either Point or operator point() const; - const point* pointer() const; + const point* pointer_() const; coord operator[](unsigned i) const; */ @@ -218,7 +215,7 @@ typedef mln_point(E) point; typedef mln_dpoint(E) dpoint; typedef mln_coord(E) coord; - const point* (E::*m1)() const = & E::pointer; + const point* (E::*m1)() const = & E::pointer_; m1 = 0; coord (E::*m2)(unsigned i) const = & E::operator[]; m2 = 0; Index: mln/core/concept/psite.hh --- mln/core/concept/psite.hh (revision 1024) +++ mln/core/concept/psite.hh (working copy) @@ -65,7 +65,7 @@ public GenPoint<E> { /* - const point* pointer() const + const point* pointer_() const { return & (exact(this)->operator point()); } Index: mln/core/concept/point.hh --- mln/core/concept/point.hh (revision 1024) +++ mln/core/concept/point.hh (working copy) @@ -55,7 +55,7 @@ /*! \brief The pointer is this point address. */ - const P* pointer() const; + const P* pointer_() const; protected: Point(); @@ -105,7 +105,7 @@ } template <typename P> - const P* Point::pointer() const + const P* Point::pointer_() const { return exact(this); } Index: mln/core/concept/pixel_iterator.hh --- mln/core/concept/pixel_iterator.hh (revision 1024) +++ mln/core/concept/pixel_iterator.hh (working copy) @@ -28,53 +28,53 @@ #ifndef MLN_CORE_CONCEPT_PIXEL_ITERATOR_HH # define MLN_CORE_CONCEPT_PIXEL_ITERATOR_HH -/*! \file mln/core/concept/fast_iterator.hh +/*! \file mln/core/concept/pixel_iterator.hh + * * \brief Definition of the concept of mln::Pixel_Iterator. */ # include <mln/core/concept/iterator.hh> # include <mln/core/concept/genpixel.hh> + namespace mln { - /*! \brief Pixel Iterator concept class. + /*! \brief Base class for the implementation of pixel iterator + * classes. * - * \see FIXME + * An iterator on pixels is an iterator that is bound to a + * particular image and that browses over a set of image pixels. + * + * \see mln::doc::Pixel_Iterator for a complete documentation of + * this class contents. */ template <typename E> - struct Pixel_Iterator : - public Iterator<E>, + struct Pixel_Iterator : public Iterator<E>, public GenPixel<E> { /* - typedef rvalue - typedef lvalue + typedef lvalue; + lvalue operator* (); - rvalue operator* () const; */ protected: Pixel_Iterator(); }; + #ifndef MLN_INCLUDE_ONLY template <typename E> Pixel_Iterator<E>::Pixel_Iterator() { - typedef mln_value(E) value; - typedef mln_rvalue(E) rvalue; typedef mln_lvalue(E) lvalue; - - lvalue (E::*m1)() = & E::operator*; - m1 = 0; - rvalue (E::*m2)() const = & E::operator*; - m2 = 0; + lvalue (E::*m)() = & E::operator*; + m = 0; } #endif // ! MLN_INCLUDE_ONLY - } // end of namespace mln Index: mln/core/concept/fast_image.hh --- mln/core/concept/fast_image.hh (revision 0) +++ mln/core/concept/fast_image.hh (revision 0) @@ -0,0 +1,96 @@ +// Copyright (C) 2007 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_CORE_CONCEPT_FAST_IMAGE_HH +# define MLN_CORE_CONCEPT_FAST_IMAGE_HH + +/*! \file mln/core/concept/fast_image.hh + * \brief Definition of the concept of mln::Fast_Image. + */ + +# include <mln/core/concept/image.hh> +# include <mln/core/trait/pixter.hh> + + +namespace mln +{ + + /*! \brief Base class for implementation of fast image classes. + * + * \see mln::doc::Fast_Image for a complete documentation of this + * class contents. + */ + template <typename E> + struct Fast_Image : public Image<E> + { + /* + unsigned border(); + + int offset(const dpoint& dp); + unsigned offset(const point& p); + point point_at_offset(unsigned o) const; + + const value* buffer() const; + */ + + protected: + Fast_Image(); + }; + + +# ifndef MLN_INCLUDE_ONLY + + template <typename E> + Fast_Image<E>::Fast_Image() + { + typedef mln_point(E) point; + typedef mln_dpoint(E) dpoint; + typedef mln_value(E) value; + + typedef mln_fwd_pixter(E) fwd_pixter; + typedef mln_bkd_pixter(E) bkd_pixter; + + int (E::*m1)(const dpoint&) = & E::offset; + m1 = 0; + unsigned (E::*m2)(const point&) = & E::offset; + m2 = 0; + point (E::*m3)(unsigned) const = & E::point_at_offset; + m3 = 0; + const value* (E::*m4)() const = & E::buffer; + m4 = 0; + unsigned (E::*m5)() const = & E::border; + m5 = 0; + + // FIXME: how to check that qixter are defined when W is unknown! + } + +# endif // ! MLN_INCLUDE_ONLY + +} // end of namespace mln + + +#endif // ! MLN_CORE_CONCEPT_FAST_IMAGE_HH Index: mln/core/concept/point_set.hh --- mln/core/concept/point_set.hh (revision 1024) +++ mln/core/concept/point_set.hh (working copy) @@ -29,6 +29,7 @@ # define MLN_CORE_CONCEPT_POINT_SET_HH /*! \file mln/core/concept/point_set.hh + * * \brief Definition of the concept of mln::Point_Set. */ @@ -50,7 +51,7 @@ /* typedef point; typedef psite; - typedef piter; + typedef fwd_piter; typedef bkd_piter; @@ -127,7 +128,6 @@ typedef mln_point(E) point; typedef mln_psite(E) psite; - typedef mln_piter(E) piter; typedef mln_fwd_piter(E) fwd_piter; typedef mln_bkd_piter(E) bkd_piter; Index: mln/core/concept/object.hh --- mln/core/concept/object.hh (revision 1024) +++ mln/core/concept/object.hh (working copy) @@ -38,6 +38,7 @@ # include <mln/core/macros.hh> # include <mln/core/contract.hh> +# include <mln/core/internal/fixme.hh> /*! \namespace mln Index: mln/core/concept/doc/image.hh --- mln/core/concept/doc/image.hh (revision 1024) +++ mln/core/concept/doc/image.hh (working copy) @@ -121,10 +121,16 @@ */ typedef void psite; - /*! \brief Piter associated type. - * \invariant This type has to derive from mln::Piter. + + /*! \brief Coordinate associated type. + */ + typedef void coord; + + /*! \brief Dpoint associated type. + * \invariant This type has to derive from mln::Dpoint. */ - typedef void piter; + typedef void dpoint; + /*! \brief Forward point iterator associated type. * \invariant This type has to derive from mln::Piter. @@ -136,6 +142,7 @@ */ typedef void bkd_piter; + /*! \brief Test if \p p belongs to the image domain. * * \param[in] p A point site. @@ -162,10 +169,6 @@ std::size_t npoints() const; }; - template <typename E> - struct Fast_Image : public Image<E> - { - }; } // end of namespace mln::doc Index: mln/core/concept/doc/genpoint.hh --- mln/core/concept/doc/genpoint.hh (revision 1024) +++ mln/core/concept/doc/genpoint.hh (working copy) @@ -79,7 +79,7 @@ * * \return A point address. */ - const point* pointer() const; + const point* pointer_() const; /*! \brief Read-only access to the \p i-th coordinate value. * Index: mln/core/concept/doc/pixel_iterator.hh --- mln/core/concept/doc/pixel_iterator.hh (revision 1024) +++ mln/core/concept/doc/pixel_iterator.hh (working copy) @@ -25,8 +25,9 @@ // reasons why the executable file might be covered by the GNU General // Public License. -/*! \file mln/core/concept/doc/fast_iterator.hh - * \brief This file documents the concept of mln::Fast_Iterator. +/*! \file mln/core/concept/doc/pixel_iterator.hh + * + * \brief This file documents the concept of mln::Pixel_Iterator. */ namespace mln @@ -34,24 +35,29 @@ namespace doc { + /*! \brief Documentation class for mln::Iterator. - * \see mln::Fast_Iterator + * \see mln::Pixel_Iterator */ template <typename E> - struct Fast_Iterator : public Iterator<E> + struct Pixel_Iterator : public Iterator<E>, + public GenPixel<E> { - /*! \brief Dereference operator. - * return the current pixel + /*! \brief Type returned by the read-write dereference operator. */ - rvalue operator* (); + typedef void lvalue; /*! \brief Dereference operator. - * return the current pixel + * + * \return The current pixel value; this value cannot be + * modified. */ lvalue operator* () const; + protected: - Fast_Iterator(); + Pixel_Iterator(); }; + } // end of namespace doc } // end of namespace mln Index: mln/core/concept/doc/fast_image.hh --- mln/core/concept/doc/fast_image.hh (revision 0) +++ mln/core/concept/doc/fast_image.hh (revision 0) @@ -0,0 +1,92 @@ +// Copyright (C) 2007 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. + +/*! \file mln/core/concept/doc/fast_image.hh + * \brief This file documents the concept of mln::Fast_Image. + */ + +namespace mln +{ + + namespace doc + { + + /*! \brief Documentation class for mln::Fast_Image. + * \see mln::Fast_Image + */ + template <typename E> + struct Fast_Image : public Image<E> + { + + /*! \brief Give the border thickness. + * + * \pre The image has to be initialized. + */ + unsigned border(); + + /*! \brief Give the offset corresponding to the delta-point \p + * dp. + * + * \param[in] p A delta-point. + * + * \pre The image has to be initialized. + */ + int offset(const dpoint& dp); + + + /*! \brief Give the offset of the point \p p. + * + * \param[in] p A point. + * + * \pre The image has to own the point \p p. + */ + unsigned offset(const point& p); + + + /*! \brief Give the point at offset \p o. + * + * \param[in] p An offset. + * + * \pre The image has to be initialized. + * \pre o < ncells() + */ + point point_at_offset(unsigned o) const; + + + /*! \brief Give a hook to the value buffer. + * + * \internal + * + * \pre The image has to be initialized. + */ + const value* buffer() const; + + }; + + } // end of namespace mln::doc + +} // end of namespace mln Index: mln/core/concept/doc/point_set.hh --- mln/core/concept/doc/point_set.hh (revision 1024) +++ mln/core/concept/doc/point_set.hh (working copy) @@ -50,10 +50,6 @@ */ typedef void psite; - /*! \brief Piter associated type. - */ - typedef void piter; - /*! \brief Forward Piter associated type. */ typedef void fwd_piter; Index: mln/core/concept/doc/genpixel.hh --- mln/core/concept/doc/genpixel.hh (revision 0) +++ mln/core/concept/doc/genpixel.hh (revision 0) @@ -0,0 +1,78 @@ +// Copyright (C) 2007 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_CORE_CONCEPT_DOC_GENPIXEL_HH +# define MLN_CORE_CONCEPT_DOC_GENPIXEL_HH + +/*! \file mln/core/concept/doc/genpixel.hh + * \brief This file documents the concept of mln::GenPixel. + */ + + +namespace mln +{ + + namespace doc + { + + /*! \brief Documentation class for mln::GenPixel. + * + * \see mln::GenPixel + */ + template <typename E> + struct GenPixel + { + + /// Value associated type. + typedef void value; + + /// Read-only value associated type. + typedef void rvalue; + + + /*! \brief Give the value of this generalized pixel. + * + * \return A read-only value. + */ + rvalue operator*() const; + + /*! \brief Give a hook to the value address. + * + * This method allows for iterators to refer to a generalized + * pixel. + * + * \return A pointer to the value address. + */ + value** address_() const; + + protected: + GenPixel(); + }; + + } // end of namespace mln::doc + +} // end of namespace mln Index: mln/core/concept/genpixel.hh --- mln/core/concept/genpixel.hh (revision 1024) +++ mln/core/concept/genpixel.hh (working copy) @@ -29,13 +29,11 @@ # define MLN_CORE_CONCEPT_GENPIXEL_HH /*! \file mln/core/concept/genpixel.hh + * * \brief Definition of the concept of mln::Genpixel. */ # include <mln/core/concept/object.hh> - -# include <mln/core/macros.hh> -# include <mln/core/contract.hh> # include <mln/core/internal/force_exact.hh> @@ -58,75 +56,29 @@ template <typename E> struct GenPixel { - /* - typedef rvalue; typedef value; - // FIXME: lvalue? - - - // FIXME: Inactivated: - // typedef ima; // not const! - // typedef psite; - // const ima& image() const; - // const psite& site() const; // FIXME ou cpy - // value* address(); + typedef rvalue; rvalue operator*() const; - lvalue operator*(); - - const value*& address() const; + value** address_() const; */ - protected: GenPixel(); - }; - /*! \brief Print a generalized pixel \p p into the output stream \p - * ostr. - * - * \param[in,out] ostr An output stream. - * \param[in] p A generalized pixel. - * - * \return The modified output stream \p ostr. - * - * \relates mln::GenPixel - */ - template <typename P> - std::ostream& operator<<(std::ostream& ostr, const GenPixel& pxl); - - # ifndef MLN_INCLUDE_ONLY - template <typename E> GenPixel<E>::GenPixel() { typedef mln_value(E) value; typedef mln_rvalue(E) rvalue; - - // FIXME: Inactivated: - -// typedef mln_ima(E) ima; -// typedef mln_psite(E) psite; -// const ima& (E::*m1)() const = & E::image; -// m1 = 0; -// const psite& (E::*m2)() const = & E::site; -// m2 = 0; - - rvalue (E::*m3)() const = & E::operator*; - m3 = 0; -// const value *const & (E::*m4)() const = & E::address; -// m4 = 0; - } - - - template <typename P> - std::ostream& operator<<(std::ostream& ostr, const GenPixel& pxl) - { - return ostr << pxl.psite() << '@' << & pxl.image(); + rvalue (E::*m1)() const = & E::operator*; + m1 = 0; + value** (E::*m2)() const = & E::address_; + m2 = 0; } # endif // ! MLN_INCLUDE_ONLY Index: mln/core/tags.hh --- mln/core/tags.hh (revision 0) +++ mln/core/tags.hh (revision 0) @@ -0,0 +1,50 @@ +// Copyright (C) 2007 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_CORE_TAGS_HH +# define MLN_CORE_TAGS_HH + +/*! \file mln/core/tags.hh + * + * \brief Define a set of tags for images types. + */ + + +namespace mln +{ + + namespace tag + { + + // FIXME: ... + + } // end of namespace mln::core + +} // end of namespace mln + + +#endif // ! MLN_CORE_TAGS_HH Index: mln/core/vline2d.hh --- mln/core/vline2d.hh (revision 1024) +++ mln/core/vline2d.hh (working copy) @@ -70,7 +70,7 @@ /*! \brief Piter type to browse a vline such as: "for each row * (decreasing), for each column (decreasing)." */ - typedef dpoints_fwd_piter<dpoint2d> bkd_qiter; // FIXME: bkd!!! + typedef internal::fixme bkd_qiter; /*! \brief Same as fwd_qiter. */ Index: mln/core/image2d_b.hh --- mln/core/image2d_b.hh (revision 1024) +++ mln/core/image2d_b.hh (working copy) @@ -39,8 +39,12 @@ # include <mln/border/thickness.hh> # include <mln/fun/all.hh> -# include <mln/core/pixter2d_b.hh> -# include <mln/core/dpoints_pixter.hh> + + +// FIXME: + +// # include <mln/core/pixter2d_b.hh> +// # include <mln/core/dpoints_pixter.hh> namespace mln @@ -64,18 +68,19 @@ typedef dpoint2d dpoint; typedef mln_fwd_piter(box2d) fwd_piter; typedef mln_bkd_piter(box2d) bkd_piter; - typedef fwd_piter piter; // end of warning - /// Forward pixel iterator associated to image2d - typedef fwd_pixter2d_b<T> fwd_pixter; - /// Foward pixel iterator on dpoints assoicated to image 2d - typedef dpoints_pixter< image2d_b<T> > fwd_qixter; + // FIXME: +// /// Forward pixel iterator associated to image2d +// typedef fwd_pixter2d_b<T> fwd_pixter; - typedef fwd_pixter pixter; - typedef fwd_qixter qixter; +// /// Foward pixel iterator on dpoints assoicated to image 2d +// typedef dpoints_pixter< image2d_b<T> > fwd_qixter; + +// typedef fwd_pixter pixter; +// typedef fwd_qixter qixter; /// Value associated type. @@ -153,14 +158,17 @@ /// Fast Image method - /// Return the offset corresponding to the dpoint \p dp. + /// Give the offset corresponding to the delta-point \p dp. int offset(const dpoint2d& dp); - /// Return the offset corresponding to the dpoint \p p. - int offset(const point2d& p); + /// Give the offset corresponding to the point \p p. + unsigned offset(const point2d& p); + + /// Give the point corresponding to the offset \p o. + point2d point_at_offset(unsigned o) const; - /// Return the point corresponding to the offset \p offset. - point2d point_at_offset(int offset); + /// Give a hook to the value buffer. + const T* buffer() const; private: @@ -340,6 +348,46 @@ deallocate_(); } + template <typename T> + const T* + image2d_b<T>::buffer() const + { + mln_precondition(this->has_data()); + return buffer_; + } + + template <typename T> + int + image2d_b<T>::offset(const dpoint2d& dp) + { + mln_precondition(this->has_data()); + int o = dp[0] * vb_.len(1) + dp[1]; + return o; + } + + template <typename T> + unsigned + image2d_b<T>::offset(const point2d& p) + { + mln_precondition(this->owns_(p)); + unsigned o = & this->operator()(p) - this->buffer_; + mln_postcondition(p = point_at_offset(o)); + return o; + } + + template <typename T> + point2d + image2d_b<T>::point_at_offset(unsigned o) const + { + mln_precondition(o < ncells()); + mln_precondition(this->has_data()); + point2d p = make::point2d(o / vb_.len(1) + vb_.min_row(), + o % vb_.len(1) + vb_.min_col()); + mln_postcondition(& this->operator()(p) = this->buffer_ + o); + return p; + } + + // private template <typename T> @@ -388,28 +436,52 @@ } } - template <typename T> - int - image2d_b<T>::offset(const dpoint2d& dp) +# endif // ! MLN_INCLUDE_ONLY + +} // end of namespace mln + + + +# include <mln/core/trait/pixter.hh> +# include <mln/core/dpoints_pixter.hh> +# include <mln/core/pixter2d_b.hh> + + +namespace mln +{ + + namespace trait { - return dp[0] * vb_.len(1) + dp[1]; - } + + // pixter template <typename T> - int - image2d_b<T>::offset(const point2d& p) + struct fwd_pixter< image2d_b<T> > { - return p[0] * vb_.len(1) + p[1]; - } + typedef fwd_pixter2d_b< image2d_b<T> > ret; + }; template <typename T> - point2d - image2d_b<T>::point_at_offset(int offset) + struct bkd_pixter< image2d_b<T> > { - return point2d(offset / (vb_.len(1) + 1), offset % (vb_.len(1) + 1)); - } + typedef internal::fixme ret; + }; -# endif // ! MLN_INCLUDE_ONLY + // qixter + + template <typename T, typename W> + struct fwd_qixter< image2d_b<T>, W > + { + typedef dpoints_fwd_pixter< image2d_b<T> > ret; + }; + + template <typename T, typename W> + struct bkd_qixter< image2d_b<T>, W > + { + typedef internal::fixme ret; + }; + + } // end of namespace mln::trait } // end of namespace mln Index: mln/metal/none.hh --- mln/metal/none.hh (revision 0) +++ mln/metal/none.hh (revision 0) @@ -0,0 +1,57 @@ +// Copyright (C) 2007 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_CORE_METAL_NONE_HH +# define MLN_CORE_METAL_NONE_HH + +/*! \file mln/core/metal/none.hh + * + * \brief Definition of a type that means "none". + */ + + +namespace mln +{ + + namespace metal + { + + + /*! \brief Type meaning "none". + * + * This type is for use when a type is required and none can be + * provided. + */ + struct none; + + + } // end of namespace mln::metal + +} // end of namespace mln + + +#endif // ! MLN_CORE_METAL_NONE_HH Index: mln/level/fast_median.hh --- mln/level/fast_median.hh (revision 1024) +++ mln/level/fast_median.hh (working copy) @@ -36,7 +36,6 @@ # include <mln/core/concept/image.hh> # include <mln/core/window2d.hh> # include <mln/accu/median.hh> -# include <ostream> namespace mln @@ -85,24 +84,20 @@ win_bot = win - (win + up), win_top = (win + up) - win; - - - //Fixme bug - point2d p; - accu::median_on<mln_value(I)> med; // initialization - p = input.domain().pmin() + up; - mln_qixter(I) - q_fp(win_fwd_plus, p, input), q_fm(win_fwd_minus, p, input), - q_bp(win_bkd_plus, p, input), q_bm(win_bkd_minus, p, input), - q_top(win_top, p, input), q_bot(win_bot, p, input); + point2d p = input.domain().pmin() + up; + + mln_qixter(I, window2d) + q_fp(input, win_fwd_plus, p), q_fm(input, win_fwd_minus, p), + q_bp(input, win_bkd_plus, p), q_bm(input, win_bkd_minus, p), + q_top(input, win_top, p), q_bot(input, win_bot, p); med.init(); { - mln_qixter(I) q(win, p, input); + mln_qixter(I, W) q(input, win, p); for_all(q) med.take(*q); } Index: mln/level/median.hh --- mln/level/median.hh (revision 1024) +++ mln/level/median.hh (working copy) @@ -160,21 +160,22 @@ template <typename I, typename O> void median(const I& input, const hline2d& win, O& output) { - const int + typedef mln_coord(I) coord; + const coord max_row = input.max_row(), min_col = input.min_col(), max_col = input.max_col(); - const unsigned half = win.length() / 2; + const coord half = win.length() / 2; point2d p; - int& row = p.row(); - int& col = p.col(); + coord& row = p.row(); + coord& col = p.col(); accu::median_on<mln_value(I)> med; for (row = input.min_row(); row <= max_row; ++row) { - int ct, cu; + coord ct, cu; // initialization (before first point of the row) med.init(); Index: mln/canvas/sbrowsing.hh --- mln/canvas/sbrowsing.hh (revision 1024) +++ mln/canvas/sbrowsing.hh (working copy) @@ -45,12 +45,18 @@ /*! FIXME: Doc! * * - * F should be - * { - * input; p; --- attributes - * void init(); - * void down(); void fwd(); void bkd(); --- methods - * } + * F shall feature: \n + * { \n + * --- as attributes: \n + * input; \n + * p; \n + * --- as methods: \n + * void init(); \n + * void down(); \n + * void fwd(); \n + * void bkd(); \n + * } \n + * */ template <typename F> void sbrowsing(F& f);

1 0

1024: Repair lineary pixel iterator.
by Nicolas Ballas 19 Jul '07

19 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Nicolas Ballas <ballas(a)lrde.epita.fr> Repair lineary pixel iterator. * mln/core/pixter2d_b.hh: Adapt pixter2d_b to the new pixel iterator hierachy. * mln/core/internal/pixel_iterator_base.hh: Fix dead code. * mln/core/internal/lineary_pixel_iterator_base.hh: New factorization class. internal/lineary_pixel_iterator_base.hh | 109 ++++++++++++++++++++++++++++++++ internal/pixel_iterator_base.hh | 101 +---------------------------- pixter2d_b.hh | 31 +++++---- 3 files changed, 131 insertions(+), 110 deletions(-) Index: mln/core/pixter2d_b.hh --- mln/core/pixter2d_b.hh (revision 1023) +++ mln/core/pixter2d_b.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_CORE_PIXTER2D_B_HH # define MLN_CORE_PIXTER2D_B_HH -# include <mln/core/internal/pixel_iterator_base.hh> +# include <mln/core/internal/lineary_pixel_iterator_base.hh> # include <mln/core/point2d.hh> # include <iostream> @@ -49,8 +49,9 @@ template <typename T> class fwd_pixter2d_b : - public internal::pixel_iterator_base_< fwd_pixter2d_b<T>, image2d_b<T> > + public internal::lineary_pixel_iterator_base_<image2d_b<T>, fwd_pixter2d_b<T> > { + typedef internal::lineary_pixel_iterator_base_<image2d_b<T>, fwd_pixter2d_b<T> > super; public: /// Image pixel value type. typedef mln_value(image2d_b<T>) value; @@ -59,37 +60,41 @@ * * \param[in] ima Image to iterate. */ - fwd_pixter2d_b(image2d_b<T>& ima); - /// Move the iterator to the next elements + fwd_pixter2d_b(image2d_b<T>& image); + /// Move the iterator on the next element. void next_(); private: + /// Size of the image border. unsigned border_size_; + /// Row offset. unsigned row_offset_; + /// End of a row. value *eor_; }; #ifndef MLN_INCLUDE_ONLY template <typename T> - fwd_pixter2d_b<T>::fwd_pixter2d_b(image2d_b<T>& ima) : - border_size_(ima.border()), - row_offset_((ima.domain().pmax()[1] - ima.domain().pmin()[1]) + fwd_pixter2d_b<T>::fwd_pixter2d_b(image2d_b<T>& image) : + super(image), + border_size_(image.border()), + row_offset_((image.domain().pmax()[1] - image.domain().pmin()[1]) + 2 * border_size_ + 1) { - this->start_ = &ima(ima.domain().pmin()); - this->eor_ = &ima(make::point2d(ima.domain().pmin()[0], ima.domain().pmax()[1])) + 1; - this->eoi_ = &ima(ima.domain().pmax()) + 1; + this->start_ = &image(image.domain().pmin()); + this->eor_ = &image(make::point2d(image.domain().pmin()[0], image.domain().pmax()[1])) + 1; + this->eoi_ = &image(image.domain().pmax()) + 1; } template <typename T> void fwd_pixter2d_b<T>::next_() { - ++(this->current_); + ++(this->value_ptr_); - if (this->current_ == this->eor_ && this->current_ != this->eoi_) + if (this->value_ptr_ == this->eor_ && this->value_ptr_ != this->eoi_) { - this->current_ += 2 * this->border_size_; + this->value_ptr_ += 2 * this->border_size_; this->eor_ += this->row_offset_; } } Index: mln/core/internal/pixel_iterator_base.hh --- mln/core/internal/pixel_iterator_base.hh (revision 1023) +++ mln/core/internal/pixel_iterator_base.hh (working copy) @@ -93,12 +93,12 @@ protected: - /// Current pixel / value - value* value_ptr_; - /// Image associated to the iterator I& image_; + /// Current pixel / value + value* value_ptr_; + // FIXME: Inactivated: // /// Psite of the pixel @@ -174,101 +174,8 @@ // return p_; // } - - - - -// FIXME: Dead code -// ################# - - - -// /*! \brief pixel_iterator_base_ class -// */ -// template <typename E, typename I> -// class pixel_iterator_base_ : public Pixel_Iterator<E> -// { -// public: -// /// I pixel value type. -// typedef mln_value(I) value; -// /// I pixel rvalue type. -// typedef mln_value(I)& rvalue; -// /// I pixel lvalue type -// typedef mln_value(I) lvalue; - -// // Go to the beginning of the image. -// void start(); -// // Go on the next element. -// void next_(); -// // Invalidate the iterator. -// void invalidate(); -// // Is the iterator referencing a correct position in the image? -// bool is_valid() const; - - -// // Return the current pixel. -// rvalue operator* (); -// lvalue operator* () const; - -// protected: -// // beginning of the image -// value* start_; -// // end of the image -// value* eoi_; -// // current position in the image -// value* current_; -// // End of Factoriasable - -// pixel_iterator_base_(); -// }; - -// #ifndef MLN_INCLUDE_ONLY - -// template <typename E, typename I> -// pixel_iterator_base_<E,I>::pixel_iterator_base_() -// { -// } - -// template <typename E, typename I> -// void pixel_iterator_base_<E,I>::start() -// { -// current_ = start_; -// } - -// template <typename E, typename I> -// void pixel_iterator_base_<E,I>::next_() -// { -// ++current_; -// } - -// template <typename E, typename I> -// void pixel_iterator_base_<E,I>::invalidate() -// { -// current_ = eoi_; -// } - -// template <typename E, typename I> -// bool pixel_iterator_base_<E,I>::is_valid() const -// { -// return (current_ != eoi_); -// } - -// template <typename E, typename I> -// typename pixel_iterator_base_<E,I>::rvalue -// pixel_iterator_base_<E,I>::operator*() -// { -// return *current_; -// } - - -// template <typename E, typename I> -// typename pixel_iterator_base_<E,I>::lvalue -// pixel_iterator_base_<E,I>::operator*() const -// { -// return *current_; -// } - #endif // ! MLN_INCLUDE_ONLY + } // end of namespace internal } // end of namespace mln Index: mln/core/internal/lineary_pixel_iterator_base.hh --- mln/core/internal/lineary_pixel_iterator_base.hh (revision 0) +++ mln/core/internal/lineary_pixel_iterator_base.hh (revision 0) @@ -0,0 +1,109 @@ +// Copyright (C) 2007 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_CORE_INTERNAL_LINEARY_PIXEL_ITERATOR_BASE_HH +# define MLN_CORE_INTERNAL_LINEARY_PIXEL_ITERATOR_BASE_HH + +/*! \file mln/core/internal/lineary_pixel_iterator_base.hh + * + * \brief Code factorization for lineary pixel iterator classes. + */ + +# include <mln/core/internal/pixel_iterator_base.hh> + +namespace mln +{ + + namespace internal + { + + /*! \brief lineary_pixel_iterator_base_ class + * + */ + template <typename I, typename E> + class lineary_pixel_iterator_base_ : public pixel_iterator_base_<I, E> + { + typedef pixel_iterator_base_<I, E> super; + public: + /// Image pixel value type. + typedef mln_value(I) value; + + /// Go to the beginning of the image. + void start(); + /// Invalidate the iterator. + void invalidate(); + /// Is the iterator referencing a correct position in the image? + bool is_valid() const; + + + protected: + /// beginning of the image + value* start_; + /// end of the image + value* eoi_; + + /// Constructor + lineary_pixel_iterator_base_(I& image); + }; + +#ifndef MLN_INCLUDE_ONLY + + template <typename I, typename E> + lineary_pixel_iterator_base_<I, E>::lineary_pixel_iterator_base_(I& image) + : super(image) + { + } + + template <typename I, typename E> + void + lineary_pixel_iterator_base_<I, E>::start() + { + this->value_ptr_ = start_; + } + + template <typename I, typename E> + void + lineary_pixel_iterator_base_<I, E>::invalidate() + { + this->value_ptr_ = eoi_; + } + + template <typename I, typename E> + bool + lineary_pixel_iterator_base_<I, E>::is_valid() const + { + return (this->value_ptr_ != eoi_); + } + +#endif // ! MLN_INCLUDE_ONLY + + } // end of namespace internal + +} // end of namespace mln + + +#endif // ! MLN_CORE_INTERNAL_LINEARY_PIXEL_ITERATOR_BASE_HH

2 1

1023: Fix pixel iterators on sets of dpoints.
by Thierry Geraud 19 Jul '07

19 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Thierry Geraud <thierry.geraud(a)lrde.epita.fr> Fix pixel iterators on sets of dpoints. * tests/fast_median.cc (test): New. * tests/naive_median.cc: Update. * mln/debug/iota.hh: New. * mln/core/pixel.hh (value_ptr_): New. * mln/core/dpoints_piter.hh (center_point): New. (update_p_): Rename as... (update): ...this. * mln/core/concept/genpixel.hh: Inactivate some code. * mln/core/dpoints_pixter.hh: Update. * mln/core/internal/pixel_iterator_base.hh: Clean-up. * mln/level/was.hmedian.hh (hmedian): Move... * mln/level/was.median.hh: ...here. * mln/level/was.hmedian.hh: Remove. * mln/level/fast_median.hh: Clean-up. Contents mln/core/concept/genpixel.hh | 36 +++-- mln/core/dpoints_piter.hh | 40 ++++- mln/core/dpoints_pixter.hh | 120 +++++++++++++--- mln/core/internal/pixel_iterator_base.hh | 223 +++++++++++++++++-------------- mln/core/pixel.hh | 20 -- mln/debug/iota.hh | 69 +++++++++ mln/level/fast_median.hh | 15 -- mln/level/was.median.hh | 61 ++++++++ tests/fast_median.cc | 47 ++++++ tests/naive_median.cc | 4 10 files changed, 468 insertions(+), 167 deletions(-) Index: tests/fast_median.cc --- tests/fast_median.cc (revision 1022) +++ tests/fast_median.cc (working copy) @@ -37,9 +37,40 @@ #include <mln/io/save_pgm.hh> #include <mln/value/int_u.hh> +#include <mln/debug/iota.hh> +#include <mln/debug/println.hh> #include <mln/level/fast_median.hh> -#include <mln/level/approx/median.hh> +#include <mln/core/dpoints_pixter.hh> +#include <mln/core/pixel.hh> + + +namespace mln +{ + + template <typename I, typename W> + void test(I& input, const W& win) + { + mln_point(I) p; + p.row() = p.col() = 1; + + { + mln_qixter(I) qix(win, p, input); + for_all(qix) + std::cout << *qix << ' '; + std::cout << " : " << qix.center_value() << std::endl; + } + + { + pixel pix(input, p); + mln_qixter(I) qix(win, pix, input); + for_all(qix) + std::cout << *qix << ' '; + std::cout << " : " << qix.center_value() << std::endl; + } + } + +} @@ -48,6 +79,17 @@ using namespace mln; using value::int_u8; +// { +// rectangle2d rect(3, 3); +// border::thickness = 4; +// image2d_b<int_u8> ima(3, 3); +// debug::iota(ima); +// debug::println(ima); +// test(ima, rect); +// } + + + { rectangle2d rect(51, 51); border::thickness = 52; @@ -57,7 +99,6 @@ level::fast_median(lena, rect, out); io::save_pgm(out, "out.pgm"); + } -// level::approx::median(lena, rec, out); -// io::save_pgm(out, "outa.pgm"); } Index: tests/naive_median.cc --- tests/naive_median.cc (revision 1022) +++ tests/naive_median.cc (working copy) @@ -46,8 +46,8 @@ int main() { - rectangle2d rec(8, 8); - border::thickness = 9; + rectangle2d rec(51, 51); + border::thickness = 52; image2d_b<int_u8> lena = io::load_pgm("../img/lena.pgm"), Index: mln/debug/iota.hh --- mln/debug/iota.hh (revision 0) +++ mln/debug/iota.hh (revision 0) @@ -0,0 +1,69 @@ +// Copyright (C) 2007 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_DEBUG_IOTA_HH +# define MLN_DEBUG_IOTA_HH + +/*! \file mln/debug/iota.hh + * + * \brief FIXME + */ + +# include <mln/core/concept/image.hh> + + +namespace mln +{ + + namespace debug + { + + /// FIXME + template <typename I> + void iota(Image& input); + + +# ifndef MLN_INCLUDE_ONLY + + template <typename I> + void iota(Image& input_) + { + unsigned i = 0; + I& input = exact(input_); + mln_piter(I) p(input.domain()); + for_all(p) + input(p) = ++i; + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::debug + +} // end of namespace mln + + +#endif // ! MLN_DEBUG_IOTA_HH Index: mln/core/pixel.hh --- mln/core/pixel.hh (revision 1022) +++ mln/core/pixel.hh (working copy) @@ -76,15 +76,15 @@ psite& site(); mln_rvalue(I) operator*() const; - typename pixel_lvalue::ret operator*(); + mln_lvalue(I) operator*(); - const value* address() const; - value* address(); + value** address() const; protected: I& ima_; psite p_; + value* value_ptr_; }; @@ -101,6 +101,7 @@ : ima_(image), p_(p) { + value_ptr_ = & ima_(p); } template <typename I> @@ -132,24 +133,17 @@ } template <typename I> - typename pixel_lvalue::ret + mln_lvalue(I) pixel::operator*() { return ima_(p_); } template <typename I> - const mln_value(I)* + mln_value(I)** pixel::address() const { - return & ima_(p_); - } - - template <typename I> - mln_value(I)* - pixel::address() - { - return & ima_(p_); + return (mln_value(I)**)(& this->value_ptr_); } # endif // ! MLN_INCLUDE_ONLY Index: mln/core/dpoints_piter.hh --- mln/core/dpoints_piter.hh (revision 1022) +++ mln/core/dpoints_piter.hh (working copy) @@ -71,13 +71,13 @@ * \param[in] p_ref Center point to iterate around. */ template <typename Dps, typename Pref> - dpoints_fwd_piter(const Dps& dps, + dpoints_fwd_piter(const Dps& dps, // FIXME: explicitly set_of_<D>? const GenPoint<Pref>& p_ref); /// Convertion to point. operator mln_point(D) () const; - /// Address of the point. + /// Address of the point this iterator designates. const point* pointer() const; /// Test the iterator validity. @@ -95,13 +95,20 @@ /// Give the i-th coordinate. coord operator[](unsigned i) const; + /// The point around which this iterator moves. + const point& center_point() const; + + /// Force this iterator to update its location to take into + /// account that its center point may have moved. + void update(); + private: const std::vector<D>& dps_; - const point& p_ref_; + const point& p_ref_; // reference point (or "center point") unsigned i_; - point p_; - void update_p_(); + point p_; // location of this iterator; p_ makes this iterator be + // itself a potential center point (Cf. the pointer() method). }; @@ -120,7 +127,7 @@ template <typename D> dpoints_fwd_piter<D>::operator mln_point(D) () const { - assert(is_valid()); + mln_precondition(is_valid()); return p_; } @@ -150,7 +157,7 @@ dpoints_fwd_piter<D>::start() { i_ = 0; - update_p_(); + update(); } template <typename D> @@ -158,12 +165,19 @@ dpoints_fwd_piter<D>::next_() { ++i_; - update_p_(); + update(); + } + + template <typename D> + const mln_point(D)& + dpoints_fwd_piter<D>::center_point() const + { + return p_ref_; } template <typename D> void - dpoints_fwd_piter<D>::update_p_() + dpoints_fwd_piter<D>::update() { if (is_valid()) p_ = p_ref_ + dps_[i_]; @@ -173,6 +187,14 @@ mln_coord(D) dpoints_fwd_piter<D>::operator[](unsigned i) const { + mln_precondition(is_valid()); + + // below we test that no update is required + // meaning that p_ref_ has not moved or that + // the user has explicitly called update() + mln_precondition(p_ref_ + dps_[i_] = p_); + // FIXME: Explain this issue in the class documentation... + return p_[i]; } Index: mln/core/concept/genpixel.hh --- mln/core/concept/genpixel.hh (revision 1022) +++ mln/core/concept/genpixel.hh (working copy) @@ -60,17 +60,22 @@ { /* - typedef ima; // not const! - typede psite; typedef rvalue; typedef value; + // FIXME: lvalue? - const ima& image() const; - const psite& site() const; // FIXME ou cpy + + // FIXME: Inactivated: + // typedef ima; // not const! + // typedef psite; + // const ima& image() const; + // const psite& site() const; // FIXME ou cpy + // value* address(); rvalue operator*() const; + lvalue operator*(); - const value* address() const; + const value*& address() const; */ protected: @@ -89,7 +94,7 @@ * \relates mln::GenPixel */ template <typename P> - std::ostream& operator<<(std::ostream& ostr, const GenPixel& p); + std::ostream& operator<<(std::ostream& ostr, const GenPixel& pxl); @@ -101,17 +106,20 @@ { typedef mln_value(E) value; typedef mln_rvalue(E) rvalue; - typedef mln_ima(E) ima; - typedef mln_psite(E) psite; - const ima& (E::*m1)() const = & E::image; - m1 = 0; - const psite& (E::*m2)() const = & E::site; - m2 = 0; + // FIXME: Inactivated: + +// typedef mln_ima(E) ima; +// typedef mln_psite(E) psite; +// const ima& (E::*m1)() const = & E::image; +// m1 = 0; +// const psite& (E::*m2)() const = & E::site; +// m2 = 0; + rvalue (E::*m3)() const = & E::operator*; m3 = 0; - const value* (E::*m4)() const = & E::address; - m4 = 0; +// const value *const & (E::*m4)() const = & E::address; +// m4 = 0; } Index: mln/core/dpoints_pixter.hh --- mln/core/dpoints_pixter.hh (revision 1022) +++ mln/core/dpoints_pixter.hh (working copy) @@ -35,6 +35,7 @@ namespace mln { + /*! \brief A generic forward iterator on image pixel of a windows and of * neighborhoods. * @@ -45,6 +46,7 @@ { typedef typename internal::pixel_iterator_base_< I, dpoints_pixter > super; public: + /// Image pixel value typedef mln_value(I) value; /// Image pixel rvalue @@ -58,73 +60,149 @@ * \param[in] p_ref Center point to iterate around. * \param[in] ima Image to iterate. */ - template <typename Dps> - dpoints_pixter(const Dps& dps, const typename I::psite& p_ref,I& ima); + template <typename Dps, typename Pref> + dpoints_pixter(const Dps& dps, + const GenPoint<Pref>& p_ref, + I& image); + + template <typename Dps, typename Pref> + dpoints_pixter(const Dps& dps, + const GenPixel<Pref>& p_ref, + I& image); - /// Set the iterator at the start. + /// Start an iteration. void start(); + /// Go to the next pixel. void next_(); + /// Invalidate the iterator. void invalidate(); - /// Is the iterator valid? + + /// Test the iterator validity. bool is_valid() const; + /// Force this iterator to update its location to take into + /// account that its center point may have moved. + void update(); + + /// The value around which this iterator moves. + const value& center_value() const; + private: + /// offset of each dpoints std::vector<int> offset_; + /// current offset - std::vector<int>::iterator i_; - /// reference pixel in the image - value* pixref_; + int i_; + + /// reference pixel / point in the image + value** value_ref_; + // or: + const mln_point(I)* p_ref_; + + + template <typename Dps> + void init_(const Dps& dps); }; #ifndef MLN_INCLUDE_ONLY template <typename I> + template <typename Dps, typename Pref> + dpoints_pixter::dpoints_pixter(const Dps& dps, + const GenPoint<Pref>& p_ref, + I& image) + : super(image) + { + p_ref_ = internal::force_exact<Pref>(p_ref).pointer(); + value_ref_ = 0; + init_(dps); + } + + + template <typename I> + template <typename Dps, typename Pref> + dpoints_pixter::dpoints_pixter(const Dps& dps, + const GenPixel<Pref>& p_ref, + I& image) + : super(image) + { + p_ref_ = 0; + value_ref_ = internal::force_exact<Pref>(p_ref).address(); + init_(dps); + } + + template <typename I> + const mln_value(I)& + dpoints_pixter::center_value() const + { + mln_invariant(value_ref_ != 0 || p_ref_ != 0); + if (p_ref_) + return image_(*p_ref_); + else + return ** value_ref_; + } + + template <typename I> template <typename Dps> - dpoints_pixter::dpoints_pixter(const Dps& dps, const typename I::psite& ref, I& ima) : super(ima) + void + dpoints_pixter::init_(const Dps& dps) { - this->p_ = ref; - pixref_ = &ima(ref); - for (typename std::vector<typename I::dpoint>::const_iterator it = dps.vec().begin(); - it != dps.vec().end(); - ++it) - offset_.push_back(ima.offset(*it)); + for (unsigned i = 0; i < dps.nelements(); ++i) + offset_.push_back(this->image_.offset(dps.element(i))); + + if (dps.nelements() > 1) + for (unsigned i = dps.nelements() - 1; i > 0; --i) + offset_[i] -= offset_[i - 1]; + + // offset_[0] is absolute; other offsets are relative. + invalidate(); } + template <typename I> + void dpoints_pixter::update() + { + if (is_valid()) + { + if (p_ref_) + this->value_ptr_ = & image_(*p_ref_) + offset_[i_]; + else + this->value_ptr_ = * value_ref_ + offset_[i_]; + } + } template <typename I> void dpoints_pixter::start() { - pixref_ = &ima_(this->p_); - i_ = offset_.begin(); - this->current_value_ = pixref_ + *i_; + i_ = 0; + update(); } template <typename I> void dpoints_pixter::next_() { ++i_; - this->current_value_ = pixref_ + *i_; + this->value_ptr_ += offset_[i_]; } template <typename I> bool dpoints_pixter::is_valid() const { - return i_ != offset_.end(); + return i_ != offset_.size(); } template <typename I> void dpoints_pixter::invalidate() { - i_ = offset_.end(); + i_ = offset_.size(); } #endif // ! MLN_INCLUDE_ONLY -} +} // end of namespace mln #endif // MLN_CORE_DPOINTS_PIXTER_HH Index: mln/core/internal/pixel_iterator_base.hh --- mln/core/internal/pixel_iterator_base.hh (revision 1022) +++ mln/core/internal/pixel_iterator_base.hh (working copy) @@ -28,10 +28,9 @@ #ifndef MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH # define MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH -/*! \file mln/core/internal/fast_iterator_base.hh +/*! \file mln/core/internal/pixel_iterator_base.hh * - * \brief Base class for Fast_Iterator concept implementation. - * Usefull for code factorisation + * \brief Base class for Pixel_Iterator concept implementation classes. */ # include <mln/core/concept/pixel_iterator.hh> @@ -44,126 +43,158 @@ { /*! \brief pixel_iterator_base_ class + * */ - template <typename Image, typename Exact> - class pixel_iterator_base_ : public Pixel_Iterator<Exact> + template <typename I, typename E> + class pixel_iterator_base_ : public Pixel_Iterator<E> { public: - /// Image type. - typedef Image ima; + /// Image value type. - typedef mln_value(Image) value; + typedef mln_value(I) value; + /// Image lvalue type. - typedef mln_lvalue(Image) lvalue; + typedef mln_lvalue(I) lvalue; + /// Image rvalue type. - typedef mln_rvalue(Image) rvalue; - /// Image psite type. - typedef mln_psite(Image) psite; + typedef mln_rvalue(I) rvalue; + - /// Get the image associated to the current pixel iterator. - const ima& image() const; /// pixel iterator value. lvalue operator* (); + /// Get the pixel iterator value. rvalue operator* () const; + + /// Address of the current iterator value/pixel. - const value* address() const; - /// Address of the current iterator value/pixel. - value* address(); + value** address() const; + + + // FIXME: Inactivated: - /// psite associated to the iterator. - const psite& site() const; - /// psite associated to the iterator. - psite& site(); +// /// I type. +// typedef I ima; + +// /// Image psite type. +// typedef mln_psite(I) psite; + +// /// Get the image associated to the current pixel iterator. +// const ima& image() const; + +// /// psite associated to the iterator. +// const psite& site() const; + +// /// psite associated to the iterator. +// psite& site(); + +// /// Address of the current iterator value/pixel. +// value* address(); protected: - /// Current pixel value - value* current_value_; + + /// Current pixel / value + value* value_ptr_; + /// Image associated to the iterator - ima& ima_; - /// Psite of the pixel - psite p_; + I& image_; + + // FIXME: Inactivated: - pixel_iterator_base_(ima& source); +// /// Psite of the pixel +// psite p_; + + pixel_iterator_base_(I& image); }; #ifndef MLN_INCLUDE_ONLY - template <typename Image, typename Exact> - pixel_iterator_base_<Image, Exact>::pixel_iterator_base_(ima& source) : - ima_(source) + template <typename I, typename E> + pixel_iterator_base_<I, E>::pixel_iterator_base_(I& image) : + image_(image), + value_ptr_(0) { } - template <typename Image, typename Exact> - const typename pixel_iterator_base_<Image, Exact>::ima& - pixel_iterator_base_<Image, Exact>::image() const - { - return ima_; - } +// template <typename I, typename E> +// const typename pixel_iterator_base_<I, E>::ima& +// pixel_iterator_base_<I, E>::image() const +// { +// return ima_; +// } - template <typename Image, typename Exact> - typename pixel_iterator_base_<Image, Exact>::lvalue - pixel_iterator_base_<Image, Exact>::operator* () + template <typename I, typename E> + mln_lvalue(I) + pixel_iterator_base_<I, E>::operator* () { - assert(exact(this)->is_valid()); - return *current_value_; + mln_precondition(exact(this)->is_valid()); + mln_precondition(value_ptr_ != 0); + return *value_ptr_; } - template <typename Image, typename Exact> - typename pixel_iterator_base_<Image, Exact>::rvalue - pixel_iterator_base_<Image, Exact>::operator* () const + template <typename I, typename E> + mln_rvalue(I) + pixel_iterator_base_<I, E>::operator* () const { - assert(exact(this)->is_valid()); - return *current_value_; + mln_precondition(exact(this)->is_valid()); + mln_precondition(value_ptr_ != 0); + return *value_ptr_; } - template <typename Image, typename Exact> - const typename pixel_iterator_base_<Image, Exact>::value* - pixel_iterator_base_<Image, Exact>::address() const + template <typename I, typename E> + mln_value(I) ** + pixel_iterator_base_<I, E>::address() const { - assert(exact(this)->is_valid()); - return current_value_; + mln_precondition(exact(this)->is_valid()); + mln_precondition(value_ptr_ != 0); + return & (mln_value(I)*)(value_ptr_); } - template <typename Image, typename Exact> - typename pixel_iterator_base_<Image, Exact>::value* - pixel_iterator_base_<Image, Exact>::address() - { - assert(exact(this)->is_valid()); - return current_value_; - } +// template <typename I, typename E> +// mln_value(I)* +// pixel_iterator_base_<I, E>::address() +// { +// mln_precondition(exact(this)->is_valid()); +// return value_ptr_; +// } + +// template <typename I, typename E> +// const typename pixel_iterator_base_<I, E>::psite& +// pixel_iterator_base_<I, E>::site() const +// { +// //FIXME: update psite +// return p_; +// } + +// template <typename I, typename E> +// typename pixel_iterator_base_<I, E>::psite& +// pixel_iterator_base_<I, E>::site() +// { +// //FIXME: update psite +// return p_; +// } + - template <typename Image, typename Exact> - const typename pixel_iterator_base_<Image, Exact>::psite& - pixel_iterator_base_<Image, Exact>::site() const - { - //FIXME: update psite - return p_; - } - template <typename Image, typename Exact> - typename pixel_iterator_base_<Image, Exact>::psite& - pixel_iterator_base_<Image, Exact>::site() - { - //FIXME: update psite - return p_; - } // FIXME: Dead code +// ################# + + + // /*! \brief pixel_iterator_base_ class // */ -// template <typename Exact, typename Image> -// class pixel_iterator_base_ : public Pixel_Iterator<Exact> +// template <typename E, typename I> +// class pixel_iterator_base_ : public Pixel_Iterator<E> // { // public: -// /// Image pixel value type. -// typedef mln_value(Image) value; -// /// Image pixel rvalue type. -// typedef mln_value(Image)& rvalue; -// /// Image pixel lvalue type -// typedef mln_value(Image) lvalue; +// /// I pixel value type. +// typedef mln_value(I) value; +// /// I pixel rvalue type. +// typedef mln_value(I)& rvalue; +// /// I pixel lvalue type +// typedef mln_value(I) lvalue; // // Go to the beginning of the image. // void start(); @@ -193,46 +224,46 @@ // #ifndef MLN_INCLUDE_ONLY -// template <typename Exact, typename Image> -// pixel_iterator_base_<Exact, Image>::pixel_iterator_base_() +// template <typename E, typename I> +// pixel_iterator_base_<E,I>::pixel_iterator_base_() // { // } -// template <typename Exact, typename Image> -// void pixel_iterator_base_<Exact, Image>::start() +// template <typename E, typename I> +// void pixel_iterator_base_<E,I>::start() // { // current_ = start_; // } -// template <typename Exact, typename Image> -// void pixel_iterator_base_<Exact, Image>::next_() +// template <typename E, typename I> +// void pixel_iterator_base_<E,I>::next_() // { // ++current_; // } -// template <typename Exact, typename Image> -// void pixel_iterator_base_<Exact, Image>::invalidate() +// template <typename E, typename I> +// void pixel_iterator_base_<E,I>::invalidate() // { // current_ = eoi_; // } -// template <typename Exact, typename Image> -// bool pixel_iterator_base_<Exact, Image>::is_valid() const +// template <typename E, typename I> +// bool pixel_iterator_base_<E,I>::is_valid() const // { // return (current_ != eoi_); // } -// template <typename Exact, typename Image> -// typename pixel_iterator_base_<Exact, Image>::rvalue -// pixel_iterator_base_<Exact, Image>::operator*() +// template <typename E, typename I> +// typename pixel_iterator_base_<E,I>::rvalue +// pixel_iterator_base_<E,I>::operator*() // { // return *current_; // } -// template <typename Exact, typename Image> -// typename pixel_iterator_base_<Exact, Image>::lvalue -// pixel_iterator_base_<Exact, Image>::operator*() const +// template <typename E, typename I> +// typename pixel_iterator_base_<E,I>::lvalue +// pixel_iterator_base_<E,I>::operator*() const // { // return *current_; // } Index: mln/level/was.median.hh --- mln/level/was.median.hh (revision 1022) +++ mln/level/was.median.hh (working copy) @@ -30,7 +30,7 @@ /*! \file mln/level/was.median.hh * - * \brief Obsolete routine for the median filter. + * \brief Obsolete routines for median filtering. */ @@ -43,6 +43,9 @@ namespace impl { + + // prefer using a canvas + template <typename I, typename W, typename O> void median_as_procedure(const I& input, const W& win, @@ -126,6 +129,62 @@ } + + // horizontal median + + template <typename I, typename O> + void hmedian(const I& input, const hline2d& win, O& output) + { + + const int + max_row = input.max_row(), + min_col = input.min_col(), + max_col = input.max_col(); + const unsigned half = win.length() / 2; + + point2d p; + int& row = p.row(); + int& col = p.col(); + + accu::median_on<mln_value(I)> med; + + for (row = input.min_row(); row <= max_row; ++row) + { + int ct, cu; + + // initialization (before first point of the row) + med.init(); + for (ct = min_col; ct < min_col + half; ++ct) + med.take(input.at(row, ct)); + + // left columns (just take new points) + for (col = min_col; col <= min_col + half; ++col, ++ct) + { + med.take(input.at(row, ct)); + output(p) = med; + } + + // middle columns (both take and untake) + cu = min_col; + for (; col <= max_col - half; ++cu, ++col, ++ct) + { + med.take(input.at(row, ct)); + med.untake(input.at(row, cu)); + output(p) = med; + } + + // right columns (now just untake old points) + for (; col <= max_col; ++cu, ++col) + { + med.untake(input.at(row, cu)); + output(p) = med; + } + } + + } // end of hmedian + + + } // end of namespace mln::level::impl } // end of namespace mln::level Index: mln/level/fast_median.hh --- mln/level/fast_median.hh (revision 1022) +++ mln/level/fast_median.hh (working copy) @@ -103,7 +103,7 @@ med.init(); { mln_qixter(I) q(win, p, input); - for_all(q) //if (input.has(q)) + for_all(q) med.take(*q); } @@ -118,10 +118,10 @@ { // "go down" - for_all(q_top) //if (input.has(q_top)) + for_all(q_top) med.untake(*q_top); - for_all(q_bot) //if (input.has(q_bot)) + for_all(q_bot) med.take(*q_bot); output(p) = med; @@ -130,11 +130,10 @@ // browse line fwd while (col < max_col) { - std::cout << "debug: psite" << q_fm.psite() << "\n"; ++col; - for_all(q_fm) //if (input.has(q_fm)) + for_all(q_fm) med.untake(*q_fm); - for_all(q_fp) //if (input.has(q_fp)) + for_all(q_fp) med.take(*q_fp); output(p) = med; } @@ -143,9 +142,9 @@ while (col > min_col) { --col; - for_all(q_bm) //if (input.has(q_bm)) + for_all(q_bm) med.untake(*q_bm); - for_all(q_bp) //if (input.has(q_bp)) + for_all(q_bp) med.take(*q_bp); output(p) = med; }

1 0

1022: Change fast_iterator hierarchy.
by Nicolas Ballas 18 Jul '07

18 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Nicolas Ballas <ballas(a)lrde.epita.fr> Change fast_iterator hierarchy. * tests/fast_median.cc: New test. * tests/pixter_point2d.cc: Update test. * mln/core/pixel.hh: Add typedefs in the class pixel_. * mln/core/macros.hh: Add new macros. * mln/core/concept/pixel_iterator.hh: New, rename Fast_Iterator concep as Pixel Iterator. * mln/core/concept/fast_iterator.hh: Remove. * mln/core/concept/genpixel.hh: Constructor fixed. * mln/core/pixter2d_b.hh: Broken!. * mln/core/image2d_b.hh: Add new method point_at_offset. * mln/core/dpoints_pixter.hh: Adapt dpoint_pixter to the new hierarchy. * mln/core/internal/pixel_iterator_base.hh: New, rename fast_iterator_base_ as pixel_iterator_base_. * mln/core/internal/fast_iterator_base.hh: Remove. * mln/level/fast_median.hh: New median using qixter. mln/core/concept/genpixel.hh | 33 ++-- mln/core/concept/pixel_iterator.hh | 81 ++++++++++ mln/core/dpoints_pixter.hh | 35 +--- mln/core/image2d_b.hh | 22 ++ mln/core/internal/pixel_iterator_base.hh | 246 +++++++++++++++++++++++++++++++ mln/core/macros.hh | 9 + mln/core/pixel.hh | 2 mln/core/pixter2d_b.hh | 5 mln/level/fast_median.hh | 176 ++++++++++++++++++++++ tests/fast_median.cc | 63 +++++++ tests/pixter_point2d.cc | 4 11 files changed, 630 insertions(+), 46 deletions(-) Index: tests/fast_median.cc --- tests/fast_median.cc (revision 0) +++ tests/fast_median.cc (revision 0) @@ -0,0 +1,63 @@ +// Copyright (C) 2007 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. + +/*! \file tests/fast_median.cc + * + * \brief Test on mln::level::fast_median. + */ + +#include <mln/core/image2d_b.hh> +#include <mln/core/rectangle2d.hh> + +#include <mln/io/load_pgm.hh> +#include <mln/io/save_pgm.hh> + +#include <mln/value/int_u.hh> +#include <mln/level/fast_median.hh> +#include <mln/level/approx/median.hh> + + + + +int main() +{ + using namespace mln; + using value::int_u8; + + rectangle2d rect(51, 51); + border::thickness = 52; + + image2d_b<int_u8> + lena = io::load_pgm("../img/lena.pgm"), + out(lena.domain()); + + level::fast_median(lena, rect, out); + io::save_pgm(out, "out.pgm"); + +// level::approx::median(lena, rec, out); +// io::save_pgm(out, "outa.pgm"); +} Index: tests/pixter_point2d.cc --- tests/pixter_point2d.cc (revision 1021) +++ tests/pixter_point2d.cc (working copy) @@ -27,7 +27,6 @@ #include <cassert> #include <mln/core/image2d_b.hh> -#include <iostream> int @@ -35,7 +34,7 @@ { using namespace mln; - image2d_b<int> ima(2, 2); + image2d_b<int> ima(20, 20); image2d_b<int>::piter p(ima.domain()); image2d_b<int>::pixter f(ima); @@ -49,7 +48,6 @@ for_all(f) { - std::cout << *f << std::endl; assert(*f == i ++); *f = 5; } Index: mln/core/pixel.hh --- mln/core/pixel.hh (revision 1021) +++ mln/core/pixel.hh (working copy) @@ -64,6 +64,8 @@ { typedef mln_psite(I) psite; typedef mln_value(I) value; + typedef mln_rvalue(I) rvalue; + typedef I ima; pixel(I& ima); pixel(I& ima, const psite& p); Index: mln/core/macros.hh --- mln/core/macros.hh (revision 1021) +++ mln/core/macros.hh (working copy) @@ -96,6 +96,15 @@ /// Shortcut to access the fwd_viter type associated to T. # define mln_fwd_viter(T) typename T::fwd_viter +// i + +/// Shortcut to access the fwd_viter type associated to T. +# define mln_image(T) typename T::image + +/// Shortcut to access the fwd_viter type associated to T. +# define mln_ima(T) typename T::ima + + // l Index: mln/core/concept/pixel_iterator.hh --- mln/core/concept/pixel_iterator.hh (revision 0) +++ mln/core/concept/pixel_iterator.hh (revision 0) @@ -0,0 +1,81 @@ +// Copyright (C) 2007 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_CORE_CONCEPT_PIXEL_ITERATOR_HH +# define MLN_CORE_CONCEPT_PIXEL_ITERATOR_HH + +/*! \file mln/core/concept/fast_iterator.hh + * \brief Definition of the concept of mln::Pixel_Iterator. + */ + +# include <mln/core/concept/iterator.hh> +# include <mln/core/concept/genpixel.hh> + +namespace mln +{ + + /*! \brief Pixel Iterator concept class. + * + * \see FIXME + */ + template <typename E> + struct Pixel_Iterator : + public Iterator<E>, + public GenPixel<E> + { + /* + typedef rvalue + typedef lvalue + lvalue operator* (); + rvalue operator* () const; + */ + protected: + Pixel_Iterator(); + }; + +#ifndef MLN_INCLUDE_ONLY + + template <typename E> + Pixel_Iterator<E>::Pixel_Iterator() + { + typedef mln_value(E) value; + typedef mln_rvalue(E) rvalue; + typedef mln_lvalue(E) lvalue; + + lvalue (E::*m1)() = & E::operator*; + m1 = 0; + rvalue (E::*m2)() const = & E::operator*; + m2 = 0; + } + +#endif // ! MLN_INCLUDE_ONLY + + +} // end of namespace mln + + +#endif // ! MLN_CORE_CONCEPT_PIXEL_ITERATOR_HH Index: mln/core/concept/genpixel.hh --- mln/core/concept/genpixel.hh (revision 1021) +++ mln/core/concept/genpixel.hh (working copy) @@ -60,12 +60,13 @@ { /* - typedef image; // not const! - typedef psite; + typedef ima; // not const! + typede psite; typedef rvalue; + typedef value; - const image& image() const; - const psite& psite() const; // FIXME ou cpy + const ima& image() const; + const psite& site() const; // FIXME ou cpy rvalue operator*() const; @@ -98,17 +99,19 @@ template <typename E> GenPixel<E>::GenPixel() { - // FIXME -// int dim = E::dim; -// mln_invariant(dim > 0); -// dim = 0; -// typedef mln_point(E) point; -// typedef mln_dpoint(E) dpoint; -// typedef mln_coord(E) coord; -// const point* (E::*m1)() const = & E::pointer; -// m1 = 0; -// coord (E::*m2)(unsigned i) const = & E::operator[]; -// m2 = 0; + typedef mln_value(E) value; + typedef mln_rvalue(E) rvalue; + typedef mln_ima(E) ima; + typedef mln_psite(E) psite; + + const ima& (E::*m1)() const = & E::image; + m1 = 0; + const psite& (E::*m2)() const = & E::site; + m2 = 0; + rvalue (E::*m3)() const = & E::operator*; + m3 = 0; + const value* (E::*m4)() const = & E::address; + m4 = 0; } Index: mln/core/pixter2d_b.hh --- mln/core/pixter2d_b.hh (revision 1021) +++ mln/core/pixter2d_b.hh (working copy) @@ -28,7 +28,7 @@ #ifndef MLN_CORE_PIXTER2D_B_HH # define MLN_CORE_PIXTER2D_B_HH -# include <mln/core/internal/fast_iterator_base.hh> +# include <mln/core/internal/pixel_iterator_base.hh> # include <mln/core/point2d.hh> # include <iostream> @@ -39,7 +39,6 @@ * \brief Pixel iterator class on a image 2d with border. */ -# include <mln/core/concept/fast_iterator.hh> namespace mln @@ -50,7 +49,7 @@ template <typename T> class fwd_pixter2d_b : - public internal::fast_iterator_base_< fwd_pixter2d_b<T>, image2d_b<T> > + public internal::pixel_iterator_base_< fwd_pixter2d_b<T>, image2d_b<T> > { public: /// Image pixel value type. Index: mln/core/image2d_b.hh --- mln/core/image2d_b.hh (revision 1021) +++ mln/core/image2d_b.hh (working copy) @@ -151,9 +151,17 @@ T& at(int row, int col); + /// Fast Image method + /// Return the offset corresponding to the dpoint \p dp. int offset(const dpoint2d& dp); + /// Return the offset corresponding to the dpoint \p p. + int offset(const point2d& p); + + /// Return the point corresponding to the offset \p offset. + point2d point_at_offset(int offset); + private: @@ -387,6 +395,20 @@ return dp[0] * vb_.len(1) + dp[1]; } + template <typename T> + int + image2d_b<T>::offset(const point2d& p) + { + return p[0] * vb_.len(1) + p[1]; + } + + template <typename T> + point2d + image2d_b<T>::point_at_offset(int offset) + { + return point2d(offset / (vb_.len(1) + 1), offset % (vb_.len(1) + 1)); + } + # endif // ! MLN_INCLUDE_ONLY } // end of namespace mln Index: mln/core/dpoints_pixter.hh --- mln/core/dpoints_pixter.hh (revision 1021) +++ mln/core/dpoints_pixter.hh (working copy) @@ -30,7 +30,7 @@ # include <cassert> # include <vector> -# include <mln/core/concept/fast_iterator.hh> +# include <mln/core/internal/pixel_iterator_base.hh> namespace mln @@ -41,8 +41,9 @@ * The parameter \c I is the image type. */ template <typename I> - class dpoints_pixter : public Fast_Iterator< dpoints_pixter > + class dpoints_pixter : public internal::pixel_iterator_base_< I, dpoints_pixter > { + typedef typename internal::pixel_iterator_base_< I, dpoints_pixter > super; public: /// Image pixel value typedef mln_value(I) value; @@ -51,14 +52,14 @@ /// Image pixel lvalue typedef mln_lvalue(I) lvalue; - /*! \brief Constructo. + /*! \brief Constructor. * * \param[in] dps Object that can provide an array of delta-points. * \param[in] p_ref Center point to iterate around. * \param[in] ima Image to iterate. */ template <typename Dps> - dpoints_pixter(const Dps& dps, const typename I::point& p_ref,I& ima); + dpoints_pixter(const Dps& dps, const typename I::psite& p_ref,I& ima); /// Set the iterator at the start. void start(); @@ -69,10 +70,6 @@ /// Is the iterator valid? bool is_valid() const; - /// Get the iterator value - rvalue operator* (); - lvalue operator* () const; - private: /// offset of each dpoints std::vector<int> offset_; @@ -87,8 +84,9 @@ template <typename I> template <typename Dps> - dpoints_pixter::dpoints_pixter(const Dps& dps, const typename I::point& ref, I& ima) + dpoints_pixter::dpoints_pixter(const Dps& dps, const typename I::psite& ref, I& ima) : super(ima) { + this->p_ = ref; pixref_ = &ima(ref); for (typename std::vector<typename I::dpoint>::const_iterator it = dps.vec().begin(); it != dps.vec().end(); @@ -100,13 +98,16 @@ template <typename I> void dpoints_pixter::start() { + pixref_ = &ima_(this->p_); i_ = offset_.begin(); + this->current_value_ = pixref_ + *i_; } template <typename I> void dpoints_pixter::next_() { ++i_; + this->current_value_ = pixref_ + *i_; } template <typename I> @@ -121,22 +122,6 @@ i_ = offset_.end(); } - template <typename I> - typename dpoints_pixter::rvalue - dpoints_pixter::operator* () - { - assert(this->is_valid()); - return *(pixref_ + *i_); - } - - template <typename I> - typename dpoints_pixter::lvalue - dpoints_pixter::operator* () const - { - assert(this->is_valid()); - return *(pixref_ + *i_); - } - #endif // ! MLN_INCLUDE_ONLY } Index: mln/core/internal/pixel_iterator_base.hh --- mln/core/internal/pixel_iterator_base.hh (revision 0) +++ mln/core/internal/pixel_iterator_base.hh (revision 0) @@ -0,0 +1,246 @@ +// Copyright (C) 2007 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_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH +# define MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH + +/*! \file mln/core/internal/fast_iterator_base.hh + * + * \brief Base class for Fast_Iterator concept implementation. + * Usefull for code factorisation + */ + +# include <mln/core/concept/pixel_iterator.hh> + + +namespace mln +{ + + namespace internal + { + + /*! \brief pixel_iterator_base_ class + */ + template <typename Image, typename Exact> + class pixel_iterator_base_ : public Pixel_Iterator<Exact> + { + public: + /// Image type. + typedef Image ima; + /// Image value type. + typedef mln_value(Image) value; + /// Image lvalue type. + typedef mln_lvalue(Image) lvalue; + /// Image rvalue type. + typedef mln_rvalue(Image) rvalue; + /// Image psite type. + typedef mln_psite(Image) psite; + + /// Get the image associated to the current pixel iterator. + const ima& image() const; + /// pixel iterator value. + lvalue operator* (); + /// Get the pixel iterator value. + rvalue operator* () const; + /// Address of the current iterator value/pixel. + const value* address() const; + /// Address of the current iterator value/pixel. + value* address(); + + /// psite associated to the iterator. + const psite& site() const; + /// psite associated to the iterator. + psite& site(); + + protected: + /// Current pixel value + value* current_value_; + /// Image associated to the iterator + ima& ima_; + /// Psite of the pixel + psite p_; + + pixel_iterator_base_(ima& source); + }; + +#ifndef MLN_INCLUDE_ONLY + + template <typename Image, typename Exact> + pixel_iterator_base_<Image, Exact>::pixel_iterator_base_(ima& source) : + ima_(source) + { + } + + template <typename Image, typename Exact> + const typename pixel_iterator_base_<Image, Exact>::ima& + pixel_iterator_base_<Image, Exact>::image() const + { + return ima_; + } + + template <typename Image, typename Exact> + typename pixel_iterator_base_<Image, Exact>::lvalue + pixel_iterator_base_<Image, Exact>::operator* () + { + assert(exact(this)->is_valid()); + return *current_value_; + } + + template <typename Image, typename Exact> + typename pixel_iterator_base_<Image, Exact>::rvalue + pixel_iterator_base_<Image, Exact>::operator* () const + { + assert(exact(this)->is_valid()); + return *current_value_; + } + + template <typename Image, typename Exact> + const typename pixel_iterator_base_<Image, Exact>::value* + pixel_iterator_base_<Image, Exact>::address() const + { + assert(exact(this)->is_valid()); + return current_value_; + } + + template <typename Image, typename Exact> + typename pixel_iterator_base_<Image, Exact>::value* + pixel_iterator_base_<Image, Exact>::address() + { + assert(exact(this)->is_valid()); + return current_value_; + } + + template <typename Image, typename Exact> + const typename pixel_iterator_base_<Image, Exact>::psite& + pixel_iterator_base_<Image, Exact>::site() const + { + //FIXME: update psite + return p_; + } + + template <typename Image, typename Exact> + typename pixel_iterator_base_<Image, Exact>::psite& + pixel_iterator_base_<Image, Exact>::site() + { + //FIXME: update psite + return p_; + } + + +// FIXME: Dead code +// /*! \brief pixel_iterator_base_ class +// */ +// template <typename Exact, typename Image> +// class pixel_iterator_base_ : public Pixel_Iterator<Exact> +// { +// public: +// /// Image pixel value type. +// typedef mln_value(Image) value; +// /// Image pixel rvalue type. +// typedef mln_value(Image)& rvalue; +// /// Image pixel lvalue type +// typedef mln_value(Image) lvalue; + +// // Go to the beginning of the image. +// void start(); +// // Go on the next element. +// void next_(); +// // Invalidate the iterator. +// void invalidate(); +// // Is the iterator referencing a correct position in the image? +// bool is_valid() const; + + +// // Return the current pixel. +// rvalue operator* (); +// lvalue operator* () const; + +// protected: +// // beginning of the image +// value* start_; +// // end of the image +// value* eoi_; +// // current position in the image +// value* current_; +// // End of Factoriasable + +// pixel_iterator_base_(); +// }; + +// #ifndef MLN_INCLUDE_ONLY + +// template <typename Exact, typename Image> +// pixel_iterator_base_<Exact, Image>::pixel_iterator_base_() +// { +// } + +// template <typename Exact, typename Image> +// void pixel_iterator_base_<Exact, Image>::start() +// { +// current_ = start_; +// } + +// template <typename Exact, typename Image> +// void pixel_iterator_base_<Exact, Image>::next_() +// { +// ++current_; +// } + +// template <typename Exact, typename Image> +// void pixel_iterator_base_<Exact, Image>::invalidate() +// { +// current_ = eoi_; +// } + +// template <typename Exact, typename Image> +// bool pixel_iterator_base_<Exact, Image>::is_valid() const +// { +// return (current_ != eoi_); +// } + +// template <typename Exact, typename Image> +// typename pixel_iterator_base_<Exact, Image>::rvalue +// pixel_iterator_base_<Exact, Image>::operator*() +// { +// return *current_; +// } + + +// template <typename Exact, typename Image> +// typename pixel_iterator_base_<Exact, Image>::lvalue +// pixel_iterator_base_<Exact, Image>::operator*() const +// { +// return *current_; +// } + +#endif // ! MLN_INCLUDE_ONLY + } // end of namespace internal + +} // end of namespace mln + + +#endif // ! MLN_CORE_INTERNAL_PIXEL_ITERATOR_BASE_HH Index: mln/level/fast_median.hh --- mln/level/fast_median.hh (revision 0) +++ mln/level/fast_median.hh (revision 0) @@ -0,0 +1,176 @@ +// Copyright (C) 2007 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_LEVEL_FAST_MEDIAN_HH +# define MLN_LEVEL_FAST_MEDIAN_HH + +/*! \file mln/level/median.hh + * + * \brief Fast Median filtering of an image. + */ + +# include <mln/core/concept/image.hh> +# include <mln/core/window2d.hh> +# include <mln/accu/median.hh> +# include <ostream> + + +namespace mln +{ + + namespace level + { + + /*! Compute in \p output the median filter of image \p input by + * the window \p win. + * + * \param[in] input The image to be filtered. + * \param[in] win The window. + * \param[in,out] output The output image. + * + * \pre \p input and \p output have to be initialized. + */ + template <typename I, typename W, typename O> + void fast_median(Image& input, const Window<W>& win, + Image<O>& output); + + +# ifndef MLN_INCLUDE_ONLY + + + namespace impl + { + + template <typename I, typename W, typename O> + void fast_median(I& input, + const W& win, + O& output) + { + mln_precondition(input.has_data()); + mln_precondition(output.has_data()); + + int + min_row = input.min_row(), max_row = input.max_row(), + min_col = input.min_col(), max_col = input.max_col(); + + window2d + win_fwd_plus = win - (win + left), + win_fwd_minus = (win + left) - win, + win_bkd_plus = win - (win + right), + win_bkd_minus = (win + right) - win, + win_bot = win - (win + up), + win_top = (win + up) - win; + + + + //Fixme bug + point2d p; + + accu::median_on<mln_value(I)> med; + + // initialization + + p = input.domain().pmin() + up; + mln_qixter(I) + q_fp(win_fwd_plus, p, input), q_fm(win_fwd_minus, p, input), + q_bp(win_bkd_plus, p, input), q_bm(win_bkd_minus, p, input), + q_top(win_top, p, input), q_bot(win_bot, p, input); + + med.init(); + { + mln_qixter(I) q(win, p, input); + for_all(q) //if (input.has(q)) + med.take(*q); + } + + int& row = p.row(); + int& col = p.col(); + bool fwd = true; + + mln_assertion(p.col() == min_col); + mln_assertion(p.row() == min_row - 1); + + for (row = min_row; row <= max_row; ++row) + { + + // "go down" + for_all(q_top) //if (input.has(q_top)) + med.untake(*q_top); + + for_all(q_bot) //if (input.has(q_bot)) + med.take(*q_bot); + + output(p) = med; + + if (fwd) + // browse line fwd + while (col < max_col) + { + std::cout << "debug: psite" << q_fm.psite() << "\n"; + ++col; + for_all(q_fm) //if (input.has(q_fm)) + med.untake(*q_fm); + for_all(q_fp) //if (input.has(q_fp)) + med.take(*q_fp); + output(p) = med; + } + else + // browse line bkd + while (col > min_col) + { + --col; + for_all(q_bm) //if (input.has(q_bm)) + med.untake(*q_bm); + for_all(q_bp) //if (input.has(q_bp)) + med.take(*q_bp); + output(p) = med; + } + // change browsing + fwd = ! fwd; + } + } + + } // end of namespace mln::level::impl + + + // facade + + template <typename I, typename W, typename O> + void fast_median(Image& input, const Window<W>& win, + Image<O>& output) + { + impl::fast_median(exact(input), exact(win), exact(output)); + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::level + +} // end of namespace mln + + +#endif // ! MLN_LEVEL_FAST_MEDIAN_HH

1 0

1021: New 2D window types in mln.
by Thierry Geraud 17 Jul '07

17 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/milena Index: ChangeLog from Thierry Geraud <thierry.geraud(a)lrde.epita.fr> New 2D window types in mln. * tests/median.cc: Clean up. * tests/hmedian.cc: New. * tests/erosion.cc: Clean up. * tests/rectangle2d.cc: Update. * mln/core/hline2d.hh: New. * mln/core/rectangle2d.hh (half_width_, half_height_): Change to... (width_, height_): ...these. Update. * mln/core/vline2d.hh: New. * mln/morpho/erosion.hh (erosion_wrt_win): New overload. * mln/level/was.median.hh: New. * mln/level/approx: New. * mln/level/approx/median.hh: New. * mln/level/median.hh (median_as_procedure): Move into... * mln/level/was.hmedian.hh: ...this new file. mln/core/hline2d.hh | 176 ++++++++++++++++++++++++++++++++++++++++++++ mln/core/rectangle2d.hh | 48 ++++++------ mln/core/vline2d.hh | 178 +++++++++++++++++++++++++++++++++++++++++++++ mln/level/approx/median.hh | 90 ++++++++++++++++++++++ mln/level/median.hh | 157 +++++++++++++++++---------------------- mln/level/was.hmedian.hh | 102 +++++++++++++++++++++++++ mln/level/was.median.hh | 136 ++++++++++++++++++++++++++++++++++ mln/morpho/erosion.hh | 31 ++++++- tests/erosion.cc | 8 +- tests/hmedian.cc | 64 ++++++++++++++++ tests/median.cc | 15 ++- tests/rectangle2d.cc | 4 - 12 files changed, 883 insertions(+), 126 deletions(-) Index: tests/median.cc --- tests/median.cc (revision 1020) +++ tests/median.cc (working copy) @@ -38,21 +38,26 @@ #include <mln/value/int_u.hh> #include <mln/level/median.hh> +#include <mln/level/approx/median.hh> -using namespace mln; -using namespace mln::value; int main() { - rectangle2d rec(64, 64); - border::thickness = 66; + using namespace mln; + using value::int_u8; + + rectangle2d rect(51, 51); + border::thickness = 52; image2d_b<int_u8> lena = io::load_pgm("../img/lena.pgm"), out(lena.domain()); - level::median(lena, rec, out); + level::median(lena, rect, out); io::save_pgm(out, "out.pgm"); + +// level::approx::median(lena, rec, out); +// io::save_pgm(out, "outa.pgm"); } Index: tests/hmedian.cc --- tests/hmedian.cc (revision 0) +++ tests/hmedian.cc (revision 0) @@ -0,0 +1,64 @@ +// Copyright (C) 2007 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. + +/*! \file tests/hmedian.cc + * + * \brief Test on the hline2d version of mln::level::median. + */ + +#include <mln/core/image2d_b.hh> +#include <mln/core/rectangle2d.hh> + +#include <mln/io/load_pgm.hh> +#include <mln/io/save_pgm.hh> + +#include <mln/value/int_u.hh> +#include <mln/level/median.hh> +#include <mln/level/compare.hh> + + + + +int main() +{ + using namespace mln; + using value::int_u8; + + border::thickness = 0; + + image2d_b<int_u8> + lena = io::load_pgm("../img/lena.pgm"), + out(lena.domain()), + ref(lena.domain()); + + level::median(lena, rectangle2d(1, 101), ref); + + level::median(lena, hline2d(101), out); + io::save_pgm(out, "out.pgm"); + + // FIXME: mln_assertion(out = ref); +} Index: tests/erosion.cc --- tests/erosion.cc (revision 1020) +++ tests/erosion.cc (working copy) @@ -45,13 +45,13 @@ #include <mln/fun/ops.hh> -using namespace mln; -using namespace mln::value; - int main() { - rectangle2d rec(4, 4);//64, 64); + using namespace mln; + using value::int_u8; + + rectangle2d rec(21, 21); border::thickness = 66; image2d_b<int_u8> Index: tests/rectangle2d.cc --- tests/rectangle2d.cc (revision 1020) +++ tests/rectangle2d.cc (working copy) @@ -38,12 +38,12 @@ { using namespace mln; - const unsigned h = 3, w = 4; + const unsigned h = 3, w = 5; rectangle2d rec(h, w); mln_assertion(rec.is_centered()); mln_assertion(rec.is_symmetric()); mln_assertion(rec = -rec); - mln_assertion(rec.nelements() = (2*h+1) * (2*w+1)); + mln_assertion(rec.nelements() = h * w); } Index: mln/core/hline2d.hh --- mln/core/hline2d.hh (revision 0) +++ mln/core/hline2d.hh (revision 0) @@ -0,0 +1,176 @@ +// Copyright (C) 2007 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_CORE_HLINE2D_HH +# define MLN_CORE_HLINE2D_HH + +/*! \file mln/core/hline2d.hh + * + * \brief Definition of the mln::hline2d window. + */ + +# include <mln/core/concept/window.hh> +# include <mln/core/internal/set_of.hh> +# include <mln/core/dpoint2d.hh> +# include <mln/core/dpoints_piter.hh> + + +namespace mln +{ + + /*! \brief Horizontal line window defined on the 2D square grid. + * + * An hline2d is centered and symmetrical; so its height is 1 and + * its width (length) is odd. + * + * For instance: \n + * o o x o o \n + * is defined with length = 5. + */ + struct hline2d : public Window< hline2d >, + public internal::set_of_<dpoint2d> + { + /// Point associated type. + typedef point2d point; + + /// Dpoint associated type. + typedef dpoint2d dpoint; + + /*! \brief Piter type to browse a hline such as: "for each row + * (increasing), for each column (increasing)." + */ + typedef dpoints_fwd_piter<dpoint2d> fwd_qiter; + + /*! \brief Piter type to browse a hline such as: "for each row + * (decreasing), for each column (decreasing)." + */ + typedef dpoints_fwd_piter<dpoint2d> bkd_qiter; // FIXME: bkd!!! + + /*! \brief Same as fwd_qiter. + */ + typedef fwd_qiter qiter; + + /*! \brief Constructor. + * + * \param[in] length Length, thus width, of the horizontal line. + * + * \pre \p length is odd. + */ + hline2d(unsigned length); + + /*! \brief Test if the window is centered. + * + * \return True. + */ + bool is_centered() const; + + /*! \brief Test if the window is symmetric. + * + * \return true. + */ + bool is_symmetric() const; + + /*! \brief Give the hline length, that is, its width. + */ + unsigned length() const; + + /*! \brief Give the maximum coordinate gap between the window + * center and a window point. + */ + unsigned delta() const; + + /// Get the symmetrical window. + hline2d sym_() const; + + protected: + unsigned length_; + }; + + + /*! \brief Print an horizontal 2D line window \p win into the output + * stream \p ostr. + * + * \param[in,out] ostr An output stream. + * \param[in] win An horizontal 2D line window. + * + * \return The modified output stream \p ostr. + * + * \relates mln::hline2d + */ + std::ostream& operator<<(std::ostream& ostr, const hline2d& win); + + + +# ifndef MLN_INCLUDE_ONLY + + hline2d::hline2d(unsigned length) + : length_(length) + { + mln_precondition(length % 2 = 1); + const int dcol = length / 2; + for (int col = - dcol; col <= dcol; ++col) + insert(make::dpoint2d(0, col)); + } + + bool hline2d::is_centered() const + { + return true; + } + + bool hline2d::is_symmetric() const + { + return true; + } + + unsigned hline2d::length() const + { + return length_; + } + + unsigned hline2d::delta() const + { + return length_ / 2; + } + + hline2d hline2d::sym_() const + { + return *this; + } + + std::ostream& operator<<(std::ostream& ostr, const hline2d& win) + { + ostr << "[line2d: length=" << win.length() << ']'; + return ostr; + } + +# endif // ! MLN_INCLUDE_ONLY + +} // end of namespace mln + + + +#endif // ! MLN_CORE_HLINE2D_HH Index: mln/core/rectangle2d.hh --- mln/core/rectangle2d.hh (revision 1020) +++ mln/core/rectangle2d.hh (working copy) @@ -47,13 +47,11 @@ * A rectangle2d is a 2D window with rectangular shape. It is * centered and symmetrical. * - * For instance: - * o o o o o - * o o x o o - * o o o o o - * is defined with half_height = 1 and half_width = 2. - * - * \todo Consider width instead of half_width (same for height). + * For instance: \n + * o o o o o \n + * o o x o o \n + * o o o o o \n + * is defined with height = 3 and width = 5. */ struct rectangle2d : public Window< rectangle2d >, public internal::set_of_<dpoint2d> @@ -80,12 +78,12 @@ /*! \brief Constructor. * - * \param[in] half_height sic - * \param[in] half_width sic + * \param[in] height sic + * \param[in] width sic * - * \pre half_height != 0 and half_width != 0 + * \pre Height and width are odd. */ - rectangle2d(unsigned half_height, unsigned half_width); + rectangle2d(unsigned height, unsigned width); /*! \brief Test if the window is centered. * @@ -116,7 +114,7 @@ rectangle2d sym_() const; protected: - unsigned half_height_, half_width_; + unsigned height_, width_; }; @@ -136,12 +134,12 @@ # ifndef MLN_INCLUDE_ONLY - rectangle2d::rectangle2d(unsigned half_height, unsigned half_width) - : half_height_(half_height), - half_width_(half_width) + rectangle2d::rectangle2d(unsigned height, unsigned width) + : height_(height), + width_(width) { - mln_precondition(half_height != 0 && half_width != 0); - const int drow = half_height, dcol = half_width; + mln_precondition(height % 2 = 1 && width % 2 = 1); + const int drow = height / 2, dcol = width / 2; for (int row = - drow; row <= drow; ++row) for (int col = - dcol; col <= dcol; ++col) insert(make::dpoint2d(row, col)); @@ -159,17 +157,17 @@ unsigned rectangle2d::height() const { - return 2 * half_height_ + 1; + return height_; } unsigned rectangle2d::width() const { - return 2 * half_width_ + 1; + return width_; } unsigned rectangle2d::delta() const { - return half_width_ > half_height_ ? half_width_ : half_height_; + return width_ > height_ ? width_ / 2 : height_ / 2; } rectangle2d rectangle2d::sym_() const @@ -177,10 +175,9 @@ return *this; } - std::ostream& operator<<(std::ostream& ostr, - const rectangle2d& win) + std::ostream& operator<<(std::ostream& ostr, const rectangle2d& win) { - ostr << "[width=" << win.width() << ", height=" << win.height() << ']'; + ostr << "[rectangle2d: width=" << win.width() << ", height=" << win.height() << ']'; return ostr; } @@ -190,4 +187,9 @@ +// when rectangle2d is involved, one surely also wants: +# include <mln/core/hline2d.hh> +# include <mln/core/vline2d.hh> + + #endif // ! MLN_CORE_RECTANGLE2D_HH Index: mln/core/vline2d.hh --- mln/core/vline2d.hh (revision 0) +++ mln/core/vline2d.hh (revision 0) @@ -0,0 +1,178 @@ +// Copyright (C) 2007 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_CORE_VLINE2D_HH +# define MLN_CORE_VLINE2D_HH + +/*! \file mln/core/vline2d.hh + * + * \brief Definition of the mln::vline2d window. + */ + +# include <mln/core/concept/window.hh> +# include <mln/core/internal/set_of.hh> +# include <mln/core/dpoint2d.hh> +# include <mln/core/dpoints_piter.hh> + + +namespace mln +{ + + /*! \brief Vertical line window defined on the 2D square grid. + * + * An vline2d is centered and symmetrical; so its width is 1 and + * its height (length) is odd. + * + * For instance: \n + * o \n + * x \n + * o \n + * is defined with length = 5. + */ + struct vline2d : public Window< vline2d >, + public internal::set_of_<dpoint2d> + { + /// Point associated type. + typedef point2d point; + + /// Dpoint associated type. + typedef dpoint2d dpoint; + + /*! \brief Piter type to browse a vline such as: "for each row + * (increasing), for each column (increasing)." + */ + typedef dpoints_fwd_piter<dpoint2d> fwd_qiter; + + /*! \brief Piter type to browse a vline such as: "for each row + * (decreasing), for each column (decreasing)." + */ + typedef dpoints_fwd_piter<dpoint2d> bkd_qiter; // FIXME: bkd!!! + + /*! \brief Same as fwd_qiter. + */ + typedef fwd_qiter qiter; + + /*! \brief Constructor. + * + * \param[in] length Length, thus height, of the vertical line. + * + * \pre \p length is odd. + */ + vline2d(unsigned length); + + /*! \brief Test if the window is centered. + * + * \return True. + */ + bool is_centered() const; + + /*! \brief Test if the window is symmetric. + * + * \return true. + */ + bool is_symmetric() const; + + /*! \brief Give the vline length, that is, its height. + */ + unsigned length() const; + + /*! \brief Give the maximum coordinate gap between the window + * center and a window point. + */ + unsigned delta() const; + + /// Get the symmetrical window. + vline2d sym_() const; + + protected: + unsigned length_; + }; + + + /*! \brief Print a vertical 2D line window \p win into the output + * stream \p ostr. + * + * \param[in,out] ostr An output stream. + * \param[in] win A vertical 2D line window. + * + * \return The modified output stream \p ostr. + * + * \relates mln::vline2d + */ + std::ostream& operator<<(std::ostream& ostr, const vline2d& win); + + + +# ifndef MLN_INCLUDE_ONLY + + vline2d::vline2d(unsigned length) + : length_(length) + { + mln_precondition(length % 2 = 1); + const int drow = length / 2; + for (int row = - drow; row <= drow; ++row) + insert(make::dpoint2d(row, 0)); + } + + bool vline2d::is_centered() const + { + return true; + } + + bool vline2d::is_symmetric() const + { + return true; + } + + unsigned vline2d::length() const + { + return length_; + } + + unsigned vline2d::delta() const + { + return length_ / 2; + } + + vline2d vline2d::sym_() const + { + return *this; + } + + std::ostream& operator<<(std::ostream& ostr, const vline2d& win) + { + ostr << "[line2d: length=" << win.length() << ']'; + return ostr; + } + +# endif // ! MLN_INCLUDE_ONLY + +} // end of namespace mln + + + +#endif // ! MLN_CORE_VLINE2D_HH Index: mln/morpho/erosion.hh --- mln/morpho/erosion.hh (revision 1020) +++ mln/morpho/erosion.hh (working copy) @@ -106,6 +106,7 @@ // FIXME: stage 3: dispatch w.r.t. fast property + // stage 2: dispatch w.r.t. the value kind template <typename I, typename W, typename O> @@ -123,19 +124,34 @@ } + // stage 1: dispatch w.r.t. the window type + // | + // V - template <typename I, typename W, typename O> // general case + template <typename I, typename W, typename O> void erosion_wrt_win(const Image& input, const Window<W>& win, Image<O>& output) { erosion_wrt_value(mln_kind(I)(), exact(input), exact(win), output); + // | + // --> call stage 2: dispatch w.r.t. the value kind } -// template <typename I, typename O> // rectangle2d -// void erosion_wrt_win(const Image& input, const rectangle2d& win, Image<O>& output) -// { -// return FIXME; -// } +# ifdef MLN_CORE_RECTANGLE2D_HH + + template <typename I, typename O> + void erosion_wrt_win(const Image& input, const rectangle2d& win, Image<O>& output) + { + O tmp(exact(output).domain()); + morpho::erosion(input, hline2d(win.width()), tmp); + morpho::erosion(tmp, vline2d(win.height()), output); + } + +# endif // MLN_CORE_RECTANGLE2D_HH + + // ^ + // | + // end of stage1 (dispatch w.r.t. the window type) } // end of namespace mln::morpho::impl @@ -146,6 +162,9 @@ template <typename I, typename W, typename O> void erosion(const Image& input, const Window<W>& win, Image<O>& output) { + mln_precondition(exact(output).domain() = exact(input).domain()); + mln_precondition(! exact(win).is_empty()); + impl::erosion_wrt_win(input, exact(win), output); if (exact(win).is_centered()) mln_postcondition(output <= input); Index: mln/level/was.median.hh --- mln/level/was.median.hh (revision 0) +++ mln/level/was.median.hh (revision 0) @@ -0,0 +1,136 @@ +// Copyright (C) 2007 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_LEVEL_WAS_MEDIAN_HH +# define MLN_LEVEL_WAS_MEDIAN_HH + +/*! \file mln/level/was.median.hh + * + * \brief Obsolete routine for the median filter. + */ + + +namespace mln +{ + + namespace level + { + + namespace impl + { + + template <typename I, typename W, typename O> + void median_as_procedure(const I& input, + const W& win, + O& output) + { + mln_precondition(input.has_data()); + mln_precondition(output.has_data()); + + int + min_row = input.min_row(), max_row = input.max_row(), + min_col = input.min_col(), max_col = input.max_col(); + + window2d + win_fwd_plus = win - (win + left), + win_fwd_minus = (win + left) - win, + win_bkd_plus = win - (win + right), + win_bkd_minus = (win + right) - win, + win_bot = win - (win + up), + win_top = (win + up) - win; + + point2d p; + mln_qiter(W) + q_fp(win_fwd_plus, p), q_fm(win_fwd_minus, p), + q_bp(win_bkd_plus, p), q_bm(win_bkd_minus, p), + q_top(win_top, p), q_bot(win_bot, p); + + accu::median_on<mln_value(I)> med; + + // initialization + + p = input.domain().pmin() + up; + med.init(); + { + mln_qiter(W) q(win, p); + for_all(q) if (input.has(q)) + med.take(input(q)); + } + + int& row = p.row(); + int& col = p.col(); + bool fwd = true; + + mln_assertion(p.col() = min_col); + mln_assertion(p.row() = min_row - 1); + + for (row = min_row; row <= max_row; ++row) + { + // "go down" + for_all(q_top) if (input.has(q_top)) + med.untake(input(q_top)); + for_all(q_bot) if (input.has(q_bot)) + med.take(input(q_bot)); + output(p) = med; + + if (fwd) + // browse line fwd + while (col < max_col) + { + ++col; + for_all(q_fm) if (input.has(q_fm)) + med.untake(input(q_fm)); + for_all(q_fp) if (input.has(q_fp)) + med.take(input(q_fp)); + output(p) = med; + } + else + // browse line bkd + while (col > min_col) + { + --col; + for_all(q_bm) if (input.has(q_bm)) + med.untake(input(q_bm)); + for_all(q_bp) if (input.has(q_bp)) + med.take(input(q_bp)); + output(p) = med; + } + + // change browsing + fwd = ! fwd; + } + } + + + } // end of namespace mln::level::impl + + } // end of namespace mln::level + +} // end of namespace mln + + +#endif // ! MLN_LEVEL_WAS_MEDIAN_HH Index: mln/level/approx/median.hh --- mln/level/approx/median.hh (revision 0) +++ mln/level/approx/median.hh (revision 0) @@ -0,0 +1,90 @@ +// Copyright (C) 2007 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_LEVEL_APPROX_MEDIAN_HH +# define MLN_LEVEL_APPROX_MEDIAN_HH + +/*! \file mln/level/median.hh + * + * \brief Approximates of some median filters of an image. + */ + +# include <mln/level/median.hh> +# include <mln/core/rectangle2d.hh> + + +namespace mln +{ + + namespace level + { + + namespace approx + { + + /*! Compute in \p output an approximate of the median filter of + * image \p input by the 2D rectangle \p win. + * + * \param[in] input The image to be filtered. + * \param[in] win The rectangle. + * \param[in,out] output The output image. + * + * The approximation is based on a vertical median ran after + * an horizontal median. + * + * \pre \p input and \p output have to be initialized. + */ + template <typename I, typename O> + void median(const Image& input, const rectangle2d& win, + Image<O>& output); + + +# ifndef MLN_INCLUDE_ONLY + + template <typename I, typename O> + void median(const Image& input_, const rectangle2d& win, + Image<O>& output_) + { + const I& input = exact(input_); + O& output = exact(output_); + mln_assertion(output.domain() = input.domain()); + + O tmp(output.domain()); + level::median(input, hline2d(win.width()), tmp); + level::median(tmp, vline2d(win.height()), output); + } + +# endif // ! MLN_INCLUDE_ONLY + + } // end of namespace mln::level::approx + + } // end of namespace mln::level + +} // end of namespace mln + + +#endif // ! MLN_LEVEL_APPROX_MEDIAN_HH Index: mln/level/median.hh --- mln/level/median.hh (revision 1020) +++ mln/level/median.hh (working copy) @@ -34,9 +34,11 @@ */ # include <mln/core/concept/image.hh> + # include <mln/core/window2d.hh> -# include <mln/accu/median.hh> +# include <mln/core/hline2d.hh> +# include <mln/accu/median.hh> # include <mln/canvas/sbrowsing.hh> @@ -66,88 +68,6 @@ namespace impl { - template <typename I, typename W, typename O> - void median_as_procedure(const I& input, - const W& win, - O& output) - { - mln_precondition(input.has_data()); - mln_precondition(output.has_data()); - - int - min_row = input.min_row(), max_row = input.max_row(), - min_col = input.min_col(), max_col = input.max_col(); - - window2d - win_fwd_plus = win - (win + left), - win_fwd_minus = (win + left) - win, - win_bkd_plus = win - (win + right), - win_bkd_minus = (win + right) - win, - win_bot = win - (win + up), - win_top = (win + up) - win; - - point2d p; - mln_qiter(W) - q_fp(win_fwd_plus, p), q_fm(win_fwd_minus, p), - q_bp(win_bkd_plus, p), q_bm(win_bkd_minus, p), - q_top(win_top, p), q_bot(win_bot, p); - - accu::median_on<mln_value(I)> med; - - // initialization - - p = input.domain().pmin() + up; - med.init(); - { - mln_qiter(W) q(win, p); - for_all(q) if (input.has(q)) - med.take(input(q)); - } - - int& row = p.row(); - int& col = p.col(); - bool fwd = true; - - mln_assertion(p.col() = min_col); - mln_assertion(p.row() = min_row - 1); - - for (row = min_row; row <= max_row; ++row) - { - // "go down" - for_all(q_top) if (input.has(q_top)) - med.untake(input(q_top)); - for_all(q_bot) if (input.has(q_bot)) - med.take(input(q_bot)); - output(p) = med; - - if (fwd) - // browse line fwd - while (col < max_col) - { - ++col; - for_all(q_fm) if (input.has(q_fm)) - med.untake(input(q_fm)); - for_all(q_fp) if (input.has(q_fp)) - med.take(input(q_fp)); - output(p) = med; - } - else - // browse line bkd - while (col > min_col) - { - --col; - for_all(q_bm) if (input.has(q_bm)) - med.untake(input(q_bm)); - for_all(q_bp) if (input.has(q_bp)) - med.take(input(q_bp)); - output(p) = med; - } - - // change browsing - fwd = ! fwd; - } - } - template <typename I, typename W, typename O> struct median_functor @@ -167,7 +87,7 @@ // ctor - median_functor(I& input_, const W& win_, O& output_) + median_functor(const I& input_, const W& win_, O& output_) : // i/o input(exact(input_)), @@ -224,15 +144,79 @@ }; // end of median_functor + + template <typename I, typename W, typename O> - void median(I& input, const W& win, O& output) + void median(const I& input, const Window<W>& win, O& output) { // FIXME: resize border! - impl::median_functor<I,W,O> f(input, win, output); + impl::median_functor<I,W,O> f(input, exact(win), output); canvas::sbrowsing(f); } + + + template <typename I, typename O> + void median(const I& input, const hline2d& win, O& output) + { + const int + max_row = input.max_row(), + min_col = input.min_col(), + max_col = input.max_col(); + const unsigned half = win.length() / 2; + + point2d p; + int& row = p.row(); + int& col = p.col(); + + accu::median_on<mln_value(I)> med; + + for (row = input.min_row(); row <= max_row; ++row) + { + int ct, cu; + + // initialization (before first point of the row) + med.init(); + for (ct = min_col; ct < min_col + half; ++ct) + med.take(input.at(row, ct)); + + // left columns (just take new points) + for (col = min_col; col <= min_col + half; ++col, ++ct) + { + med.take(input.at(row, ct)); + output(p) = med; + } + + // middle columns (both take and untake) + cu = min_col; + for (; col <= max_col - half; ++cu, ++col, ++ct) + { + med.take(input.at(row, ct)); + med.untake(input.at(row, cu)); + output(p) = med; + } + + // right columns (now just untake old points) + for (; col <= max_col; ++cu, ++col) + { + med.untake(input.at(row, cu)); + output(p) = med; + } + } + } + + + // FIXME: Use transpose. + +// template <typename I, typename O> +// void median(const I& input, const vline2d& win, O& output) +// { + +// median(, hline2d(win.length()), output); +// } + + } // end of namespace mln::level::impl @@ -242,6 +226,7 @@ void median(const Image& input, const Window<W>& win, Image<O>& output) { + mln_assertion(exact(output).domain() = exact(input).domain()); impl::median(exact(input), exact(win), exact(output)); } Index: mln/level/was.hmedian.hh --- mln/level/was.hmedian.hh (revision 0) +++ mln/level/was.hmedian.hh (revision 0) @@ -0,0 +1,102 @@ +// Copyright (C) 2007 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_LEVEL_WAS_HMEDIAN_HH +# define MLN_LEVEL_WAS_HMEDIAN_HH + +/*! \file mln/level/was.hmedian.hh + * + * \brief Median filtering of an image. + */ + + +namespace mln +{ + + namespace level + { + + namespace impl + { + + template <typename I, typename O> + void hmedian(const I& input, const hline2d& win, O& output) + { + const int + max_row = input.max_row(), + min_col = input.min_col(), + max_col = input.max_col(); + const unsigned half = win.length() / 2; + + point2d p; + int& row = p.row(); + int& col = p.col(); + + accu::median_on<mln_value(I)> med; + + for (row = input.min_row(); row <= max_row; ++row) + { + int ct, cu; + + // initialization (before first point of the row) + med.init(); + for (ct = min_col; ct < min_col + half; ++ct) + med.take(input.at(row, ct)); + + // left columns (just take new points) + for (col = min_col; col <= min_col + half; ++col, ++ct) + { + med.take(input.at(row, ct)); + output(p) = med; + } + + // middle columns (both take and untake) + cu = min_col; + for (; col <= max_col - half; ++cu, ++col, ++ct) + { + med.take(input.at(row, ct)); + med.untake(input.at(row, cu)); + output(p) = med; + } + + // right columns (now just untake old points) + for (; col <= max_col; ++cu, ++col) + { + med.untake(input.at(row, cu)); + output(p) = med; + } + } + } + + } // end of namespace mln::level::impl + + } // end of namespace mln::level + +} // end of namespace mln + + +#endif // ! MLN_LEVEL_WAS_HMEDIAN_HH

1 0

1020: Use SCOOP mechanism from static/stc/, not olena/oln/stc/.
by Roland Levillain 17 Jul '07

17 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/olena Index: ChangeLog from Roland Levillain <roland(a)lrde.epita.fr> Use SCOOP mechanism from static/stc/, not olena/oln/stc/. * oln/core/equipment.hh: Include stc/scoop.hh (resp. stc/scoop.hxx) instead of oln/stc/scoop.hh (resp. oln/stc/scoop.hxx). * oln/stc/scoop.hh, oln/stc/scoop.hxx: Remove. * oln/stc: Remove directory. equipment.hh | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) Index: oln/core/equipment.hh --- oln/core/equipment.hh (revision 1019) +++ oln/core/equipment.hh (working copy) @@ -30,7 +30,7 @@ # include <mlc/contract.hh> # ifndef OLENA_USE_SCOOP_ALT -# include <oln/stc/scoop.hh> // FIXME: Remove "oln/" later. +# include <stc/scoop.hh> # else # include <stc/scoop-alt.hh> # endif @@ -46,7 +46,7 @@ { # ifndef OLENA_USE_SCOOP_ALT -# include <oln/stc/scoop.hxx> // FIXME: Remove "oln/" later. +# include <stc/scoop.hxx> # else # include <stc/scoop-alt.inc> # endif

1 0

1019: Update stc/scoop.{hh, hxx} with the latest version from Olena.
by Roland Levillain 17 Jul '07

17 Jul '07

https://svn.lrde.epita.fr/svn/oln/trunk/static Index: ChangeLog from Roland Levillain <roland(a)lrde.epita.fr> Update stc/scoop.{hh,hxx} with the latest version from Olena. * stc/was.scoop.hh: Rename as... * stc/scoop-old-2006.hh: ...this. * stc/scoop.hh: Rename as... * stc/scoop-old-2007.hh: ...this. * stc/scoop.hh, stc/scoop.hxx: New. Imported from olena/stc/. 0 files changed

1 0