DGtal 1.4.0
Loading...
Searching...
No Matches
testSegmentation.cpp
Go to the documentation of this file.
1
16//LICENSE-END
34#include <cstdio>
35#include <cmath>
36#include <fstream>
37#include <vector>
38#include <iostream>
39#include <iterator>
40
41
42
43
44#include "DGtal/base/Common.h"
45#include "DGtal/base/Exceptions.h"
46
47#include "DGtal/io/boards/Board2D.h"
48#include "DGtal/io/Color.h"
49
50#include "DGtal/base/Circulator.h"
51
52#include "DGtal/geometry/curves/ArithmeticalDSSComputer.h"
53#include "DGtal/geometry/curves/FreemanChain.h"
54#include "DGtal/geometry/curves/GreedySegmentation.h"
55#include "DGtal/geometry/curves/SaturatedSegmentation.h"
56
57
58#include "ConfigTest.h"
59
60
61using namespace DGtal;
62using namespace std;
63using namespace LibBoard;
64
66// Functions for testing class GreedySegmentation
68
69
70
74template <typename Iterator, typename Board>
75void draw(const Iterator& itb, const Iterator& ite, Board& aBoard)
76{
77
78 for (Iterator i(itb); i != ite; ++i) {
79
80 typename Iterator::SegmentComputer segment(*i);
81
82 aBoard << SetMode(segment.primitive().className(), "BoundingBox" )
83 << segment.primitive(); // draw bounding box
84
85 }
86
87}
88
92template <typename Iterator, typename Board>
93void drawCCP(const Iterator& itb, const Iterator& ite, Board& aBoard)
94{
95
96 typedef typename Iterator::SegmentComputer::ConstIterator PointIterator;
97
98 for (Iterator i(itb); i != ite; ++i) {
99
100 //choose pen color
101 CustomPenColor* aPenColor;
102
103 if ( !(i.intersectNext() && i.intersectPrevious()) ) {
104
105 aPenColor = new CustomPenColor( Color::Black );
106
107 } else {
108
109 //end points
110
111 PointIterator begin = i->begin(); --begin;
112 PointIterator end = i->end();
113
114 //parameters
115 int mu = i->mu();
116 int omega = i->omega();
117
118 //configurations
119 if ( (i->remainder(begin)<=mu-1)&&
120 (i->remainder(end)<=mu-1) ) { //concave
121 aPenColor = new CustomPenColor( Color::Green);
122 } else if ( (i->remainder(begin)>=mu+omega)&&
123 (i->remainder(end)>=mu+omega) ) { //convex
124 aPenColor = new CustomPenColor( Color::Blue );
125 } else if ( (i->remainder(begin)>=mu+omega)&&
126 (i->remainder(end)<=mu-1) ) { //convex to concave
127 aPenColor = new CustomPenColor( Color::Yellow );
128 } else if ( (i->remainder(begin)<=mu-1)&&
129 (i->remainder(end)>=mu+omega) ) { //concave to convex
130 aPenColor = new CustomPenColor( Color::Yellow );
131 } else { //pb
132 aPenColor = new CustomPenColor( Color::Red );
133 }
134
135 }
136
137 // draw each segment
138 aBoard << SetMode(i->primitive().className(), "BoundingBox" )
139 << CustomStyle( i->primitive().className(), aPenColor )
140 << i->primitive();
141
142 }
143
144}
145
149template <typename Iterator, typename Board>
150void segmentationIntoDSSs(const Iterator& itb, const Iterator& ite,
151 const Iterator& sitb, const Iterator& site,
152 const string& aMode, Board& aBoard)
153{
154 typedef typename IteratorCirculatorTraits<Iterator>::Value::Coordinate Coordinate;
155 typedef ArithmeticalDSSComputer<Iterator,Coordinate,4> RecognitionAlgorithm;
156 typedef GreedySegmentation<RecognitionAlgorithm> Segmentation;
157
158 RecognitionAlgorithm algo;
159 Segmentation s(itb,ite,algo);
160 s.setSubRange(sitb,site);
161 s.setMode(aMode);
162
163 typename Segmentation::SegmentComputerIterator i = s.begin();
164 typename Segmentation::SegmentComputerIterator end = s.end();
165
166 draw<typename Segmentation::SegmentComputerIterator, Board>
167 (i,end,aBoard);
168
169}
170
174template <typename Iterator, typename Board>
175void segmentationIntoMaximalDSSs(const Iterator& itb, const Iterator& ite,
176 const Iterator& sitb, const Iterator& site,
177 const string& aMode, Board& aBoard)
178{
179 typedef typename IteratorCirculatorTraits<Iterator>::Value::Coordinate Coordinate;
180 typedef ArithmeticalDSSComputer<Iterator,Coordinate,4> RecognitionAlgorithm;
181 typedef SaturatedSegmentation<RecognitionAlgorithm> Segmentation;
182
183 RecognitionAlgorithm algo;
184 Segmentation s(itb,ite,algo);
185 s.setSubRange(sitb,site);
186 s.setMode(aMode);
187
188 typename Segmentation::SegmentComputerIterator i = s.begin();
189 typename Segmentation::SegmentComputerIterator end = s.end();
190
191 drawCCP<typename Segmentation::SegmentComputerIterator, Board>
192 (i,end,aBoard);
193
194}
195
196
200bool greedySegmentationVisualTest()
201{
202
203 typedef int Coordinate;
204 typedef FreemanChain<Coordinate> FC;
205 typedef FreemanChain<Coordinate>::ConstIterator ConstIterator;
206
207
208
209 std::string filename = testPath + "samples/manche.fc";
210
211 std::fstream fst;
212 fst.open (filename.c_str(), std::ios::in);
213 FC fc(fst);
214
216// open digital curve
217
218 trace.beginBlock("Segmentation of a whole range (mode1)");
219{
220 Board2D aBoard;
221 aBoard << SetMode("PointVector", "Grid") << fc;
222
223 segmentationIntoDSSs<ConstIterator,Board2D>
224 (fc.begin(),fc.end(),
225 fc.begin(),fc.end(),
226 "Truncate",aBoard);
227
228 aBoard.saveEPS("WholeOpenCurveWithItMode1.eps");
229}
230 trace.endBlock();
231
232
233 trace.beginBlock("Segmentation of a whole range (mode3)");
234{
235 Board2D aBoard;
236 aBoard << SetMode("PointVector", "Grid") << fc;
237
238 segmentationIntoDSSs<ConstIterator,Board2D>
239 (fc.begin(),fc.end(),
240 fc.begin(),fc.end(),
241 "DoNotTruncate",aBoard);
242
243 aBoard.saveEPS("WholeOpenCurveWithItMode3.eps");
244}
245 trace.endBlock();
246
247 trace.beginBlock("Segmentation of a whole range (mode2)");
248{
249 Board2D aBoard;
250 aBoard << SetMode("PointVector", "Grid") << fc;
251
252 segmentationIntoDSSs<ConstIterator,Board2D>
253 (fc.begin(),fc.end(),
254 fc.begin(),fc.end(),
255 "Truncate+1",aBoard);
256
257 aBoard.saveEPS("WholeOpenCurveWithItMode2.eps");
258}
259 trace.endBlock();
260
261
263// subrange
264
265 typedef vector<PointVector<2,Coordinate> > Curve;
266 typedef Curve::const_iterator RAConstIterator;
267
268 Curve vPts;
269 vPts.assign ( fc.begin(), fc.end() );
270
271 RAConstIterator start = vPts.begin()+15;
272 RAConstIterator stop = vPts.begin()+200;
273
274trace.info() << *start << " " << *stop << endl;
275
276 trace.beginBlock("Segmentation of a subrange (mode1)");
277{
278 Board2D aBoard;
279 aBoard << SetMode("PointVector", "Grid") << fc;
280 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
281
282 segmentationIntoDSSs<RAConstIterator,Board2D>
283 (vPts.begin(),vPts.end(),
284 start,stop,
285 "Truncate",aBoard);
286
287 aBoard.saveEPS("PartOpenCurveWithItMode1.eps");
288}
289 trace.endBlock();
290
291 trace.beginBlock("Segmentation of a subrange (mode2)");
292{
293 Board2D aBoard;
294 aBoard << SetMode("PointVector", "Grid") << fc;
295 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
296
297 segmentationIntoDSSs<RAConstIterator,Board2D>
298 (vPts.begin(),vPts.end(),
299 start,stop,
300 "Truncate+1",aBoard);
301
302 aBoard.saveEPS("PartOpenCurveWithItMode2.eps");
303}
304 trace.endBlock();
305
306
307 trace.beginBlock("Segmentation of a subrange (mode3)");
308{
309 Board2D aBoard;
310 aBoard << SetMode("PointVector", "Grid") << fc;
311 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
312
313 segmentationIntoDSSs<RAConstIterator,Board2D>
314 (vPts.begin(),vPts.end(),
315 start,stop,
316 "DoNotTruncate",aBoard);
317
318 aBoard.saveEPS("PartOpenCurveWithItMode3.eps");
319}
320 trace.endBlock();
321
322
324// using circulators
325
326 typedef Circulator<RAConstIterator> ConstCirculator;
327 ConstCirculator c(vPts.begin(),vPts.begin(),vPts.end());
328
329 trace.beginBlock("Segmentation of a whole range (mode1) with circulators");
330{
331 Board2D aBoard;
332 aBoard << SetMode("PointVector", "Grid") << fc;
333
334 segmentationIntoDSSs<ConstCirculator,Board2D>
335 (c,c,c,c,
336 "Truncate",aBoard);
337
338 aBoard.saveEPS("WholeOpenCurveWithCircMode1.eps");
339}
340 trace.endBlock();
341
342 trace.beginBlock("Segmentation of a whole range (mode2) with circulators");
343{
344 Board2D aBoard;
345 aBoard << SetMode("PointVector", "Grid") << fc;
346
347 segmentationIntoDSSs<ConstCirculator,Board2D>
348 (c,c,c,c,
349 "Truncate+1",aBoard);
350
351 aBoard.saveEPS("WholeOpenCurveWithCircMode2.eps");
352}
353 trace.endBlock();
354
355
356 trace.beginBlock("Segmentation of a whole range (mode3) with circulators");
357{
358 Board2D aBoard;
359 aBoard << SetMode("PointVector", "Grid") << fc;
360
361 segmentationIntoDSSs<ConstCirculator,Board2D>
362 (c,c,c,c,
363 "DoNotTruncate",aBoard);
364
365 aBoard.saveEPS("WholeOpenCurveWithCircMode3.eps");
366}
367 trace.endBlock();
368
370// subrange with circulators
371
372 ConstCirculator cstart(start,vPts.begin(),vPts.end());
373 ConstCirculator cstop(stop,vPts.begin(),vPts.end());
374
375 trace.beginBlock("Segmentation of a subrange (mode1) with circulators");
376{
377 Board2D aBoard;
378 aBoard << SetMode("PointVector", "Grid") << fc;
379 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
380
381 segmentationIntoDSSs<ConstCirculator,Board2D>
382 (c,c,
383 cstart,cstop,
384 "Truncate",aBoard);
385
386 aBoard.saveEPS("PartOpenCurveWithCircMode1.eps");
387}
388 trace.endBlock();
389
390 trace.beginBlock("Segmentation of a subrange (mode2) with circulators");
391{
392 Board2D aBoard;
393 aBoard << SetMode("PointVector", "Grid") << fc;
394 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
395
396 segmentationIntoDSSs<ConstCirculator,Board2D>
397 (c,c,
398 cstart,cstop,
399 "Truncate+1",aBoard);
400
401 aBoard.saveEPS("PartOpenCurveWithCircMode2.eps");
402}
403 trace.endBlock();
404
405
406 trace.beginBlock("Segmentation of a subrange (mode3) with circulators");
407{
408 Board2D aBoard;
409 aBoard << SetMode("PointVector", "Grid") << fc;
410 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
411
412 segmentationIntoDSSs<ConstCirculator,Board2D>
413 (c,c,
414 cstart,cstop,
415 "DoNotTruncate",aBoard);
416
417 aBoard.saveEPS("PartOpenCurveWithCircMode3.eps");
418}
419 trace.endBlock();
420
422// closed digital curve
423
424 std::stringstream ss(stringstream::in | stringstream::out);
425 ss << "31 16 1112121212121221212121221212212222232232323332333333332333332330333033003030000010001001001000100010101010111" << endl;
426
427 // Construct the Freeman chain
428 FC fc2(ss);
429
430 Curve vPts2;
431 vPts2.assign( fc2.begin(), fc2.end() );
432 vPts2.insert( vPts2.end(), fc2.begin(), fc2.end() );
433 vPts2.insert( vPts2.end(), fc2.begin(), fc2.end() );
434
435//copy ( vPts2.begin(), vPts2.end(), ostream_iterator<PointVector<2,Coordinate> >(cout,"\n") ) ;
436
437 RAConstIterator start2 = vPts2.begin() + fc2.size()+1;
438 RAConstIterator stop2 = start2 + fc2.size()+1;
439
440trace.info() << *start2 << " " << *stop2 << endl;
441
442 trace.beginBlock("Segmentation of a subrange of a duplicated range (mode1)");
443{
444 Board2D aBoard;
445 aBoard << SetMode("PointVector", "Grid") << fc2;
446 aBoard << SetMode("PointVector", "Paving") << *start2 << *stop2;
447
448 segmentationIntoDSSs<RAConstIterator,Board2D>
449 (vPts2.begin(),vPts2.end(),
450 start2,stop2,
451 "Truncate",aBoard);
452
453 aBoard.saveEPS("DuplicatedCurveWithItMode1.eps");
454}
455 trace.endBlock();
456
457 trace.beginBlock("Segmentation of a subrange of a duplicated range (mode2)");
458{
459 Board2D aBoard;
460 aBoard << SetMode("PointVector", "Grid") << fc2;
461 aBoard << SetMode("PointVector", "Paving") << *start2 << *stop2;
462
463 segmentationIntoDSSs<RAConstIterator,Board2D>
464 (vPts2.begin(),vPts2.end(),
465 start2,stop2,
466 "Truncate+1",aBoard);
467
468 aBoard.saveEPS("DuplicatedCurveWithItMode2.eps");
469}
470 trace.endBlock();
471
472
473 trace.beginBlock("Segmentation of a subrange of a duplicated range (mode3)");
474{
475 Board2D aBoard;
476 aBoard << SetMode("PointVector", "Grid") << fc2;
477 aBoard << SetMode("PointVector", "Paving") << *start2 << *stop2;
478
479 segmentationIntoDSSs<RAConstIterator,Board2D>
480 (vPts2.begin(),vPts2.end(),
481 start2,stop2,
482 "DoNotTruncate",aBoard);
483
484 aBoard.saveEPS("DuplicatedCurveWithItMode3.eps");
485}
486 trace.endBlock();
487
488
490// closed digital curve with circulators
491
492 ConstCirculator c2(start2,start2,stop2);
493
494 trace.beginBlock("Segmentation of a whole range (mode1) with circulators");
495{
496 Board2D aBoard;
497 aBoard << SetMode("PointVector", "Grid") << fc2;
498
499 segmentationIntoDSSs<ConstCirculator,Board2D>
500 (c2,c2,c2,c2,
501 "Truncate",aBoard);
502
503 aBoard.saveEPS("ClosedCurveWithCircMode1.eps");
504}
505 trace.endBlock();
506
507 trace.beginBlock("Segmentation of a whole range (mode2) with circulators");
508{
509 Board2D aBoard;
510 aBoard << SetMode("PointVector", "Grid") << fc2;
511
512 segmentationIntoDSSs<ConstCirculator,Board2D>
513 (c2,c2,c2,c2,
514 "Truncate+1",aBoard);
515
516 aBoard.saveEPS("ClosedCurveWithCircMode2.eps");
517}
518 trace.endBlock();
519
520
521 trace.beginBlock("Segmentation of a whole range (mode3) with circulators");
522{
523 Board2D aBoard;
524 aBoard << SetMode("PointVector", "Grid") << fc2;
525
526 segmentationIntoDSSs<ConstCirculator,Board2D>
527 (c2,c2,c2,c2,
528 "DoNotTruncate",aBoard);
529
530 aBoard.saveEPS("ClosedCurveWithCircMode3.eps");
531}
532 trace.endBlock();
533
534
536 return true;
537}
538
539
541// Functions for testing class SaturatedSegmentation
543
544
545
549bool SaturatedSegmentationVisualTest()
550{
551
552 typedef int Coordinate;
553 typedef FreemanChain<Coordinate> FC;
554 typedef FreemanChain<Coordinate>::ConstIterator ConstIterator;
555
556 std::string filename = testPath + "samples/manche.fc";
557
558 std::fstream fst;
559 fst.open (filename.c_str(), std::ios::in);
560 FC fc(fst);
561
563// whole open curve
564
565 trace.beginBlock("Simple saturated Segmentation (mode=First)");
566{
567 Board2D aBoard;
568 aBoard << SetMode("PointVector", "Grid") << fc;
569
570 segmentationIntoMaximalDSSs<ConstIterator,Board2D>
571 (fc.begin(),fc.end(),
572 fc.begin(),fc.end(),
573 "First",aBoard);
574
575 aBoard.saveEPS("MSOpenCurve2.eps");
576}
577 trace.endBlock();
578
579 trace.beginBlock("Simple saturated Segmentation (mode=MostCentered)");
580{
581 Board2D aBoard;
582 aBoard << SetMode("PointVector", "Grid") << fc;
583
584 segmentationIntoMaximalDSSs<ConstIterator,Board2D>
585 (fc.begin(),fc.end(),
586 fc.begin(),fc.end(),
587 "MostCentered",aBoard);
588
589 aBoard.saveEPS("MSOpenCurve3.eps");
590}
591 trace.endBlock();
592
593 trace.beginBlock("Simple saturated Segmentation (mode=Last)");
594{
595 Board2D aBoard;
596 aBoard << SetMode("PointVector", "Grid") << fc;
597
598 segmentationIntoMaximalDSSs<ConstIterator,Board2D>
599 (fc.begin(),fc.end(),
600 fc.begin(),fc.end(),
601 "Last",aBoard);
602
603 aBoard.saveEPS("MSOpenCurve4.eps");
604}
605 trace.endBlock();
606
608// subrange
609
610 typedef vector<PointVector<2,Coordinate> > Curve;
611 typedef Curve::const_iterator RAConstIterator;
612
613 Curve vPts;
614 vPts.assign ( fc.begin(), fc.end() );
615
616 RAConstIterator start = vPts.begin()+190;
617 RAConstIterator stop = vPts.begin()+400;
618
619trace.info() << *start << " " << *stop << endl;
620
621
622 trace.beginBlock("saturated Segmentation of a subrange (mode=First)");
623{
624 Board2D aBoard;
625 aBoard << SetMode("PointVector", "Grid") << fc;
626 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
627
628 segmentationIntoMaximalDSSs<RAConstIterator,Board2D>
629 (vPts.begin(),vPts.end(),
630 start,stop,
631 "First",aBoard);
632
633 aBoard.saveEPS("MSOpenCurvePart2.eps");
634}
635 trace.endBlock();
636
637 trace.beginBlock("saturated Segmentation of a subrange (mode=MostCentered)");
638{
639 Board2D aBoard;
640 aBoard << SetMode("PointVector", "Grid") << fc;
641 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
642
643 segmentationIntoMaximalDSSs<RAConstIterator,Board2D>
644 (vPts.begin(),vPts.end(),
645 start,stop,
646 "MostCentered",aBoard);
647
648 aBoard.saveEPS("MSOpenCurvePart3.eps");
649}
650 trace.endBlock();
651
652 trace.beginBlock("saturated Segmentation of a subrange (mode=Last)");
653{
654 Board2D aBoard;
655 aBoard << SetMode("PointVector", "Grid") << fc;
656 aBoard << SetMode("PointVector", "Paving") << *start << *stop;
657
658 segmentationIntoMaximalDSSs<RAConstIterator,Board2D>
659 (vPts.begin(),vPts.end(),
660 start,stop,
661 "Last",aBoard);
662
663 aBoard.saveEPS("MSOpenCurvePart4.eps");
664}
665 trace.endBlock();
666
667
669// using circulators
670
671 typedef Circulator<RAConstIterator> ConstCirculator;
672 ConstCirculator c(vPts.begin(),vPts.begin(),vPts.end());
673
674
675 trace.beginBlock("Segmentation of a range bounded by circulators (mode=First)");
676{
677 Board2D aBoard;
678 aBoard << SetMode("PointVector", "Grid") << fc;
679
680 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
681 (c,c,c,c,
682 "First",aBoard);
683
684 aBoard.saveEPS("MSOpenCurveWithCirc2.eps");
685}
686 trace.endBlock();
687
688 trace.beginBlock("Segmentation of a range bounded by circulators (mode=MostCentered)");
689{
690 Board2D aBoard;
691 aBoard << SetMode("PointVector", "Grid") << fc;
692
693 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
694 (c,c,c,c,
695 "MostCentered",aBoard);
696
697 aBoard.saveEPS("MSOpenCurveWithCirc3.eps");
698}
699 trace.endBlock();
700
701 trace.beginBlock("Segmentation of a range bounded by circulators (mode=Last)");
702{
703 Board2D aBoard;
704 aBoard << SetMode("PointVector", "Grid") << fc;
705
706 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
707 (c,c,c,c,
708 "Last",aBoard);
709
710 aBoard.saveEPS("MSOpenCurveWithCirc4.eps");
711}
712 trace.endBlock();
713
714
716// subrange with circulators
717
718 ConstCirculator cstart(start,vPts.begin(),vPts.end());
719 ConstCirculator cstop(stop,vPts.begin(),vPts.end());
720
721
722 trace.beginBlock("Segmentation of a subrange with circulators (mode=First)");
723{
724 Board2D aBoard;
725 aBoard << SetMode("PointVector", "Grid") << fc;
726 aBoard << SetMode("PointVector", "Paving") << *cstart << *cstop;
727
728 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
729 (c,c,cstart,cstop,
730 "First",aBoard);
731
732 aBoard.saveEPS("MSOpenCurvePartWithCirc2.eps");
733}
734 trace.endBlock();
735
736 trace.beginBlock("Segmentation of a subrange with circulators (mode=MostCentered)");
737{
738 Board2D aBoard;
739 aBoard << SetMode("PointVector", "Grid") << fc;
740 aBoard << SetMode("PointVector", "Paving") << *cstart << *cstop;
741
742 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
743 (c,c,cstart,cstop,
744 "MostCentered",aBoard);
745
746 aBoard.saveEPS("MSOpenCurvePartWithCirc3.eps");
747}
748 trace.endBlock();
749
750 trace.beginBlock("Segmentation of a subrange with circulators (mode=Last)");
751{
752 Board2D aBoard;
753 aBoard << SetMode("PointVector", "Grid") << fc;
754 aBoard << SetMode("PointVector", "Paving") << *cstart << *cstop;
755
756 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
757 (c,c,cstart,cstop,
758 "Last",aBoard);
759
760 aBoard.saveEPS("MSOpenCurvePartWithCirc4.eps");
761}
762 trace.endBlock();
763
765// closed digital curve with circulators
766
767 std::stringstream ss(stringstream::in | stringstream::out);
768 ss << "31 16 1112121212121221212121221212212222232232323332333333332333332330333033003030000010001001001000100010101010111" << endl;
769
770 // Construct the Freeman chain
771 FC fc2(ss);
772
773 Curve vPts2;
774 vPts2.assign( fc2.begin(), fc2.end() );
775 vPts2.insert( vPts2.end(), fc2.begin(), fc2.end() );
776 vPts2.insert( vPts2.end(), fc2.begin(), fc2.end() );
777
778 RAConstIterator start2 = vPts2.begin() + fc2.size()+1;
779 RAConstIterator stop2 = start2 + fc2.size()+1;
780
781 trace.info() << *start2 << " " << *stop2 << endl;
782
783 trace.beginBlock("Satured Segmentation of a duplicated range (mode2)");
784{
785 Board2D aBoard;
786 aBoard << SetMode("PointVector", "Grid") << fc2;
787 aBoard << SetMode("PointVector", "Paving") << *start2 << *stop2;
788
789 segmentationIntoMaximalDSSs<RAConstIterator,Board2D>
790 (vPts2.begin(),vPts2.end(),
791 start2,stop2,
792 "First",aBoard);
793
794 aBoard.saveEPS("MSClosedCurveWithIt2.eps");
795}
796 trace.endBlock();
797
798 trace.beginBlock("Satured Segmentation of a duplicated range (mode3)");
799{
800 Board2D aBoard;
801 aBoard << SetMode("PointVector", "Grid") << fc2;
802 aBoard << SetMode("PointVector", "Paving") << *start2 << *stop2;
803
804 segmentationIntoMaximalDSSs<RAConstIterator,Board2D>
805 (vPts2.begin(),vPts2.end(),
806 start2,stop2,
807 "MostCentered",aBoard);
808
809 aBoard.saveEPS("MSClosedCurveWithIt3.eps");
810}
811 trace.endBlock();
812
813
814 trace.beginBlock("Satured Segmentation of a duplicated range (mode4)");
815{
816 Board2D aBoard;
817 aBoard << SetMode("PointVector", "Grid") << fc2;
818 aBoard << SetMode("PointVector", "Paving") << *start2 << *stop2;
819
820 segmentationIntoMaximalDSSs<RAConstIterator,Board2D>
821 (vPts2.begin(),vPts2.end(),
822 start2,stop2,
823 "Last",aBoard);
824
825 aBoard.saveEPS("MSClosedCurveWithIt4.eps");
826}
827 trace.endBlock();
828
829
832
833 Curve vPts3;
834 vPts3.assign( fc2.begin(), fc2.end() );
835
836 ConstCirculator c2(vPts3.begin(),vPts3.begin(),vPts3.end());
837
838 trace.beginBlock("saturated Segmentation of a closed digital curve (mode=First)");
839{
840 Board2D aBoard;
841 aBoard << SetMode("PointVector", "Grid") << fc2;
842 aBoard << SetMode("PointVector", "Paving") << *c2;
843
844 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
845 (c2,c2,c2,c2,
846 "First",aBoard);
847
848 aBoard.saveEPS("MSClosedCurveWithCirc2.eps");
849}
850 trace.endBlock();
851
852 trace.beginBlock("saturated Segmentation of a closed digital curve (mode=MostCentered)");
853{
854 Board2D aBoard;
855 aBoard << SetMode("PointVector", "Grid") << fc2;
856 aBoard << SetMode("PointVector", "Paving") << *c2;
857
858 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
859 (c2,c2,c2,c2,
860 "MostCentered",aBoard);
861
862 aBoard.saveEPS("MSClosedCurveWithCirc3.eps");
863}
864 trace.endBlock();
865
866 trace.beginBlock("saturated Segmentation of a closed digital curve (mode=Last)");
867{
868 Board2D aBoard;
869 aBoard << SetMode("PointVector", "Grid") << fc2;
870 aBoard << SetMode("PointVector", "Paving") << *c2;
871
872 segmentationIntoMaximalDSSs<ConstCirculator,Board2D>
873 (c2,c2,c2,c2,
874 "Last",aBoard);
875
876 aBoard.saveEPS("MSClosedCurveWithCirc4.eps");
877}
878 trace.endBlock();
879
881
882 return true;
883
884}
885
886
890bool SaturatedSegmentationTest()
891{
892
893 typedef int Coordinate;
894 typedef FreemanChain<Coordinate> FC;
895
896 std::string filename = testPath + "samples/BigBall2.fc";
897
898 std::fstream fst;
899 fst.open (filename.c_str(), std::ios::in);
900 FC fc(fst);
901
902 typedef PointVector<2,Coordinate> Point;
903
904 vector<Point> vPts;
905 vPts.assign(fc.begin(),fc.end());
906
907 typedef vector<Point>::const_iterator ConstIterator;
908 typedef Circulator<ConstIterator> ConstCirculator;
909
910 Circulator<ConstIterator> c(vPts.begin(), vPts.begin(), vPts.end() );
911
912 typedef ArithmeticalDSSComputer<ConstCirculator,Coordinate,4> RecognitionAlgorithm;
913 typedef SaturatedSegmentation<RecognitionAlgorithm> Segmentation;
914
915 trace.beginBlock("saturated Segmentation");
916 trace.info() << filename << endl;
917
918 RecognitionAlgorithm algo;
919 Segmentation s(c,c,algo);
920 s.setMode("First");
921
922 Segmentation::SegmentComputerIterator begin = s.begin();
923 Segmentation::SegmentComputerIterator end = s.end();
924
925 unsigned int compteur = 0;
926 for (Segmentation::SegmentComputerIterator i(begin); i != end; ++i) {
927 ++compteur;
928 }
929
930 trace.info() << "# nbpts nbsegments " << endl;
931 trace.info() << fc.size()+1 << " " << compteur << endl;
932
933 trace.endBlock();
934
935 return (compteur == 4295);
936}
937
941
942int main(int argc, char **argv)
943{
944
945 trace.beginBlock ( "Testing class GreedyDecomposition and SaturatedSegmentation" );
946 trace.info() << "Args:";
947 for ( int i = 0; i < argc; ++i )
948 trace.info() << " " << argv[ i ];
949 trace.info() << endl;
950
951 bool res = greedySegmentationVisualTest()
952&& SaturatedSegmentationVisualTest()
953&& SaturatedSegmentationTest()
954;
955
956 trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
957 trace.endBlock();
958
959 return res ? 0 : 1;
960
961}
DGtal::ArithmeticalDSSComputer
Aim: This class is a wrapper around ArithmeticalDSS that is devoted to the dynamic recognition of dig...
Definition ArithmeticalDSSComputer.h:87
DGtal::Board2D
Aim: This class specializes a 'Board' class so as to display DGtal objects more naturally (with <<)....
Definition Board2D.h:71
DGtal::Circulator
Aim: Provides an adapter for classical iterators that can iterate through the underlying data structu...
Definition Circulator.h:86
DGtal::Color::Yellow
static const Color Yellow
Definition Color.h:422
DGtal::Color::Green
static const Color Green
Definition Color.h:417
DGtal::Color::Red
static const Color Red
Definition Color.h:416
DGtal::Color::Blue
static const Color Blue
Definition Color.h:419
DGtal::Color::Black
static const Color Black
Definition Color.h:413
DGtal::GreedySegmentation
Aim: Computes the greedy segmentation of a range given by a pair of ConstIterators....
Definition GreedySegmentation.h:150
DGtal::PointVector
Aim: Implements basic operations that will be used in Point and Vector classes.
Definition PointVector.h:593
DGtal::SaturatedSegmentation::SegmentComputerIterator
Aim: Specific iterator to visit all the maximal segments of a saturated segmentation.
Definition SaturatedSegmentation.h:180
DGtal::SaturatedSegmentation
Aim: Computes the saturated segmentation, that is the whole set of maximal segments within a range gi...
Definition SaturatedSegmentation.h:154
DGtal::Trace::beginBlock
void beginBlock(const std::string &keyword="")
DGtal::Trace::emphase
std::ostream & emphase()
DGtal::Trace::info
std::ostream & info()
DGtal::Trace::endBlock
double endBlock()
LibBoard::Board
Class for EPS, FIG or SVG drawings.
Definition Board.h:35
LibBoard::Board::saveEPS
void saveEPS(const char *filename, PageSize size=Board::BoundingBox, double margin=10.0) const
Definition Board.cpp:804
ConstIterator
MyDigitalSurface::ConstIterator ConstIterator
Definition greedy-plane-segmentation-ex2.cpp:93
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
std
STL namespace.
DGtal::CustomPenColor
Custom style class redefining the pen color. You may use Board2D::Color::None for transparent color.
Definition Board2D.h:313
DGtal::IteratorCirculatorTraits
Aim: Provides nested types for both iterators and circulators: Type, Category, Value,...
Definition IteratorCirculatorTraits.h:294
DGtal::SetMode
Modifier class in a Board2D stream. Useful to choose your own mode for a given class....
Definition Board2D.h:247
Segmentation
SaturatedSegmentation< SegmentComputer > Segmentation
Definition testArithmeticalDSSComputerOnSurfels.cpp:56
main
int main()
Definition testBits.cpp:56
Point
MyPointD Point
Definition testClone2.cpp:383
greedySegmentationVisualTest
bool greedySegmentationVisualTest()
Definition testSegmentation.cpp:200
segmentationIntoMaximalDSSs
void segmentationIntoMaximalDSSs(const Iterator &itb, const Iterator &ite, const Iterator &sitb, const Iterator &site, const string &aMode, Board &aBoard)
Definition testSegmentation.cpp:175
SaturatedSegmentationVisualTest
bool SaturatedSegmentationVisualTest()
Definition testSegmentation.cpp:549
SaturatedSegmentationTest
bool SaturatedSegmentationTest()
Definition testSegmentation.cpp:890
drawCCP
void drawCCP(const Iterator &itb, const Iterator &ite, Board &aBoard)
Definition testSegmentation.cpp:93
draw
void draw(const Iterator &itb, const Iterator &ite, Board &aBoard)
Definition testSegmentation.cpp:75
segmentationIntoDSSs
void segmentationIntoDSSs(const Iterator &itb, const Iterator &ite, const Iterator &sitb, const Iterator &site, const string &aMode, Board &aBoard)
Definition testSegmentation.cpp:150