Merge pull request #1320 from vissarion/add_test_cases

[test] Add test case for goegrphic intersection for polygons touching North Pole

test/algorithms/intersects/intersects_sph_geo.cpp

@@ -79,3 +79,29 @@ int test_main( int , char* [] )
 
     return 0;
 }
+
+//https://github.com/boostorg/geometry/issues/480
+namespace spherical {
+  using point = boost::geometry::model::d2::point_xy<double, boost::geometry::cs::spherical_equatorial<boost::geometry::degree> >;
+  using linestring = boost::geometry::model::linestring<point>;
+  using polygon = boost::geometry::model::polygon<point>;
+  using multipolygon = boost::geometry::model::multi_polygon<polygon>;
+  using mpoint = boost::geometry::model::multi_point<point>;
+}
+
+BOOST_AUTO_TEST_CASE(usaIntersectBigTile90) {
+  auto const big_tile_wkt = "POLYGON((-182.25 -2.25,-182.25 85,-90 85,4.5 85,4.5 -2.25,-90 -2.25,-182.25 -2.25))";
+  auto const big_tile_wkt_on90 = "POLYGON((-182.25 -2.25,-182.25 90,-90 90,4.5 90,4.5 -2.25,-90 -2.25,-182.25 -2.25))";
+  auto const usa_wkt = "POLYGON((-125.0840938165823 48.18451089135078,-123.2484442951568 48.28402780487058,-123.3222396381058 49.00207166780857,-106.2294164816153 48.99936106028161,-94.95738881549684 49.37019439405368,-91.56749988130491 48.04377777876205,-88.36986182651864 48.30606296052222,-82.51861102256366 45.33863888434129,-82.67972207213995 41.67655556472526,-74.97255548445517 44.98347221572623,-71.50101930291031 45.01344579588383,-69.22446981284281 47.4598397141629,-67.79019834014535 47.06722478688557,-67.04184250949656 44.36199518537477,-70.23005912341451 43.2124078816444,-69.75267183371977 41.13341339615248,-73.71086681700677 40.38246941502378,-75.73925419455345 36.92336928786909,-75.32570262625269 35.10289000177789,-81.00038138644489 31.36986317745486,-80.21558018606714 24.87617100882401,-82.31642199133823 24.40399007010678,-81.37277484337314 25.35110770766313,-82.98991779576328 28.87424738341057,-86.63551567900944 30.19587268651689,-88.60142307258893 30.01998991835538,-89.40197910608261 28.705405956028,-94.4542152059084 29.30489831467678,-96.8117767445914 27.80056771982362,-97.42440989188675 25.84041996728249,-99.0845302237628 26.39797732094523,-101.4011198718317 29.77120000013631,-103.2862798604784 28.97691995351909,-106.4506798904261 31.76422996252971,-111.074824919583 31.33223902311842,-114.7199604399525 32.71865525370477,-117.4636899326359 32.58936795833381,-118.1537433047147 33.49758670948682,-118.5200689113933 32.61376470726387,-119.9333407371885 33.34100477748434,-118.6596572250774 33.83057945961541,-120.7240783642898 33.9956452157932,-124.7155563780111 40.37995067155807,-125.0840938165823 48.18451089135078))";
+
+  spherical::polygon big_tile, big_tile_on90, usa;
+  boost::geometry::read_wkt(big_tile_wkt, big_tile);
+  boost::geometry::read_wkt(big_tile_wkt_on90, big_tile_on90);
+  boost::geometry::read_wkt(usa_wkt, usa);
+
+  spherical::multipolygon out1, out2;
+  boost::geometry::intersection(big_tile_on90, usa, out1);
+  boost::geometry::intersection(big_tile, usa, out2);
+  BOOST_CHECK(out1.size() == 1);
+  BOOST_CHECK(boost::geometry::area(out1) == boost::geometry::area(out2));
+}