123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Copyright (c) 2014-2015, Oracle and/or its affiliates.
- // Licensed under the Boost Software License version 1.0.
- // http://www.boost.org/users/license.html
- // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- #ifndef BOOST_TEST_MODULE
- #define BOOST_TEST_MODULE test_disjoint_coverage
- #endif
- // unit test to test disjoint for all geometry combinations
- #include <iostream>
- #include <boost/test/included/unit_test.hpp>
- #include <boost/geometry/core/tag.hpp>
- #include <boost/geometry/core/tags.hpp>
- #include <boost/geometry/strategies/strategies.hpp>
- #include <boost/geometry/io/wkt/wkt.hpp>
- #include <boost/geometry/io/dsv/write.hpp>
- #include <boost/geometry/geometries/geometries.hpp>
- #include <boost/geometry/algorithms/disjoint.hpp>
- #include <from_wkt.hpp>
- #ifdef HAVE_TTMATH
- #include <boost/geometry/extensions/contrib/ttmath_stub.hpp>
- #endif
- namespace bg = ::boost::geometry;
- //============================================================================
- struct test_disjoint
- {
- template <typename Geometry1, typename Geometry2>
- static inline void apply(std::string const& case_id,
- Geometry1 const& geometry1,
- Geometry2 const& geometry2,
- bool expected_result)
- {
- bool result = bg::disjoint(geometry1, geometry2);
- BOOST_CHECK_MESSAGE(result == expected_result,
- "case ID: " << case_id << ", G1: " << bg::wkt(geometry1)
- << ", G2: " << bg::wkt(geometry2) << " -> Expected: "
- << expected_result << ", detected: " << result);
- result = bg::disjoint(geometry2, geometry1);
- BOOST_CHECK_MESSAGE(result == expected_result,
- "case ID: " << case_id << ", G1: " << bg::wkt(geometry2)
- << ", G2: " << bg::wkt(geometry1) << " -> Expected: "
- << expected_result << ", detected: " << result);
- #ifdef BOOST_GEOMETRY_TEST_DEBUG
- std::cout << "case ID: " << case_id << "; G1 - G2: ";
- std::cout << bg::wkt(geometry1) << " - ";
- std::cout << bg::wkt(geometry2) << std::endl;
- std::cout << std::boolalpha;
- std::cout << "expected/computed result: "
- << expected_result << " / " << result << std::endl;
- std::cout << std::endl;
- std::cout << std::noboolalpha;
- #endif
- }
- };
- //============================================================================
- // areal-areal geometries
- template <typename P>
- inline void test_box_box()
- {
- typedef bg::model::box<P> B;
- typedef test_disjoint tester;
- tester::apply("b-b-01",
- from_wkt<B>("BOX(2 2,3 3)"),
- from_wkt<B>("BOX(0 0,2 2)"),
- false);
- tester::apply("b-b-02",
- from_wkt<B>("BOX(1 1,3 3)"),
- from_wkt<B>("BOX(0 0,2 2)"),
- false);
- tester::apply("b-b-03",
- from_wkt<B>("BOX(3 3,4 4)"),
- from_wkt<B>("BOX(0 0,2 2)"),
- true);
- }
- template <typename P>
- inline void test_ring_box()
- {
- typedef bg::model::box<P> B;
- typedef bg::model::ring<P, false, false> R; // ccw, open
- typedef test_disjoint tester;
- tester::apply("r-b-01",
- from_wkt<B>("BOX(2 2,3 3)"),
- from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("r-b-02",
- from_wkt<B>("BOX(1 1,3 3)"),
- from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("r-b-03",
- from_wkt<B>("BOX(3 3,4 4)"),
- from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"),
- true);
- }
- template <typename P>
- inline void test_polygon_box()
- {
- typedef bg::model::box<P> B;
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef test_disjoint tester;
- tester::apply("pg-b-01",
- from_wkt<B>("BOX(2 2,3 3)"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("pg-b-02",
- from_wkt<B>("BOX(1 1,3 3)"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("pg-b-03",
- from_wkt<B>("BOX(3 3,4 4)"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- true);
- }
- template <typename P>
- inline void test_multipolygon_box()
- {
- typedef bg::model::box<P> B;
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef bg::model::multi_polygon<PL> MPL;
-
- typedef test_disjoint tester;
- tester::apply("mpg-b-01",
- from_wkt<B>("BOX(2 2,3 3)"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("mpg-b-02",
- from_wkt<B>("BOX(1 1,3 3)"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("mpg-b-03",
- from_wkt<B>("BOX(3 3,4 4)"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- true);
- }
- template <typename P>
- inline void test_ring_ring()
- {
- typedef bg::model::ring<P, false, false> R; // ccw, open
- typedef test_disjoint tester;
- tester::apply("r-r-01",
- from_wkt<R>("POLYGON((2 2,2 3,3 3,3 2))"),
- from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("r-r-02",
- from_wkt<R>("POLYGON((1 1,1 3,3 3,3 1))"),
- from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("r-r-03",
- from_wkt<R>("POLYGON((3 3,3 4,4 4,4 3))"),
- from_wkt<R>("POLYGON((0 0,2 0,2 2,0 2))"),
- true);
- }
- template <typename P>
- inline void test_polygon_ring()
- {
- typedef bg::model::ring<P, false, false> R; // ccw, open
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef test_disjoint tester;
- tester::apply("pg-r-01",
- from_wkt<R>("POLYGON((2 2,2 3,3 3,3 2))"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("pg-r-02",
- from_wkt<R>("POLYGON((1 1,1 3,3 3,3 1))"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("pg-r-03",
- from_wkt<R>("POLYGON((3 3,3 4,4 4,4 3))"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- true);
- }
- template <typename P>
- inline void test_multipolygon_ring()
- {
- typedef bg::model::ring<P, false, false> R; // ccw, open
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef bg::model::multi_polygon<PL> MPL;
- typedef test_disjoint tester;
- tester::apply("mpg-r-01",
- from_wkt<R>("POLYGON((2 2,2 3,3 3,3 2))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("mpg-r-02",
- from_wkt<R>("POLYGON((1 1,1 3,3 3,3 1))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("mpg-r-03",
- from_wkt<R>("POLYGON((3 3,3 4,4 4,4 3))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- true);
- }
- template <typename P>
- inline void test_polygon_polygon()
- {
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef test_disjoint tester;
- tester::apply("pg-pg-01",
- from_wkt<PL>("POLYGON((2 2,2 3,3 3,3 2))"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("pg-pg-02",
- from_wkt<PL>("POLYGON((1 1,1 3,3 3,3 1))"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- false);
- tester::apply("pg-pg-03",
- from_wkt<PL>("POLYGON((3 3,3 4,4 4,4 3))"),
- from_wkt<PL>("POLYGON((0 0,2 0,2 2,0 2))"),
- true);
- tester::apply("pg-pg-04",
- from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9))"),
- from_wkt<PL>("POLYGON((3 3,6 3,6 6,3 6))"),
- false);
- // polygon with a hole which entirely contains the other polygon
- tester::apply("pg-pg-05",
- from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(2 2,2 7,7 7,7 2))"),
- from_wkt<PL>("POLYGON((3 3,6 3,6 6,3 6))"),
- true);
- // polygon with a hole, but the inner ring intersects the other polygon
- tester::apply("pg-pg-06",
- from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(3 2,3 7,7 7,7 2))"),
- from_wkt<PL>("POLYGON((2 3,6 3,6 6,2 6))"),
- false);
- // polygon with a hole, but the other polygon is entirely contained
- // between the inner and outer rings.
- tester::apply("pg-pg-07",
- from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(6 2,6 7,7 7,7 2))"),
- from_wkt<PL>("POLYGON((3 3,5 3,5 6,3 6))"),
- false);
- // polygon with a hole and the outer ring of the other polygon lies
- // between the inner and outer, but without touching either.
- tester::apply("pg-pg-08",
- from_wkt<PL>("POLYGON((0 0,9 0,9 9,0 9),(3 3,3 6,6 6,6 3))"),
- from_wkt<PL>("POLYGON((2 2,7 2,7 7,2 7))"),
- false);
- {
- typedef bg::model::polygon<P> PL; // cw, closed
- // https://svn.boost.org/trac/boost/ticket/10647
- tester::apply("ticket-10647",
- from_wkt<PL>("POLYGON((0 0, 0 5, 5 5, 5 0, 0 0)(1 1, 4 1, 4 4, 1 4, 1 1))"),
- from_wkt<PL>("POLYGON((2 2, 2 3, 3 3, 3 2, 2 2))"),
- true);
- }
- }
- template <typename P>
- inline void test_polygon_multipolygon()
- {
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef bg::model::multi_polygon<PL> MPL;
- typedef test_disjoint tester;
- tester::apply("pg-mpg-01",
- from_wkt<PL>("POLYGON((2 2,2 3,3 3,3 2))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("pg-mpg-02",
- from_wkt<PL>("POLYGON((1 1,1 3,3 3,3 1))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("pg-mpg-03",
- from_wkt<PL>("POLYGON((3 3,3 4,4 4,4 3))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- true);
- }
- template <typename P>
- inline void test_multipolygon_multipolygon()
- {
- typedef bg::model::polygon<P, false, false> PL; // ccw, open
- typedef bg::model::multi_polygon<PL> MPL;
- typedef test_disjoint tester;
- tester::apply("mpg-mpg-01",
- from_wkt<MPL>("MULTIPOLYGON(((2 2,2 3,3 3,3 2)))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("mpg-mpg-02",
- from_wkt<MPL>("MULTIPOLYGON(((1 1,1 3,3 3,3 1)))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- false);
- tester::apply("mpg-mpg-03",
- from_wkt<MPL>("MULTIPOLYGON(((3 3,3 4,4 4,4 3)))"),
- from_wkt<MPL>("MULTIPOLYGON(((0 0,2 0,2 2,0 2)))"),
- true);
- }
- //============================================================================
- template <typename CoordinateType>
- inline void test_areal_areal()
- {
- typedef bg::model::point<CoordinateType, 2, bg::cs::cartesian> point_type;
- test_polygon_polygon<point_type>();
- test_polygon_multipolygon<point_type>();
- test_polygon_ring<point_type>();
- test_polygon_box<point_type>();
- test_multipolygon_multipolygon<point_type>();
- test_multipolygon_ring<point_type>();
- test_multipolygon_box<point_type>();
- test_ring_ring<point_type>();
- test_ring_box<point_type>();
- test_box_box<point_type>();
- }
- //============================================================================
- BOOST_AUTO_TEST_CASE( test_areal_areal_all )
- {
- test_areal_areal<double>();
- test_areal_areal<int>();
- #ifdef HAVE_TTMATH
- test_areal_areal<ttmath_big>();
- #endif
- }
|