123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153 |
- // Boost.Geometry (aka GGL, Generic Geometry Library)
- // Unit Test
- // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
- // Copyright (c) 2008-2012 Bruno Lalande, Paris, France.
- // Copyright (c) 2009-2012 Mateusz Loskot, London, UK.
- // This file was modified by Oracle on 2015, 2017.
- // Modifications copyright (c) 2015-2017 Oracle and/or its affiliates.
- // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
- // Parts of Boost.Geometry are redesigned from Geodan's Geographic Library
- // (geolib/GGL), copyright (c) 1995-2010 Geodan, Amsterdam, the Netherlands.
- // Use, modification and distribution is subject to the Boost Software License,
- // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
- // http://www.boost.org/LICENSE_1_0.txt)
- #include <geometry_test_common.hpp>
- #include <boost/concept_check.hpp>
- #include <boost/geometry/algorithms/assign.hpp>
- #include <boost/geometry/algorithms/distance.hpp>
- #include <boost/geometry/geometries/point.hpp>
- #include <boost/geometry/srs/spheroid.hpp>
- #include <boost/geometry/strategies/concepts/distance_concept.hpp>
- #include <boost/geometry/strategies/geographic/distance_thomas.hpp>
- #include <boost/geometry/strategies/geographic/side_thomas.hpp>
- #include <test_common/test_point.hpp>
- #ifdef HAVE_TTMATH
- # include <boost/geometry/extensions/contrib/ttmath_stub.hpp>
- #endif
- template <typename P1, typename P2>
- void test_distance(double lon1, double lat1, double lon2, double lat2, double expected_km)
- {
- // Set radius type, but for integer coordinates we want to have floating point radius type
- typedef typename bg::promote_floating_point
- <
- typename bg::coordinate_type<P1>::type
- >::type rtype;
- typedef bg::srs::spheroid<rtype> stype;
- typedef bg::strategy::distance::thomas<stype> thomas_type;
- typedef bg::strategy::distance::geographic<bg::strategy::thomas, stype> geographic_type;
- BOOST_CONCEPT_ASSERT
- (
- (bg::concepts::PointDistanceStrategy<thomas_type, P1, P2>)
- );
- thomas_type thomas;
- geographic_type geographic;
- typedef typename bg::strategy::distance
- ::services::return_type<thomas_type, P1, P2>::type return_type;
- P1 p1;
- P2 p2;
- bg::assign_values(p1, lon1, lat1);
- bg::assign_values(p2, lon2, lat2);
- BOOST_CHECK_CLOSE(thomas.apply(p1, p2), return_type(1000.0 * expected_km), 0.001);
- BOOST_CHECK_CLOSE(geographic.apply(p1, p2), return_type(1000.0 * expected_km), 0.001);
- BOOST_CHECK_CLOSE(bg::distance(p1, p2, thomas), return_type(1000.0 * expected_km), 0.001);
- }
- template <typename PS, typename P>
- void test_side(double lon1, double lat1,
- double lon2, double lat2,
- double lon, double lat,
- int expected_side)
- {
- // Set radius type, but for integer coordinates we want to have floating point radius type
- typedef typename bg::promote_floating_point
- <
- typename bg::coordinate_type<PS>::type
- >::type rtype;
- typedef bg::srs::spheroid<rtype> stype;
- typedef bg::strategy::side::thomas<stype> strategy_type;
- typedef bg::strategy::side::geographic<bg::strategy::thomas, stype> strategy2_type;
- strategy_type strategy;
- strategy2_type strategy2;
- PS p1, p2;
- P p;
- bg::assign_values(p1, lon1, lat1);
- bg::assign_values(p2, lon2, lat2);
- bg::assign_values(p, lon, lat);
- int side = strategy.apply(p1, p2, p);
- int side2 = strategy2.apply(p1, p2, p);
- BOOST_CHECK_EQUAL(side, expected_side);
- BOOST_CHECK_EQUAL(side2, expected_side);
- }
- template <typename P1, typename P2>
- void test_all()
- {
- test_distance<P1, P2>(0, 90, 1, 80, 1116.825795); // polar
- test_distance<P1, P2>(0, -90, 1, -80, 1116.825795); // polar
- test_distance<P1, P2>(4, 52, 4, 52, 0.0); // no point difference
- test_distance<P1, P2>(4, 52, 3, 40, 1336.025365); // normal case
- test_side<P1, P2>(0, 0, 0, 1, 0, 2, 0);
- test_side<P1, P2>(0, 0, 0, 1, 0, -2, 0);
- test_side<P1, P2>(10, 0, 10, 1, 10, 2, 0);
- test_side<P1, P2>(10, 0, 10, -1, 10, 2, 0);
- test_side<P1, P2>(10, 0, 10, 1, 0, 2, 1); // left
- test_side<P1, P2>(10, 0, 10, -1, 0, 2, -1); // right
- test_side<P1, P2>(-10, -10, 10, 10, 10, 0, -1); // right
- test_side<P1, P2>(-10, -10, 10, 10, -10, 0, 1); // left
- test_side<P1, P2>(170, -10, -170, 10, -170, 0, -1); // right
- test_side<P1, P2>(170, -10, -170, 10, 170, 0, 1); // left
- }
- template <typename P>
- void test_all()
- {
- test_all<P, P>();
- }
- int test_main(int, char* [])
- {
- //test_all<float[2]>();
- //test_all<double[2]>();
- test_all<bg::model::point<int, 2, bg::cs::geographic<bg::degree> > >();
- test_all<bg::model::point<float, 2, bg::cs::geographic<bg::degree> > >();
- test_all<bg::model::point<double, 2, bg::cs::geographic<bg::degree> > >();
- #if defined(HAVE_TTMATH)
- test_all<bg::model::point<ttmath::Big<1,4>, 2, bg::cs::geographic<bg::degree> > >();
- test_all<bg::model::point<ttmath_big, 2, bg::cs::geographic<bg::degree> > >();
- #endif
- return 0;
- }
|