CGAL 6.1 - 2D Polygon Repair

Loading...
Searching...
No Matches
Polygon_repair/repair_union_intersect_2.cpp
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Polygon_with_holes_2.h>
#include <CGAL/Multipolygon_with_holes_2.h>
#include <CGAL/draw_multipolygon_with_holes_2.h>
#include <CGAL/Polygon_repair/repair.h>
#include <CGAL/IO/WKT.h>
#include <iostream>
#include <sstream>
using K = CGAL::Exact_predicates_inexact_constructions_kernel;
using Point_2 = K::Point_2;
using Polygon_2 = CGAL::Polygon_2<K>;
using Polygon_with_holes_2 = CGAL::Polygon_with_holes_2<K>;
using Multipolygon_with_holes_2 = CGAL::Multipolygon_with_holes_2<K>;
int
main(int argc, char* argv[])
{
Multipolygon_with_holes_2 pA;
if (argc == 2) {
{
std::ifstream in(argv[1]);
CGAL::IO::read_multi_polygon_WKT(in, pA);
}
} else {
std::istringstream is("MULTIPOLYGON( ((0 0, 20 0, 20 20, 0 20), (1 1, 1 19, 19 19, 19 1) ) , (( 10 -2, 12 -2, 12 22, 10 22)) )");
//std::istringstream is("MULTIPOLYGON( ((0 0, 2 0, 2 3, 0 3) ) )"); // (0.1 0.1, 0.1 0.4, 0.4 0.1)
CGAL::IO::read_multi_polygon_WKT(is, pA);
}
Multipolygon_with_holes_2 mpwh = CGAL::Polygon_repair::repair(pA, CGAL::Polygon_repair::Union_rule());
{
std::ofstream out("union.wkt");
CGAL::IO::write_multi_polygon_WKT(out, mpwh);
#ifdef CGAL_USE_BASIC_VIEWER
CGAL::draw(mpwh);
#endif
}
mpwh = CGAL::Polygon_repair::repair(pA, CGAL::Polygon_repair::Intersection_rule());
{
std::ofstream out("intersection.wkt");
CGAL::IO::write_multi_polygon_WKT(out, mpwh);
CGAL::draw(mpwh);
}
{
Polygon_2 pB;
pB.push_back(Point_2(-1,-1));
pB.push_back(Point_2(1,-1));
pB.push_back(Point_2(1,1));
pB.push_back(Point_2(-1,1));
mpwh = CGAL::Polygon_repair::join(mpwh, pB);
std::ofstream out("joinn.wkt");
CGAL::IO::write_multi_polygon_WKT(out, mpwh);
CGAL::draw(mpwh);
}
return 0;
}
CGAL::Exact_predicates_inexact_constructions_kernel
CGAL::Multipolygon_with_holes_2
CGAL::Polygon_2
CGAL::Polygon_with_holes_2
CGAL::draw
void draw(const P &p, const GSOptions &gso)
CGAL::Polygon_repair::repair
Multipolygon_with_holes_2< Kernel, Container > repair(const Polygon_2< Kernel, Container > &p, Rule=Rule())
repairs polygon p using the given rule
Definition: repair.h:51
CGAL::Polygon_repair::join
Multipolygon_with_holes_2< Kernel, Container > join(const Multipolygon_with_holes_2< Kernel, Container > &pa)
computes the union of all polygons with holes in p
Definition: repair.h:821
CGAL::IO::write_multi_polygon_WKT
std::ostream & write_multi_polygon_WKT(std::ostream &out, MultiPolygon &polygons)
CGAL::IO::read_multi_polygon_WKT
bool read_multi_polygon_WKT(std::istream &in, MultiPolygon &polygons)
CGAL::Polygon_repair::Intersection_rule
Tag class to select the intersection rule when calling CGAL::Polygon_repair::repair().
Definition: Intersection_rule.h:30
CGAL::Polygon_repair::Union_rule
Tag class to select the union rule when calling CGAL::Polygon_repair::repair().
Definition: Union_rule.h:29