ceph/src/boost/libs/geometry/test/srs/projections_combined.cpp

   1 // Boost.Geometry (aka GGL, Generic Geometry Library)
   2 // Unit Test
   3
   4 // Copyright (c) 2015 Barend Gehrels, Amsterdam, the Netherlands.
   5
   6 // This file was modified by Oracle on 2017, 2018.
   7 // Modifications copyright (c) 2017-2018, Oracle and/or its affiliates.
   8 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
   9
  10 // Use, modification and distribution is subject to the Boost Software License,
  11 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
  12 // http://www.boost.org/LICENSE_1_0.txt)
  13
  14
  15 #if defined(_MSC_VER)
  16 #pragma warning( disable : 4305 ) // truncation double -> float
  17 #endif // defined(_MSC_VER)
  18
  19 #define BOOST_GEOMETRY_SRS_ENABLE_STATIC_PROJECTION_HYBRID_INTERFACE
  20
  21 #include <geometry_test_common.hpp>
  22
  23 #include <boost/geometry/srs/projection.hpp>
  24
  25 #include <boost/geometry/core/coordinate_type.hpp>
  26 #include <boost/geometry/algorithms/make.hpp>
  27
  28 #include <boost/geometry.hpp>
  29 #include <boost/geometry/geometries/geometries.hpp>
  30 #include <boost/geometry/geometries/point_xy.hpp>
  31
  32 #include "proj4.hpp"
  33
  34 namespace srs = bg::srs;
  35
  36 template <typename P>
  37 bool check_expected(P const& p1, P const& p2)
  38 {
  39     return bg::math::abs(bg::get<0>(p1) - bg::get<0>(p2)) <= 1.0
  40         && bg::math::abs(bg::get<1>(p1) - bg::get<1>(p2)) <= 1.0;
  41 }
  42
  43 template <typename StaticParams, typename GeoPoint, typename CartPoint>
  44 void test_forward(std::string const& id, GeoPoint const& geo_point, CartPoint const& cart_point,
  45                   StaticParams const& params, std::string const& proj4 = "")
  46 {
  47     try
  48     {
  49         srs::projection<StaticParams> prj = params;
  50
  51         CartPoint xy;
  52         prj.forward(geo_point, xy);
  53
  54         bool ok = check_expected(xy, cart_point);
  55
  56         BOOST_CHECK_MESSAGE(ok,
  57                 " id: " << id
  58                 << " point: " << bg::wkt(xy)
  59                 << " different than expected: " << bg::wkt(cart_point));
  60
  61         if (! proj4.empty())
  62         {
  63             srs::projection<> prj2 = srs::proj4(proj4);
  64
  65             CartPoint xy2;
  66             prj2.forward(geo_point, xy2);
  67
  68             bool eq2 = bg::equals(xy, xy2);
  69
  70             BOOST_CHECK_MESSAGE(eq2,
  71                 " id: " << id << " result of static: "
  72                 << bg::wkt(xy) << " different than dynamic: " << bg::wkt(xy2));
  73
  74 #ifdef TEST_WITH_PROJ4
  75             pj_projection prj3(proj4);
  76
  77             CartPoint xy3;
  78             prj3.forward(geo_point, xy3);
  79
  80             bool eq3 = bg::equals(xy, xy3);
  81
  82             BOOST_CHECK_MESSAGE(eq3,
  83                 " id: " << id << " result: "
  84                 << bg::wkt(xy) << " different than proj4: " << bg::wkt(xy3));
  85 #endif // TEST_WITH_PROJ4
  86         }
  87     }
  88     catch(bg::projection_exception const& e)
  89     {
  90         std::cout << "Exception in " << id << " : " << e.code() << std::endl;
  91     }
  92     catch(...)
  93     {
  94         std::cout << "Exception (unknown) in " << id << std::endl;
  95     }
  96 }
  97
  98 template <typename StaticParams, typename GeoPoint, typename CartPoint>
  99 void test_forward(std::string const& id, GeoPoint const& geo_point, CartPoint const& cart_point,
 100                   std::string const& proj4 = "")
 101 {
 102     test_forward(id, geo_point, cart_point, StaticParams(), proj4);
 103 }
 104
 105 template <typename T>
 106 void test_all()
 107 {
 108     typedef bg::model::point<T, 2, bg::cs::geographic<bg::degree> > geo_point_type;
 109     typedef bg::model::point<T, 2, bg::cs::cartesian> cart;
 110
 111     geo_point_type amsterdam = bg::make<geo_point_type>(4.8925, 52.3731);
 112     geo_point_type utrecht   = bg::make<geo_point_type>(5.1213, 52.0907);
 113
 114     geo_point_type anchorage = bg::make<geo_point_type>(-149.90, 61.22);
 115     geo_point_type juneau    = bg::make<geo_point_type>(-134.42, 58.30);
 116
 117     geo_point_type auckland   = bg::make<geo_point_type>(174.74, -36.84);
 118     geo_point_type wellington = bg::make<geo_point_type>(177.78, -41.29);
 119
 120     geo_point_type aspen  = bg::make<geo_point_type>(-106.84, 39.19);
 121     geo_point_type denver = bg::make<geo_point_type>(-104.88, 39.76);
 122
 123     using namespace srs::spar;
 124
 125     // IGH (internally using moll/sinu)
 126     {
 127         typedef parameters<proj_igh, ellps_sphere, units_m> params_t;
 128         std::string sparams = "+proj=igh +ellps=sphere +units=m";
 129         test_forward<params_t>("igh-am", amsterdam, cart(1489299.1509520211, 5776413.4260336142), sparams);
 130         test_forward<params_t>("igh-ut", utrecht, cart(1498750.6627020084, 5747394.3313896423), sparams);
 131         test_forward<params_t>("igh-as", aspen, cart(-11708973.126426676, 4357727.1232166551), sparams);
 132         test_forward<params_t>("igh-de", denver, cart(-11536624.264589204, 4421108.2015589233), sparams);
 133         test_forward<params_t>("igh-au", auckland, cart(18658819.353676274, -4096419.1686476548), sparams);
 134         test_forward<params_t>("igh-we", wellington, cart(18733710.557981707, -4591140.1631184481), sparams);
 135         test_forward<params_t>("igh-an", anchorage, cart(-14275110.630537530, 6648284.9393376000), sparams);
 136         test_forward<params_t>("igh-ju", juneau, cart(-13421076.123140398, 6368936.3597440729), sparams);
 137     }
 138
 139     // Using moll
 140     {
 141         typedef parameters<proj_ob_tran, ellps_wgs84, o_proj<proj_moll>, units_m, o_lat_p<>, o_lon_p<> > params_t;
 142         params_t params = params_t(proj_ob_tran(), ellps_wgs84(), o_proj<proj_moll>(), units_m(), o_lat_p<>(10), o_lon_p<>(90));
 143         std::string sparams = "+proj=ob_tran +ellps=WGS84 +o_proj=moll +units=m +o_lat_p=10 +o_lon_p=90 +o_lon_o=11.50";
 144         test_forward<params_t>("obt-m-am", amsterdam, cart(8688778.3518596273, -3348126.4518623645), params, sparams);
 145         test_forward<params_t>("obt-m-ut", utrecht, cart(8693109.5205437448, -3379708.1134765535), params, sparams);
 146         test_forward<params_t>("obt-m-as", aspen, cart(3691751.3259231807, 2371456.9674431868), params, sparams);
 147         test_forward<params_t>("obt-m-de", denver, cart(3764685.2104777521, 2185616.0182080171), params, sparams);
 148     }
 149
 150     // Using sinu
 151     {
 152         // In Proj4 >= 5.0.0 ob_tran projection doesn't overwrite the ellipsoid's parameters to make sphere in underlying projection
 153         // So in order to use spherical sinu projection WGS84-compatible sphere has to be set manually.
 154         typedef parameters<proj_ob_tran, a<>, es<>, o_proj<proj_sinu>, units_m, o_lat_p<>, o_lon_p<> > params_t;
 155         params_t params = params_t(proj_ob_tran(), a<>(6378137), es<>(0), o_proj<proj_sinu>(), units_m(), o_lat_p<>(10), o_lon_p<>(90));
 156         std::string sparams = "+proj=ob_tran +a=6378137 +es=0 +o_proj=sinu +units=m +o_lat_p=10 +o_lon_p=90 +o_lon_o=11.50";
 157         test_forward<params_t>("obt-s-am", amsterdam, cart(9220221.4221933037, -3059652.3579233172), params, sparams);
 158         test_forward<params_t>("obt-s-ut", utrecht, cart(9216281.0977674071, -3089427.4415689218), params, sparams);
 159         test_forward<params_t>("obt-s-as", aspen, cart(4010672.3356677019, 2150730.9484995930), params, sparams);
 160         test_forward<params_t>("obt-s-de", denver, cart(4103945.8062708224, 1979964.9315176210), params, sparams);
 161     }
 162 }
 163
 164 int test_main(int, char* [])
 165 {
 166     test_all<double>();
 167
 168     return 0;
 169 }