ceph/src/boost/libs/geometry/test/algorithms/correct_closure.cpp

   1 // Boost.Geometry
   2 // Unit Test
   3
   4 // Copyright (c) 2017 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 #include <sstream>
  11
  12 #include <geometry_test_common.hpp>
  13
  14 #include <boost/geometry/strategies/strategies.hpp>
  15
  16 #include <boost/geometry/geometries/point_xy.hpp>
  17 #include <boost/geometry/geometries/ring.hpp>
  18 #include <boost/geometry/geometries/polygon.hpp>
  19
  20 #include <boost/geometry/algorithms/correct_closure.hpp>
  21 #include <boost/geometry/io/wkt/read.hpp>
  22 #include <boost/geometry/io/wkt/write.hpp>
  23
  24 #include <boost/variant/variant.hpp>
  25
  26
  27 template <typename BaseGeometry, typename Geometry>
  28 void check_geometry(Geometry const& geometry, std::string const& expected)
  29 {
  30     std::ostringstream out;
  31     out << bg::wkt_manipulator<Geometry>(geometry, false);
  32     BOOST_CHECK_EQUAL(out.str(), expected);
  33 }
  34
  35 template <typename Geometry>
  36 void test_geometry(std::string const& wkt, std::string const& expected)
  37 {
  38     Geometry geometry;
  39     bg::read_wkt(wkt, geometry);
  40
  41     // Test tye type
  42     bg::correct_closure(geometry);
  43     check_geometry<Geometry>(geometry, expected);
  44
  45     // Test varianted type
  46     boost::variant<Geometry> v(geometry);
  47     bg::correct_closure(v);
  48     check_geometry<Geometry>(v, expected);
  49 }
  50
  51 template <typename P>
  52 void test_all()
  53 {
  54     typedef bg::model::ring<P, true, true> cw_closed_ring_type;
  55     typedef bg::model::ring<P, true, false> cw_open_ring_type;
  56     typedef bg::model::ring<P, false, true> ccw_closed_ring_type;
  57     typedef bg::model::ring<P, false, false> ccw_open_ring_type;
  58
  59     // Define clockwise and counter clockwise polygon
  60     std::string cw_ring       = "POLYGON((0 0,0 1,1 1,1 0,0 0))";
  61     std::string cw_open_ring  = "POLYGON((0 0,0 1,1 1,1 0))";
  62
  63     std::string ccw_ring      = "POLYGON((0 0,1 0,1 1,0 1,0 0))";
  64     std::string ccw_open_ring = "POLYGON((0 0,1 0,1 1,0 1))";
  65
  66     // Cases which should be closed or opened
  67     test_geometry<cw_closed_ring_type>(cw_open_ring, cw_ring);
  68     test_geometry<cw_open_ring_type>(cw_ring, cw_open_ring);
  69     test_geometry<ccw_closed_ring_type>(ccw_open_ring, ccw_ring);
  70     test_geometry<ccw_open_ring_type>(ccw_ring, ccw_open_ring);
  71
  72     // Cases which are incorrect but should still be closed or opened
  73     test_geometry<cw_closed_ring_type>(ccw_open_ring, ccw_ring);
  74     test_geometry<ccw_open_ring_type>(cw_ring, cw_open_ring);
  75
  76     // Cases where no action is necessary (even if order is incorrect)
  77     test_geometry<cw_closed_ring_type>(cw_ring, cw_ring);
  78     test_geometry<cw_closed_ring_type>(ccw_ring, ccw_ring);
  79     test_geometry<cw_open_ring_type>(cw_open_ring, cw_open_ring);
  80     test_geometry<cw_open_ring_type>(ccw_open_ring, ccw_open_ring);
  81     test_geometry<ccw_closed_ring_type>(cw_ring, cw_ring);
  82     test_geometry<ccw_closed_ring_type>(ccw_ring, ccw_ring);
  83     test_geometry<ccw_open_ring_type>(cw_open_ring, cw_open_ring);
  84     test_geometry<ccw_open_ring_type>(ccw_open_ring, ccw_open_ring);
  85
  86     // Polygon cases
  87     std::string cw_polygon =
  88             "POLYGON((0 0,0 4,4 4,4 0,0 0),(1 1,2 1,2 2,1 2,1 1))";
  89
  90     std::string cw_open_polygon =
  91             "POLYGON((0 0,0 4,4 4,4 0),(1 1,2 1,2 2,1 2))";
  92
  93     typedef bg::model::polygon<P, true, true> cw_closed_polygon_type;
  94     typedef bg::model::polygon<P, true, false> cw_open_polygon_type;
  95
  96     test_geometry<cw_closed_polygon_type>(cw_open_polygon, cw_polygon);
  97     test_geometry<cw_open_polygon_type>(cw_polygon, cw_open_polygon);
  98
  99     test_geometry<cw_closed_polygon_type>(cw_polygon, cw_polygon);
 100     test_geometry<cw_open_polygon_type>(cw_open_polygon, cw_open_polygon);
 101 }
 102
 103
 104 int test_main(int, char* [])
 105 {
 106     test_all<bg::model::d2::point_xy<int> >();
 107     test_all<bg::model::d2::point_xy<float> >();
 108     test_all<bg::model::d2::point_xy<double> >();
 109
 110     test_all<bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > >();
 111     test_all<bg::model::point<double, 2, bg::cs::geographic<bg::degree> > >();
 112
 113     return 0;
 114 }