$extrastylesheet
Olena  User documentation 2.1
An Image Processing Platform
 All Classes Namespaces Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
distance_and_closest_point_geodesic.hh
1 // Copyright (C) 2008, 2009, 2012 EPITA Research and Development
2 // Laboratory (LRDE)
3 //
4 // This file is part of Olena.
5 //
6 // Olena is free software: you can redistribute it and/or modify it under
7 // the terms of the GNU General Public License as published by the Free
8 // Software Foundation, version 2 of the License.
9 //
10 // Olena is distributed in the hope that it will be useful,
11 // but WITHOUT ANY WARRANTY; without even the implied warranty of
12 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 // General Public License for more details.
14 //
15 // You should have received a copy of the GNU General Public License
16 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
17 //
18 // As a special exception, you may use this file as part of a free
19 // software project without restriction. Specifically, if other files
20 // instantiate templates or use macros or inline functions from this
21 // file, or you compile this file and link it with other files to produce
22 // an executable, this file does not by itself cause the resulting
23 // executable to be covered by the GNU General Public License. This
24 // exception does not however invalidate any other reasons why the
25 // executable file might be covered by the GNU General Public License.
26 
27 #ifndef MLN_TRANSFORM_DISTANCE_AND_CLOSEST_POINT_GEODESIC_HH
28 # define MLN_TRANSFORM_DISTANCE_AND_CLOSEST_POINT_GEODESIC_HH
29 
33 
34 # include <mln/canvas/distance_geodesic.hh>
35 # include <mln/transform/internal/closest_point_functor.hh>
36 
37 # include <mln/data/fill.hh>
38 # include <mln/util/couple.hh>
39 
40 
41 namespace mln
42 {
43 
44  namespace transform
45  {
46 
61  template <typename I, typename N, typename D>
62  util::couple<mln_ch_value(I,D), mln_ch_value(I,mln_psite(I))>
63  distance_and_closest_point_geodesic(const Image<I>& input,
64  const Neighborhood<N>& nbh,
65  D max);
66 
67 
83  template <typename P, typename N, typename D>
84  util::couple<mln_image_from_grid(mln_grid(P),D),
85  mln_image_from_grid(mln_grid(P),unsigned)>
86  distance_and_closest_point_geodesic(const p_array<P>& pset,
87  const box<P>& closest_point_domain,
88  const Neighborhood<N>& nbh,
89  D max);
90 
91 
92 # ifndef MLN_INCLUDE_ONLY
93 
94  template <typename I, typename N, typename D>
95  inline
96  util::couple<mln_ch_value(I,D), mln_ch_value(I, mln_psite(I))>
97  distance_and_closest_point_geodesic(const Image<I>& input,
98  const Neighborhood<N>& nbh,
99  D max)
100  {
101  mln_trace("transform::distance_closest_point_geodesic");
102 
103  mln_precondition(exact(input).is_valid());
104  mln_precondition(exact(nbh).is_valid());
105 
106  internal::closest_point_functor_with_sites<I> f;
107  mln_ch_value(I,D) dmap = mln::canvas::distance_geodesic(input, nbh,
108  max, f);
109 
110  return make::couple(dmap, f.cp_ima);
111  }
112 
113  template <typename P, typename N, typename D>
114  inline
115  util::couple<mln_image_from_grid(mln_grid(P),D), mln_image_from_grid(mln_grid(P),unsigned)>
116  distance_and_closest_point_geodesic(const p_array<P>& pset,
117  const box<P>& closest_point_domain,
118  const Neighborhood<N>& nbh,
119  D max)
120  {
121  mln_trace("transform::distance_and_closest_point_geodesic");
122 
123  mln_precondition(closest_point_domain.is_valid());
124  mln_precondition(exact(nbh).is_valid());
125  mln_precondition(geom::bbox(pset) <= closest_point_domain);
126 
127  typedef mln_image_from_grid(mln_grid(P), bool) I;
128  internal::closest_point_functor_with_indexes<I> f(pset);
129 
130  I ima(closest_point_domain);
131  data::fill(ima, false);
132  data::fill((ima | pset).rw(), true);
133 
134  mln_ch_value(I,D) dmap = mln::canvas::distance_geodesic(ima, nbh, max, f);
135 
136  return make::couple(dmap, f.cp_ima);
137  }
138 
139 # endif // ! MLN_INCLUDE_ONLY
140 
141  } // end of namespace mln::transform
142 
143 } // end of namespace mln
144 
145 
146 #endif // ! MLN_TRANSFORM_DISTANCE_AND_CLOSEST_POINT_GEODESIC_HH