Commit 480e6ac5 authored by Nicolas Neri's avatar Nicolas Neri
Browse files

Add mirroring for borders.

	* border_mirror.cc: Test file for border mirroring.
	* border_mirror.hh: Fix border mirroring.

git-svn-id: https://svn.lrde.epita.fr/svn/oln/trunk@1121 4aad255d-cdde-0310-9447-f3009e2ae8c0
parent 3d4a2797
2007-09-17 Guillaume Duhamel <guillaume.duhamel@lrde.epita.fr>
Add mirroring for borders.
* border_mirror.cc: Test file for border mirroring.
* border_mirror.hh: Fix border mirroring.
2007-09-17 Matthieu Garrigues <garrigues@lrde.epita.fr>
Improve pnm support.
......
// Copyright (C) 2007 EPITA Research and Development Laboratory
//
// This file is part of the Olena Library. This library is free
// software; you can redistribute it and/or modify it under the terms
// of the GNU General Public License version 2 as published by the
// Free Software Foundation.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this library; see the file COPYING. If not, write to
// the Free Software Foundation, 51 Franklin Street, Fifth Floor,
// Boston, MA 02111-1307, USA.
//
// As a special exception, you may use this file as part of a free
// software library without restriction. Specifically, if other files
// instantiate templates or use macros or inline functions from this
// file, or you compile this file and link it with other files to
// produce an executable, this file does not by itself cause the
// resulting executable to be covered by the GNU General Public
// License. This exception does not however invalidate any other
// reasons why the executable file might be covered by the GNU General
// Public License.
/*! \file tests/border_fill.cc
*
* \brief Tests on mln::border::fill.
*/
#include "border_mirror.hh"
#include <mln/core/image2d_b.hh>
#include <mln/value/int_u8.hh>
#include <mln/pw/all.hh>
#include <mln/debug/println_with_border.hh>
using namespace mln;
int
main (void)
{
image2d_b<value::int_u8> i1(5, 7);
// Fill with randomized value.
for (unsigned int i = 0; i < i1.ncells (); ++i)
i1[i] = i;//(i * 4452) % 10;
std::cout << "before mirror"
<< std::endl
<< std::endl;
debug::println_with_border(i1);
border::mirror (i1);
std::cout << "after mirror"
<< std::endl
<< std::endl;
debug::println_with_border(i1);
}
// Copyright (C) 2007 EPITA Research and Development Laboratory
//
// This file is part of the Olena Library. This library is free
// software; you can redistribute it and/or modify it under the terms
// of the GNU General Public License version 2 as published by the
// Free Software Foundation.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this library; see the file COPYING. If not, write to
// the Free Software Foundation, 51 Franklin Street, Fifth Floor,
// Boston, MA 02111-1307, USA.
//
// As a special exception, you may use this file as part of a free
// software library without restriction. Specifically, if other files
// instantiate templates or use macros or inline functions from this
// file, or you compile this file and link it with other files to
// produce an executable, this file does not by itself cause the
// resulting executable to be covered by the GNU General Public
// License. This exception does not however invalidate any other
// reasons why the executable file might be covered by the GNU General
// Public License.
#ifndef MLN_BORDER_MIRROR_HH
# define MLN_BORDER_MIRROR_HH
/*! \file mln/border/mirror.hh
*
* \brief FIXME.
*/
# include <mln/core/concept/image.hh>
# include <mln/core/internal/fixme.hh>
# include <mln/level/memset_.hh>
//# include <mln/core/line_piter.hh>
#include <mln/geom/nrows.hh>
#include <mln/geom/ncols.hh>
#include <mln/core/image2d_b.hh>
#include <mln/core/image2d_b.hh>
#include <mln/core/pixel.hh>
namespace mln
{
namespace border
{
/*! Mirror the virtual (outer) border of image \p ima with the
* (inner) level contents of this image.
*
* \param[in,out] ima The image whose border is to be mirrored.
*
* \pre \p ima has to be initialized.
*
* \todo Implement it + optimize with memset if possible.
*/
template <typename I>
void mirror(const Fast_Image<I>& ima);
# ifndef MLN_INCLUDE_ONLY
template <typename I>
void mirror(const Fast_Image<I>& ima_)
{
const I& ima = exact(ima_);
mln_precondition(ima.has_data());
// internal::fixme();
// FIX
std::size_t border = ima.border ();
std::size_t nbrows = geom::max_row(ima) - geom::min_row(ima);
std::size_t nbcols = geom::max_col(ima) - geom::min_col(ima);
std::size_t real_nbcols = (nbcols + 1) + 2 * border;
std::size_t start = real_nbcols * border + border;
std::size_t s = start;
// duplicate top left corner
for (std::size_t i = 0; i < border; ++i)
for (std::size_t j = 0; j < border; ++j)
const_cast<I&>(ima)[i * ((nbcols + 1) + 2 * border) + j] = ima[s];
// duplicate top left corner
s = start + nbcols;
for (std::size_t i = 0; i < border; ++i)
for (std::size_t j = 1; j <= border; ++j)
const_cast<I&>(ima)[i * ((nbcols + 1) + 2 * border) + (nbcols + border + j)] = ima[s];
// duplicate bottom left corner
s = start + (nbrows * real_nbcols);
for (std::size_t i = 1; i <= border; ++i)
for (std::size_t j = 1; j <= border; ++j)
const_cast<I&>(ima)[s - i + (j * (real_nbcols))] = ima[s];
// duplicate bottom right corner
s = start + (nbrows * real_nbcols) + nbcols;
for (std::size_t i = 1; i <= border; ++i)
for (std::size_t j = 1; j <= border; ++j)
const_cast<I&>(ima)[s + i + (j * real_nbcols)] = ima[s];
// mirror top border
s = start;
for (std::size_t i = 0; i <= nbcols; ++i)
for (std::size_t j = 1; j <= border; ++j)
const_cast<I&>(ima)[s + i - (j * real_nbcols)] = ima[s + i + ((j - 1)* real_nbcols)];
// mirror left border
s = start;
for (std::size_t i = 0; i <= nbrows; ++i)
for (std::size_t j = 0; j <= border; ++j)
const_cast<I&>(ima)[s + (i * real_nbcols) - j] = ima[s + (i * real_nbcols) + (j - 1)];
// mirror right border
s = start;
for (std::size_t i = 0; i <= nbrows; ++i)
for (std::size_t j = 1; j <= border; ++j)
const_cast<I&>(ima)[s + (i * real_nbcols + nbcols) + j] = ima[s + (i * real_nbcols + nbcols) - (j - 1)];
// mirror bottom border
s = start + (nbrows * real_nbcols);
for (std::size_t i = 0; i <= nbcols; ++i)
for (std::size_t j = 1; j <= border; ++j)
const_cast<I&>(ima)[s + i + (j * real_nbcols)] = ima[s + i - ((j - 1)* real_nbcols)];
//END FIX
}
# endif // ! MLN_INCLUDE_ONLY
} // end of namespace mln::border
} // end of namespace mln
#endif // ! MLN_BORDER_MIRROR_HH
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment