March 2005 - Olena-patches

proto-1.0 93: Add dilation and erosion
by Damien Thivolle 20 Mar '05

20 Mar '05

Index: ChangeLog from Damien Thivolle <damien(a)lrde.epita.fr> * tests/morpho/tests/dilation: New. * tests/morpho/tests/erosion: New. * tests/level/Makefile.am: Add missing inclusions. * oln/core/abstract/image_with_data.hh: Use resize_border for const images. * oln/utils/key.hh: Add md5 key implementation. * oln/utils/key.hxx: Likewise. * oln/utils/buffer.hh: Add md5 buffer implementation. * oln/utils/buffer.hxx: Likewise. * oln/utils/md5.hh: Add md5 algorithm implementation. * oln/utils/md5.hxx: Likewise. * oln/utils/clone.hh: Make it respect algorithm coding style. * oln/morpho/dilation.hh: New. Classical dilation implementation * oln/morpho/stat.hh: Replace obsolete use of ::dim by size properties. * oln/morpho/erosion.hh: New. Classical erosion implementation * oln/io/write_image_2d_pnm.hh: Fix binary file handling. * oln/io/read_image_2d_pnm.hh: Likewise. oln/core/abstract/image_with_data.hh | 4 oln/io/read_image_2d_pnm.hh | 4 oln/io/write_image_2d_pnm.hh | 1 oln/morpho/dilation.hh | 178 +++++++++++++++++++ oln/morpho/erosion.hh | 209 ++++++++++++++++++++++ oln/morpho/stat.hh | 12 - oln/utils/buffer.hh | 127 +++++++++++++ oln/utils/buffer.hxx | 163 +++++++++++++++++ oln/utils/clone.hh | 53 +++-- oln/utils/key.hh | 122 +++++++++++++ oln/utils/key.hxx | 80 ++++++++ oln/utils/md5.hh | 241 ++++++++++++++++++++++++++ oln/utils/md5.hxx | 324 +++++++++++++++++++++++++++++++++++ tests/level/Makefile.am | 5 tests/morpho/tests/dilation | 32 +++ tests/morpho/tests/erosion | 31 +++ 16 files changed, 1554 insertions(+), 32 deletions(-) Index: tests/morpho/tests/dilation --- tests/morpho/tests/dilation (revision 0) +++ tests/morpho/tests/dilation (revision 0) @@ -0,0 +1,32 @@ +#include "data.hh" +#include <oln/utils/md5.hh> + +#include <oln/io/read_image.hh> +#include <oln/basics2d.hh> +#include <oln/morpho/dilation.hh> +#include <oln/level/compare.hh> +#include <ntg/all.hh> + +bool check() +{ + oln::utils::key::value_type data_key[16] = + {0xbf, 0x67, 0x3, 0xa1, 0x2d, 0x6c, 0xac, 0xe9, + 0xf3, 0x7c, 0x19, 0x4c, 0x3a, 0x3b, 0x4d, 0xaa}; + + oln::utils::key key(data_key); + + typedef oln::image2d<ntg::bin> im_type; + im_type ima; + ima = oln::io::read(rdata("object.pbm")); + + if (oln::utils::md5(oln::morpho::dilation(ima, oln::win_c8p())) == key) + { + std::cout << "OK" << std::endl; + return false; + } + else + { + std::cout << "FAIL" << std::endl; + return true; + } +} Index: tests/morpho/tests/erosion --- tests/morpho/tests/erosion (revision 0) +++ tests/morpho/tests/erosion (revision 0) @@ -0,0 +1,31 @@ +#include "data.hh" +#include <oln/utils/md5.hh> + +#include <oln/io/read_image.hh> +#include <oln/basics2d.hh> +#include <oln/morpho/erosion.hh> +#include <oln/level/compare.hh> +#include <ntg/all.hh> + +bool check() +{ + oln::utils::key::value_type data_key[16] = + {0x7c, 0x4f, 0xf, 0xc7, 0x27, 0x1, 0x98, 0x81, + 0xfd, 0xb, 0x2e, 0xf6, 0x1c, 0x84, 0xc9, 0x1e}; + + oln::utils::key key(data_key); + + typedef oln::image2d<ntg::bin> im_type; + im_type ima; + ima = oln::io::read(rdata("object.pbm")); + if (oln::utils::md5(oln::morpho::erosion(ima, oln::win_c8p())) == key) + { + std::cout << "OK" << std::endl; + return false; + } + else + { + std::cout << "FAIL" << std::endl; + return true; + } +} Index: tests/level/Makefile.am --- tests/level/Makefile.am (revision 92) +++ tests/level/Makefile.am (working copy) @@ -1,4 +1,5 @@ include ../check/Makefile.runtests +include ../check/Makefile.check -CLEANFILES += cc 2nd_cc connected_out lena_inverted set_level \ - set_level2 threshold_128 + +check-local: check-runtests Index: oln/core/abstract/image_with_data.hh --- oln/core/abstract/image_with_data.hh (revision 92) +++ oln/core/abstract/image_with_data.hh (working copy) @@ -178,9 +178,9 @@ return data_ != 0; } - void resize_border(size_t new_border, bool copy_border = false) + void resize_border(size_t new_border, bool copy_border = false) const { - this->data_->resize_border(new_border, copy_border); + (const_cast<storage_type&>(*this->data_)).resize_border(new_border, copy_border); } Index: oln/utils/buffer.hxx --- oln/utils/buffer.hxx (revision 0) +++ oln/utils/buffer.hxx (revision 0) @@ -0,0 +1,163 @@ +// Copyright (C) 2004 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, but +// WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_OLN_UTILS_BUFFER_HXX +# define OLENA_OLN_UTILS_BUFFER_HXX + +// Ctor implementation +inline +buffer::buffer(): current_bit_(0), + data_(), + len1_(ntg_zero_val(value_type)), + len2_(ntg_zero_val(value_type)), + padding_appended_(false), + length_appended_(false) +{ + data_.reserve(length_); + data_.push_back(ntg_zero_val(value_type)); +} + +// add implementation +template <class E> +inline +void buffer::add(const E &e, bool count) +{ + ntg_storage_type(E) mask = 1 << (ntg_nbits(E) - 1); + ntg_storage_type(E) e_cast = e; + + while (mask) + { + push_back(e_cast & mask, !count); + mask >>= 1; + } +} + +// push_back implementation +inline +void buffer::push_back(bool bit, bool is_padding) +{ + precondition(current_bit_ < 32); + + current_bit_ += 1; + data_.back() *= 2; + if (bit) + data_.back() += ntg_unit_val(value_type); + if (current_bit_ == 32) + { + if (data_.size() == data_.capacity()) + data_.reserve(data_.capacity() + length_); + current_bit_ = 0; + data_.push_back(ntg_zero_val(value_type)); + } + if (!is_padding) + { + if (len1_ == ntg_max_val(value_type)) + { + if (len2_ == ntg_max_val(value_type)) + len2_ = ntg_unit_val(value_type); + else + len2_ += ntg_unit_val(value_type); + } + len1_ += ntg_unit_val(value_type); + } +} + +// reorder implementation +inline +buffer::value_type buffer::reorder(value_type x) const +{ + ntg_storage_type_(value_type) r = 0; + ntg_storage_type_(value_type) mask = 0xFF; + + for (unsigned i = 0; i < (sizeof(ntg_storage_type_(value_type))); ++i) + { + r <<= 8; + r |= (mask & ntg_storage_type_(value_type)(x)) >> (i * 8); + mask <<= 8; + } + return r; +} + +// [] operator implementation +inline +ntg::int_u32 buffer::operator[] (unsigned n) const +{ + precondition(n < data_.size()); + + return reorder(data_[n]); +} + +// append padding implementation +inline +void buffer::append_padding() +{ + push_back(true, true); + while (((data_.size() * 32) % 512) != 448) + push_back(false, true); + padding_appended_ = true; +} + +// append_length implementation +inline +void buffer::append_length() +{ + precondition(padding_appended_); + + add(reorder(len2_), false); + add(reorder(len1_), false); + length_appended_ = true; +} + +// size implementation. +inline +unsigned buffer::size() const +{ + return data_.size(); +} + +// bit len implementation. +inline +unsigned buffer::bit_len() const +{ + return data_.size() * 32; +} + +// is_append_padding implementation. +inline +bool buffer::is_append_padding() const +{ + return padding_appended_; +} + +// is append_length implementation. +inline +bool buffer::is_append_length() const +{ + return length_appended_; +} + +#endif // !OLENA_OLN_UTILS_BUFFER_HXX Index: oln/utils/md5.hh --- oln/utils/md5.hh (revision 0) +++ oln/utils/md5.hh (revision 0) @@ -0,0 +1,241 @@ +// Copyright (C) 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, but +// WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_OLN_UTILS_MD5_HH +# define OLENA_OLN_UTILS_MD5_HH + +# include <ntg/all.hh> +# include <vector> + +# include <oln/core/abstract/image_vectorialness.hh> +# include <oln/core/properties.hh> +# include <oln/core/abstract/piter.hh> + +# include <oln/utils/key.hh> +# include <oln/utils/buffer.hh> + +namespace oln { + namespace utils { + /*! + ** \brief Class used to compute a MD5 digest. + */ + class MD5 + { + public: + typedef ntg::int_u32 value_type; ///< Value to work on. + typedef ntg_storage_type_(value_type) storage_type; ///< Representation type of value type. + + /*! + ** \brief Constructor. + ** + ** Initialization from a buffer. + ** + ** \arg The buffer to use. + */ + explicit MD5(buffer b); + + /*! + ** \brief Compute the digest. + */ + key execute(); + + /*! + ** \brief Generate the key. + ** + ** \pre To avoid meaningless result, you should call execute before. + */ + key gen_key(); + + protected: + /*! + ** \brief initialization of a MD5 algorithm step. + */ + void sav(); + + /*! + ** \brief Round 1 of the algorithm. + */ + void round1(); + + /*! + ** \brief Round 2 of the algorithm. + */ + void round2(); + + /*! + ** \brief Round 3 of the algorithm. + */ + void round3(); + + /*! + ** \brief Round 4 of the algorithm. + */ + void round4(); + + /*! + ** \brief A step of the algorithm. + */ + void step(); + + /*! + ** \brief Sub step to update a register value. + */ + template <class Fun> + void sub_step(const Fun &fun, + storage_type &a, + const storage_type &b, + const storage_type &c, + const storage_type &d, + unsigned k, + unsigned s, + unsigned i); + + /*! + ** \brief Generate a sub part of the key. + ** + ** \arg v Vector of values (output). + ** \arg x Data to put in v. + ** \arg base Where to put x in v. + */ + void gen_sub_key(std::vector<key::value_type> &v, + storage_type x, + unsigned base); + + + buffer buffer_; ///< The buffer to process. + std::vector<storage_type> X; ///< Vector of 16 words. + std::vector<storage_type> T; ///< Vector of 64 words. + unsigned current_; ///< Current position in the buffer. + storage_type aa; ///< Save of a. + storage_type bb; ///< Save of b. + storage_type cc; ///< Save of c. + storage_type dd; ///< Save of d. + storage_type a; ///< A register. + storage_type b; ///< B register. + storage_type c; ///< C register. + storage_type d; ///< D register. + + /*! + ** \brief Functor for bit operations. + ** + ** F(X,Y,Z) = XY v not(X) Z + */ + struct F_ + { + storage_type operator()(const storage_type &x, + const storage_type &y, + const storage_type &z) const + { + // return z ^ (x & (y ^ z)); + return x & y | ~x & z; + } + }; + + F_ F; + + /*! + ** \brief Functor for bit operations. + ** + ** G(X,Y,Z) = XZ v Y not(Z) + */ + struct G_ + { + storage_type operator()(const storage_type &x, + const storage_type &y, + const storage_type &z) const + { + return x & z | y & ~z; + } + }; + + G_ G; + + /*! + ** \brief Functor for bit operations. + ** + ** H(X,Y,Z) = X xor Y xor Z + */ + struct H_ + { + storage_type operator()(const storage_type &x, + const storage_type &y, + const storage_type &z) const + { + return x ^ y ^ z; + } + }; + + H_ H; + + /*! + ** \brief Functor for bit operations. + ** + ** I(X,Y,Z) = Y xor (X v not(Z)) + */ + struct I_ + { + storage_type operator()(const storage_type &x, + const storage_type &y, + const storage_type &z) const + { + return y ^ (x | ~z); + } + }; + + I_ I; + + }; + +#include <oln/utils/md5.hxx> + + /*! + ** \brief Compute The Md5 value of an image. + ** + ** \param I Exact type of the image. + ** + ** \arg im Image to process. + ** + ** Non vectorial image version. + */ + template <class I> + key md5(const oln::abstract::non_vectorial_image &im); + + /*! + ** \brief Compute The Md5 value of an image. + ** + ** \param I Exact type of the image. + ** + ** \arg im Image to process. + ** + ** Vectorial image version. + */ + template <class I> + key md5(const oln::abstract::vectorial_image &im); + + } // !utils +} // !oln +#endif // !OLENA_OLN_UTILS_MD5_HH Index: oln/utils/key.hh --- oln/utils/key.hh (revision 0) +++ oln/utils/key.hh (revision 0) @@ -0,0 +1,122 @@ +// Copyright (C) 2004 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, but +// WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_OLN_UTILS_KEY_HH +# define OLENA_OLN_UTILS_KEY_HH +# include <ntg/all.hh> +# include <vector> +# include <iostream> + +namespace oln { + namespace utils { + /*! + ** \brief 16 bytes key + ** + ** Object to have a nice representation of a MD5 result. + */ + class key + { + public: + typedef ntg::int_u8 value_type; ///< Used data type. + + /*! + ** \brief Constructor used to initialize the key. + ** + ** \pre There must at least 16 elements in data. + ** \arg data Inuput data. + */ + explicit key(const std::vector<value_type> &data); + + /*! + ** \brief Constructor used to initialize the key. + ** + ** \pre There must at least 16 elements in data. + ** \arg data Inuput data. + */ + explicit key(const value_type *data); + + /*! + ** \brief [] operator. + ** + ** \return the ith byte of the key. + ** + ** \arg i Index of the byte wanted. + ** + ** This is the non const version of the operator. + */ + value_type &operator[](unsigned i); + + /*! + ** \brief [] operator. + ** + ** \return the ith byte of the key. + ** + ** \arg i Index of the byte wanted. + ** + ** This is the const version of the operator. + */ + const value_type &operator[](unsigned i) const; + + /*! + ** \brief Check equality between two keys. + ** + ** \arg k Key to be compared with. + */ + bool operator ==(const key &k); + + /*! + ** \brief Check whether two key are different. + ** + ** \arg k Key to be compared with. + */ + bool operator !=(const key &k); + + /*! + ** \brief Write of representation of the key on a stream. + ** + ** stream Stream to put the key on. + ** k Key to represent. + */ + friend std::ostream &operator<<(std::ostream &stream, const key &k) + { + stream << "{"; + for (unsigned i = 0; i < 15; ++i) + stream << "0x" << std::hex << k[i] << ", "; + stream << "0x" << std::hex << k[15] << "}"; + return stream; + } + + + protected: + value_type data_[16]; ///< Internal data. + }; + +# include <oln/utils/key.hxx> + + } // !utils +} // !oln +#endif // !OLENA_OLN_UTILS_KEY_HH Index: oln/utils/buffer.hh --- oln/utils/buffer.hh (revision 0) +++ oln/utils/buffer.hh (revision 0) @@ -0,0 +1,127 @@ +// Copyright (C) 2004 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, but +// WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + + +#ifndef OLENA_OLN_UTILS_BUFFER_HH +# define OLENA_OLN_UTILS_BUFFER_HH + +# include <ntg/int.hh> +# include <vector> + +namespace oln { + namespace utils { + /*! + ** \brief Buffer used for MD5 data type abstraction. + */ + class buffer + { + public: + typedef ntg::int_u32 value_type; ///< Used data type. + + /*! + ** \brief Initialization of data. + */ + buffer(); + + /*! + ** \brief Add every bits of e in data. + ** + ** \arg e Element to work on. + ** \arg count Tell if you want to count e size in the buffer size. + */ + template <class E> + void add(const E &e, bool count = true); + + /*! + ** \brief Push back a bit at the end of the buffer. + ** + ** true -> push back a 1\n + ** false -> push back a 0\n + ** + ** \arg bit Bit to push. + ** \arg is_padding Are you adding padding bytes ? + */ + void push_back(bool bit, bool is_padding = false); + + /*! + ** \brief Return the nth word. + */ + ntg::int_u32 operator[] (unsigned n) const; + + /*! + ** \brief Append padding bits. + */ + void append_padding(); + + /*! + ** \brief Append length. + */ + void append_length(); + + /*! + ** \brief Return the number of words. + */ + unsigned size() const; + + /*! + ** \brief Return the length in bits. + */ + unsigned bit_len() const; + + /*! + ** \brief Tell if the buffer has already been padded. + */ + bool is_append_padding() const; + + /*! + ** \brief Tell if the length has already been appended. + */ + bool is_append_length() const; + + protected: + /*! + ** \brief Change order of data in a word. + ** + ** \arg x Data to reorder. + */ + value_type reorder(value_type x) const; + + static const unsigned length_ = 100; ///< Capacity chunk. + /// Number of the current bit in the current word. + unsigned char current_bit_; + std::vector<value_type> data_; ///< The buffer data. + value_type len1_; ///< first word length + value_type len2_; ///< second word length + bool padding_appended_; ///< Status of padding. + bool length_appended_; ///< Status of length appending. + }; + +# include <oln/utils/buffer.hxx> + + } // !utils +} // !oln +#endif // !OLENA_OLN_UTILS_BUFFER_HH Index: oln/utils/md5.hxx --- oln/utils/md5.hxx (revision 0) +++ oln/utils/md5.hxx (revision 0) @@ -0,0 +1,324 @@ +// Copyright (C) 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, but +// WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_OLN_UTILS_MD5_HXX +# define OLENA_OLN_UTILS_MD5_HXX + +// Ctor. +inline +MD5::MD5(buffer b): buffer_(b), + X(16), + T(), + current_(0), + a(0x67452301), + b(0xefcdab89), + c(0x98badcfe), + d(0x10325476) + // a(0x01234567), + // b(0x89abcdef), + // c(0xfedcba98), + // d(0x76543210) +{ + precondition(buffer_.is_append_padding()); + precondition(buffer_.is_append_length()); + T.reserve(64); + T.push_back(0xd76aa478); /* 1 */ + T.push_back(0xe8c7b756); /* 2 */ + T.push_back(0x242070db); /* 3 */ + T.push_back(0xc1bdceee); /* 4 */ + T.push_back(0xf57c0faf); /* 5 */ + T.push_back(0x4787c62a); /* 6 */ + T.push_back(0xa8304613); /* 7 */ + T.push_back(0xfd469501); /* 8 */ + T.push_back(0x698098d8); /* 9 */ + T.push_back(0x8b44f7af); /* 10 */ + T.push_back(0xffff5bb1); /* 11 */ + T.push_back(0x895cd7be); /* 12 */ + T.push_back(0x6b901122); /* 13 */ + T.push_back(0xfd987193); /* 14 */ + T.push_back(0xa679438e); /* 15 */ + T.push_back(0x49b40821); /* 16 */ + T.push_back(0xf61e2562); /* 17 */ + T.push_back(0xc040b340); /* 18 */ + T.push_back(0x265e5a51); /* 19 */ + T.push_back(0xe9b6c7aa); /* 20 */ + T.push_back(0xd62f105d); /* 21 */ + T.push_back(0x2441453); /* 22 */ + T.push_back(0xd8a1e681); /* 23 */ + T.push_back(0xe7d3fbc8); /* 24 */ + T.push_back(0x21e1cde6); /* 25 */ + T.push_back(0xc33707d6); /* 26 */ + T.push_back(0xf4d50d87); /* 27 */ + T.push_back(0x455a14ed); /* 28 */ + T.push_back(0xa9e3e905); /* 29 */ + T.push_back(0xfcefa3f8); /* 30 */ + T.push_back(0x676f02d9); /* 31 */ + T.push_back(0x8d2a4c8a); /* 32 */ + T.push_back(0xfffa3942); /* 33 */ + T.push_back(0x8771f681); /* 34 */ + T.push_back(0x6d9d6122); /* 35 */ + T.push_back(0xfde5380c); /* 36 */ + T.push_back(0xa4beea44); /* 37 */ + T.push_back(0x4bdecfa9); /* 38 */ + T.push_back(0xf6bb4b60); /* 39 */ + T.push_back(0xbebfbc70); /* 40 */ + T.push_back(0x289b7ec6); /* 41 */ + T.push_back(0xeaa127fa); /* 42 */ + T.push_back(0xd4ef3085); /* 43 */ + T.push_back(0x4881d05); /* 44 */ + T.push_back(0xd9d4d039); /* 45 */ + T.push_back(0xe6db99e5); /* 46 */ + T.push_back(0x1fa27cf8); /* 47 */ + T.push_back(0xc4ac5665); /* 48 */ + T.push_back(0xf4292244); /* 49 */ + T.push_back(0x432aff97); /* 50 */ + T.push_back(0xab9423a7); /* 51 */ + T.push_back(0xfc93a039); /* 52 */ + T.push_back(0x655b59c3); /* 53 */ + T.push_back(0x8f0ccc92); /* 54 */ + T.push_back(0xffeff47d); /* 55 */ + T.push_back(0x85845dd1); /* 56 */ + T.push_back(0x6fa87e4f); /* 57 */ + T.push_back(0xfe2ce6e0); /* 58 */ + T.push_back(0xa3014314); /* 59 */ + T.push_back(0x4e0811a1); /* 60 */ + T.push_back(0xf7537e82); /* 61 */ + T.push_back(0xbd3af235); /* 62 */ + T.push_back(0x2ad7d2bb); /* 63 */ + T.push_back(0xeb86d391); /* 64 */ +} + +// sav() implementation. +inline +void MD5::sav() +{ + for (unsigned i = 0; i < 16; ++i) + X[i] = buffer_[current_ + i]; + current_ += 16; + aa = a; + bb = b; + cc = c; + dd = d; +} + +// round1() implementation. +inline +void MD5::round1() +{ + sub_step(F, a, b, c, d, 0, 7, 1); + sub_step(F, d, a, b, c, 1, 12, 2); + sub_step(F, c, d, a, b, 2, 17, 3); + sub_step(F, b, c, d, a, 3, 22, 4); + sub_step(F, a, b, c, d, 4, 7, 5); + sub_step(F, d, a, b, c, 5, 12, 6); + sub_step(F, c, d, a, b, 6, 17, 7); + sub_step(F, b, c, d, a, 7, 22, 8); + sub_step(F, a, b, c, d, 8, 7, 9); + sub_step(F, d, a, b, c, 9, 12, 10); + sub_step(F, c, d, a, b, 10, 17, 11); + sub_step(F, b, c, d, a, 11, 22, 12); + sub_step(F, a, b, c, d, 12, 7, 13); + sub_step(F, d, a, b, c, 13, 12, 14); + sub_step(F, c, d, a, b, 14, 17, 15); + sub_step(F, b, c, d, a, 15, 22, 16); +} + +// round2() implementation. +inline +void MD5::round2() +{ + sub_step(G, a, b, c, d, 1, 5, 17); + sub_step(G, d, a, b, c, 6, 9, 18); + sub_step(G, c, d, a, b, 11, 14, 19); + sub_step(G, b, c, d, a, 0, 20, 20); + sub_step(G, a, b, c, d, 5, 5, 21); + sub_step(G, d, a, b, c, 10, 9, 22); + sub_step(G, c, d, a, b, 15, 14, 23); + sub_step(G, b, c, d, a, 4, 20, 24); + sub_step(G, a, b, c, d, 9, 5, 25); + sub_step(G, d, a, b, c, 14, 9, 26); + sub_step(G, c, d, a, b, 3, 14, 27); + sub_step(G, b, c, d, a, 8, 20, 28); + sub_step(G, a, b, c, d, 13, 5, 29); + sub_step(G, d, a, b, c, 2, 9, 30); + sub_step(G, c, d, a, b, 7, 14, 31); + sub_step(G, b, c, d, a, 12, 20, 32); +} + +// round3() implementation. +inline +void MD5::round3() +{ + sub_step(H, a, b, c, d, 5, 4, 33); + sub_step(H, d, a, b, c, 8, 11, 34); + sub_step(H, c, d, a, b, 11, 16, 35); + sub_step(H, b, c, d, a, 14, 23, 36); + sub_step(H, a, b, c, d, 1, 4, 37); + sub_step(H, d, a, b, c, 4, 11, 38); + sub_step(H, c, d, a, b, 7, 16, 39); + sub_step(H, b, c, d, a, 10, 23, 40); + sub_step(H, a, b, c, d, 13, 4, 41); + sub_step(H, d, a, b, c, 0, 11, 42); + sub_step(H, c, d, a, b, 3, 16, 43); + sub_step(H, b, c, d, a, 6, 23, 44); + sub_step(H, a, b, c, d, 9, 4, 45); + sub_step(H, d, a, b, c, 12, 11, 46); + sub_step(H, c, d, a, b, 15, 16, 47); + sub_step(H, b, c, d, a, 2, 23, 48); +} + +// round 4() implementation. +inline +void MD5::round4() +{ + sub_step(I, a, b, c, d, 0, 6, 49); + sub_step(I, d, a, b, c, 7, 10, 50); + sub_step(I, c, d, a, b, 14, 15, 51); + sub_step(I, b, c, d, a, 5, 21, 52); + sub_step(I, a, b, c, d, 12, 6, 53); + sub_step(I, d, a, b, c, 3, 10, 54); + sub_step(I, c, d, a, b, 10, 15, 55); + sub_step(I, b, c, d, a, 1, 21, 56); + sub_step(I, a, b, c, d, 8, 6, 57); + sub_step(I, d, a, b, c, 15, 10, 58); + sub_step(I, c, d, a, b, 6, 15, 59); + sub_step(I, b, c, d, a, 13, 21, 60); + sub_step(I, a, b, c, d, 4, 6, 61); + sub_step(I, d, a, b, c, 11, 10, 62); + sub_step(I, c, d, a, b, 2, 15, 63); + sub_step(I, b, c, d, a, 9, 21, 64); +} + +// step implementation. +inline +void MD5::step() +{ + sav(); + round1(); + round2(); + round3(); + round4(); + a += aa; + b += bb; + c += cc; + d += dd; +} + +// sub step implementation. +template <class Fun> +inline +void MD5::sub_step(const Fun &fun, + storage_type &a, + const storage_type &b, + const storage_type &c, + const storage_type &d, + unsigned k, + unsigned s, + unsigned i) +{ + a += fun(b, c, d) + (X[k] + T[i - 1]); + a = a<<s | a>>(32-s); + a += b; + //a = b + ((a + fun(b, c, d) + X[k] + T[i - 1]) << s); +} + +// gen_sub_key implementation. +inline +void MD5::gen_sub_key(std::vector<key::value_type> &v, + storage_type x, + unsigned base) +{ + base *= sizeof(value_type); + + for (unsigned i = 0; i < sizeof(value_type); ++i) + v[base + i] = (x & (0xFF << (i * 8))) >> (i * 8); +} + +// gen_key implementation. + +inline +key MD5::gen_key() +{ + std::vector<key::value_type> k; + + k.reserve(16); + gen_sub_key(k, a, 0); + gen_sub_key(k, b, 1); + gen_sub_key(k, c, 2); + gen_sub_key(k, d, 3); + return key(k); +} + +// execute implementation. +inline +key MD5::execute() +{ + while (buffer_.size() > current_) + step(); + return gen_key(); +} + + +// md5 implementation. +template <class I> +inline +key md5(const oln::abstract::non_vectorial_image &im) +{ + buffer b; + oln_type_of(I, fwd_piter) p(im.size()); + + for_all(p) + b.add(im[p].value()); + b.append_padding(); + b.append_length(); + + oln::utils::MD5 md5(b); + return md5.execute(); +} + +// md5 implementation. +template <class I> +inline +key md5(const oln::abstract::vectorial_image &im) +{ + oln_type_of(I, fwd_piter) p(im.size()); + buffer b; + // unsigned nb = + + for_all(p) + { + for (unsigned i = 0; i < ntg_nb_comp(oln_type_of(I, value)); ++i) + b.add(im[p].value()[i]); + } + b.append_padding(); + b.append_length(); + + oln::utils::MD5 md5(b); + return md5.execute(); +} + +#endif // !OLENA_OLN_UTILS_MD5_HXX Index: oln/utils/clone.hh --- oln/utils/clone.hh (revision 92) +++ oln/utils/clone.hh (working copy) @@ -35,56 +35,63 @@ // fwd decl namespace utils { + + namespace impl { template <typename I> struct clone_type; + } + } // category template <typename I> - struct set_category < utils::clone_type > + struct set_category < utils::impl::clone_type > { typedef category::image ret; }; // super_type template <typename I> - struct set_super_type < utils::clone_type > + struct set_super_type < utils::impl::clone_type > { - typedef abstract::op<I, utils::clone_type > ret; + typedef abstract::op<I, utils::impl::clone_type > ret; }; namespace utils { - template <typename I> - struct clone_type : abstract::op<I, clone_type > - { - typedef abstract::op<I, clone_type > super_type; + namespace impl { - mlc::box input_; + template <typename I> + struct clone_type : abstract::op<I, clone_type > + { + typedef abstract::op<I, clone_type > super_type; - clone_type(I& input) : input_(input) - {} + mlc::box input_; - void impl_run() - { - I ima(input_->size()); - oln_type_of(I, fwd_piter) it(input_->size()); + clone_type(I& input) : input_(input) + {} - for_all(it) - { - ima[it] = (*input_)[it]; - } + void impl_run() + { + I ima(input_->size()); + oln_type_of(I, fwd_piter) it(input_->size()); - *this->image_ = ima; - } - }; + for_all(it) + { + ima[it] = (*input_)[it]; + } + *this->image_ = ima; + } + }; + } + template <typename I> - clone_type clone(abstract::image& ima) + impl::clone_type clone(abstract::image& ima) { - clone_type tmp(ima.exact()); + impl::clone_type tmp(ima.exact()); tmp.run(); return tmp; } Index: oln/utils/key.hxx --- oln/utils/key.hxx (revision 0) +++ oln/utils/key.hxx (revision 0) @@ -0,0 +1,80 @@ +// Copyright (C) 2004 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, but +// WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_OLN_UTILS_KEY_HXX +# define OLENA_OLN_UTILS_KEY_HXX + +// Ctor implementation. +inline +key::key(const std::vector<value_type> &data) +{ + for (unsigned i = 0; i < 16; ++i) + data_[i] = data[i]; +} + +// Ctor implementation. +inline +key::key(const value_type *data) +{ + for (unsigned i = 0; i < 16; ++i) + data_[i] = data[i]; +} + +// [] operator implementation +inline +key::value_type &key::operator[](unsigned i) +{ + precondition(i < 16); + return data_[i]; +} + +// [] operator +inline +const key::value_type &key::operator[](unsigned i) const +{ + precondition(i < 16); + return data_[i]; +} + +// == operator +inline +bool key::operator ==(const key &k) +{ + for (unsigned i = 0; i < 16; ++i) + if (data_[i] != k[i]) + return false; + return true; +} + +// != operator +inline +bool key::operator !=(const key &k) +{ + return !(*this == k); +} + +#endif // !OLENA_OLN_UTILS_KEY_HXX Index: oln/morpho/dilation.hh --- oln/morpho/dilation.hh (revision 0) +++ oln/morpho/dilation.hh (revision 0) @@ -0,0 +1,178 @@ +// Copyright (C) 2001, 2002, 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this filek as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_MORPHO_DILATION_HH +# define OLENA_MORPHO_DILATION_HH + +# include <oln/basics.hh> +# include <oln/morpho/stat.hh> +// # include <oln/morpho/fast_morpho.hh> +# include <mlc/cmp.hh> + +namespace oln { + + namespace morpho { + + /*! + ** \brief Processing dilation. + ** + ** Compute the morphological dilation of input using se + ** as structural element.\n + ** + ** On grey-scale images, each point is replaced by the maximum value + ** of its neighbors, as indicated by se. On binary images, + ** a logical or is performed between neighbors.\n + ** + ** The morpho::fast version of this function use a different + ** algorithm: This algorithm is described in + ** Implementation of morphological operations from: + ** M. Van Droogenbroeck and H. Talbot. + ** "Fast computation of morphological operations with arbitrary + ** structuring elements". Pattern Recognition Letters, + ** 17(14):1451-1460, 1996.\n + ** + ** An histogram of the value of the neighborhood indicated by + ** se is updated while iterating over all point of the + ** image. Doing so is more efficient when the + ** structural element is large.\n + ** + ** \param I Exact type of the input image. + ** \param E Exact type of the neighborhood. + ** + ** \arg input The input image. + ** \arg se Structuring element to use. + ** + ** \code + ** #include <oln/basics2d.hh> + ** #include <oln/morpho/dilation.hh> + ** #include <oln/level/compare.hh> + ** #include <ntg/all.hh> + ** int main() + ** { + ** typedef oln::image2d<ntg::bin> im_type; + ** + ** im_type im1(oln::load(IMG_IN "object.pbm")); + ** save(oln::morpho::dilation(im1, oln::win_c8p()), + ** IMG_OUT "oln_morpho_dilation.pbm"); + ** } + ** \endcode + ** + ** \image html object_pbm.png + ** \image latex object_pbm.png + ** => + ** \image html oln_morpho_dilation.png + ** \image latex oln_morpho_dilation.png + */ + + // fwd decl + namespace impl { + template <typename I, typename E> struct dilation_type; + } + + } + + // category + template <typename I, typename E> + struct set_category< morpho::impl::dilation_type<I,E> > { typedef category::image ret; }; + + // super_type + template <typename I, typename E> + struct set_super_type< morpho::impl::dilation_type<I,E> > + { + typedef abstract::op<I, morpho::impl::dilation_type<I, E> > ret; + }; + + namespace morpho { + + namespace impl { + + template <class I, class E> + struct dilation_type : abstract::op<I, dilation_type<I, E> > + { + mlc::box<const I> input_; + mlc::box<const E> se_; + + dilation_type(const abstract::non_vectorial_image& input, + const abstract::struct_elt<E>& se) : + input_(input.exact()), + se_(se.exact()) + {} + + }; + + template <class I, class E> + struct dilation_type_classic : dilation_type<I, E> + { + typedef dilation_type<I, E> super_type; + + dilation_type_classic(const abstract::non_vectorial_image& input, + const abstract::struct_elt<E>& se) : + super_type(input, se) + {} + + void impl_run() + { + mlc::is_true<mlc::type::eq<oln_type_of(I, size), + oln_type_of(E, size)>::ret>::ensure(); + I output(this->input_->size()); + this->input_->resize_border(this->se_->get_delta()); + oln_type_of(I, fwd_piter) p(this->input_->size()); + + for_all (p) + output[p] = morpho::max(*this->input_, p, *this->se_); + *this->image_ = output; + } + }; + + } + + template<class I, class E> + impl::dilation_type<I, E> + dilation(const abstract::non_vectorial_image& input, + const abstract::struct_elt<E>& se) + { + impl::dilation_type_classic<I, E> tmp(input, se); + tmp.run(); + return tmp; + } + + +// namespace fast { +// template<class I, class E> +// oln_concrete_type(I) +// dilation(const abstract::non_vectorial_image& input, +// const abstract::struct_elt<E>& se) +// { +// return fast_morpho<I, E, utils::histogram_max<oln_value_type(I)> > +// (input, se); +// } +// } + } // end of morpho + +} // end of oln + +#endif // ! OLENA_MORPHO_DILATION_HH Index: oln/morpho/stat.hh --- oln/morpho/stat.hh (revision 92) +++ oln/morpho/stat.hh (working copy) @@ -63,7 +63,9 @@ static V max(const I& input, const oln_type_of(I, point)& p, const E& se) { - // FIXME: test dim I == dim E + // FIXME: test dim I == dim + mlc::is_true<mlc::type::eq<oln_type_of(I, size), + oln_type_of(E, size)>::ret>::ensure(); oln_type_of(E, fwd_witer) dp(se); dp.start(); @@ -88,6 +90,8 @@ min(const I& input, const oln_type_of(I, point)& p, const E& se) { // FIXME: test dim I == dim E + mlc::is_true<mlc::type::eq<oln_type_of(I, size), + oln_type_of(E, size)>::ret>::ensure(); oln_type_of(E, fwd_witer) dp(se); dp.start(); @@ -108,7 +112,8 @@ static ntg::bin max(const I& input, const oln_type_of(I, point)& p, const E& se) { - mlc::eq<I::dim, E::dim>::ensure(); + mlc::is_true<mlc::type::eq<oln_type_of(I, size), + oln_type_of(E, size)>::ret>::ensure(); oln_type_of(E, fwd_witer) dp(se); for_all (dp) if (input[p + dp] == true) @@ -119,7 +124,8 @@ static ntg::bin min(const I& input, const oln_type_of(I, point)& p, const E& se) { - mlc::eq<I::dim, E::dim>::ensure(); + mlc::is_true<mlc::type::eq<oln_type_of(I, size), + oln_type_of(E, size)>::ret>::ensure(); oln_type_of(E, fwd_witer) dp(se); for_all (dp) if (input[p + dp] == false) Index: oln/morpho/erosion.hh --- oln/morpho/erosion.hh (revision 0) +++ oln/morpho/erosion.hh (revision 0) @@ -0,0 +1,209 @@ +// Copyright (C) 2001, 2002, 2003, 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_MORPHO_EROSION_HH +# define OLENA_MORPHO_EROSION_HH + +# include <oln/basics.hh> +# include <oln/morpho/stat.hh> +// # include <oln/morpho/fast_morpho.hh> +# include <oln/core/abstract/op.hh> +# include <mlc/cmp.hh> + +namespace oln { + + namespace morpho { + /*! + ** \brief Perform a morphological erosion. + ** + ** Compute the morphological erosion of input using se + ** as structuring element. + ** + ** On grey-scale images, each point is replaced by the minimum + ** value of its neighbors, as indicated by se. On binary + ** images, a logical and is performed between neighbors. The + ** morpho::fast version of this function use a different + ** algorithm: an histogram of the value of the neighborhood + ** indicated by se is updated while iterating over all point of + ** the image. Doing so is more efficient when the structuring + ** element is large. + ** + ** \param I Exact type of the input image. + ** \param E Exact type of the structuring element. + ** + ** + ** \arg input Input image. + ** \arg se Structuring element to use. + ** + ** \code + ** #include <oln/basics2d.hh> + ** #include <oln/morpho/erosion.hh> + ** #include <oln/level/compare.hh> + ** #include <ntg/all.hh> + ** int main() + ** { + ** typedef oln::image2d<ntg::bin> im_type; + ** + ** im_type im1(oln::load(IMG_IN "object.pbm")); + ** save(oln::morpho::erosion(im1, oln::win_c8p()), + ** IMG_OUT "oln_morpho_erosion.pbm"); + ** } + ** \endcode + ** + ** \image html object_pbm.png + ** \image latex object_pbm.png + ** => + ** \image html oln_morpho_erosion.png + ** \image latex oln_morpho_erosion.png + ** + ** \see oln::morpho::fast::erosion() + */ + + // fwd decl + namespace impl { + template <typename I, typename E> struct erosion_type; + } + + } + + // category + template <typename I, typename E> + struct set_category< morpho::impl::erosion_type<I,E> > { typedef category::image ret; }; + + // super_type + template <typename I, typename E> + struct set_super_type< morpho::impl::erosion_type<I,E> > + { + typedef abstract::op<I, morpho::impl::erosion_type<I, E> > ret; + }; + + namespace morpho { + + namespace impl { + + template <class I, class E> + struct erosion_type : abstract::op<I, erosion_type<I, E> > + { + mlc::box<const I> input_; + mlc::box<const E> se_; + + erosion_type(const abstract::non_vectorial_image& input, + const abstract::struct_elt<E>& se) : + input_(input.exact()), + se_(se.exact()) + {} + + }; + + template <class I, class E> + struct erosion_type_classic : erosion_type<I, E> + { + typedef erosion_type<I, E> super_type; + + erosion_type_classic(const abstract::non_vectorial_image& input, + const abstract::struct_elt<E>& se) : + super_type(input, se) + {} + + void impl_run() + { + mlc::is_true<mlc::type::eq<oln_type_of(I, size), + oln_type_of(E, size)>::ret>::ensure(); + I output(this->input_->size()); + this->input_->resize_border(this->se_->get_delta()); + oln_type_of(I, fwd_piter) p(this->input_->size()); + + for_all (p) + output[p] = morpho::min(*this->input_, p, *this->se_); + *this->image_ = output; + } + }; + + } + + + template<class I, class E> + impl::erosion_type<I, E> + erosion(const abstract::non_vectorial_image& input, + const abstract::struct_elt<E>& se) + { + impl::erosion_type_classic<I, E> tmp(input, se); + tmp.run(); + return tmp; + } + + namespace fast { + /*! + ** \brief Perform a morphological erosion. + ** + ** Compute the morphological erosion of input using se + ** as structuring element. + ** + ** \param I Exact type of the input image. + ** \param E Exact type of the structuring element. + ** + ** + ** \arg input Input image. + ** \arg se Structuring element to use. + ** + ** \code + ** #include <oln/basics2d.hh> + ** #include <oln/morpho/erosion.hh> + ** #include <oln/level/compare.hh> + ** #include <ntg/all.hh> + ** int main() + ** { + ** typedef oln::image2d<ntg::bin> im_type; + ** + ** im_type im1(oln::load(IMG_IN "object.pbm")); + ** save(oln::morpho::fast::erosion(im1, oln::win_c8p()), + ** IMG_OUT "oln_morpho_fast_erosion.pbm"); + ** } + ** \endcode + ** + ** \image html object_pbm.png + ** \image latex object_pbm.png + ** => + ** \image html oln_morpho_fast_erosion.png + ** \image latex oln_morpho_fast_erosion.png + */ +// template<class I, class E> +// oln_concrete_type(I) +// erosion(const abstract::non_vectorial_image& input, +// const abstract::struct_elt<E>& se) +// { +// return fast_morpho<I, E, utils::histogram_min<oln_value_type(I)> > +// (input, se); +// } + } // end of fast + + } // end of morpho + +} // end of oln + + +#endif // ! OLENA_MORPHO_EROSION_HH Index: oln/io/write_image_2d_pnm.hh --- oln/io/write_image_2d_pnm.hh (revision 92) +++ oln/io/write_image_2d_pnm.hh (working copy) @@ -116,6 +116,7 @@ if (bin_offset == -1) { + bin_v = ~bin_v; ostr_.write(&bin_v, 1); bin_offset = 7; bin_v = 0; Index: oln/io/read_image_2d_pnm.hh --- oln/io/read_image_2d_pnm.hh (revision 92) +++ oln/io/read_image_2d_pnm.hh (working copy) @@ -172,9 +172,9 @@ offset = 7; } if ((int)(v & (1<<offset--)) == 0) + c = 1; + else c = 0; - else - c = 1; } };

1 0

proto-1.0 92: Add test for morpho min and max
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * tests/morpho/tests/stat: New. Test morpho min and max. stat | 19 +++++++++++++++++++ 1 files changed, 19 insertions(+) Index: tests/morpho/tests/stat --- tests/morpho/tests/stat (revision 0) +++ tests/morpho/tests/stat (revision 0) @@ -0,0 +1,19 @@ +#include <ntg/all.hh> +#include <oln/basics2d.hh> +#include <oln/level/fill.hh> +#include <oln/fancy/iota.hh> +#include <oln/morpho/stat.hh> + +#include "check.hh" +#include "data.hh" + +bool check() +{ + oln::window2d w2 = oln::win_c8p(); + oln::image2d<ntg::int_u8> im(10, 10); + oln::fancy::iota(im); + if (oln::morpho::max(im, oln::point2d(5, 5), w2) != 67 || + oln::morpho::min(im, oln::point2d(5, 5), w2) != 45) + return true; + return false; +}

1 0

proto-1.0 91: Add arith tests suite to configure.ac
by Damien Thivolle 20 Mar '05

20 Mar '05

Index: ChangeLog from Damien Thivolle <damien(a)lrde.epita.fr> * configure.ac: Add arith tests to olena tests suite. configure.ac | 1 + 1 files changed, 1 insertion(+) Index: configure.ac --- configure.ac (revision 90) +++ configure.ac (working copy) @@ -99,6 +99,7 @@ # OLN_TESTS defined in oln-local.m4 OLN_TESTS([sanity]) + OLN_TESTS([arith]) OLN_TESTS([convert]) OLN_TESTS([core]) OLN_TESTS([io])

1 0

proto-1.0 90: Add arith::min and arith::max
by Damien Thivolle 20 Mar '05

20 Mar '05

Index: ChangeLog from Damien Thivolle <damien(a)lrde.epita.fr> * tests/arith: New. * tests/arith/tests: New. * tests/arith/tests/min: New. * tests/arith/tests/max: New. * tests/arith/Makefile.am: New. * oln/basics.hh: Include oln/core/abstract/op.hh. * oln/arith: New. * oln/arith/min.hh: Compute the pointwise minimum between two images. * oln/arith/max.hh: Compute the pointwise maximum between two images. oln/arith/max.hh | 104 ++++++++++++++++++++++++++++++++++++++++++++++++ oln/arith/min.hh | 104 ++++++++++++++++++++++++++++++++++++++++++++++++ oln/basics.hh | 1 oln/makefile.src | 2 tests/arith/Makefile.am | 5 ++ tests/arith/tests/max | 25 +++++++++++ tests/arith/tests/min | 27 ++++++++++++ 7 files changed, 268 insertions(+) Index: tests/arith/tests/min --- tests/arith/tests/min (revision 0) +++ tests/arith/tests/min (revision 0) @@ -0,0 +1,27 @@ +#include <oln/basics2d.hh> +#include <oln/arith/min.hh> +#include <oln/level/fill.hh> +#include <oln/level/compare.hh> +#include <ntg/int.hh> + + + +bool check() +{ + oln::image2d<ntg::int_u8> ima1(10, 10); + oln::image2d<ntg::int_u8> ima2(10, 10); + + oln::level::fill(ima1, 10); + oln::level::fill(ima2, 20); + + oln::image2d<ntg::int_u8> ima; + + ima = oln::arith::min(ima1, ima2); + + if (oln::level::is_equal(ima, ima1)) + return false; + return true; +} + + + Index: tests/arith/tests/max --- tests/arith/tests/max (revision 0) +++ tests/arith/tests/max (revision 0) @@ -0,0 +1,25 @@ +#include <oln/basics2d.hh> +#include <oln/arith/max.hh> +#include <oln/level/fill.hh> +#include <oln/level/compare.hh> +#include <ntg/int.hh> + + + +bool check() +{ + oln::image2d<ntg::int_u8> ima1(10, 10); + oln::image2d<ntg::int_u8> ima2(10, 10); + + oln::level::fill(ima1, 10); + oln::level::fill(ima2, 20); + + oln::image2d<ntg::int_u8> ima; + + ima = oln::arith::max(ima1, ima2); + + if (oln::level::is_equal(ima, ima2)) + return false; + return true; +} + Index: tests/arith/Makefile.am --- tests/arith/Makefile.am (revision 0) +++ tests/arith/Makefile.am (revision 0) @@ -0,0 +1,5 @@ +include ../check/Makefile.runtests +include ../check/Makefile.check + + +check-local: check-runtests Index: oln/basics.hh --- oln/basics.hh (revision 89) +++ oln/basics.hh (working copy) @@ -58,6 +58,7 @@ # include <oln/core/abstract/point.hh> # include <oln/core/abstract/images.hh> # include <oln/core/abstract/entry.hh> +# include <oln/core/abstract/op.hh> # include <oln/core/abstract/piter.hh> # include <oln/core/abstract/witer.hh> Index: oln/makefile.src --- oln/makefile.src (revision 89) +++ oln/makefile.src (working copy) @@ -5,6 +5,8 @@ OLN_DEP = \ all.hh \ + arith/max.hh \ + arith/min.hh \ basics.hh \ basics1d.hh \ basics2d.hh \ Index: oln/arith/min.hh --- oln/arith/min.hh (revision 0) +++ oln/arith/min.hh (revision 0) @@ -0,0 +1,104 @@ +// Copyright (C) 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_ARITH_MIN_HH +# define OLENA_ARITH_MIN_HH + +# include <oln/basics.hh> +# include <ntg/all.hh> + +namespace oln { + + namespace arith { + // fwd decl + namespace impl { + template <typename I> struct min_type; + } + + } + + // category + template <typename I> + struct set_category< arith::impl::min_type > + { + typedef category::image ret; + }; + + // super_type + template <typename I> + struct set_super_type< arith::impl::min_type > + { + typedef abstract::op<I, arith::impl::min_type > ret; + }; + + namespace arith { + + namespace impl { + + template <class I> + struct min_type : abstract::op<I, min_type > + { + mlc::box<const I> input1_; + mlc::box<const I> input2_; + + min_type(const abstract::non_vectorial_image& input1, + const abstract::non_vectorial_image& input2) : + input1_(input1.exact()), + input2_(input2.exact()) + {} + + void impl_run() + { + precondition(input1_->size() == input2_->size()); + I output(input1_->size()); + oln_type_of(I, fwd_piter) p(input1_->size()); + + for_all(p) + output[p] = ntg::min((*input1_)[p].value(), (*input2_)[p].value()); + + *this->image_ = output; + } + + }; + + } + + template <typename I> + impl::min_type min(const abstract::non_vectorial_image& input1, + const abstract::non_vectorial_image& input2) + { + impl::min_type tmp(input1, input2); + tmp.run(); + return tmp; + } + + } + +} // end of namespace oln + + +#endif // ! OLENA_ARITH_MIN_HH Index: oln/arith/max.hh --- oln/arith/max.hh (revision 0) +++ oln/arith/max.hh (revision 0) @@ -0,0 +1,104 @@ +// Copyright (C) 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_ARITH_MAX_HH +# define OLENA_ARITH_MAX_HH + +# include <oln/basics.hh> +# include <ntg/all.hh> + +namespace oln { + + namespace arith { + // fwd decl + namespace impl { + template <typename I> struct max_type; + } + + } + + // category + template <typename I> + struct set_category< arith::impl::max_type > + { + typedef category::image ret; + }; + + // super_type + template <typename I> + struct set_super_type< arith::impl::max_type > + { + typedef abstract::op<I, arith::impl::max_type > ret; + }; + + namespace arith { + + namespace impl { + + template <class I> + struct max_type : abstract::op<I, max_type > + { + mlc::box<const I> input1_; + mlc::box<const I> input2_; + + max_type(const abstract::non_vectorial_image& input1, + const abstract::non_vectorial_image& input2) : + input1_(input1.exact()), + input2_(input2.exact()) + {} + + void impl_run() + { + precondition(input1_->size() == input2_->size()); + I output(input1_->size()); + oln_type_of(I, fwd_piter) p(input1_->size()); + + for_all(p) + output[p] = ntg::max((*input1_)[p].value(), (*input2_)[p].value()); + + *this->image_ = output; + } + + }; + + } + + template <typename I> + impl::max_type max(const abstract::non_vectorial_image& input1, + const abstract::non_vectorial_image& input2) + { + impl::max_type tmp(input1, input2); + tmp.run(); + return tmp; + } + + } + +} // end of namespace oln + + +#endif // ! OLENA_ARITH_MAX_HH

1 0

proto-1.0 89: Add delta to structuring element and neighborhood
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * oln/core/abstract/struct_elt.hh: Add delta to stucturing element.. * oln/core/2d/window2d.hh: Likewise. * oln/core/abstract/neighborhood.hh: Add delta to neighborhood. * oln/core/2d/neighborhood2d.hh: Likewise. * oln/core/accum.hh: New. Accumulator needed for delta computation. * oln/makefile.src: Add missing files. core/2d/neighborhood2d.hh | 25 ++++++++++------ core/2d/window2d.hh | 10 +++++- core/abstract/neighborhood.hh | 27 +++++++++++++++-- core/abstract/struct_elt.hh | 25 ++++++++++++++-- core/accum.hh | 65 ++++++++++++++++++++++++++++++++++++++++++ makefile.src | 4 ++ 6 files changed, 140 insertions(+), 16 deletions(-) Index: oln/core/abstract/struct_elt.hh --- oln/core/abstract/struct_elt.hh (revision 88) +++ oln/core/abstract/struct_elt.hh (working copy) @@ -32,6 +32,7 @@ # include <oln/core/properties.hh> # include <oln/core/coord.hh> # include <oln/core/abstract/dpoint.hh> +# include <oln/core/accum.hh> # include <vector> # include <utility> # include <iostream> @@ -132,6 +133,18 @@ return this->exact().impl_at(i); } + coord_t + get_delta() const + { + return this->exact().impl_get_delta(); + } + + coord_t + delta_update(const dpoint_type& dp) + { + return this->exact().impl_delta_update(dp); + } + protected: bool @@ -145,9 +158,16 @@ { if (!(impl_has(dp))) this->dp_.push_back(dp); + this->delta_update(dp); return this->exact(); } + coord_t + impl_get_delta() const + { + return delta_; + } + unsigned impl_card() const { @@ -161,15 +181,16 @@ return dp_[i]; } - struct_elt() : dp_() + struct_elt() : dp_(), delta_(0) {}; - struct_elt(unsigned size) : dp_() + struct_elt(unsigned size) : dp_(), delta_(0) { dp_.reserve(size); }; std::vector<dpoint_type> dp_; + max_accumulator<coord_t> delta_; }; } // end of abstract Index: oln/core/abstract/neighborhood.hh --- oln/core/abstract/neighborhood.hh (revision 88) +++ oln/core/abstract/neighborhood.hh (working copy) @@ -32,6 +32,7 @@ # include <oln/core/properties.hh> # include <oln/core/coord.hh> # include <oln/core/abstract/dpoint.hh> +# include <oln/core/accum.hh> # include <vector> # include <utility> # include <iostream> @@ -129,6 +130,18 @@ return this->exact().impl_at(i); } + coord_t + get_delta() const + { + return this->exact().impl_get_delta(); + } + + coord_t + delta_update(const dpoint_type& dp) + { + return this->exact().impl_delta_update(dp); + } + protected: bool @@ -145,9 +158,16 @@ this->dp_.push_back(dp); this->dp_.push_back(-dp); } + this->delta_update(dp); return this->exact(); } - + + coord_t + impl_get_delta() const + { + return delta_; + } + unsigned impl_card() const { @@ -161,15 +181,16 @@ return dp_[i]; } - neighborhood() : dp_() + neighborhood() : dp_(), delta_(0) {}; - neighborhood(unsigned size) : dp_() + neighborhood(unsigned size) : dp_(), delta_(0) { dp_.reserve(size); }; std::vector<dpoint_type> dp_; + max_accumulator<coord_t> delta_; }; } // end of abstract Index: oln/core/2d/neighborhood2d.hh --- oln/core/2d/neighborhood2d.hh (revision 88) +++ oln/core/2d/neighborhood2d.hh (working copy) @@ -31,7 +31,7 @@ # include <oln/core/2d/dpoint2d.hh> # include <oln/core/2d/size2d.hh> # include <oln/core/coord.hh> -# include <oln/core/abstract/struct_elt.hh> +# include <oln/core/abstract/neighborhood.hh> namespace oln { @@ -39,29 +39,30 @@ // category template <> - struct set_category< neighborhood2d > { typedef category::struct_elt ret; }; + struct set_category< neighborhood2d > + { typedef category::neighborhood ret; }; // super_type template <> struct set_super_type< neighborhood2d > { - typedef abstract::struct_elt< neighborhood2d > ret; + typedef abstract::neighborhood< neighborhood2d > ret; }; template <> - struct set_props< category::struct_elt, neighborhood2d > - : public props_of<category::struct_elt> + struct set_props< category::neighborhood, neighborhood2d > + : public props_of<category::neighborhood> { typedef dpoint2d dpoint_type; typedef size2d size_type; }; - class neighborhood2d : public abstract::struct_elt<neighborhood2d> + class neighborhood2d : public abstract::neighborhood<neighborhood2d> { public: - typedef abstract::struct_elt<neighborhood2d> super_type; + typedef abstract::neighborhood<neighborhood2d> super_type; /*! ** \brief Construct a neighborhood of 2 dimensions. @@ -108,8 +109,14 @@ return std::string("neighborhood2d"); } - protected: - + coord_t + impl_delta_update(const dpoint_type& dp) + { + delta_(abs(dp.row())); + delta_(abs(dp.col())); + return delta_; + } + }; inline const neighborhood2d& Index: oln/core/2d/window2d.hh --- oln/core/2d/window2d.hh (revision 88) +++ oln/core/2d/window2d.hh (working copy) @@ -110,8 +110,14 @@ return std::string("window2d"); } - protected: - + coord_t + impl_delta_update(const dpoint_type& dp) + { + delta_(abs(dp.row())); + delta_(abs(dp.col())); + return delta_; + } + }; Index: oln/core/accum.hh --- oln/core/accum.hh (revision 0) +++ oln/core/accum.hh (revision 0) @@ -0,0 +1,65 @@ +// Copyright (C) 2001, 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_CORE_ACCUM_HH +# define OLENA_CORE_ACCUM_HH + + +/*! \class max_accumulator +** +** This is a \a functor. It saves the maximum T value +** that has been passed as an argument of its operator() +** method. To retrieve the value saved, just use the +** max_accumulator as a T instance. +*/ + +template <class T> +struct max_accumulator +{ + + max_accumulator (T t) : acc_(t) + {} + + void + operator()(T t) + { + if (t > acc_) + acc_ = t; + } + + operator T() const + { + return acc_; + } + +private: + + T acc_; + +}; + +#endif // OLENA_CORE_ACCUM_HH Index: oln/makefile.src --- oln/makefile.src (revision 88) +++ oln/makefile.src (working copy) @@ -11,6 +11,9 @@ basics3d.hh \ config/pconf.hh \ config/system.hh \ + convert/value_to_point.hh \ + morpho/stat.hh \ + morpho/erosion.hh \ core/1d/array1d.hh \ core/1d/dpoint1d.hh \ core/1d/fwd_piter1d.hh \ @@ -50,6 +53,7 @@ core/abstract/size.hh \ core/abstract/neighborhood.hh \ core/abstract/struct_elt.hh \ + core/accum.hh \ core/coord.hh \ core/compose.hh \ core/properties.hh \

1 0

proto-1.0 88: Activate testsuite for convert and morpho
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * configure.ac: Activate test suite for morpho and convert. configure.ac | 2 ++ 1 files changed, 2 insertions(+) Index: configure.ac --- configure.ac (revision 86) +++ configure.ac (working copy) @@ -99,9 +99,11 @@ # OLN_TESTS defined in oln-local.m4 OLN_TESTS([sanity]) + OLN_TESTS([convert]) OLN_TESTS([core]) OLN_TESTS([io]) OLN_TESTS([level]) + OLN_TESTS([morpho]) OLN_TESTS([utils]) AC_SUBST([OLN_TESTS_SUBDIRS])

1 0

proto-1.0 87: Add temporary conversion of value to point
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * oln/convert/value_to_point.hh: New. Conversion of value to point. This is temporary: waiting for a coherent `convert' in this prototype. * tests/convert/tests/value_to_point: New. * tests/convert: New. * tests/convert/tests: New. * tests/convert/Makefile.am: New. * oln/convert: New. oln/convert/value_to_point.hh | 106 +++++++++++++++++++++++++++++++++++++ tests/convert/Makefile.am | 5 + tests/convert/tests/value_to_point | 27 +++++++++ 3 files changed, 138 insertions(+) Index: tests/convert/tests/value_to_point --- tests/convert/tests/value_to_point (revision 0) +++ tests/convert/tests/value_to_point (revision 0) @@ -0,0 +1,27 @@ + +#include <iostream> + +#include <ntg/basics.hh> +#include <ntg/color.hh> +#include <oln/convert/value_to_point.hh> + +#include "check.hh" +#include "data.hh" + +bool check() +{ + oln::point1d p = oln::convert::value_to_point(42u); + oln::point3d p2 = oln::convert::value_to_point(ntg::rgb_8(1,6,64)); + std::cout << p.index() << std::endl; + std::cout << p2.row() << std::endl; + std::cout << p2.col() << std::endl; + std::cout << p2.slice() << std::endl; + std::cout << p2 << std::endl; +fdas + if (p.index() != 42 || + p2.row() != 6 || + p2.col() != 64 || + p2.slice() != 1) + return true; + return false; +} Index: tests/convert/Makefile.am --- tests/convert/Makefile.am (revision 0) +++ tests/convert/Makefile.am (revision 0) @@ -0,0 +1,5 @@ +include ../check/Makefile.runtests +include ../check/Makefile.check + + +check-local: check-runtests Index: oln/convert/value_to_point.hh --- oln/convert/value_to_point.hh (revision 0) +++ oln/convert/value_to_point.hh (revision 0) @@ -0,0 +1,106 @@ +// Copyright (C) 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_VALUE_TO_POINT +# define OLENA_VALUE_TO_POINT + +# include <ntg/basics.hh> + +# include <oln/core/1d/point1d.hh> +# include <oln/core/3d/point3d.hh> + +// FIXME: waiting for a coherent convert. + +namespace oln { + + namespace convert { + /*! Convert a value of pixel to a point. + ** + ** For example, transform an RGB color to a 3D point + ** (ntg::rgb_8 => oln::point3d). + ** This function is useful to build the histogram. \n + ** Example: + ** \verbatim + ** f(oln::convert::value_to_point<ntg::rgb_8>()(ntg::rgb_8(1,6,64))); + ** // is equivalent to: + ** f(oln::point3d(1, 6, 64)); + ** \endverbatim + */ + + namespace internal { + + template <typename ValueType> + struct value_to_point + { + typedef ValueType value_type; + typedef point1d point_type; + + value_to_point(const value_type& v) + : v_(v) + {} + + operator point_type() + { + return point1d(v_ - ntg_min_val(value_type)); + } + + const value_type& v_; + }; + + template <unsigned Qbits, template <unsigned> class S> + struct value_to_point< ntg::color<3, Qbits, S> > + { + typedef ntg::color<3, Qbits, S> value_type; + typedef point3d point_type; + + value_to_point(const value_type& v) + : v_(v) + {} + + operator point_type() + { + return point_type(coord_t(v_[0]), coord_t(v_[1]), coord_t(v_[2])); + } + + const value_type& v_; + }; + + + } + + template <typename ValueType> + typename internal::value_to_point<ValueType>::point_type + value_to_point(const ValueType& v) + { + return internal::value_to_point<ValueType>(v); + } + + } // convert + +} // oln + +#endif

1 0

proto-1.0 86: Add morpho min and max
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * oln/morpho: New. * oln/morpho/stat.hh: New. Add morphological min and max. * oln/core/abstract/piter.hh: Add `for_all_remaining' macro. * oln/core/abstract/point.hh: Fix a return type (was abstract). * tests/morpho: New. Add test suite for morpho. * tests/morpho/tests: New. Likewise. * tests/morpho/Makefile.am: New. Likewise. * oln/basics.hh: Add headers. oln/basics.hh | 2 oln/core/abstract/piter.hh | 6 + oln/core/abstract/point.hh | 4 oln/morpho/stat.hh | 185 +++++++++++++++++++++++++++++++++++++++++++++ tests/morpho/Makefile.am | 1 5 files changed, 196 insertions(+), 2 deletions(-) Index: tests/morpho/Makefile.am --- tests/morpho/Makefile.am (revision 0) +++ tests/morpho/Makefile.am (revision 0) @@ -0,0 +1 @@ +include ../check/Makefile.runtests Index: oln/basics.hh --- oln/basics.hh (revision 85) +++ oln/basics.hh (working copy) @@ -57,8 +57,10 @@ # include <oln/core/abstract/size.hh> # include <oln/core/abstract/point.hh> # include <oln/core/abstract/images.hh> +# include <oln/core/abstract/entry.hh> # include <oln/core/abstract/piter.hh> +# include <oln/core/abstract/witer.hh> // # include <oln/core/abstract/niter.hh> # include <oln/core/abstract/struct_elt.hh> Index: oln/morpho/stat.hh --- oln/morpho/stat.hh (revision 0) +++ oln/morpho/stat.hh (revision 0) @@ -0,0 +1,185 @@ +// Copyright (C) 2001, 2002, 2003, 2004, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_MORPHO_STAT_HH +# define OLENA_MORPHO_STAT_HH + +# include <oln/basics.hh> +# include <ntg/bin.hh> + +namespace oln { + + namespace morpho { + + namespace internal { + + /*! + ** \brief Min and Max on a structuring element. + ** + ** We need to use this inner definition in order to specialize + ** max and min on binary images. + ** + ** \param I Image exact type. + ** \param E Structuring element type. + ** \param V Associated value type. + */ + template <class I, class E, class V =oln_type_of(I, value)> + struct stat_ + { + /*! + ** \brief Maximum of a structuring element. + ** + ** Look for the maximum in the structuring element se disposed + ** on the image input, at the point p. + ** + ** \arg input Input image. + ** \arg p Point of the image to move the structuring element on. + ** \arg se The structuring element to use. + */ + static V + max(const I& input, const oln_type_of(I, point)& p, const E& se) + { + // FIXME: test dim I == dim E + + oln_type_of(E, fwd_witer) dp(se); + dp.start(); + V val = input[p + dp]; + for_all_remaining (dp) + if (val < input[p + dp].value()) + val = input[p + dp].value(); + return val; + } + + /*! + ** \brief Minimum of a structuring element. + ** + ** Look for the minimum in the structuring element se disposed + ** on the image input, at the point p. + ** + ** \arg input Input image. + ** \arg p Point of the image to move the structuring element on. + ** \arg se The structuring element to use. + */ + static V + min(const I& input, const oln_type_of(I, point)& p, const E& se) + { + // FIXME: test dim I == dim E + + oln_type_of(E, fwd_witer) dp(se); + dp.start(); + V val = input[p + dp]; + for_all_remaining (dp) + if (val > input[p + dp].value()) + val = input[p + dp].value(); + return val; + } + + }; + + /* Binary specialization. */ + + template <class I, class E> + struct stat_<I, E, ntg::bin> + { + static ntg::bin + max(const I& input, const oln_type_of(I, point)& p, const E& se) + { + mlc::eq<I::dim, E::dim>::ensure(); + oln_type_of(E, fwd_witer) dp(se); + for_all (dp) + if (input[p + dp] == true) + return true; + return false; + } + + static ntg::bin + min(const I& input, const oln_type_of(I, point)& p, const E& se) + { + mlc::eq<I::dim, E::dim>::ensure(); + oln_type_of(E, fwd_witer) dp(se); + for_all (dp) + if (input[p + dp] == false) + return false; + return true; + } + + }; + + } // internal + + /*! + ** \brief Maximum of a structuring element. + ** + ** Look for the maximum in the structuring element se disposed + ** on the image input, at the point p. + ** + ** \param I Image exact type. + ** \param E Structuring element type. + ** + ** \arg input Input image. + ** \arg p Point of the image to move the structuring element on. + ** \arg se The structuring element to use. + */ + template<class I, class E> + oln_type_of(I, value) + max(const abstract::non_vectorial_image& input, + const oln_type_of(I, point)& p, + const abstract::struct_elt<E>& se) + { + // FIXME: test dim I == dim E + return internal::stat_<I, E>::max(input.exact(), p, se.exact()); + } + + /*! ** \brief Minimum of a structuring element. + ** + ** Look for the minimum in the structuring element se disposed + ** on the image input, at the point p. + ** + ** \param I Image exact type. + ** \param E Structuring element type. + ** + ** \arg input Input image. + ** \arg p Point of the image to move the structuring element on. + ** \arg se The structuring element to use. + */ + template<class I, class E> + oln_type_of(I, value) + min(const abstract::non_vectorial_image& input, + const oln_type_of(I, point)& p, + // const mlc_exact_type(I)::iter_type& p, + const abstract::struct_elt<E>& se) + { + // FIXME: test dim I == dim E + return internal::stat_<I, E>::min(input.exact(), p, se.exact()); + } + + } // end of namespace morpho + +} // end of namespace oln + + +#endif // ! OLENA_MORPHO_STAT_HH Index: oln/core/abstract/piter.hh --- oln/core/abstract/piter.hh (revision 85) +++ oln/core/abstract/piter.hh (working copy) @@ -37,9 +37,13 @@ -# define for_all(p) for(p.start(); p.is_valid(); p.next()) +# define for_all(p) \ + for(p.start(); p.is_valid(); p.next()) +# define for_all_remaining(p) \ + for(; p.is_valid(); p.next()) + namespace oln { // fwd decl Index: oln/core/abstract/point.hh --- oln/core/abstract/point.hh (revision 85) +++ oln/core/abstract/point.hh (working copy) @@ -87,6 +87,8 @@ struct point : public mlc::any__best_memory<E> { + typedef E exact_type; + /*! \brief Test equality of two points. Nota bene: this method ** is abstract-like. ** @@ -111,7 +113,7 @@ typedef oln_type_of(E, dpoint) dpoint_type; - const point operator+(const dpoint_type& dp) const + const exact_type operator+(const dpoint_type& dp) const { return this->exact().impl_plus(dp); }

1 0

proto-1.0 85: Add window iterator (missing file)
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * oln/core/2d/fwd_witer2d.hh: New. Forward window iterator in 2d. fwd_witer2d.hh | 68 +++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 files changed, 68 insertions(+) Index: oln/core/2d/fwd_witer2d.hh --- oln/core/2d/fwd_witer2d.hh (revision 0) +++ oln/core/2d/fwd_witer2d.hh (revision 0) @@ -0,0 +1,68 @@ +// Copyright (C) 2001, 2003, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_CORE_FWD_WITER2D_HH +# define OLENA_CORE_FWD_WITER2D_HH + +# include <oln/core/abstract/witer.hh> +# include <oln/core/2d/window2d.hh> + +# include <string> + +namespace oln { + + // fwd decl + struct fwd_witer2d; + + // category + template <> + struct set_category<fwd_witer2d> { typedef category::witer ret; }; + + // props + template <> + struct set_props < category::witer, fwd_witer2d > : public props_of<category::witer> + { + typedef window2d se_type; + }; + + struct fwd_witer2d : public abstract::witer< fwd_witer2d > + { + + typedef abstract::witer<fwd_witer2d> super_type; + + fwd_witer2d(const se_type& se) : + super_type(se) + { + this->exact_ptr = this; + this->invalidate(); + } + + }; + +} // oln + +#endif // OLENA_CORE_FWD_WITER2D_HH

1 0

proto-1.0 84: Add windows iterator
by Simon Odou 20 Mar '05

20 Mar '05

Index: ChangeLog from Simon Odou <simon(a)lrde.epita.fr> * tests/core/tests/fwd_witer2d: New. Test for forward window iterator in 2d. * oln/makefile.src: Add new headers. * oln/core/abstract/struct_elt.hh: Fix a bad prototype. * oln/core/abstract/neighborhood.hh: Likewise. * oln/core/abstract/witer.hh: New. Add window iterator. * oln/core/2d/window2d.hh: Add `fwd_witer' property.. * oln/core/properties.hh: Likewise. * oln/basics2d.hh: Add new headers. oln/basics2d.hh | 1 oln/core/2d/window2d.hh | 2 oln/core/abstract/neighborhood.hh | 2 oln/core/abstract/struct_elt.hh | 5 + oln/core/abstract/witer.hh | 142 ++++++++++++++++++++++++++++++++++++++ oln/core/properties.hh | 3 oln/makefile.src | 2 tests/core/tests/fwd_witer2d | 25 ++++++ 8 files changed, 180 insertions(+), 2 deletions(-) Index: tests/core/tests/fwd_witer2d --- tests/core/tests/fwd_witer2d (revision 0) +++ tests/core/tests/fwd_witer2d (revision 0) @@ -0,0 +1,25 @@ +#include <oln/core/2d/fwd_witer2d.hh> +#include <oln/core/abstract/piter.hh> + +#include <ntg/all.hh> + +#include <cassert> + +using namespace oln; +using namespace ntg; + +template <typename E> +bool foo(const E& w) +{ + oln_type_of(E, fwd_witer) i(w); + unsigned cpt = 0; + for_all(i) + ++cpt; + return not (cpt == w.card()); +} + +int main() +{ + window2d w = win_c8p(); + return foo(w); +} Index: oln/makefile.src --- oln/makefile.src (revision 83) +++ oln/makefile.src (working copy) @@ -25,6 +25,7 @@ core/2d/point2d.hh \ core/2d/size2d.hh \ core/2d/neighborhood2d.hh \ + core/2d/fwd_witer2d.hh \ core/2d/window2d.hh \ core/3d/array3d.hh \ core/3d/fwd_piter3d.hh \ @@ -44,6 +45,7 @@ core/abstract/internal/image_impl.hh \ core/abstract/op.hh \ core/abstract/piter.hh \ + core/abstract/witer.hh \ core/abstract/point.hh \ core/abstract/size.hh \ core/abstract/neighborhood.hh \ Index: oln/core/abstract/struct_elt.hh --- oln/core/abstract/struct_elt.hh (revision 83) +++ oln/core/abstract/struct_elt.hh (working copy) @@ -60,12 +60,14 @@ mlc_decl_prop(category::struct_elt, dpoint_type); mlc_decl_prop(category::struct_elt, size_type); + mlc_decl_prop(category::struct_elt, fwd_witer_type); static void echo(std::ostream& ostr) { ostr << "props_of( category::struct_elt, " << typeid(type).name() << ") = {" << " size_type = " << typeid(size_type).name() << "," << std::endl + << " fwd_witer_type = " << typeid(fwd_witer_type).name() << "," << std::endl << " dpoint_type = " << typeid(dpoint_type).name() << " }" << std::endl; } @@ -73,6 +75,7 @@ mlc_register_prop(category::struct_elt, dpoint_type); mlc_register_prop(category::struct_elt, size_type); + mlc_register_prop(category::struct_elt, fwd_witer_type); namespace abstract { @@ -100,7 +103,7 @@ } bool - has(const abstract::dpoint<dpoint_type>& dp) const + has(const dpoint_type& dp) const { return this->exact().impl_has(dp.exact()); } Index: oln/core/abstract/witer.hh --- oln/core/abstract/witer.hh (revision 0) +++ oln/core/abstract/witer.hh (revision 0) @@ -0,0 +1,142 @@ +// Copyright (C) 2001, 2003, 2005 EPITA Research and Development Laboratory +// +// This file is part of the Olena Library. This library is free +// software; you can redistribute it and/or modify it under the terms +// of the GNU General Public License version 2 as published by the +// Free Software Foundation. +// +// This library is distributed in the hope that it will be useful, +// but WITHOUT ANY WARRANTY; without even the implied warranty of +// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU +// General Public License for more details. +// +// You should have received a copy of the GNU General Public License +// along with this library; see the file COPYING. If not, write to +// the Free Software Foundation, 59 Temple Place - Suite 330, Boston, +// MA 02111-1307, USA. +// +// As a special exception, you may use this file as part of a free +// software library without restriction. Specifically, if other files +// instantiate templates or use macros or inline functions from this +// file, or you compile this file and link it with other files to +// produce an executable, this file does not by itself cause the +// resulting executable to be covered by the GNU General Public +// License. This exception does not however invalidate any other +// reasons why the executable file might be covered by the GNU General +// Public License. + +#ifndef OLENA_CORE_ABSTRACT_WITER_HH +# define OLENA_CORE_ABSTRACT_WITER_HH + +# include <mlc/any.hh> +# include <mlc/types.hh> + +# include <oln/core/properties.hh> + +# include <string> + +namespace oln { + + // fwd decl + namespace abstract { + template <typename E> struct witer; + } + + // category + template <typename E> + struct set_category< abstract::witer<E> > { typedef category::witer ret; }; + + /// properties of any type in category::witer + // + template <typename type> + struct props_of < category::witer, type > + { + typedef mlc::true_type user_defined_; + + mlc_decl_prop(category::witer, se_type); + + static void echo(std::ostream& ostr) + { + ostr << "props_of( category::witer, " + << typeid(type).name() << ") = {" + << " se_type = " << typeid(se_type).name() + << " }" << std::endl; + } + }; + + mlc_register_prop(category::witer, se_type); + + namespace abstract { + + template <typename E> + struct witer : public mlc::any__best_speed<E> + { + + typedef witer<E> self_type; + + typedef oln_type_of(E, se) se_type; + typedef oln_type_of(se_type, dpoint) dpoint_type; + + void start() + { + this->exact().impl_start(); + } + + void next() + { + precondition(this->is_valid()); + this->exact().impl_next(); + } + + bool is_valid() const + { + return this->exact().impl_is_valid(); + } + + operator dpoint_type() const + { + precondition(this->is_valid()); + return this->se_[pos_]; + } + + void invalidate() + { + this->exact().impl_invalidate(); + postcondition(! this->is_valid()); + } + + protected: + + void impl_start() + { + pos_ = 0; + } + + void impl_next() + { + ++pos_; + } + + bool impl_is_valid() const + { + return pos_ != se_.card(); + } + + void impl_invalidate() + { + pos_ = se_.card(); + } + + witer(const se_type& se) + : se_(se), pos_(0) + {} + + const se_type& se_; + unsigned pos_; + }; + + } // abstract + +} // oln + +#endif // OLENA_CORE_ABSTRACT_WITER_HH Index: oln/core/abstract/neighborhood.hh --- oln/core/abstract/neighborhood.hh (revision 83) +++ oln/core/abstract/neighborhood.hh (working copy) @@ -100,7 +100,7 @@ } bool - has(const abstract::dpoint<dpoint_type>& dp) const + has(const dpoint_type& dp) const { return this->exact().impl_has(dp.exact()); } Index: oln/core/2d/window2d.hh --- oln/core/2d/window2d.hh (revision 83) +++ oln/core/2d/window2d.hh (working copy) @@ -36,6 +36,7 @@ namespace oln { class window2d; // forward declaration + class fwd_witer2d; // category template <> @@ -54,6 +55,7 @@ { typedef dpoint2d dpoint_type; typedef size2d size_type; + typedef fwd_witer2d fwd_witer_type; }; class window2d : public abstract::struct_elt<window2d> Index: oln/core/properties.hh --- oln/core/properties.hh (revision 83) +++ oln/core/properties.hh (working copy) @@ -46,6 +46,7 @@ struct point; struct size; struct piter; + struct witer; struct niter; // FIXME: ... } @@ -60,9 +61,11 @@ struct dpoint_type; struct fwd_piter_type; struct bkd_piter_type; + struct fwd_witer_type; struct iter_type; struct delegated_type; struct size_type; + struct se_type; struct image_constness_type; struct image_dimension_type; Index: oln/basics2d.hh --- oln/basics2d.hh (revision 83) +++ oln/basics2d.hh (working copy) @@ -35,6 +35,7 @@ # include <oln/core/2d/size2d.hh> # include <oln/core/2d/point2d.hh> # include <oln/core/2d/image2d.hh> +# include <oln/core/2d/fwd_witer2d.hh> # include <oln/core/2d/window2d.hh> # include <oln/core/2d/neighborhood2d.hh>

1 0