154 lines
5.0 KiB
C++
154 lines
5.0 KiB
C++
// 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;
|
|
}
|