DGtal 1.4.2
Loading...
Searching...
No Matches
PlaneProbingRNeighborhood.ih
1/**
2 * This program is free software: you can redistribute it and/or modify
3 * it under the terms of the GNU Lesser General Public License as
4 * published by the Free Software Foundation, either version 3 of the
5 * License, or (at your option) any later version.
6 *
7 * This program is distributed in the hope that it will be useful,
8 * but WITHOUT ANY WARRANTY; without even the implied warranty of
9 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
10 * GNU General Public License for more details.
11 *
12 * You should have received a copy of the GNU General Public License
13 * along with this program. If not, see <http://www.gnu.org/licenses/>.
14 *
15 **/
16
17/**
18 * @file
19 * @author Jocelyn Meyron (\c jocelyn.meyron@liris.cnrs.fr )
20 * Laboratoire d'InfoRmatique en Image et Systemes d'information - LIRIS (CNRS, UMR 5205), CNRS, France
21 *
22 * @date 2020/12/04
23 *
24 * Implementation of inline methods defined in PlaneProbingRNeighborhood.h
25 *
26 * This file is part of the DGtal library.
27 */
28
29
30//////////////////////////////////////////////////////////////////////////////
31#include <cstdlib>
32//////////////////////////////////////////////////////////////////////////////
33
34///////////////////////////////////////////////////////////////////////////////
35// IMPLEMENTATION of inline methods.
36///////////////////////////////////////////////////////////////////////////////
37
38///////////////////////////////////////////////////////////////////////////////
39// ----------------------- Standard services ------------------------------
40
41// ------------------------------------------------------------------------
42template < typename TPredicate >
43inline
44DGtal::PlaneProbingRNeighborhood<TPredicate>::
45PlaneProbingRNeighborhood(Predicate const& aPredicate, Point const& aQ, Triangle const& aM)
46 : DGtal::PlaneProbingNeighborhood<TPredicate>(aPredicate, aQ, aM)
47{}
48
49// ------------------------------------------------------------------------
50template < typename TPredicate >
51inline
52DGtal::PlaneProbingRNeighborhood<TPredicate>::
53~PlaneProbingRNeighborhood()
54{}
55
56///////////////////////////////////////////////////////////////////////////////
57// ----------------------- Plane Probing services ------------------------------
58
59// ------------------------------------------------------------------------
60template < typename TPredicate >
61inline
62typename DGtal::PlaneProbingRNeighborhood<TPredicate>::HexagonState
63DGtal::PlaneProbingRNeighborhood<TPredicate>::hexagonState ()
64{
65 this->myCandidates.clear();
66
67 std::array< bool, 6 > state = { false, false, false, false, false, false };
68 for (int i = 0; i < 6; ++i)
69 {
70 PointOnProbingRay r = this->myNeighborhood[i].getBase();
71
72 if (this->isNeighbor(r))
73 {
74 state[i] = this->myPredicate(this->absolutePoint(r));
75
76 if (state[i])
77 {
78 PointOnProbingRay pt = closestPointOnRayLogWithPredicate(r.getBase());
79
80 assert(pt == closestPointOnRayLinearWithPredicate(r.getBase()));
81 assert(this->myPredicate(this->absolutePoint(pt)));
82
83 this->myCandidates.push_back(pt);
84 }
85 }
86 }
87
88 return this->classify(state);
89}
90
91// ------------------------------------------------------------------------
92template < typename TPredicate >
93template < typename TPointAdapter >
94inline
95TPointAdapter
96DGtal::PlaneProbingRNeighborhood<TPredicate>::closestPointOnRayLogWithPredicate (TPointAdapter const& aRay) const
97{
98 assert(this->myPredicate(this->absolutePoint(aRay)));
99
100 // Exponential march
101 TPointAdapter Xk = aRay, Xl = aRay.next(1);
102 while (this->myPredicate(this->absolutePoint(Xk)) &&
103 this->isSmallest(this->relativePoint(Xk), this->relativePoint(Xl))) {
104 Integer d = Xl.position() - Xk.position();
105 Xk = Xl;
106 Xl = Xl.next(2 * d);
107 }
108 Xk = Xk.previous(Integer((Xl.position() - Xk.position()) / 2));
109
110 // Binary search
111 Integer d = Xl.position() - Xk.position();
112 while (d > 4) {
113 assert(this->myPredicate(this->absolutePoint(Xk)));
114
115 TPointAdapter Xalpha = Xk.next(Integer(d / 4)),
116 Xbeta = Xk.next(Integer(d / 2)),
117 Xgamma = Xk.next(Integer(3*d/4));
118
119 assert(Xk.position() < Xalpha.position() && Xalpha.position() < Xbeta.position() &&
120 Xbeta.position() < Xgamma.position() && Xgamma.position() < Xl.position());
121
122 if (this->myPredicate(this->absolutePoint(Xbeta)) &&
123 this->isSmallest(this->relativePoint(Xbeta), this->relativePoint(Xgamma))) {
124 Xk = Xbeta;
125 } else if (! this->myPredicate(this->absolutePoint(Xalpha)) ||
126 this->isSmallest(this->relativePoint(Xbeta), this->relativePoint(Xalpha))) {
127 Xl = Xbeta;
128 } else {
129 Xk = Xalpha;
130 Xl = Xgamma;
131 }
132
133 d = Xl.position() - Xk.position();
134 }
135
136 return closestPointOnRayLinearWithPredicate(Xk);
137}
138
139// ------------------------------------------------------------------------
140template < typename TPredicate >
141template < typename TPointAdapter >
142inline
143TPointAdapter
144DGtal::PlaneProbingRNeighborhood<TPredicate>::closestPointOnRayLinearWithPredicate (TPointAdapter const& aRay) const
145{
146 assert(this->myPredicate(this->absolutePoint(aRay)));
147
148 TPointAdapter previousX = aRay, currentX = previousX.next(1);
149 while (this->myPredicate(this->absolutePoint(currentX)) &&
150 this->isSmallest(this->relativePoint(previousX), this->relativePoint(currentX))) {
151 previousX = currentX;
152 currentX = previousX.next(1);
153 }
154
155 assert(this->myPredicate(this->absolutePoint(previousX)));
156
157 return previousX;
158}
159
160///////////////////////////////////////////////////////////////////////////////
161// Interface - public :
162
163/**
164 * Writes/Displays the object on an output stream.
165 * @param out the output stream where the object is written.
166 */
167template <typename TPredicate>
168inline
169void
170DGtal::PlaneProbingRNeighborhood<TPredicate>::selfDisplay ( std::ostream & out ) const
171{
172 out << "[PlaneProbingRNeighborhood]";
173}
174
175/**
176 * Checks the validity/consistency of the object.
177 * @return 'true' if the object is valid, 'false' otherwise.
178 */
179template <typename TPredicate>
180inline
181bool
182DGtal::PlaneProbingRNeighborhood<TPredicate>::isValid() const
183{
184 return true;
185}
186
187
188
189///////////////////////////////////////////////////////////////////////////////
190// Implementation of inline functions //
191
192template <typename TPredicate>
193inline
194std::ostream&
195DGtal::operator<< ( std::ostream & out,
196 const PlaneProbingRNeighborhood<TPredicate> & object )
197{
198 object.selfDisplay( out );
199 return out;
200}
201
202// //
203///////////////////////////////////////////////////////////////////////////////
204
205