border_duplicate.hh 4.63 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
// 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_DUPLICATE_HH
# define MLN_BORDER_DUPLICATE_HH

/*! \file mln/border/duplicate.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
  {

    /*! Assign the virtual (outer) border of image \p ima with the
     *  dupplicate of the inner border of this image.
     *
     * \param[in,out] ima The image whose border is to be duplicated.
     *
     * \pre \p ima has to be initialized.
     *
     * \todo Implement it + optimize with memcpy if possible.
     */
    template <typename I>
    void duplicate(const Fast_Image<I>& ima);


# ifndef MLN_INCLUDE_ONLY

    template <typename I>
    void duplicate(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 + 1; ++i)
	for (std::size_t j = 0; j < border + 1; ++j)
	  const_cast<I&>(ima)[i * ((nbcols + 1) + 2 * border) + j] = ima[s];

      // duplicate top border
      s = start;
      for (std::size_t i = 1; i <= nbcols - 1; ++i)
	for (std::size_t j = 0; j <= border; ++j)
	  const_cast<I&>(ima)[s + i - (j * real_nbcols)] = ima[s + i];

      // duplicate top left corner
      s = start + nbcols;
      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) + (nbcols + border + j)] = ima[s];

      // duplicate left border
      s = start;
      for (std::size_t i = 1; i <= nbrows - 1; ++i)
	for (std::size_t j = 1; j <= border; ++j)
	  const_cast<I&>(ima)[s - j + (i * real_nbcols)] = ima[s + (i * real_nbcols)];

      // duplicate right border
      s = start;
      for (std::size_t i = 1; i <= nbrows - 1; ++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)];


      // duplicate bottom left corner
      s = start + (nbrows * real_nbcols);
      for (std::size_t i = 0; i <= border; ++i)
	for (std::size_t j = 0; j <= border; ++j)
	  const_cast<I&>(ima)[s - i + (j * (real_nbcols))] = ima[s];

      // duplicate bottom border
      s = start + (nbrows * real_nbcols);
      for (std::size_t i = 1; i <= nbcols - 1; ++i)
	for (std::size_t j = 0; j <= border; ++j)
	  const_cast<I&>(ima)[s + i + (j * real_nbcols)] = ima[s + i];

      // duplicate bottom right corner
      s = start + (nbrows * real_nbcols) + nbcols;
      for (std::size_t i = 0; i <= border; ++i)
	for (std::size_t j = 0; j <= border; ++j)
	  const_cast<I&>(ima)[s + i + (j * real_nbcols)] = ima[s];



      //END FIX
    }

# endif // ! MLN_INCLUDE_ONLY

  } // end of namespace mln::border

} // end of namespace mln


#endif // ! MLN_BORDER_DUPLICATE_HH