/Src/Dependencies/Boost/libs/geometry/doc/src/examples/algorithms/difference_inserter.cpp

http://hadesmem.googlecode.com/ · C++ · 88 lines · 36 code · 18 blank · 34 comment · 0 complexity · 281b06e4933cc10fe05b780cfe6b82af MD5 · raw file

  1. // Boost.Geometry (aka GGL, Generic Geometry Library)
  2. // QuickBook Example
  3. // Copyright (c) 2011 Barend Gehrels, Amsterdam, the Netherlands.
  4. // Use, modification and distribution is subject to the Boost Software License,
  5. // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
  6. // http://www.boost.org/LICENSE_1_0.txt)
  7. //[difference_inserter
  8. //` Shows how the difference_inserter function can be used
  9. #include <iostream>
  10. #include <vector>
  11. #include <boost/geometry.hpp>
  12. #include <boost/geometry/geometries/point_xy.hpp>
  13. #include <boost/geometry/geometries/polygon.hpp>
  14. #include <boost/geometry/domains/gis/io/wkt/wkt.hpp>
  15. #include <boost/foreach.hpp>
  16. /*<-*/ #include "create_svg_overlay.hpp" /*->*/
  17. int main()
  18. {
  19. typedef boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > polygon;
  20. polygon green, blue;
  21. boost::geometry::read_wkt(
  22. "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)"
  23. "(4.0 2.0, 4.2 1.4, 4.8 1.9, 4.4 2.2, 4.0 2.0))", green);
  24. boost::geometry::read_wkt(
  25. "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);
  26. std::vector<polygon> output;
  27. // Note that this sample simulates the symmetric difference,
  28. // which is also available as a separate algorithm.
  29. // It chains the output iterator returned by the function to the second instance.
  30. boost::geometry::difference_inserter<polygon>
  31. (
  32. green, blue,
  33. boost::geometry::difference_inserter<polygon>
  34. (
  35. blue, green,
  36. std::back_inserter(output)
  37. )
  38. );
  39. int i = 0;
  40. std::cout << "(blue \ green) u (green \ blue):" << std::endl;
  41. BOOST_FOREACH(polygon const& p, output)
  42. {
  43. std::cout << i++ << ": " << boost::geometry::area(p) << std::endl;
  44. }
  45. /*<-*/ create_svg("difference_inserter.svg", green, blue, output); /*->*/
  46. return 0;
  47. }
  48. //]
  49. //[difference_inserter_output
  50. /*`
  51. Output:
  52. [pre
  53. (blue \\ green) u (green \\ blue):
  54. 0: 0.525154
  55. 1: 0.015
  56. 2: 0.181136
  57. 3: 0.128798
  58. 4: 0.340083
  59. 5: 0.307778
  60. 6: 0.02375
  61. 7: 0.542951
  62. 8: 0.0149697
  63. 9: 0.226855
  64. 10: 0.839424
  65. [$img/algorithms/sym_difference.png]
  66. ]
  67. */
  68. //]