| // Boost.Geometry (aka GGL, Generic Geometry Library) |
| // Unit Test |
| |
| // Copyright (c) 2014-2015, Oracle and/or its affiliates. |
| |
| // Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle |
| |
| // Licensed under the Boost Software License version 1.0. |
| // http://www.boost.org/users/license.html |
| |
| #include <iostream> |
| |
| #ifndef BOOST_TEST_MODULE |
| #define BOOST_TEST_MODULE test_distance_areal_areal |
| #endif |
| |
| #include <boost/test/included/unit_test.hpp> |
| |
| #include "test_distance_common.hpp" |
| |
| |
| typedef bg::model::point<int,2,bg::cs::cartesian> int_point_type; |
| typedef bg::model::point<double,2,bg::cs::cartesian> point_type; |
| typedef bg::model::polygon<point_type, false> polygon_type; |
| typedef bg::model::multi_polygon<polygon_type> multi_polygon_type; |
| typedef bg::model::ring<point_type, false> ring_type; |
| typedef bg::model::box<int_point_type> int_box_type; |
| typedef bg::model::box<point_type> box_type; |
| |
| namespace services = bg::strategy::distance::services; |
| typedef bg::default_distance_result<point_type>::type return_type; |
| |
| typedef bg::strategy::distance::projected_point<> point_segment_strategy; |
| typedef bg::strategy::distance::pythagoras_box_box<> box_box_strategy; |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_polygon_polygon(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "polygon/polygon distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<polygon_type, polygon_type> tester; |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,5 20,5 25,-5 25,-5 20))", |
| 10, 100, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,-5 5,5 5,5 20,-5 20))", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,-5 -20,5 -20,5 20,-5 20))", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10),\ |
| (-5 -5,-5 5,5 5,5 -5,-5 -5))", |
| "polygon((-1 -1,0 0,-1 0,-1 -1))", |
| 4, 16, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_polygon_multipolygon(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "polygon/multipolygon distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries |
| < |
| polygon_type, multi_polygon_type |
| > tester; |
| |
| tester::apply("polygon((12 0,14 0,19 0,19.9 -1,12 0))", |
| "multipolygon(((-10 -10,10 -10,10 10,-10 10,-10 -10)),\ |
| ((20 -1,21 2,30 -10,20 -1)))", |
| 0.1, 0.01, strategy); |
| |
| tester::apply("polygon((19 0,19.9 -1,12 0,20.5 0.5,19 0))", |
| "multipolygon(((-10 -10,10 -10,10 10,-10 10,-10 -10)),\ |
| ((20 -1,21 2,30 -10,20 -1)))", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_polygon_ring(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "polygon/ring distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<polygon_type, ring_type> tester; |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,5 20,5 25,-5 25,-5 20))", |
| 10, 100, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,-5 5,5 5,5 20,-5 20))", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,-5 -20,5 -20,5 20,-5 20))", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_multipolygon_multipolygon(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "multipolygon/multipolygon distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries |
| < |
| multi_polygon_type, multi_polygon_type |
| > tester; |
| |
| tester::apply("multipolygon(((12 0,14 0,14 1,12 0)),\ |
| ((18 0,19 0,19.9 -1,18 0)))", |
| "multipolygon(((-10 -10,10 -10,10 10,-10 10,-10 -10)),\ |
| ((20 -1,21 2,30 -10,20 -1)))", |
| 0.1, 0.01, strategy); |
| |
| tester::apply("multipolygon(((18 0,19 0,19.9 -1,18 0)),\ |
| ((12 0,14 0,20.5 0.5,12 0)))", |
| "multipolygon(((-10 -10,10 -10,10 10,-10 10,-10 -10)),\ |
| ((20 -1,21 2,30 -10,20 -1)))", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_multipolygon_ring(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "multipolygon/ring distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries |
| < |
| multi_polygon_type, ring_type |
| > tester; |
| |
| tester::apply("multipolygon(((-10 -10,10 -10,10 10,-10 10,-10 -10)),\ |
| ((20 -1,21 2,30 -10,20 -1)))", |
| "polygon((12 0,14 0,19 0,19.9 -1,12 0))", |
| 0.1, 0.01, strategy); |
| |
| tester::apply("multipolygon(((-10 -10,10 -10,10 10,-10 10,-10 -10)),\ |
| ((20 -1,21 2,30 -10,20 -1)))", |
| "polygon((19 0,19.9 -1,12 0,20.5 0.5,19 0))", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_ring_ring(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "ring/ring distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<ring_type, ring_type> tester; |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,5 20,5 25,-5 25,-5 20))", |
| 10, 100, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,-5 5,5 5,5 20,-5 20))", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((-10 -10,10 -10,10 10,-10 10,-10 -10))", |
| "polygon((-5 20,-5 -20,5 -20,5 20,-5 20))", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_box_box(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "box/box distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<int_box_type, int_box_type> int_tester; |
| typedef test_distance_of_geometries<box_type, box_type> tester; |
| |
| int_tester::apply("box(5 5,10 10)", |
| "box(0 0,1 1)", |
| sqrt(32.0), 32, strategy); |
| |
| tester::apply("box(5 5,10 10)", |
| "box(0 0,1 1)", |
| sqrt(32.0), 32, strategy); |
| |
| tester::apply("box(3 8,13 18)", |
| "box(0 0,5 5)", |
| 3, 9, strategy); |
| |
| tester::apply("box(5 5,10 10)", |
| "box(0 0,5 5)", |
| 0, 0, strategy); |
| |
| tester::apply("box(5 5,10 10)", |
| "box(0 0,6 6)", |
| 0, 0, strategy); |
| |
| tester::apply("box(3 5,13 15)", |
| "box(0 0,5 5)", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_polygon_box(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "polygon/box distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<polygon_type, box_type> tester; |
| |
| tester::apply("polygon((10 10,10 5,5 5,5 10,10 10))", |
| "box(0 0,1 1)", |
| sqrt(32.0), 32, strategy); |
| |
| tester::apply("polygon((10 10,10 5,5 5,5 10,10 10))", |
| "box(0 0,5 5)", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((10 10,10 5,5 5,5 10,10 10))", |
| "box(0 0,6 6)", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((10 10,15 5,10 0,5 5,10 10))", |
| "box(5 0,7.5 2.5)", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((10 10,15 5,10 0,5 5,10 10))", |
| "box(5 0,6 1)", |
| sqrt(4.5), 4.5, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_multipolygon_box(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "multipolygon/box distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<multi_polygon_type, box_type> tester; |
| |
| tester::apply("multipolygon(((-10 -10,-10 -9,-9 -9,-9 -10,-10 -10)),\ |
| ((2 2,2 3,3 3,3 2,2 2)))", |
| "box(0 0,1 1)", |
| sqrt(2.0), 2, strategy); |
| |
| tester::apply("multipolygon(((-10 -10,-10 -9,-9 -9,-9 -10,-10 -10)),\ |
| ((2 2,2 3,3 3,3 2,2 2)))", |
| "box(0 0,2 2)", |
| 0, 0, strategy); |
| |
| tester::apply("multipolygon(((-10 -10,-10 -9,-9 -9,-9 -10,-10 -10)),\ |
| ((2 2,2 3,3 3,3 2,2 2)))", |
| "box(0 0,2.5 2)", |
| 0, 0, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Strategy> |
| void test_distance_ring_box(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "ring/box distance tests" << std::endl; |
| #endif |
| typedef test_distance_of_geometries<ring_type, box_type> tester; |
| |
| tester::apply("polygon((10 10,10 5,5 5,5 10,10 10))", |
| "box(0 0,1 1)", |
| sqrt(32.0), 32, strategy); |
| |
| tester::apply("polygon((10 10,10 5,5 5,5 10,10 10))", |
| "box(0 0,5 5)", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((10 10,10 5,5 5,5 10,10 10))", |
| "box(0 0,6 6)", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((10 10,15 5,10 0,5 5,10 10))", |
| "box(5 0,7.5 2.5)", |
| 0, 0, strategy); |
| |
| tester::apply("polygon((10 10,15 5,10 0,5 5,10 10))", |
| "box(5 0,6 1)", |
| sqrt(4.5), 4.5, strategy); |
| } |
| |
| //=========================================================================== |
| |
| template <typename Point, typename Strategy> |
| void test_more_empty_input_areal_areal(Strategy const& strategy) |
| { |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << std::endl; |
| std::cout << "testing on empty inputs... " << std::flush; |
| #endif |
| bg::model::polygon<Point> polygon_empty; |
| bg::model::multi_polygon<bg::model::polygon<Point> > multipolygon_empty; |
| bg::model::ring<Point> ring_empty; |
| |
| bg::model::polygon<Point> polygon = |
| from_wkt<bg::model::polygon<Point> >("polygon((0 0,1 0,0 1))"); |
| |
| bg::model::multi_polygon<bg::model::polygon<Point> > multipolygon = |
| from_wkt |
| < |
| bg::model::multi_polygon<bg::model::polygon<Point> > |
| >("multipolygon(((0 0,1 0,0 1)))"); |
| |
| bg::model::ring<Point> ring = |
| from_wkt<bg::model::ring<Point> >("polygon((0 0,1 0,0 1))"); |
| |
| // 1st geometry is empty |
| test_empty_input(polygon_empty, polygon, strategy); |
| test_empty_input(polygon_empty, multipolygon, strategy); |
| test_empty_input(polygon_empty, ring, strategy); |
| test_empty_input(multipolygon_empty, polygon, strategy); |
| test_empty_input(multipolygon_empty, multipolygon, strategy); |
| test_empty_input(multipolygon_empty, ring, strategy); |
| test_empty_input(ring_empty, polygon, strategy); |
| test_empty_input(ring_empty, multipolygon, strategy); |
| test_empty_input(ring_empty, ring, strategy); |
| |
| // 2nd geometry is empty |
| test_empty_input(polygon, polygon_empty, strategy); |
| test_empty_input(polygon, multipolygon_empty, strategy); |
| test_empty_input(polygon, ring_empty, strategy); |
| test_empty_input(multipolygon, polygon_empty, strategy); |
| test_empty_input(multipolygon, multipolygon_empty, strategy); |
| test_empty_input(multipolygon, ring_empty, strategy); |
| test_empty_input(ring, polygon_empty, strategy); |
| test_empty_input(ring, multipolygon_empty, strategy); |
| test_empty_input(ring, ring_empty, strategy); |
| |
| // both geometries are empty |
| test_empty_input(polygon_empty, polygon_empty, strategy); |
| test_empty_input(polygon_empty, multipolygon_empty, strategy); |
| test_empty_input(polygon_empty, ring_empty, strategy); |
| test_empty_input(multipolygon_empty, polygon_empty, strategy); |
| test_empty_input(multipolygon_empty, multipolygon_empty, strategy); |
| test_empty_input(multipolygon_empty, ring_empty, strategy); |
| test_empty_input(ring_empty, polygon_empty, strategy); |
| test_empty_input(ring_empty, multipolygon_empty, strategy); |
| test_empty_input(ring_empty, ring_empty, strategy); |
| |
| #ifdef BOOST_GEOMETRY_TEST_DEBUG |
| std::cout << "done!" << std::endl; |
| #endif |
| } |
| |
| //=========================================================================== |
| |
| BOOST_AUTO_TEST_CASE( test_all_polygon_polygon ) |
| { |
| test_distance_polygon_polygon(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_polygon_multipolygon ) |
| { |
| test_distance_polygon_multipolygon(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_polygon_ring ) |
| { |
| test_distance_polygon_ring(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_multipolygon_multipolygon ) |
| { |
| test_distance_multipolygon_multipolygon(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_multipolygon_ring ) |
| { |
| test_distance_multipolygon_ring(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_ring_ring ) |
| { |
| test_distance_ring_ring(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_box_box ) |
| { |
| test_distance_box_box(box_box_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_polygon_box ) |
| { |
| test_distance_polygon_box(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_multipolygon_box ) |
| { |
| test_distance_multipolygon_box(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_ring_box ) |
| { |
| test_distance_ring_box(point_segment_strategy()); |
| } |
| |
| BOOST_AUTO_TEST_CASE( test_all_empty_input_areal_areal ) |
| { |
| test_more_empty_input_areal_areal<point_type>(point_segment_strategy()); |
| } |