libs/geometry/test/strategies/distance_cross_track.cpp

   1 // Boost.Geometry
   2 // Unit Test
   3
   4 // Copyright (c) 2019 Oracle and/or its affiliates.
   5
   6 // Contributed and/or modified by Vissarion Fysikopoulos, on behalf of Oracle
   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
  14 #include <sstream>
  15
  16 #include "../formulas/test_formula.hpp"
  17 #include "distance_cross_track_cases.hpp"
  18
  19 #include <boost/geometry/strategies/strategies.hpp>
  20 #include <boost/geometry/srs/spheroid.hpp>
  21
  22 struct error{
  23     double long_distance;
  24     double short_distance;
  25     double very_short_distance;
  26     double very_very_short_distance;
  27 };
  28
  29 void check_result(double const& result, double const& expected,
  30                   double const& reference, error const& reference_error)
  31 {
  32     BOOST_GEOMETRY_CHECK_CLOSE(result, expected, 0.1,
  33     std::setprecision(20) << "result {" << result
  34                           << "} different than expected {" << expected << "}.");
  35
  36     double reference_error_value = result > 2000 ? reference_error.long_distance
  37                                  : result > 100  ? reference_error.short_distance
  38                                  : result > 20   ? reference_error.very_short_distance
  39                                  : reference_error.very_very_short_distance;
  40
  41     BOOST_GEOMETRY_CHECK_CLOSE(result, reference, reference_error_value,
  42         std::setprecision(20) << "result {" << result
  43                               << "} different than reference {"
  44                               << reference << "}.");
  45 }
  46
  47 template <typename Point>
  48 void test_all(expected_results const& results)
  49 {
  50     double const d2r = bg::math::d2r<double>();
  51
  52     double lon1r = results.p1.lon * d2r;
  53     double lat1r = results.p1.lat * d2r;
  54     double lon2r = results.p2.lon * d2r;
  55     double lat2r = results.p2.lat * d2r;
  56     double lon3r = results.p3.lon * d2r;
  57     double lat3r = results.p3.lat * d2r;
  58
  59     typedef bg::srs::spheroid<double> Spheroid;
  60
  61     // WGS84
  62     Spheroid spheroid(6378137.0, 6356752.3142451793);
  63
  64     error errors [] =
  65     {
  66         {0.00000001, 0.00000001, 0.00000001, 0.000001}, //vincenty
  67         {0.0002, 0.002, 0.01, 0.2}, //thomas
  68         {0.002, 0.4, 15, 25}, //andoyer
  69         {1, 6, 15, 200} //spherical
  70     };
  71
  72     //vincenty
  73     double distance = bg::strategy::distance::detail::geographic_cross_track<bg::strategy::vincenty, Spheroid, double, true>(spheroid)
  74             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
  75     check_result(distance, results.vincenty_bisection, results.reference, errors[0]);
  76
  77     distance = bg::strategy::distance::geographic_cross_track<bg::strategy::vincenty, Spheroid, double>(spheroid)
  78             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
  79     check_result(distance, results.vincenty, results.reference, errors[0]);
  80
  81     //thomas
  82     distance = bg::strategy::distance::detail::geographic_cross_track<bg::strategy::thomas, Spheroid, double, true>(spheroid)
  83             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
  84     check_result(distance, results.thomas_bisection, results.reference, errors[1]);
  85
  86     distance = bg::strategy::distance::geographic_cross_track<bg::strategy::thomas, Spheroid, double>(spheroid)
  87             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
  88     check_result(distance, results.thomas, results.reference, errors[1]);
  89
  90     //andoyer
  91     distance = bg::strategy::distance::detail::geographic_cross_track<bg::strategy::andoyer, Spheroid, double, true>(spheroid)
  92             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
  93     check_result(distance, results.andoyer_bisection, results.reference, errors[2]);
  94
  95     distance = bg::strategy::distance::geographic_cross_track<bg::strategy::andoyer, Spheroid, double>(spheroid)
  96             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
  97     check_result(distance, results.andoyer, results.reference, errors[2]);
  98
  99     //spherical
 100     distance = bg::strategy::distance::cross_track<>(bg::formula::mean_radius<double>(spheroid))
 101             .apply(Point(lon3r, lat3r), Point(lon1r, lat1r), Point(lon2r, lat2r));
 102     check_result(distance, results.spherical, results.reference, errors[3]);
 103
 104 }
 105
 106 int test_main(int, char*[])
 107 {
 108     typedef bg::model::point<double, 2, bg::cs::geographic<bg::radian> > point;
 109
 110     for (size_t i = 0; i < expected_size; ++i)
 111     {
 112         test_all<point>(expected[i]);
 113     }
 114
 115     return 0;
 116 }