ceph/src/boost/libs/geometry/test/algorithms/set_operations/union/union_aa_sph.cpp

   1 // Boost.Geometry
   2 // Unit Test
   3
   4 // Copyright (c) 2017 Adam Wulkiewicz, Lodz, Poland.
   5
   6 // Copyright (c) 2017-2019, Oracle and/or its affiliates.
   7 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
   8
   9 // Use, modification and distribution is subject to the Boost Software License,
  10 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
  11 // http://www.boost.org/LICENSE_1_0.txt)
  12
  13 #include <iostream>
  14 #include <string>
  15
  16 #include "test_union.hpp"
  17
  18
  19 struct exterior_points_counter
  20 {
  21     exterior_points_counter() : count(0) {}
  22
  23     template <typename Polygon>
  24     void operator()(Polygon const& poly)
  25     {
  26         count += boost::size(bg::exterior_ring(poly));
  27     }
  28
  29     std::size_t count;
  30 };
  31
  32 struct interiors_counter
  33     : exterior_points_counter
  34 {
  35     template <typename Polygon>
  36     void operator()(Polygon const& poly)
  37     {
  38         count += boost::size(bg::interior_rings(poly));
  39     }
  40 };
  41
  42 void test_spherical_one(std::string const& wkt1, std::string const& wkt2,
  43                         std::size_t count, std::size_t exterior_points_count, std::size_t interiors_count,
  44                         double expected_area)
  45 {
  46     typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > point;
  47     typedef bg::model::polygon<point> polygon;
  48     typedef bg::model::multi_polygon<polygon> multipolygon;
  49
  50     polygon p1, p2;
  51
  52     boost::geometry::read_wkt(wkt1, p1);
  53     boost::geometry::read_wkt(wkt2, p2);
  54
  55     multipolygon result;
  56
  57     boost::geometry::union_(p1, p2, result);
  58
  59     double result_area = bg::area(result);
  60
  61     std::size_t result_count = boost::size(result);
  62     std::size_t result_exterior_points = std::for_each(boost::begin(result),
  63                                                        boost::end(result),
  64                                                        exterior_points_counter()).count;
  65     std::size_t result_interiors = std::for_each(boost::begin(result),
  66                                                  boost::end(result),
  67                                                  interiors_counter()).count;
  68     if (expected_area >= 0)
  69     {
  70         BOOST_CHECK_EQUAL(result_count, count);
  71         BOOST_CHECK_EQUAL(result_exterior_points, exterior_points_count);
  72         BOOST_CHECK_EQUAL(result_interiors, interiors_count);
  73         BOOST_CHECK_CLOSE(result_area, expected_area, 0.001);
  74     }
  75     else
  76     {
  77         BOOST_CHECK_EQUAL(result_count, 0u);
  78         BOOST_CHECK_EQUAL(result_area, 0.0);
  79     }
  80 }
  81
  82
  83 void test_spherical()
  84 {
  85     // https://github.com/boostorg/geometry/issues/475
  86     test_spherical_one("POLYGON((-78.4072265625001 43.06652924482626,-78.4072265625 43.06740063068311,-78.4063141178686 43.06653210403569,-78.4072265625001 43.06652924482626))",
  87                        "POLYGON((-78.55968743491499 43.06594969590624,-78.55036227331367 43.07380195109801,-78.53503704605811 43.08248347074284,-78.51769210872999 43.08880392487917,-78.49899564953199 43.09251971058174,-78.47966844278045 43.09348761253013,-78.46045580120891 43.09167037120638,-78.44209853911326 43.08713812460473,-78.42530412309867 43.08006566649393,-78.41071917768537 43.07072563376782,-78.40631359930124 43.06653210565861,-78.55968743491499 43.06594969590624))",
  88                        1, 12, 0, 0.00000064324358632259458);
  89     test_spherical_one("POLYGON((-121.5731370931394 37.95996093777254,-121.5731374919342 37.9599609375,-121.5774655572748 37.96661505459264,-121.5814229252572 37.9775390625,-121.5814227715605 37.97753906261685,-121.5824154607857 37.98027887997713,-121.5839036001453 37.9944445317401,-121.5835741960922 37.99673349281881,-121.4641571044922 37.9967043028667,-121.4641571044935 37.97756945550759,-121.4641571044935 37.95998526681451,-121.4641571044935 37.92462240375453,-121.4758884395752 37.92196920341666,-121.493771910322 37.92080607027856,-121.5115987448441 37.92242636171255,-121.5286857982694 37.92676802819236,-121.5443780170238 37.93366470898578,-121.5580733707665 37.94285205174562,-121.5692458339493 37.95397776814401,-121.5731370931394 37.95996093777254))",
  90                        "POLYGON((-121.6413116455078 37.9967043028667,-121.5810851310329 37.99673411655077,-121.5697730327269 37.99673629028738,-121.5608282121158 37.99281120211074,-121.5471234211378 37.98362406015312,-121.5359500167068 37.97249627286875,-121.527736882307 37.95985622934639,-121.5227984361816 37.94619029268395,-121.5213227340147 37.93202401606796,-121.5233645060913 37.91790188887406,-121.5288433558941 37.90436639967157,-121.537547143277 37.89193722240091,-121.5491403777151 37.8810913202222,-121.5572385386645 37.87598800167662,-121.5490395877234 37.86781333310999,-121.5420190642174 37.85699822330511,-121.6413116455078 37.85696553311346,-121.6413116455089 37.8645849275209,-121.6413116455089 37.874858067486,-121.6413116455089 37.89983322075715,-121.6413116455089 37.95996381547175,-121.6413116455089 37.97754168382039,-121.6413116455089 37.97754250238781,-121.6413116455078 37.9967043028667))",
  91                        1, 24, 0, 0.0000047757848548868592);
  92     test_spherical_one("POLYGON((-121.3246465532687 35.88962804614829,-121.2883758544922 35.88960686771812,-121.2883758544928 35.85058599290677,-121.290026580334 35.8505859375,-121.3246465532687 35.88962804614829))",
  93                        "POLYGON((-121.2883758544928 35.85055605452366,-121.29 35.850556,-121.3056335092293 35.8681640625,-121.2883758544928 35.86816464188531,-121.2883758544928 35.85055605452366))",
  94                        1, 8, 0, 0.00000018266345129866598);
  95     test_spherical_one("POLYGON((-120.4666783969804 37.2997584488572,-120.466678397 37.3000076100001,-120.466678399 37.3166742810001,-120.450011727 37.3166742830001,-120.4500117249686 37.2997528818001,-120.4666783969804 37.2997584488572))",
  96                        "POLYGON((-120.466678399 37.3166742810001,-120.466678396991 37.29975844885721,-120.566678432986 37.29974277357464,-120.477336505269 37.35026345984599,-120.471699839074 37.3480665431537,-120.466678402489 37.345749649536,-120.466269200472 37.345560843499,-120.461071086478 37.3427585963354,-120.456130850491 37.3396734837438,-120.45147257817 37.3363205671631,-120.447118969699 37.3327162133727,-120.443091229047 37.3288780141012,-120.439408960746 37.3248246996654,-120.4360900747 37.3205760470716,-120.4334972505575 37.31667428399087,-120.4333450539725 37.31644525379688,-120.433150699472 37.3161527830495,-120.43060510448 37.3115764824821,-120.428465631443 37.306869462763,-120.426742635401 37.3020546745901,-120.4261294725058 37.29974083085035,-120.4333450519679 37.29974497771207,-120.4500117249686 37.2997528818001,-120.450011727 37.3166742830001,-120.466678399 37.3166742810001))",
  97                        1, 20, 0, 0.00000096333365336696321);
  98 }
  99
 100
 101 int test_main(int, char* [])
 102 {
 103     test_spherical();
 104
 105     return 0;
 106 }