524 lines
18 KiB
C++
524 lines
18 KiB
C++
// Boost.Geometry (aka GGL, Generic Geometry Library)
|
|
// Unit Test
|
|
|
|
// Copyright (c) 2007-2014 Barend Gehrels, Amsterdam, the Netherlands.
|
|
// Copyright (c) 2008-2014 Bruno Lalande, Paris, France.
|
|
// Copyright (c) 2009-2014 Mateusz Loskot, London, UK.
|
|
// Copyright (c) 2013-2014 Adam Wulkiewicz, Lodz, Poland.
|
|
|
|
// This file was modified by Oracle on 2014, 2017.
|
|
// Modifications copyright (c) 2014-2017, Oracle and/or its affiliates.
|
|
|
|
// Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
|
|
// 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 <sstream>
|
|
|
|
#include <boost/mpl/if.hpp>
|
|
#include <boost/type_traits/is_integral.hpp>
|
|
#include <boost/type_traits/is_same.hpp>
|
|
#include <geometry_test_common.hpp>
|
|
|
|
#include <boost/geometry/algorithms/comparable_distance.hpp>
|
|
#include <boost/geometry/algorithms/make.hpp>
|
|
|
|
#include <boost/geometry/geometries/geometries.hpp>
|
|
#include <boost/geometry/geometries/point_xy.hpp>
|
|
#include <boost/geometry/io/wkt/read.hpp>
|
|
#include <boost/geometry/strategies/strategies.hpp>
|
|
|
|
template <typename P>
|
|
void test_distance_result()
|
|
{
|
|
typedef typename bg::default_distance_result<P, P>::type distance_type;
|
|
|
|
P p1 = bg::make<P>(0, 0);
|
|
P p2 = bg::make<P>(3, 0);
|
|
P p3 = bg::make<P>(0, 4);
|
|
|
|
distance_type dr12 = bg::comparable_distance(p1, p2);
|
|
distance_type dr13 = bg::comparable_distance(p1, p3);
|
|
distance_type dr23 = bg::comparable_distance(p2, p3);
|
|
|
|
BOOST_CHECK_CLOSE(dr12, 9.000, 0.001);
|
|
BOOST_CHECK_CLOSE(dr13, 16.000, 0.001);
|
|
BOOST_CHECK_CLOSE(dr23, 25.000, 0.001);
|
|
|
|
}
|
|
|
|
template <typename P>
|
|
void test_distance_point()
|
|
{
|
|
P p1;
|
|
bg::set<0>(p1, 1);
|
|
bg::set<1>(p1, 1);
|
|
|
|
P p2;
|
|
bg::set<0>(p2, 2);
|
|
bg::set<1>(p2, 2);
|
|
|
|
typename bg::coordinate_type<P>::type d = bg::comparable_distance(p1, p2);
|
|
BOOST_CHECK_CLOSE(d, 2.0, 0.001);
|
|
}
|
|
|
|
template <typename P>
|
|
void test_distance_segment()
|
|
{
|
|
typedef typename bg::coordinate_type<P>::type coordinate_type;
|
|
|
|
P s1 = bg::make<P>(2, 2);
|
|
P s2 = bg::make<P>(3, 3);
|
|
|
|
// Check points left, right, projected-left, projected-right, on segment
|
|
P p1 = bg::make<P>(0, 0);
|
|
P p2 = bg::make<P>(4, 4);
|
|
P p3 = bg::make<P>(2.4, 2.6);
|
|
P p4 = bg::make<P>(2.6, 2.4);
|
|
P p5 = bg::make<P>(2.5, 2.5);
|
|
|
|
bg::model::referring_segment<P const> const seg(s1, s2);
|
|
|
|
coordinate_type d1 = bg::comparable_distance(p1, seg); BOOST_CHECK_CLOSE(d1, 8.0, 0.001);
|
|
coordinate_type d2 = bg::comparable_distance(p2, seg); BOOST_CHECK_CLOSE(d2, 2.0, 0.001);
|
|
coordinate_type d3 = bg::comparable_distance(p3, seg); BOOST_CHECK_CLOSE(d3, 0.02, 0.001);
|
|
coordinate_type d4 = bg::comparable_distance(p4, seg); BOOST_CHECK_CLOSE(d4, 0.02, 0.001);
|
|
coordinate_type d5 = bg::comparable_distance(p5, seg); BOOST_CHECK_CLOSE(d5, 0.0, 0.001);
|
|
|
|
// Reverse case
|
|
coordinate_type dr1 = bg::comparable_distance(seg, p1); BOOST_CHECK_CLOSE(dr1, d1, 0.001);
|
|
coordinate_type dr2 = bg::comparable_distance(seg, p2); BOOST_CHECK_CLOSE(dr2, d2, 0.001);
|
|
}
|
|
|
|
template <typename P>
|
|
void test_distance_linestring()
|
|
{
|
|
bg::model::linestring<P> points;
|
|
points.push_back(bg::make<P>(1, 1));
|
|
points.push_back(bg::make<P>(3, 3));
|
|
|
|
P p = bg::make<P>(2, 1);
|
|
|
|
typename bg::coordinate_type<P>::type d = bg::comparable_distance(p, points);
|
|
BOOST_CHECK_CLOSE(d, 0.5, 0.001);
|
|
|
|
p = bg::make<P>(5, 5);
|
|
d = bg::comparable_distance(p, points);
|
|
BOOST_CHECK_CLOSE(d, 8.0, 0.001);
|
|
|
|
|
|
bg::model::linestring<P> line;
|
|
line.push_back(bg::make<P>(1,1));
|
|
line.push_back(bg::make<P>(2,2));
|
|
line.push_back(bg::make<P>(3,3));
|
|
|
|
p = bg::make<P>(5, 5);
|
|
|
|
d = bg::comparable_distance(p, line);
|
|
BOOST_CHECK_CLOSE(d, 8.0, 0.001);
|
|
|
|
// Reverse case
|
|
d = bg::comparable_distance(line, p);
|
|
BOOST_CHECK_CLOSE(d, 8.0, 0.001);
|
|
}
|
|
|
|
template <typename P>
|
|
void test_all()
|
|
{
|
|
test_distance_result<P>();
|
|
test_distance_point<P>();
|
|
test_distance_segment<P>();
|
|
test_distance_linestring<P>();
|
|
}
|
|
|
|
template <typename T>
|
|
void test_double_result_from_integer()
|
|
{
|
|
typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
|
|
|
|
point_type point;
|
|
|
|
// Check linestring
|
|
bg::model::linestring<point_type> linestring;
|
|
bg::read_wkt("POINT(2 2)", point);
|
|
bg::read_wkt("LINESTRING(4 1,1 4)", linestring);
|
|
|
|
double normal_distance = bg::distance(point, linestring);
|
|
double comparable_distance = bg::comparable_distance(point, linestring);
|
|
|
|
BOOST_CHECK_CLOSE(normal_distance, std::sqrt(0.5), 0.001);
|
|
BOOST_CHECK_CLOSE(comparable_distance, 0.5, 0.001);
|
|
|
|
// Check polygon
|
|
bg::model::polygon<point_type> polygon;
|
|
bg::read_wkt("POLYGON((0 0,1 9,8 1,0 0),(1 1,4 1,1 4,1 1))", polygon);
|
|
|
|
normal_distance = bg::distance(point, polygon);
|
|
comparable_distance = bg::comparable_distance(point, polygon);
|
|
|
|
BOOST_CHECK_CLOSE(normal_distance, std::sqrt(0.5), 0.001);
|
|
BOOST_CHECK_CLOSE(comparable_distance, 0.5, 0.001);
|
|
}
|
|
|
|
template <typename T>
|
|
struct test_variant_different_default_strategy
|
|
{
|
|
static inline void apply()
|
|
{
|
|
typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
|
|
typedef bg::model::segment<point_type> segment_type;
|
|
typedef bg::model::box<point_type> box_type;
|
|
typedef boost::variant<point_type, segment_type, box_type> variant_type;
|
|
|
|
point_type point;
|
|
bg::read_wkt("POINT(1 3)", point);
|
|
|
|
segment_type seg;
|
|
bg::read_wkt("LINESTRING(1 1,4 4)", seg);
|
|
|
|
box_type box;
|
|
bg::read_wkt("BOX(-1 -1,0 0)", box);
|
|
|
|
variant_type v1, v2;
|
|
|
|
BOOST_MPL_ASSERT((
|
|
boost::is_same
|
|
<
|
|
typename bg::comparable_distance_result
|
|
<
|
|
variant_type, variant_type, bg::default_strategy
|
|
>::type,
|
|
typename bg::comparable_distance_result
|
|
<
|
|
point_type, point_type, bg::default_strategy
|
|
>::type
|
|
>
|
|
));
|
|
|
|
// Default strategy
|
|
v1 = point;
|
|
v2 = point;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(point, point),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, point),
|
|
bg::comparable_distance(point, point),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
|
|
bg::comparable_distance(point, point),
|
|
0.0001);
|
|
v1 = point;
|
|
v2 = seg;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(point, seg),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, seg),
|
|
bg::comparable_distance(point, seg),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
|
|
bg::comparable_distance(point, seg), 0.0001);
|
|
v1 = point;
|
|
v2 = box;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(point, box),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, box),
|
|
bg::comparable_distance(point, box),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
|
|
bg::comparable_distance(point, box), 0.0001);
|
|
}
|
|
};
|
|
|
|
template <typename T, typename ExpectedResultType = double>
|
|
struct test_variant_same_default_strategy
|
|
{
|
|
static inline void apply()
|
|
{
|
|
typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
|
|
typedef bg::model::segment<point_type> segment_type;
|
|
typedef bg::model::linestring<point_type> linestring_type;
|
|
typedef boost::variant
|
|
<
|
|
point_type, segment_type, linestring_type
|
|
> variant_type;
|
|
|
|
point_type point;
|
|
bg::read_wkt("POINT(1 3)", point);
|
|
|
|
segment_type seg;
|
|
bg::read_wkt("LINESTRING(1 1,4 4)", seg);
|
|
|
|
linestring_type linestring;
|
|
bg::read_wkt("LINESTRING(-1 -1,-1 0,0 0,0 -1,-1 -1)", linestring);
|
|
|
|
variant_type v1, v2;
|
|
|
|
BOOST_MPL_ASSERT((
|
|
boost::is_same
|
|
<
|
|
typename bg::comparable_distance_result
|
|
<
|
|
variant_type, variant_type, bg::default_strategy
|
|
>::type,
|
|
ExpectedResultType
|
|
>
|
|
));
|
|
|
|
BOOST_MPL_ASSERT((
|
|
boost::is_same
|
|
<
|
|
typename bg::comparable_distance_result
|
|
<
|
|
point_type, point_type, bg::default_strategy
|
|
>::type,
|
|
ExpectedResultType
|
|
>
|
|
));
|
|
|
|
// Default strategy
|
|
v1 = point;
|
|
v2 = point;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(point, point),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, point),
|
|
bg::comparable_distance(point, point),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
|
|
bg::comparable_distance(point, point),
|
|
0.0001);
|
|
v1 = point;
|
|
v2 = seg;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(point, seg),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, seg),
|
|
bg::comparable_distance(point, seg),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
|
|
bg::comparable_distance(point, seg),
|
|
0.0001);
|
|
v1 = point;
|
|
v2 = linestring;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(point, linestring),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, linestring),
|
|
bg::comparable_distance(point, linestring),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(point, v2),
|
|
bg::comparable_distance(point, linestring),
|
|
0.0001);
|
|
}
|
|
};
|
|
|
|
template <typename T, typename ExpectedResultType = T>
|
|
struct test_variant_with_strategy
|
|
{
|
|
static inline void apply()
|
|
{
|
|
typedef bg::strategy::distance::projected_point<T> strategy_type;
|
|
|
|
typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
|
|
typedef bg::model::segment<point_type> segment_type;
|
|
typedef bg::model::linestring<point_type> linestring_type;
|
|
typedef bg::model::multi_linestring
|
|
<
|
|
linestring_type
|
|
> multi_linestring_type;
|
|
typedef boost::variant
|
|
<
|
|
segment_type, linestring_type, multi_linestring_type
|
|
> variant_type;
|
|
|
|
segment_type seg;
|
|
bg::read_wkt("LINESTRING(1 1,4 4)", seg);
|
|
|
|
linestring_type ls;
|
|
bg::read_wkt("LINESTRING(-1 -1,-1 0,0 0,0 -1,-1 -1)", ls);
|
|
|
|
multi_linestring_type mls;
|
|
bg::read_wkt("MULTILINESTRING((10 0,20 0),(30 0,40 0))", mls);
|
|
|
|
variant_type v1, v2;
|
|
|
|
strategy_type strategy;
|
|
|
|
BOOST_MPL_ASSERT((
|
|
boost::is_same
|
|
<
|
|
typename bg::comparable_distance_result
|
|
<
|
|
variant_type, variant_type, strategy_type
|
|
>::type,
|
|
ExpectedResultType
|
|
>
|
|
));
|
|
|
|
BOOST_MPL_ASSERT((
|
|
boost::is_same
|
|
<
|
|
typename bg::comparable_distance_result
|
|
<
|
|
segment_type, linestring_type, strategy_type
|
|
>::type,
|
|
ExpectedResultType
|
|
>
|
|
));
|
|
|
|
// Passed strategy
|
|
v1 = seg;
|
|
v2 = seg;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
|
|
bg::comparable_distance(seg, seg, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, seg, strategy),
|
|
bg::comparable_distance(seg, seg, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(seg, v2, strategy),
|
|
bg::comparable_distance(seg, seg, strategy),
|
|
0.0001);
|
|
v1 = seg;
|
|
v2 = ls;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
|
|
bg::comparable_distance(seg, ls, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, ls, strategy),
|
|
bg::comparable_distance(seg, ls, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(seg, v2, strategy),
|
|
bg::comparable_distance(seg, ls, strategy),
|
|
0.0001);
|
|
v1 = seg;
|
|
v2 = mls;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
|
|
bg::comparable_distance(seg, mls, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, mls, strategy),
|
|
bg::comparable_distance(seg, mls, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(seg, v2, strategy),
|
|
bg::comparable_distance(seg, mls, strategy),
|
|
0.0001);
|
|
v1 = ls;
|
|
v2 = mls;
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, v2, strategy),
|
|
bg::comparable_distance(ls, mls, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(v1, mls, strategy),
|
|
bg::comparable_distance(ls, mls, strategy),
|
|
0.0001);
|
|
BOOST_CHECK_CLOSE(bg::comparable_distance(ls, v2, strategy),
|
|
bg::comparable_distance(ls, mls, strategy),
|
|
0.0001);
|
|
}
|
|
};
|
|
|
|
template <typename T, bool IsIntergral = boost::is_integral<T>::value>
|
|
struct check_result
|
|
{
|
|
template <typename ExpectedResult>
|
|
static inline void apply(T const& value,
|
|
ExpectedResult const& expected_value)
|
|
{
|
|
BOOST_CHECK_EQUAL(value, expected_value);
|
|
}
|
|
};
|
|
|
|
template <typename T>
|
|
struct check_result<T, false>
|
|
{
|
|
template <typename ExpectedResult>
|
|
static inline void apply(T const& value,
|
|
ExpectedResult const& expected_value)
|
|
{
|
|
BOOST_CHECK_CLOSE(value, expected_value, 0.0001);
|
|
}
|
|
};
|
|
|
|
template <typename T>
|
|
struct test_variant_boxes
|
|
{
|
|
static inline void apply()
|
|
{
|
|
typedef bg::model::point<T, 2, bg::cs::cartesian> point_type;
|
|
typedef bg::model::box<point_type> box_type;
|
|
typedef boost::variant<box_type> variant_type;
|
|
|
|
box_type box1, box2;
|
|
bg::read_wkt("BOX(-1 -1,0 0)", box1);
|
|
bg::read_wkt("BOX(1 1,2 2)", box2);
|
|
|
|
variant_type v1 = box1, v2 = box2;
|
|
|
|
typedef typename boost::mpl::if_c
|
|
<
|
|
boost::is_float<T>::value,
|
|
double,
|
|
typename bg::util::detail::default_integral::type
|
|
>::type expected_result_type;
|
|
|
|
BOOST_MPL_ASSERT((
|
|
boost::is_same
|
|
<
|
|
typename bg::comparable_distance_result
|
|
<
|
|
variant_type, variant_type, bg::default_strategy
|
|
>::type,
|
|
expected_result_type
|
|
>
|
|
));
|
|
|
|
// Default strategy
|
|
check_result<T>::apply(bg::comparable_distance(v1, v2),
|
|
bg::comparable_distance(box1, box2));
|
|
check_result<T>::apply(bg::comparable_distance(v1, box2),
|
|
bg::comparable_distance(box1, box2));
|
|
check_result<T>::apply(bg::comparable_distance(box1, v2),
|
|
bg::comparable_distance(box1, box2));
|
|
}
|
|
};
|
|
|
|
|
|
int test_main(int, char* [])
|
|
{
|
|
test_double_result_from_integer<int>();
|
|
test_double_result_from_integer<boost::long_long_type>();
|
|
|
|
test_all<bg::model::d2::point_xy<float> >();
|
|
test_all<bg::model::d2::point_xy<double> >();
|
|
|
|
#ifdef HAVE_TTMATH
|
|
test_all<bg::model::d2::point_xy<ttmath_big> >();
|
|
#endif
|
|
|
|
// test variant support
|
|
test_variant_different_default_strategy<double>::apply();
|
|
|
|
test_variant_same_default_strategy<double>::apply();
|
|
test_variant_same_default_strategy<int>::apply();
|
|
test_variant_same_default_strategy<long>::apply();
|
|
|
|
test_variant_with_strategy<double>::apply();
|
|
test_variant_with_strategy<float>::apply();
|
|
test_variant_with_strategy<long double>::apply();
|
|
test_variant_with_strategy<int, double>::apply();
|
|
#ifdef HAVE_TTMATH
|
|
test_variant_with_strategy<ttmath_big>::apply();
|
|
#endif
|
|
|
|
test_variant_boxes<double>::apply();
|
|
test_variant_boxes<int>::apply();
|
|
test_variant_boxes<long>::apply();
|
|
|
|
return 0;
|
|
}
|