• Home
  • Line#
  • Scopes#
  • Navigate#
  • Raw
  • Download
1 // Boost.Geometry (aka GGL, Generic Geometry Library)
2 // QuickBook Example
3 
4 // Copyright (c) 2011-2012 Barend Gehrels, Amsterdam, the Netherlands.
5 
6 // Use, modification and distribution is subject to the Boost Software License,
7 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
8 // http://www.boost.org/LICENSE_1_0.txt)
9 
10 //[union
11 //` Shows how to get a united geometry of two polygons
12 
13 #include <iostream>
14 #include <vector>
15 
16 #include <boost/geometry.hpp>
17 #include <boost/geometry/geometries/point_xy.hpp>
18 #include <boost/geometry/geometries/polygon.hpp>
19 
20 #include <boost/foreach.hpp>
21 /*<-*/ #include "create_svg_overlay.hpp" /*->*/
22 
main()23 int main()
24 {
25     typedef boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > polygon;
26 
27     polygon green, blue;
28 
29     boost::geometry::read_wkt(
30         "POLYGON((2 1.3,2.4 1.7,2.8 1.8,3.4 1.2,3.7 1.6,3.4 2,4.1 3,5.3 2.6,5.4 1.2,4.9 0.8,2.9 0.7,2 1.3)"
31             "(4.0 2.0, 4.2 1.4, 4.8 1.9, 4.4 2.2, 4.0 2.0))", green);
32 
33     boost::geometry::read_wkt(
34         "POLYGON((4.0 -0.5 , 3.5 1.0 , 2.0 1.5 , 3.5 2.0 , 4.0 3.5 , 4.5 2.0 , 6.0 1.5 , 4.5 1.0 , 4.0 -0.5))", blue);
35 
36     std::vector<polygon> output;
37     boost::geometry::union_(green, blue, output);
38 
39     int i = 0;
40     std::cout << "green || blue:" << std::endl;
41     BOOST_FOREACH(polygon const& p, output)
42     {
43         std::cout << i++ << ": " << boost::geometry::area(p) << std::endl;
44     }
45 
46     /*<-*/ create_svg("union.svg", green, blue, output); /*->*/
47     return 0;
48 }
49 
50 //]
51 
52 
53 //[union_output
54 /*`
55 Output:
56 [pre
57 green || blue:
58 0: 5.64795
59 
60 [$img/algorithms/union.png]
61 
62 ]
63 */
64 //]
65 
66