pixter3d.hh

// Copyright (C) 2007, 2008, 2009 EPITA Research and Development Laboratory (LRDE)
//
// This file is part of Olena.
//
// Olena is free software: you can redistribute it and/or modify it under
// the terms of the GNU General Public License as published by the Free
// Software Foundation, version 2 of the License.
//
// Olena 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 Olena.  If not, see <http://www.gnu.org/licenses/>.
//
// As a special exception, you may use this file as part of a free
// software project 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_PIXTER3D_HH
# define MLN_CORE_PIXTER3D_HH


# include <mln/core/internal/pixel_iterator_base.hh>
# include <mln/core/alias/point3d.hh>
# include <mln/geom/size3d.hh>
# include <mln/opt/at.hh>


namespace mln
{

  /*------------------.
  | fwd_pixter3d<I>.  |
  `------------------*/

  template <typename I>
  class fwd_pixter3d
    : public internal::forward_pixel_iterator_base_< I, fwd_pixter3d<I> >
  {
    typedef internal::forward_pixel_iterator_base_< I, fwd_pixter3d<I> > super_;

  public:
    typedef I image;

    fwd_pixter3d(I& image);

    void next_();

    void start_();

  private:

    const unsigned border_x2_;

    const unsigned row_offset_;

    mln_qlf_value(I)* eor_;

    const unsigned next_sli_offset_;

    const unsigned next_srow_offset_;

    const unsigned sli_offset_;

    mln_qlf_value(I)* eos_;

    using super_::image_;
  };


  /*------------------.
  | bkd_pixter3d<I>.  |
  `------------------*/

  template <typename I>
  class bkd_pixter3d
    : public internal::backward_pixel_iterator_base_< I, bkd_pixter3d<I> >
  {
    typedef internal::backward_pixel_iterator_base_< I, bkd_pixter3d<I> > super_;

  public:
    typedef I image;

    bkd_pixter3d(I& image);

    void next_();

    void start_();

  private:

    const unsigned border_x2_;

    const unsigned row_offset_;

    mln_qlf_value(I)* bor_;

    const unsigned next_sli_offset_;

    const unsigned next_srow_offset_;

    const unsigned sli_offset_;

    mln_qlf_value(I)* bos_;

    using super_::image_;
  };


#ifndef MLN_INCLUDE_ONLY

  /*------------------.
  | fwd_pixter3d<I>.  |
  `------------------*/

  template <typename I>
  inline
  fwd_pixter3d<I>::fwd_pixter3d(I& image)
    : super_(image),
      border_x2_(2 * image.border()),
      row_offset_(image.bbox().ncols() + border_x2_),
      next_sli_offset_(row_offset_ * border_x2_ + border_x2_),
      next_srow_offset_(next_sli_offset_ + image.bbox().ncols()),
      sli_offset_((image.bbox().ncols() + border_x2_) *
                  (image.bbox().nrows() + border_x2_))
  {
    mln_precondition(image.is_valid());
  }

  template <typename I>
  inline
  void
  fwd_pixter3d<I>::start_()
  {
    eor_ = & opt::at(image_,
                     geom::min_sli(image_),
                     geom::min_row(image_),
                     geom::max_col(image_)) + 1;
    eos_ = & opt::at(image_,
                     geom::min_sli(image_),
                     geom::max_row(image_),
                     geom::max_col(image_)) + 1;
  }

  template <typename I>
  inline
  void
  fwd_pixter3d<I>::next_()
  {
    ++this->value_ptr_;
    if (this->value_ptr_ == eos_ && this->value_ptr_ != this->eoi_)
    {
      this->value_ptr_ += next_sli_offset_;
      eos_ += sli_offset_;
      eor_ += next_srow_offset_;
    }
    else if (this->value_ptr_ == eor_ && this->value_ptr_ != this->eoi_)
    {
      this->value_ptr_ += border_x2_;
      eor_ += row_offset_;
    }
  }


  /*------------------.
  | bkd_pixter3d<I>.  |
  `------------------*/

  template <typename I>
  inline
  bkd_pixter3d<I>::bkd_pixter3d(I& image)
    : super_(image),
      border_x2_(2 * image.border()),
      row_offset_(image.bbox().ncols() + border_x2_),
      next_sli_offset_(row_offset_ * border_x2_ + border_x2_),
      next_srow_offset_(next_sli_offset_ + image.bbox().ncols()),
      sli_offset_((image.bbox().ncols() + border_x2_) *
                  (image.bbox().nrows() + border_x2_))
  {
    mln_precondition(image.is_valid());
  }

  template <typename I>
  inline
  void
  bkd_pixter3d<I>::start_()
  {
    bor_ = & opt::at(image_,
                     geom::max_sli(image_),
                     geom::max_row(image_),
                     geom::min_col(image_)) - 1;
    bos_ = & opt::at(image_,
                     geom::max_sli(image_),
                     geom::min_row(image_),
                     geom::min_col(image_)) - 1;
  }

  template <typename I>
  inline
  void
  bkd_pixter3d<I>::next_()
  {
    --this->value_ptr_;
    if (this->value_ptr_ == bos_ && this->value_ptr_ != this->boi_)
    {
      this->value_ptr_ -= next_sli_offset_;
      bos_ -= sli_offset_;
      bor_ -= next_srow_offset_;
    }
    else if (this->value_ptr_ == bor_ && this->value_ptr_ != this->boi_)
    {
      this->value_ptr_ -= border_x2_;
      bor_ -= row_offset_;
    }
  }

#endif // ! MLN_INCLUDE_ONLY

} // end of namespace mln


#endif // ! MLN_CORE_PIXTER3D_HH