DGtal 1.4.0
Loading...
Searching...
No Matches
testChamferDT.cpp
Go to the documentation of this file.
1
31#include <iostream>
32#include <set>
33#include <algorithm>
34#include "DGtal/base/Common.h"
35#include "ConfigTest.h"
36#include "DGtal/helpers/StdDefs.h"
37#include "DGtal/geometry/volumes/distance/ChamferNorm2D.h"
38#include "DGtal/geometry/volumes/distance/CMetricSpace.h"
39#include "DGtal/geometry/volumes/distance/CSeparableMetric.h"
41
42using namespace DGtal;
43using namespace Z2i;
44
46// Functions for testing class ChamferDT.
48
49
50
51bool checkCMetricConcept()
52{
53 BOOST_CONCEPT_ASSERT(( concepts::CMetricSpace<experimental::ChamferNorm2D<Z2i::Space > > ));
54 BOOST_CONCEPT_ASSERT(( concepts::CSeparableMetric<experimental::ChamferNorm2D<Z2i::Space > > ));
55 return true;
56}
57
58
59bool testChamferSimple()
60{
61 unsigned int nbok = 0;
62 unsigned int nb = 0;
63 trace.beginBlock ( "Testing simple chamfer mask localization...");
64
65 experimental::ChamferNorm2D<Z2i::Space>::LessThanAngular lthan;
66 experimental::ChamferNorm2D<Z2i::Space>::LessOrEqThanAngular lqthan;
67
68 experimental::ChamferNorm2D<Space>::Directions dirs34;
69 experimental::ChamferNorm2D<Space>::Directions normals34;
70
71 //3-4 mask
72 dirs34.push_back(Z2i::Vector(0,-1));
73 dirs34.push_back(Z2i::Vector(1,-1));
74 dirs34.push_back(Z2i::Vector(1,0));
75 dirs34.push_back(Z2i::Vector(1,1));
76 normals34.push_back(Z2i::Vector(1,-3));
77 normals34.push_back(Z2i::Vector(3,-1));
78 normals34.push_back(Z2i::Vector(3,1));
79 normals34.push_back(Z2i::Vector(1,3));
80
81 experimental::ChamferNorm2D<Space> aMask(dirs34,normals34);
82
83
84
85
86 //Testing cone discovery
87 Vector d(101,100);
88 trace.info() << "Direction "<<d<<std::endl;
89 trace.info() << " -> cone "<< *aMask.getCone(d)
90 << " -- " << *(aMask.getCone(d)+1) <<std::endl;
91 nbok += ( lthan(d, *(aMask.getCone(d)+1))) ? 1 : 0;
92 nb++;
93 trace.info() << "(" << nbok << "/" << nb << ") "
94 << "(d < cone+1) returns true" << std::endl;
95 nbok += (lqthan(*(aMask.getCone(d)),d)) ? 1 : 0;
96 nb++;
97 trace.info() << "(" << nbok << "/" << nb << ") "
98 << "(cone <= d) returns true" << std::endl;
99
100 //testing canonical
101 Vector d3(-101,100);
102 trace.info() << "Direction "<<d3<<std::endl;
103 trace.info() << " -> cone "<< *aMask.getCone(d3)
104 << " -- " << *(aMask.getCone(d3)+1) <<std::endl;
105 nbok += ( lthan(aMask.canonicalRay(d3), *(aMask.getCone(d3)+1))) ? 1 : 0;
106 nb++;
107 trace.info() << "(" << nbok << "/" << nb << ") "
108 << "(d < cone+1) returns true" << std::endl;
109 nbok += (lqthan(*(aMask.getCone(d3)),aMask.canonicalRay(d3))) ? 1 : 0;
110 nb++;
111 trace.info() << "(" << nbok << "/" << nb << ") "
112 << "(cone <= d) returns true" << std::endl;
113
114
115 //aligned
116 Vector d2(8,8);
117 trace.info() << "Direction "<<d2<<std::endl;
118 trace.info() << " -> cone "<< *aMask.getCone(d2)
119 << " -- " << *(aMask.getCone(d2)+1) <<std::endl;
120 if ((aMask.getCone(d2)+1) != aMask.end())
121 {
122 nbok += ( lthan(d2, *(aMask.getCone(d2)+1))) ? 1 : 0;
123 nb++;
124 trace.info() << "(" << nbok << "/" << nb << ") "
125 << "(d2 < cone+1) returns true" << std::endl;
126 }
127 nbok += (lqthan(*(aMask.getCone(d2)),d2) ) ? 1 : 0;
128 nb++;
129 trace.info() << "(" << nbok << "/" << nb << ") "
130 << "(cone <= d2) returns true" << std::endl;
131
132
133 //Distance test
134 Point p(1,1), q(102,101);
135 trace.info()<< "Distance between "<<p<<" and "<< q << " = "<< aMask(p,q)<<std::endl;
136 trace.info()<< "Distance between "<<q<<" and "<< p << " = "<< aMask(q,p)<<std::endl;
137 nbok += ( aMask(p,q) == aMask(q,p)) ? 1 : 0;
138 nb++;
139 trace.info() << "(" << nbok << "/" << nb << ") "
140 << "d(p,q)==d(q,p)" << std::endl;
141
142
143 //Testing Last cone:
144 Vector d4(1,10);
145 trace.info() << "Direction "<<d4<<std::endl;
146 trace.info() << " -> cone "<< *aMask.getCone(d4)
147 << " -- " << *(aMask.getCone(d4)+1) <<std::endl;
148 nbok += ( *aMask.getCone(d4) == Point(1,1)) ? 1 : 0;
149 nb++;
150 trace.info() << "(" << nbok << "/" << nb << ") "
151 << "Last cone 1" << std::endl;
152 trace.info() << " -> normal "<< aMask.getNormalFromCone(aMask.getCone(d4)) <<std::endl;;
153 nbok += (aMask.getNormalFromCone(aMask.getCone(d4)) == Point(1,3)) ? 1 : 0;
154 nb++;
155 trace.info() << "(" << nbok << "/" << nb << ") "
156 << "Last cone 2" << std::endl;
157
158 //Testing Last cone:
159 Vector d5(0,10);
160 trace.info() << "Direction "<<d5<<std::endl;
161 trace.info() << " -> cone "<< *aMask.getCone(d5)
162 << " -- " << *(aMask.getCone(d5)+1) <<std::endl;
163 nbok += ( *aMask.getCone(d5) == Point(1,1)) ? 1 : 0;
164 nb++;
165 trace.info() << "(" << nbok << "/" << nb << ") "
166 << "Last cone 1" << std::endl;
167 trace.info() << " -> normal "<< aMask.getNormalFromCone(aMask.getCone(d5)) <<std::endl;;
168 nbok += (aMask.getNormalFromCone(aMask.getCone(d5)) == Point(1,3)) ? 1 : 0;
169 nb++;
170 trace.info() << "(" << nbok << "/" << nb << ") "
171 << "Last cone 2" << std::endl;
172
173 trace.info() << " upper="<<*std::upper_bound(aMask.begin(), aMask.end(), d4, lthan)<<std::endl;
174
175 trace.endBlock();
176 return nbok == nb;
177}
178
179bool testBasicMasks()
180{
181 unsigned int nbok = 0;
182 unsigned int nb = 0;
183 trace.beginBlock ( "Testing basic chamfer masks...");
184
185 experimental::ChamferNorm2D<Space>::Directions dirs34;
186 experimental::ChamferNorm2D<Space>::Directions normals34;
187
188 //3-4 mask
189 dirs34.push_back(Z2i::Vector(0,-1));
190 dirs34.push_back(Z2i::Vector(1,-1));
191 dirs34.push_back(Z2i::Vector(1,0));
192 dirs34.push_back(Z2i::Vector(1,1));
193 normals34.push_back(Z2i::Vector(1,-3));
194 normals34.push_back(Z2i::Vector(3,-1));
195 normals34.push_back(Z2i::Vector(3,1));
196 normals34.push_back(Z2i::Vector(1,3));
197
198 experimental::ChamferNorm2D<Space> mask34(dirs34,normals34);
199
200 Point orig(0,0);
201 Point p(3,0);
202 trace.info() <<"Distance "<<p<<" = "<< mask34(orig,p)<<std::endl;;
203 trace.info() <<"Cone "<< *(mask34.getCone(p))<< " "<< *(mask34.getCone(p)+1)<< std::endl;
204
205 Point q(3,1);
206 trace.info() <<"Distance "<<q<<" = "<< mask34(orig,q)<<std::endl;;
207 trace.info() <<"Cone "<< *(mask34.getCone(q))<< " "<< *(mask34.getCone(q)+1)<< std::endl;
208
209
210 trace.endBlock();
211 return nbok == nb;
212}
213
214bool testIntersection()
215{
216 unsigned int nbok = 0;
217 unsigned int nb = 0;
218 trace.beginBlock ( "Testing intersection...");
219
220 Point P(0,0),Q(7,3),Qb(7,-3), Q2(2, -10), Q3(2,10);
221 Point Lmin(10,-10), Lmax(10,10);
222 experimental::ChamferNorm2D<Space> mask(1);
223
224 trace.info() << "Intersection "<<P<<" "<<Q<<" = "
225 << mask.getLowerRayIntersection(P,Q,Lmin,Lmax, 1)<<std::endl;
226 nbok += ( mask.getLowerRayIntersection(P,Q,Lmin,Lmax, 1) == 4) ? 1 : 0;
227 nb++;
228 trace.info() << "(" << nbok << "/" << nb << ") "
229 << "inter==4" << std::endl;
230
231 trace.info() << "Intersection "<<P<<" "<<Qb<<" = "
232 << mask.getLowerRayIntersection(P,Qb,Lmin,Lmax, 1)<<std::endl;
233 nbok += ( mask.getLowerRayIntersection(P,Qb,Lmin,Lmax, 1) == -5) ? 1 : 0;
234 nb++;
235 trace.info() << "(" << nbok << "/" << nb << ") "
236 << "inter==-5" << std::endl;
237
238 trace.info() << "Intersection "<<P<<" "<<Q2<<" = "
239 << mask.getLowerRayIntersection(P,Q2,Lmin,Lmax, 1)<<std::endl;
240 nbok += ( mask.getLowerRayIntersection(P,Q2,Lmin,Lmax, 1) <= Lmin[1]) ? 1 : 0;
241 nb++;
242 trace.info() << "(" << nbok << "/" << nb << ") "
243 << "inter==Lmin" << std::endl;
244
245 trace.info() << "Intersection "<<P<<" "<<Q3<<" = "
246 << mask.getLowerRayIntersection(P,Q3,Lmin,Lmax, 1)<<std::endl;
247 nbok += ( mask.getLowerRayIntersection(P,Q3,Lmin,Lmax, 1) >= Lmax[1]) ? 1 : 0;
248 nb++;
249 trace.info() << "(" << nbok << "/" << nb << ") "
250 << "inter==Lmax" << std::endl;
251
252 trace.info()<<std::endl;
253
254 Point PP(0,0),QQ(3,7),QQb(-3,7), QQ2(-10, 2), QQ3(10, 2);
255 Point LLmin(-10,10), LLmax(10,10);
256
257 trace.info() << "Intersection "<<PP<<" "<<QQ<<" = "
258 << mask.getLowerRayIntersection(PP,QQ,LLmin,LLmax, 0)<<std::endl;
259 nbok += ( mask.getLowerRayIntersection(PP,QQ,LLmin,LLmax, 0) == 4) ? 1 : 0;
260 nb++;
261 trace.info() << "(" << nbok << "/" << nb << ") "
262 << "inter==4" << std::endl;
263
264 trace.info() << "Intersection "<<PP<<" "<<QQb<<" = "
265 << mask.getLowerRayIntersection(PP,QQb,LLmin,LLmax, 0)<<std::endl;
266 nbok += ( mask.getLowerRayIntersection(PP,QQb,LLmin,LLmax, 0) == -5) ? 1 : 0;
267 nb++;
268 trace.info() << "(" << nbok << "/" << nb << ") "
269 << "inter==-5" << std::endl;
270
271 trace.info() << "Intersection "<<PP<<" "<<QQ2<<" = "
272 << mask.getLowerRayIntersection(PP,QQ2,LLmin,LLmax, 0)<<std::endl;
273 nbok += ( mask.getLowerRayIntersection(PP,QQ2,LLmin,LLmax, 0) <= LLmin[0]) ? 1 : 0;
274 nb++;
275 trace.info() << "(" << nbok << "/" << nb << ") "
276 << "inter==Lmin" << std::endl;
277
278 trace.info() << "Intersection "<<PP<<" "<<QQ3<<" = "
279 << mask.getLowerRayIntersection(PP,QQ3,LLmin,LLmax, 0)<<std::endl;
280 nbok += ( mask.getLowerRayIntersection(PP,QQ3,LLmin,LLmax, 0) >= LLmax[0]) ? 1 : 0;
281 nb++;
282 trace.info() << "(" << nbok << "/" << nb << ") "
283 << "inter==Lmax" << std::endl;
284
285 return nbok == nb;
286}
287
288
289bool testShrink()
290{
291 unsigned int nbok = 0;
292 unsigned int nb = 0;
293 trace.beginBlock ( "Testing shrinking...");
294
295 //5-7-11 metic
296 typedef experimental::ChamferNorm2D<Z2i::Space> Metric;
297 Metric::Directions dirs5711;
298 Metric::Directions normals5711;
299 //5-7-11 mask
300 dirs5711.push_back(Z2i::Vector(0,-1));
301 dirs5711.push_back(Z2i::Vector(1,-2));
302 dirs5711.push_back(Z2i::Vector(1,-1));
303 dirs5711.push_back(Z2i::Vector(2,-1));
304 dirs5711.push_back(Z2i::Vector(1,0));
305 dirs5711.push_back(Z2i::Vector(2,1));
306 dirs5711.push_back(Z2i::Vector(1,1));
307 dirs5711.push_back(Z2i::Vector(1,2));
308
309 normals5711.push_back(Z2i::Vector(1,-5));
310 normals5711.push_back(Z2i::Vector(3,-4));
311 normals5711.push_back(Z2i::Vector(4,-3));
312 normals5711.push_back(Z2i::Vector(5,-1));
313 normals5711.push_back(Z2i::Vector(5,1));
314 normals5711.push_back(Z2i::Vector(4,3));
315 normals5711.push_back(Z2i::Vector(3,4));
316 normals5711.push_back(Z2i::Vector(1,5));
317
318 Metric mask5711(dirs5711,normals5711);
319
320 //Setting:
321 Point P(0,0);
322 Point Q(1,-2), QQ(1,2);
323 Point Lmin(10,0);
324 Point Lmax(10,10);
325 Metric::ConstIterator itBeg = mask5711.begin();
326 Metric::ConstIterator itEnd = mask5711.end();
327
328 Point midPoint, nextMidPoint;
329
330 Metric::ConstIterator cone = mask5711.shrinkPSubMask(itBeg, itEnd, P, Q, Lmin, Lmax, 1, midPoint, nextMidPoint);
331 trace.info() <<"Shrink returns the cone "<< *cone<<" " <<*(cone+1)<<std::endl;
332 trace.info() <<"MidPoint "<< midPoint<<" " <<nextMidPoint<<std::endl<<std::endl;
333 nbok += ( midPoint == Point(10,0)) ? 1 : 0;
334 nb++;
335 trace.info() << "(" << nbok << "/" << nb << ") "
336 << "midPoint = (10,0)" << std::endl;
337 nbok += ( nextMidPoint == Point(10,5)) ? 1 : 0;
338 nb++;
339 trace.info() << "(" << nbok << "/" << nb << ") "
340 << "nextMidPoint = (10,5)" << std::endl;
341
342
343 trace.warning()<<" Shrinking with QQ(1,2)"<<std::endl;
344
345 Metric::ConstIterator cone2 = mask5711.shrinkPSubMask(itBeg, itEnd, P, QQ, Lmin, Lmax, 1, midPoint, nextMidPoint);
346 trace.info() <<"Shrink returns the cone2 "<< *cone2<<" " <<*(cone2+1)<<std::endl;
347 trace.info() <<"MidPoint "<< midPoint<<" " <<nextMidPoint<<std::endl<<std::endl;
348 nbok += ( midPoint == Point(10,-5)) ? 1 : 0;
349 nb++;
350 trace.info() << "(" << nbok << "/" << nb << ") "
351 << "midPoint = (10,0)" << std::endl;
352 nbok += ( nextMidPoint == Point(10,0)) ? 1 : 0;
353 nb++;
354 trace.info() << "(" << nbok << "/" << nb << ") "
355 << "nextMidPoint = (10,5)" << std::endl;
356
357
358 trace.beginBlock("Testing Symmetry");
359 Point LLmin(-10,0);
360 Point LLmax(-10,10);
361 Point QQQ(-1,-2);
362 Metric::ConstIterator cone3 = mask5711.shrinkP(itBeg, itEnd, P, QQQ, LLmin, LLmax, 1, midPoint, nextMidPoint);
363 trace.info() <<"Shrink returns the cone3 "<< *cone3<<" " <<*(cone3+1)<<std::endl;
364 trace.info() <<"MidPoint "<< midPoint<<" " <<nextMidPoint<<std::endl<<std::endl;
365 nbok += ( midPoint == Point(-10,0)) ? 1 : 0;
366 nb++;
367 trace.info() << "(" << nbok << "/" << nb << ") "
368 << "midPoint = (-10,0)" << std::endl;
369 nbok += ( nextMidPoint == Point(-10,5)) ? 1 : 0;
370 nb++;
371 trace.info() << "(" << nbok << "/" << nb << ") "
372 << "nextMidPoint = (-10,5)" << std::endl;
373 trace.endBlock();
374
375 //horizontal
376 trace.beginBlock("Testing H/V symmetry");
377 Point LLLmin(0,10);
378 Point LLLmax(10,10);
379 Point Qhori(-2,1);
380 Metric::ConstIterator cone4 = mask5711.shrinkP(itBeg, itEnd, P, Qhori, LLLmin, LLLmax, 0, midPoint, nextMidPoint);
381 trace.info() <<"Shrink returns the cone3 "<< *cone4<<" " <<*(cone4+1)<<std::endl;
382 trace.info() <<"MidPoint "<< midPoint<<" " <<nextMidPoint<<std::endl<<std::endl;
383 nbok += ( *cone4 == *cone) ? 1 : 0;
384 nb++;
385 trace.info() << "(" << nbok << "/" << nb << ") "
386 << "horiz==vert" << std::endl;
387 nbok += ( midPoint == Point(0,10)) ? 1 : 0;
388 nb++;
389 trace.info() << "(" << nbok << "/" << nb << ") "
390 << "midPoint = (0,10)" << std::endl;
391 nbok += ( nextMidPoint == Point(5,10)) ? 1 : 0;
392 nb++;
393 trace.info() << "(" << nbok << "/" << nb << ") "
394 << "nextMidPoint = (5,10)" << std::endl;
395
396 Point LLLLmin(0,-10);
397 Point LLLLmax(10,-10);
398 Point QQhori(-2,-1);
399 Metric::ConstIterator cone5 = mask5711.shrinkP(itBeg, itEnd, P, QQhori, LLLLmin, LLLLmax, 0, midPoint, nextMidPoint);
400 trace.info() <<"Shrink returns the cone3 "<< *cone4<<" " <<*(cone4+1)<<std::endl;
401 trace.info() <<"MidPoint "<< midPoint<<" " <<nextMidPoint<<std::endl;
402 nbok += ( *cone5 == *cone) ? 1 : 0;
403 nb++;
404 trace.info() << "(" << nbok << "/" << nb << ") "
405 << "horiz==vert" << std::endl;
406 nbok += ( midPoint == Point(0,-10)) ? 1 : 0;
407 nb++;
408 trace.info() << "(" << nbok << "/" << nb << ") "
409 << "midPoint = (0,-10)" << std::endl;
410 nbok += ( nextMidPoint == Point(5,-10)) ? 1 : 0;
411 nb++;
412 trace.info() << "(" << nbok << "/" << nb << ") "
413 << "nextMidPoint = (5,-10)" << std::endl;
414 trace.endBlock();
415
416 trace.endBlock();
417 return nbok == nb;
418}
419bool testDoubleShrink()
420{
421 unsigned int nbok = 0;
422 unsigned int nb = 0;
423 trace.beginBlock ( "Testing Double Shrinking...");
424
425 //5-7-11 metic
426 typedef experimental::ChamferNorm2D<Z2i::Space> Metric;
427 Metric::Directions dirs5711;
428 Metric::Directions normals5711;
429 //5-7-11 mask
430 dirs5711.push_back(Z2i::Vector(0,-1));
431 dirs5711.push_back(Z2i::Vector(1,-2));
432 dirs5711.push_back(Z2i::Vector(1,-1));
433 dirs5711.push_back(Z2i::Vector(2,-1));
434 dirs5711.push_back(Z2i::Vector(1,0));
435 dirs5711.push_back(Z2i::Vector(2,1));
436 dirs5711.push_back(Z2i::Vector(1,1));
437 dirs5711.push_back(Z2i::Vector(1,2));
438
439 normals5711.push_back(Z2i::Vector(1,-5));
440 normals5711.push_back(Z2i::Vector(3,-4));
441 normals5711.push_back(Z2i::Vector(4,-3));
442 normals5711.push_back(Z2i::Vector(5,-1));
443 normals5711.push_back(Z2i::Vector(5,1));
444 normals5711.push_back(Z2i::Vector(4,3));
445 normals5711.push_back(Z2i::Vector(3,4));
446 normals5711.push_back(Z2i::Vector(1,5));
447
448 Metric mask5711(dirs5711,normals5711);
449
450 //Setting:
451 trace.beginBlock("Tessting shrinking using Q");
452 Point P(0,0);
453 Point Q(1,-2), QQ(1,2);
454 Point Lmin(10,-10);
455 Point Lmax(10,10);
456 Point midPointP,nextMidPointP;
457 Point midPointQ,nextMidPointQ;
458
459 Metric::ConstIterator itBeg = mask5711.begin();
460 Metric::ConstIterator itEnd = mask5711.end();
461
462 Metric::ConstIterator cone = mask5711.shrinkP(itBeg, itEnd, P, Q, Lmin, Lmax, 1, midPointP, nextMidPointP);
463 trace.info() <<" P - Shrink returns the cone "<< *cone<<" " <<*(cone+1)<<std::endl;
464 trace.info() <<" P - MidPoint "<< midPointP<<" " <<nextMidPointP<<std::endl<<std::endl;
465
466 Metric::ConstIterator cone2 = mask5711.shrinkP(itBeg, itEnd, Q, P, Lmin, Lmax, 1, midPointQ, nextMidPointQ);
467 trace.info() <<" Q - Shrink returns the cone "<< *cone2<<" " <<*(cone2+1)<<std::endl;
468 trace.info() <<" Q - MidPoint "<< midPointQ<<" " <<nextMidPointQ<<std::endl<<std::endl;
469 trace.info() << "Checking Voro cell" <<std::endl;
470
471 double dpmidp = mask5711(P,midPointP);
472 double dqmidp = mask5711(Q,midPointP);
473 trace.info() << ((dpmidp < dqmidp) ? "MidP closer to P" : "Mid closer to Q") << std::endl;
474 double dpnextmidp = mask5711(P,nextMidPointP);
475 double dqnextmidp = mask5711(Q,nextMidPointP);
476 trace.info() << ((dpnextmidp < dqnextmidp) ? "NextMidP closer to P" : "Mid closer to Q") << std::endl;
477 nbok += ( (dpmidp < dqmidp) != (dpnextmidp < dqnextmidp)) ? 1 : 0;
478 nb++;
479 trace.info() << "(" << nbok << "/" << nb << ") "
480 << " Shrink P contains the Voronoi Edge" << std::endl;
481
482 double dpmidq = mask5711(P,midPointQ);
483 double dqmidq = mask5711(Q,midPointQ);
484 trace.info() << ((dpmidq < dqmidq) ? "MidQ closer to P" : "Mid closer to Q") << std::endl;
485 double dpnextmidq = mask5711(P,nextMidPointQ);
486 double dqnextmidq = mask5711(Q,nextMidPointQ);
487 trace.info() << ((dpnextmidq < dqnextmidq) ? "NextMidP closer to P" : "Mid closer to Q") << std::endl;
488 nbok += ( (dpmidq < dqmidq) != (dpnextmidq < dqnextmidq)) ? 1 : 0;
489 nb++;
490 trace.info() << "(" << nbok << "/" << nb << ") "
491 << " Shrink Q contains the Voronoi Edge" << std::endl;
492
493 trace.beginBlock("Computing the lower");
494 Metric::Abscissa vorocell;
495 Metric::Vector normalP = mask5711.getNormalFromCone(cone);
496 Metric::Vector normalQ = mask5711.getNormalFromCone(cone2);
497 vorocell = static_cast<Metric::Abscissa>(floor((double) (P[1]*normalP[1] - Q[1]*normalQ[1] - (Lmin[0] - P[0])*normalP[0] + (Lmin[0] - Q[0])*normalQ[0]) /(normalP[1] - normalQ[1]) ));
498 Point voro(Lmin[0], vorocell);
499
500 trace.emphase() << " Lower Voronoi Cell Point "<< voro << " absc = "<<vorocell<<std::endl;
501 nbok += ( (voro >= midPointP) && (voro<nextMidPointP)) ? 1 : 0;
502 nb++;
503 trace.info() << "(" << nbok << "/" << nb << ") "
504 << " Voro cell in the cone of P" << std::endl;
505 nbok += ( (voro >= midPointQ) && (voro<nextMidPointQ)) ? 1 : 0;
506 nb++;
507 trace.info() << "(" << nbok << "/" << nb << ") "
508 << " Voro cell in the cone of Q" << std::endl;
509
510
511 nbok += ( (vorocell == mask5711.getLowerVoronoiEdgeAbscissa(P, Q, Lmin, Lmax, 1))) ? 1 : 0;
512 nb++;
513 trace.info() << "(" << nbok << "/" << nb << ") "
514 << " Voro cell == getLowerVoronoiEdgeAbscissa " << std::endl;
515
516 trace.endBlock();
517
518 trace.endBlock();
519
520
521 //Setting2
522 trace.beginBlock("Testing double shrinking on QQ");
523 Metric::ConstIterator coneQQ = mask5711.shrinkP(itBeg, itEnd, P, QQ, Lmin, Lmax, 1, midPointP, nextMidPointP);
524 trace.info() <<" P - Shrink returns the cone "<< *coneQQ<<" " <<*(coneQQ+1)<<std::endl;
525 trace.info() <<" P - MidPoint "<< midPointP<<" " <<nextMidPointP<<std::endl;
526 trace.info() <<" P - cone distance to P("<< mask5711(midPointP,P)<<","
527 << mask5711(nextMidPointP,P)<<")" <<std::endl;
528 trace.info() <<" P - cone distance to QQ("<< mask5711(midPointP,QQ)<<","
529 << mask5711(nextMidPointP,QQ)<<")" <<std::endl<<std::endl;
530
531 Metric::ConstIterator coneQQ2 = mask5711.shrinkP(itBeg, itEnd, QQ, P, Lmin, Lmax, 1, midPointQ, nextMidPointQ);
532 trace.info() <<" QQ - Shrink returns the cone "<< *coneQQ2<<" " <<*(coneQQ2+1)<<std::endl;
533 trace.info() <<" QQ - MidPoint "<< midPointQ<<" " <<nextMidPointQ<<std::endl;
534 trace.info() <<" QQ - cone distance to QQ("<< mask5711(midPointQ,QQ)<<","
535 << mask5711(nextMidPointQ,QQ)<<")" <<std::endl;
536 trace.info() <<" QQ - cone distance to P("<< mask5711(midPointQ,P)<<","
537 << mask5711(nextMidPointQ,P)<<")" <<std::endl;
538 trace.info() << "Checking Voro cell" <<std::endl<<std::endl;
539
540 dpmidp = mask5711(P,midPointP);
541 dqmidp = mask5711(QQ,midPointP);
542 trace.info() << ((dpmidp < dqmidp) ? "MidP closer to P" : "MidP closer to QQ") << std::endl;
543 dpnextmidp = mask5711(P,nextMidPointP);
544 dqnextmidp = mask5711(QQ,nextMidPointP);
545 trace.info() << ((dpnextmidp < dqnextmidp) ? "NextMidP closer to P" : "NextMidP closer to QQ") << std::endl;
546 nbok += ( (dpmidp < dqmidp) != (dpnextmidp < dqnextmidp)) ? 1 : 0;
547 nb++;
548 trace.info() << "(" << nbok << "/" << nb << ") "
549 << " Shrink P contains the Voronoi Edge" << std::endl;
550 trace.info()<< "Distances : ConeP<->P("<<dpmidp<<","<<dpnextmidp<<") ConeP<->Q("<<dqmidp<<","<<dqnextmidp<<")"<<std::endl;
551
552
553 dpmidq = mask5711(P,midPointQ);
554 dqmidq = mask5711(QQ,midPointQ);
555 trace.info() << ((dpmidq < dqmidq) ? "MidQ closer to P" : "MidQ closer to QQ") << std::endl;
556 dpnextmidq = mask5711(P,nextMidPointQ);
557 dqnextmidq = mask5711(QQ,nextMidPointQ);
558 trace.info() << ((dpnextmidq < dqnextmidq) ? "NextMidQ closer to P" : "NextMidQ closer to QQ") << std::endl;
559 nbok += ( (dpmidq < dqmidq) != (dpnextmidq < dqnextmidq)) ? 1 : 0;
560 nb++;
561 trace.info() << "(" << nbok << "/" << nb << ") "
562 << " Shrink QQ contains the Voronoi Edge" << std::endl;
563 trace.info()<< "Distances : ConeQ<->P("<<dpmidq<<","<<dpnextmidq<<") ConeQ<->QQ("<<dqmidq<<","<<dqnextmidq<<")"<<std::endl;
564 trace.endBlock();
565
566 trace.endBlock();
567
568 return nbok == nb;
569}
570
571
572bool testDoubleShrinkHorizontal()
573{
574 unsigned int nbok = 0;
575 unsigned int nb = 0;
576 trace.beginBlock ( "Testing Double Shrinking Horizontal...");
577
578 //5-7-11 metic
579 typedef experimental::ChamferNorm2D<Z2i::Space> Metric;
580 Metric::Directions dirs5711;
581 Metric::Directions normals5711;
582 //5-7-11 mask
583 dirs5711.push_back(Z2i::Vector(0,-1));
584 dirs5711.push_back(Z2i::Vector(1,-2));
585 dirs5711.push_back(Z2i::Vector(1,-1));
586 dirs5711.push_back(Z2i::Vector(2,-1));
587 dirs5711.push_back(Z2i::Vector(1,0));
588 dirs5711.push_back(Z2i::Vector(2,1));
589 dirs5711.push_back(Z2i::Vector(1,1));
590 dirs5711.push_back(Z2i::Vector(1,2));
591
592 normals5711.push_back(Z2i::Vector(1,-5));
593 normals5711.push_back(Z2i::Vector(3,-4));
594 normals5711.push_back(Z2i::Vector(4,-3));
595 normals5711.push_back(Z2i::Vector(5,-1));
596 normals5711.push_back(Z2i::Vector(5,1));
597 normals5711.push_back(Z2i::Vector(4,3));
598 normals5711.push_back(Z2i::Vector(3,4));
599 normals5711.push_back(Z2i::Vector(1,5));
600
601 Metric mask5711(dirs5711,normals5711);
602
603 //Setting:
604 trace.beginBlock("Testing shrinking using Q");
605 Point P(0,0);
606 Point Q(-2,1), QQ(2,1);
607 Point Lmin(-10,10);
608 Point Lmax(10,10);
609 Point midPointP,nextMidPointP;
610 Point midPointQ,nextMidPointQ;
611
612 Metric::ConstIterator itBeg = mask5711.begin();
613 Metric::ConstIterator itEnd = mask5711.end();
614
615 Metric::ConstIterator cone = mask5711.shrinkP(itBeg, itEnd, P, Q, Lmin, Lmax, 0, midPointP, nextMidPointP);
616 trace.info() <<" P - Shrink returns the cone "<< *cone<<" " <<*(cone+1)<<std::endl;
617 trace.info() <<" P - MidPoint "<< midPointP<<" " <<nextMidPointP<<std::endl<<std::endl;
618
619 Metric::ConstIterator cone2 = mask5711.shrinkP(itBeg, itEnd, Q, P, Lmin, Lmax, 0, midPointQ, nextMidPointQ);
620 trace.info() <<" Q - Shrink returns the cone "<< *cone2<<" " <<*(cone2+1)<<std::endl;
621 trace.info() <<" Q - MidPoint "<< midPointQ<<" " <<nextMidPointQ<<std::endl<<std::endl;
622 trace.info() << "Checking Voro cell" <<std::endl;
623
624
625 double dpmidp = mask5711(P,midPointP);
626 double dqmidp = mask5711(Q,midPointP);
627 trace.info() << ((dpmidp < dqmidp) ? "MidP closer to P" : "Mid closer to Q") << std::endl;
628 double dpnextmidp = mask5711(P,nextMidPointP);
629 double dqnextmidp = mask5711(Q,nextMidPointP);
630 trace.info() << ((dpnextmidp < dqnextmidp) ? "NextMidP closer to P" : "Mid closer to Q") << std::endl;
631 nbok += ( (dpmidp < dqmidp) != (dpnextmidp < dqnextmidp)) ? 1 : 0;
632 nb++;
633 trace.info() << "(" << nbok << "/" << nb << ") "
634 << " Shrink P contains the Voronoi Edge" << std::endl;
635
636
637 double dpmidq = mask5711(P,midPointQ);
638 double dqmidq = mask5711(Q,midPointQ);
639 trace.info() << ((dpmidq < dqmidq) ? "MidQ closer to P" : "Mid closer to Q") << std::endl;
640 double dpnextmidq = mask5711(P,nextMidPointQ);
641 double dqnextmidq = mask5711(Q,nextMidPointQ);
642 trace.info() << ((dpnextmidq < dqnextmidq) ? "NextMidP closer to P" : "Mid closer to Q") << std::endl;
643 nbok += ( (dpmidq < dqmidq) != (dpnextmidq < dqnextmidq)) ? 1 : 0;
644 nb++;
645 trace.info() << "(" << nbok << "/" << nb << ") "
646 << " Shrink Q contains the Voronoi Edge" << std::endl;
647 trace.endBlock();
648
649 //Setting2
650 trace.beginBlock("Testing double shrinking on QQ");
651 Metric::ConstIterator coneQQ = mask5711.shrinkP(itBeg, itEnd, P, QQ, Lmin, Lmax, 0, midPointP, nextMidPointP);
652 trace.info() <<" P - Shrink returns the cone "<< *coneQQ<<" " <<*(coneQQ+1)<<std::endl;
653 trace.info() <<" P - MidPoint "<< midPointP<<" " <<nextMidPointP<<std::endl;
654 trace.info() <<" P - cone distance to P("<< mask5711(midPointP,P)<<","
655 << mask5711(nextMidPointP,P)<<")" <<std::endl;
656 trace.info() <<" P - cone distance to QQ("<< mask5711(midPointP,QQ)<<","
657 << mask5711(nextMidPointP,QQ)<<")" <<std::endl<<std::endl;
658
659 Metric::ConstIterator coneQQ2 = mask5711.shrinkP(itBeg, itEnd, QQ, P, Lmin, Lmax, 0, midPointQ, nextMidPointQ);
660 trace.info() <<" QQ - Shrink returns the cone "<< *coneQQ2<<" " <<*(coneQQ2+1)<<std::endl;
661 trace.info() <<" QQ - MidPoint "<< midPointQ<<" " <<nextMidPointQ<<std::endl;
662 trace.info() <<" QQ - cone distance to QQ("<< mask5711(midPointQ,QQ)<<","
663 << mask5711(nextMidPointQ,QQ)<<")" <<std::endl;
664 trace.info() <<" QQ - cone distance to P("<< mask5711(midPointQ,P)<<","
665 << mask5711(nextMidPointQ,P)<<")" <<std::endl;
666 trace.info() << "Checking Voro cell" <<std::endl<<std::endl;
667
668 dpmidp = mask5711(P,midPointP);
669 dqmidp = mask5711(QQ,midPointP);
670 trace.info() << ((dpmidp < dqmidp) ? "MidP closer to P" : "MidP closer to QQ") << std::endl;
671 dpnextmidp = mask5711(P,nextMidPointP);
672 dqnextmidp = mask5711(QQ,nextMidPointP);
673 trace.info() << ((dpnextmidp < dqnextmidp) ? "NextMidP closer to P" : "NextMidP closer to QQ") << std::endl;
674 nbok += ( (dpmidp < dqmidp) != (dpnextmidp < dqnextmidp)) ? 1 : 0;
675 nb++;
676 trace.info() << "(" << nbok << "/" << nb << ") "
677 << " Shrink P contains the Voronoi Edge" << std::endl;
678 trace.info()<< "Distances : ConeP<->P("<<dpmidp<<","<<dpnextmidp<<") ConeP<->Q("<<dqmidp<<","<<dqnextmidp<<")"<<std::endl;
679
680
681 dpmidq = mask5711(P,midPointQ);
682 dqmidq = mask5711(QQ,midPointQ);
683 trace.info() << ((dpmidq < dqmidq) ? "MidQ closer to P" : "MidQ closer to QQ") << std::endl;
684 dpnextmidq = mask5711(P,nextMidPointQ);
685 dqnextmidq = mask5711(QQ,nextMidPointQ);
686 trace.info() << ((dpnextmidq < dqnextmidq) ? "NextMidQ closer to P" : "NextMidQ closer to QQ") << std::endl;
687 nbok += ( (dpmidq < dqmidq) != (dpnextmidq < dqnextmidq)) ? 1 : 0;
688 nb++;
689 trace.info() << "(" << nbok << "/" << nb << ") "
690 << " Shrink QQ contains the Voronoi Edge" << std::endl;
691 trace.info()<< "Distances : ConeQ<->P("<<dpmidq<<","<<dpnextmidq<<") ConeQ<->QQ("<<dqmidq<<","<<dqnextmidq<<")"<<std::endl;
692 trace.endBlock();
693
694 trace.endBlock();
695
696
697 return nbok == nb;
698}
699
700
701bool testHiddenBy()
702{
703 unsigned int nbok = 0;
704 unsigned int nb = 0;
705 trace.beginBlock ( "Testing HiddenBy...");
706
707 //5-7-11 metic
708 typedef experimental::ChamferNorm2D<Z2i::Space> Metric;
709 Metric::Directions dirs5711;
710 Metric::Directions normals5711;
711 //5-7-11 mask
712 dirs5711.push_back(Z2i::Vector(0,-1));
713 dirs5711.push_back(Z2i::Vector(1,-2));
714 dirs5711.push_back(Z2i::Vector(1,-1));
715 dirs5711.push_back(Z2i::Vector(2,-1));
716 dirs5711.push_back(Z2i::Vector(1,0));
717 dirs5711.push_back(Z2i::Vector(2,1));
718 dirs5711.push_back(Z2i::Vector(1,1));
719 dirs5711.push_back(Z2i::Vector(1,2));
720
721 normals5711.push_back(Z2i::Vector(1,-5));
722 normals5711.push_back(Z2i::Vector(3,-4));
723 normals5711.push_back(Z2i::Vector(4,-3));
724 normals5711.push_back(Z2i::Vector(5,-1));
725 normals5711.push_back(Z2i::Vector(5,1));
726 normals5711.push_back(Z2i::Vector(4,3));
727 normals5711.push_back(Z2i::Vector(3,4));
728 normals5711.push_back(Z2i::Vector(1,5));
729
730 Metric mask5711(dirs5711,normals5711);
731
732 //Setting:
733 Point U(0,0);
734 Point V(0,3);
735 Point VV(-10,3);
736 Point W(0,6);
737 Point Lmin(10,-10);
738 Point Lmax(10,10);
739
740 Point Usym(20,0);
741 Point Vsym(20,3);
742 Point VVsym(30,3);
743 Point Wsym(20,6);
744
745 Point LminH(-10,10);
746 Point LmaxH(10,10);
747
748 Point UH(0,0);
749 Point VH(3,0);
750 Point VVH(3,-10);
751 Point WH(6,0);
752
753
754 nbok += ( mask5711.hiddenBy(U,V, W, Lmin, Lmax, 1) == false ) ? 1 : 0;
755 nb++;
756 trace.info() << "(" << nbok << "/" << nb << ") "
757 << " H(U,V,W) is false" << std::endl;
758
759 nbok += ( mask5711.hiddenBy(U,VV, W, Lmin, Lmax, 1) == true ) ? 1 : 0;
760 nb++;
761 trace.info() << "(" << nbok << "/" << nb << ") "
762 << " H(U,VV,W) is true" << std::endl;
763
764 nbok += ( mask5711.hiddenBy(Usym,Vsym, Wsym, Lmin, Lmax, 1) == false ) ? 1 : 0;
765 nb++;
766 trace.info() << "(" << nbok << "/" << nb << ") "
767 << " H(Usym,Vsym,Wsym) is false" << std::endl;
768
769 nbok += ( mask5711.hiddenBy(Usym,VVsym, Wsym, Lmin, Lmax, 1) == true ) ? 1 : 0;
770 nb++;
771 trace.info() << "(" << nbok << "/" << nb << ") "
772 << " H(Usym,VVsym,Wsym) is true" << std::endl;
773
774 nbok += ( mask5711.hiddenBy(UH,VH, WH, LminH, LmaxH, 0) == false ) ? 1 : 0;
775 nb++;
776 trace.info() << "(" << nbok << "/" << nb << ") "
777 << " H(UH,VH,WH) is false" << std::endl;
778
779 nbok += ( mask5711.hiddenBy(UH,VVH, WH, LminH, LmaxH, 0) == true ) ? 1 : 0;
780 nb++;
781 trace.info() << "(" << nbok << "/" << nb << ") "
782 << " H(UH,VVH,WH) is true" << std::endl;
783
784 trace.endBlock();
785 return nb == nbok;
786}
787
788
789
790bool testDoubleOnSegment()
791{
792 unsigned int nbok = 0;
793 unsigned int nb = 0;
794
795 //5-7-11 metic
796 typedef experimental::ChamferNorm2D<Z2i::Space> Metric;
797 Metric::Directions dirs5711;
798 Metric::Directions normals5711;
799 //5-7-11 mask
800 dirs5711.push_back(Z2i::Vector(0,-1));
801 dirs5711.push_back(Z2i::Vector(1,-2));
802 dirs5711.push_back(Z2i::Vector(1,-1));
803 dirs5711.push_back(Z2i::Vector(2,-1));
804 dirs5711.push_back(Z2i::Vector(1,0));
805 dirs5711.push_back(Z2i::Vector(2,1));
806 dirs5711.push_back(Z2i::Vector(1,1));
807 dirs5711.push_back(Z2i::Vector(1,2));
808
809 normals5711.push_back(Z2i::Vector(1,-5));
810 normals5711.push_back(Z2i::Vector(3,-4));
811 normals5711.push_back(Z2i::Vector(4,-3));
812 normals5711.push_back(Z2i::Vector(5,-1));
813 normals5711.push_back(Z2i::Vector(5,1));
814 normals5711.push_back(Z2i::Vector(4,3));
815 normals5711.push_back(Z2i::Vector(3,4));
816 normals5711.push_back(Z2i::Vector(1,5));
817
818 Metric mask5711(dirs5711,normals5711);
819
820 trace.beginBlock ( "Testing when Q or P belongs to the segment...");
821 //Setting:
822 Point P(0,0);
823 Point Pb(0,-1);
824 Point Pup(0,10);
825 Point Q(10,5);
826 Point Lmin(10,-10);
827 Point Lmax(10,10);
828 Point midPointP,nextMidPointP;
829 Point midPointQ,nextMidPointQ;
830
831 Metric::ConstIterator itBeg = mask5711.begin();
832 Metric::ConstIterator itEnd = mask5711.end();
833
834 Metric::ConstIterator cone = mask5711.shrinkP(itBeg, itEnd, P, Q, Lmin, Lmax, 1, midPointP, nextMidPointP);
835 trace.info() <<" P - Shrink returns the cone "<< *cone<<" " <<*(cone+1)<<std::endl;
836 trace.info() <<" P - MidPoint "<< midPointP<<" " <<nextMidPointP<<std::endl<<std::endl;
837
838
839 Metric::ConstIterator cone2 = mask5711.shrinkP(itBeg, itEnd, Q, P, Lmin, Lmax, 1, midPointQ, nextMidPointQ);
840 Metric::ConstIterator cone3 = mask5711.shrinkP(itBeg, itEnd, Q, Pup, Lmin, Lmax, 1, midPointQ, nextMidPointQ);
841 trace.info() <<" Q - P Shrink returns the cone "<< *cone2<<" " <<*(cone2+1)<<std::endl;
842 trace.info() <<" Q - Pup Shrink returns the cone "<< *cone3<<" " <<*(cone3+1)<<std::endl;
843 trace.info() <<" Q - MidPoint "<< midPointQ<<" " <<nextMidPointQ<<std::endl<<std::endl;
844 nbok += ( cone2 == itBeg ) ? 1 : 0;
845 nb++;
846 trace.info() << "(" << nbok << "/" << nb << ") "
847 << " P's cone returns itBegin" << std::endl;
848
849 trace.endBlock();
850
851 Point PP(10,5);
852 Point QQ(3,4);
853 Point QQup(3,8);
854
855 Metric::Abscissa pos = mask5711.getLowerVoronoiEdgeAbscissa(PP, QQ, Lmin, Lmax, 1);
856 trace.warning() <<" Abscissa of voro call (PP,QQ) = "<< pos <<std::endl;
857 nbok += ( pos < 5 ) ? 1 : 0;
858 nb++;
859 trace.info() << "(" << nbok << "/" << nb << ") "
860 << " PP-QQ cell is lesser than 5" << std::endl;
861
862
863 Metric::Abscissa pos2 = mask5711.getLowerVoronoiEdgeAbscissa(PP, QQup, Lmin, Lmax, 1);
864 trace.warning() <<" Abscissa of voro call (PP,QQup) = "<< pos2 <<std::endl;
865 nbok += ( pos2 > 5 ) ? 1 : 0;
866 nb++;
867 trace.info() << "(" << nbok << "/" << nb << ") "
868 << " PP-QQ cell is greater than 5" << std::endl;
869
870
871 Metric::ConstIterator conebis = mask5711.shrinkP(itBeg, itEnd, P, Pb, Lmin, Lmax, 1, midPointP, nextMidPointP);
872 Metric::Abscissa posbis = mask5711.getLowerVoronoiEdgeAbscissa(P, Pb, Lmin, Lmax, 1);
873 trace.info() <<" Pbis - Shrink returns the cone "<< *conebis<<" " <<*(conebis+1)<<std::endl;
874 trace.info() <<" Pbis - Shrink returns the cone "<< *conebis<<" " <<*(conebis+1)<<std::endl;
875 trace.info() <<" Pbis - Pos "<< posbis<<std::endl;
876
877
878 return nbok == nb;
879}
880
881
882bool testSpecialCase()
883{
884 unsigned int nbok = 0;
885 unsigned int nb = 0;
886
887 //5-7-11 metic
888 typedef experimental::ChamferNorm2D<Z2i::Space> Metric;
889 Metric::Directions dirs5711;
890 Metric::Directions normals5711;
891 //5-7-11 mask
892 dirs5711.push_back(Z2i::Vector(0,-1));
893 dirs5711.push_back(Z2i::Vector(1,-2));
894 dirs5711.push_back(Z2i::Vector(1,-1));
895 dirs5711.push_back(Z2i::Vector(2,-1));
896 dirs5711.push_back(Z2i::Vector(1,0));
897 dirs5711.push_back(Z2i::Vector(2,1));
898 dirs5711.push_back(Z2i::Vector(1,1));
899 dirs5711.push_back(Z2i::Vector(1,2));
900
901 normals5711.push_back(Z2i::Vector(1,-5));
902 normals5711.push_back(Z2i::Vector(3,-4));
903 normals5711.push_back(Z2i::Vector(4,-3));
904 normals5711.push_back(Z2i::Vector(5,-1));
905 normals5711.push_back(Z2i::Vector(5,1));
906 normals5711.push_back(Z2i::Vector(4,3));
907 normals5711.push_back(Z2i::Vector(3,4));
908 normals5711.push_back(Z2i::Vector(1,5));
909
910 Metric mask5711(dirs5711,normals5711);
911
912 trace.beginBlock ( "Testing special cases...");
913 //Setting:
914 Point U(12,-12);
915 Point V(6,0);
916 Point W(-16,3);
917 Point Lmin(11,-20);
918 Point Lmax(11,20);
919 Point midPointP,nextMidPointP;
920 Point midPointQ,nextMidPointQ;
921
922 Metric::ConstIterator itBeg = mask5711.begin();
923 Metric::ConstIterator itEnd = mask5711.end();
924
925
926 Metric::ConstIterator cone2 = mask5711.shrinkP(itBeg, itEnd, V, W, Lmin, Lmax, 1, midPointQ, nextMidPointQ);
927 Metric::ConstIterator cone3 = mask5711.shrinkP(itBeg, itEnd, W, V, Lmin, Lmax, 1, midPointQ, nextMidPointQ);
928
929 trace.info() <<" V - W Shrink returns the cone "<< *cone2<<" " <<*(cone2+1)<<std::endl;
930 trace.info() <<" W - V Shrink returns the cone "<< *cone3<<" " <<*(cone3+1)<<std::endl;
931 trace.info() <<" Q - MidPoint "<< midPointQ<<" " <<nextMidPointQ<<std::endl<<std::endl;
932
933 Metric::Abscissa pos = mask5711.getLowerVoronoiEdgeAbscissa(V, W, Lmin, Lmax, 1);
934 trace.info() << " Pos = "<<pos<<std::endl;
935
936
937
938 nbok += ( mask5711.hiddenBy(U,V, W, Lmin, Lmax, 1) == false ) ? 1 : 0;
939 nb++;
940 trace.info() << "(" << nbok << "/" << nb << ") "
941 << " H(U,V,W) is false" << std::endl;
942
943
944
945 return nbok == nb;
946
947
948
949}
950
952// Standard services - public :
953
954int main( int argc, char** argv )
955{
956 trace.beginBlock ( "Testing class ChamferDT" );
957 trace.info() << "Args:";
958 for ( int i = 0; i < argc; ++i )
959 trace.info() << " " << argv[ i ];
960 trace.info() << std::endl;
961
962 bool res = testChamferSimple() && checkCMetricConcept() && testBasicMasks() && testIntersection()
963 && testShrink()
964 && testDoubleShrink()
965 && testDoubleShrinkHorizontal()
966 && testHiddenBy()
967 && testDoubleOnSegment()
968 && testSpecialCase(); // && ... other tests
969 trace.emphase() << ( res ? "Passed." : "Error." ) << std::endl;
970 trace.endBlock();
971 return res ? 0 : 1;
972}
973// //
DGtal::Trace::beginBlock
void beginBlock(const std::string &keyword="")
DGtal::Trace::warning
std::ostream & warning()
DGtal::Trace::emphase
std::ostream & emphase()
DGtal::Trace::info
std::ostream & info()
DGtal::Trace::endBlock
double endBlock()
DGtal::experimental::ChamferNorm2D
Aim: implements a model of CSeparableMetric for Chamfer and path based norms.
Definition ChamferNorm2D.h:87
DGtal::experimental::ChamferNorm2D::getNormalFromCone
Vector getNormalFromCone(ConstIterator aCone) const
DGtal::experimental::ChamferNorm2D::canonicalRay
Vector canonicalRay(const Vector &aRay) const
DGtal::experimental::ChamferNorm2D::getCone
ConstIterator getCone(const Vector &aDirection, ConstIterator aBegin, ConstIterator aEnd) const
DGtal::experimental::ChamferNorm2D::Directions
std::vector< Vector > Directions
Container for set of directions.
Definition ChamferNorm2D.h:103
DGtal::Z2i::Point
Space::Point Point
Definition StdDefs.h:95
DGtal
DGtal is the top-level namespace which contains all DGtal functions and types.
Definition ClosedIntegerHalfPlane.h:49
DGtal::trace
Trace trace
Definition Common.h:153
DGtal::concepts::CMetricSpace
Aim: defines the concept of metric spaces.
Definition CMetricSpace.h:125
DGtal::concepts::CSeparableMetric
Aim: defines the concept of separable metrics.
Definition CSeparableMetric.h:102
main
int main()
Definition testBits.cpp:56
testDoubleOnSegment
bool testDoubleOnSegment()
Definition testChamferDT.cpp:790
testDoubleShrinkHorizontal
bool testDoubleShrinkHorizontal()
Definition testChamferDT.cpp:572
testDoubleShrink
bool testDoubleShrink()
Definition testChamferDT.cpp:419
testIntersection
bool testIntersection()
Definition testChamferDT.cpp:214
testShrink
bool testShrink()
Definition testChamferDT.cpp:289
testBasicMasks
bool testBasicMasks()
Definition testChamferDT.cpp:179
testHiddenBy
bool testHiddenBy()
Definition testChamferDT.cpp:701
testChamferSimple
bool testChamferSimple()
Definition testChamferDT.cpp:59
checkCMetricConcept
bool checkCMetricConcept()
Definition testChamferDT.cpp:51
testSpecialCase
bool testSpecialCase()
Definition testChamferDT.cpp:882