Polygon_repair/repair_non_zero_polygon_2.cpp

#include <iostream>
#include <fstream>
#include <string>
 
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Polygon_repair/repair.h>
#include <CGAL/IO/WKT.h>
#ifdef CGAL_USE_BASIC_VIEWER
#include <CGAL/draw_multipolygon_with_holes_2.h>
#endif
 
using Kernel = CGAL::Exact_predicates_inexact_constructions_kernel;
using Point_2 = Kernel::Point_2;
using Polygon_2 = CGAL::Polygon_2<Kernel>;
using Polygon_with_holes_2 = CGAL::Polygon_with_holes_2<Kernel>;
using Multipolygon_with_holes_2 = CGAL::Multipolygon_with_holes_2<Kernel>;
 
int main(int argc, char* argv[]) {
 
  const char* filename = (argc > 1) ? argv[1] : "data/winding.wkt";
 
  std::ifstream in(filename);
  Polygon_with_holes_2 pin;
  CGAL::IO::read_polygon_WKT(in, pin);
 
  Multipolygon_with_holes_2 mp = CGAL::Polygon_repair::repair(pin, CGAL::Polygon_repair::Non_zero_rule());
#ifdef CGAL_USE_BASIC_VIEWER
  CGAL::draw(mp);
#endif
 
  if (mp.number_of_polygons_with_holes() > 1) {
    CGAL::IO::write_multi_polygon_WKT(std::cout, mp);
  } else {
    CGAL::IO::write_polygon_WKT(std::cout, mp.polygons_with_holes()[0]);
  }
  return 0;
}