boost_mirror/geometry
1
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2025-05-09 23:24:02 +00:00
Code Issues Packages Projects Releases Wiki Activity
geometry/test/strategies/segment_intersection_sph.cpp
Adam Wulkiewicz bc77e1255e [test] Adapt tests to new interfaces of various utilities in the library. 
default_strategy services, partition, get_turns, strategies.
2017-01-28 04:09:49 +01:00

410 lines
16 KiB
C++
Raw Blame History

// Boost.Geometry
// Unit Test
// Copyright (c) 2016-2017, Oracle and/or its affiliates.
// Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
// 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/geometry/algorithms/equals.hpp>
#include <boost/geometry/geometries/point.hpp>
#include <boost/geometry/geometries/segment.hpp>
#include <boost/geometry/io/wkt/read.hpp>
#include <boost/geometry/io/wkt/write.hpp>
#include <boost/geometry/policies/relate/direction.hpp>
#include <boost/geometry/policies/relate/intersection_points.hpp>
#include <boost/geometry/policies/relate/tupled.hpp>
#include <boost/geometry/strategies/spherical/intersection.hpp>
template <typename P1, typename P2, typename T>
bool equals_relaxed(P1 const& p1, P2 const& p2, T const& eps_scale)
{
typedef typename bg::select_coordinate_type<P1, P2>::type calc_t;
calc_t c1 = 1;
calc_t p10 = bg::get<0>(p1);
calc_t p11 = bg::get<1>(p1);
calc_t p20 = bg::get<0>(p2);
calc_t p21 = bg::get<1>(p2);
calc_t relaxed_eps = std::numeric_limits<calc_t>::epsilon()
* bg::math::detail::greatest(c1, bg::math::abs(p10), bg::math::abs(p11), bg::math::abs(p20), bg::math::abs(p21))
* eps_scale;
calc_t diff0 = p10 - p20;
// needed e.g. for -179.999999 - 180.0
if (diff0 < -180)
diff0 += 360;
return bg::math::abs(diff0) <= relaxed_eps
&& bg::math::abs(p11 - p21) <= relaxed_eps;
}
template <typename S1, typename S2, typename P>
void test_spherical_strategy_one(S1 const& s1, S2 const& s2, char m, std::size_t expected_count, P const& ip0 = P(), P const& ip1 = P())
{
typedef typename bg::coordinate_type<P>::type coord_t;
typedef bg::policies::relate::segments_tupled
<
bg::policies::relate::segments_intersection_points
<
bg::segment_intersection_points<P, bg::segment_ratio<coord_t> >
>,
bg::policies::relate::segments_direction
> policy_t;
typedef bg::strategy::intersection::relate_spherical_segments<> strategy_t;
typedef typename policy_t::return_type return_type;
// NOTE: robust policy is currently ignored
return_type res = strategy_t::apply(s1, s2, policy_t(), 0);
size_t const res_count = boost::get<0>(res).count;
char const res_method = boost::get<1>(res).how;
BOOST_CHECK_MESSAGE(res_method == m,
"IP method: " << res_method << " different than expected: " << m
<< " for " << bg::wkt(s1) << " and " << bg::wkt(s2));
BOOST_CHECK_MESSAGE(res_count == expected_count,
"IP count: " << res_count << " different than expected: " << expected_count
<< " for " << bg::wkt(s1) << " and " << bg::wkt(s2));
// The EPS is scaled because during the conversion various angles may be not converted
// to cartesian 3d the same way which results in a different intersection point
// For more information read the info in spherical intersection strategy file.
int eps_scale = res_method != 'i' ? 1 : 1000;
if (res_count > 0 && expected_count > 0)
{
P const& res_i0 = boost::get<0>(res).intersections[0];
BOOST_CHECK_MESSAGE(equals_relaxed(res_i0, ip0, eps_scale),
"IP0: " << std::setprecision(16) << bg::wkt(res_i0) << " different than expected: " << bg::wkt(ip0)
<< " for " << bg::wkt(s1) << " and " << bg::wkt(s2));
}
if (res_count > 1 && expected_count > 1)
{
P const& res_i1 = boost::get<0>(res).intersections[1];
BOOST_CHECK_MESSAGE(equals_relaxed(res_i1, ip1, eps_scale),
"IP1: " << std::setprecision(16) << bg::wkt(res_i1) << " different than expected: " << bg::wkt(ip1)
<< " for " << bg::wkt(s1) << " and " << bg::wkt(s2));
}
}
template <typename T>
T translated(T v, double t)
{
v += T(t);
if (v > 180)
v -= 360;
return v;
}
template <typename S1, typename S2, typename P>
void test_spherical_strategy(S1 const& s1, S2 const& s2, char m, std::size_t expected_count, P const& ip0 = P(), P const& ip1 = P())
{
S1 s1t = s1;
S2 s2t = s2;
P ip0t = ip0;
P ip1t = ip1;
double t = 0;
for (int i = 0; i < 4; ++i, t += 90)
{
bg::set<0, 0>(s1t, translated(bg::get<0, 0>(s1), t));
bg::set<1, 0>(s1t, translated(bg::get<1, 0>(s1), t));
bg::set<0, 0>(s2t, translated(bg::get<0, 0>(s2), t));
bg::set<1, 0>(s2t, translated(bg::get<1, 0>(s2), t));
if (expected_count > 0)
bg::set<0>(ip0t, translated(bg::get<0>(ip0), t));
if (expected_count > 1)
bg::set<0>(ip1t, translated(bg::get<0>(ip1), t));
test_spherical_strategy_one(s1t, s2t, m, expected_count, ip0t, ip1t);
}
}
template <typename S1, typename S2, typename P>
void test_spherical_strategy(std::string const& s1_wkt, std::string const& s2_wkt,
char m, std::size_t expected_count,
std::string const& ip0_wkt = "", std::string const& ip1_wkt = "")
{
S1 s1;
S2 s2;
P ip0, ip1;
bg::read_wkt(s1_wkt, s1);
bg::read_wkt(s2_wkt, s2);
if (! ip0_wkt.empty())
bg::read_wkt(ip0_wkt, ip0);
if (!ip1_wkt.empty())
bg::read_wkt(ip1_wkt, ip1);
test_spherical_strategy(s1, s2, m, expected_count, ip0, ip1);
}
template <typename S, typename P>
void test_spherical_strategy(std::string const& s1_wkt, std::string const& s2_wkt,
char m, std::size_t expected_count,
std::string const& ip0_wkt = "", std::string const& ip1_wkt = "")
{
test_spherical_strategy<S, S, P>(s1_wkt, s2_wkt, m, expected_count, ip0_wkt, ip1_wkt);
}
template <typename T>
void test_spherical()
{
typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > point_t;
typedef bg::model::segment<point_t> segment_t;
// crossing X
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-45 -45, 45 45)", "SEGMENT(-45 45, 45 -45)", 'i', 1, "POINT(0 0)");
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-45 -45, 45 45)", "SEGMENT(45 -45, -45 45)", 'i', 1, "POINT(0 0)");
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(45 45, -45 -45)", "SEGMENT(-45 45, 45 -45)", 'i', 1, "POINT(0 0)");
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(45 45, -45 -45)", "SEGMENT(45 -45, -45 45)", 'i', 1, "POINT(0 0)");
// crossing X
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-1 1, 1 -1)", 'i', 1, "POINT(0 0)");
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(1 -1, -1 1)", 'i', 1, "POINT(0 0)");
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(1 1, -1 -1)", "SEGMENT(-1 1, 1 -1)", 'i', 1, "POINT(0 0)");
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(1 1, -1 -1)", "SEGMENT(1 -1, -1 1)", 'i', 1, "POINT(0 0)");
// equal
// //
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-45 -45, 45 45)", "SEGMENT(-45 -45, 45 45)", 'e', 2, "POINT(-45 -45)", "POINT(45 45)");
// //
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-45 -45, 45 45)", "SEGMENT(45 45, -45 -45)", 'e', 2, "POINT(-45 -45)", "POINT(45 45)");
// starting outside s1
// /
// |
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-2 -2, -1 -1)", 'a', 1, "POINT(-1 -1)");
// /
// /|
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-2 -2, 0 0)", 'm', 1, "POINT(0 0)");
// /|
// / |
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-2 -2, 1 1)", 't', 1, "POINT(1 1)");
// |/
// /|
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-2 -2, 2 2)", 'i', 1, "POINT(0 0)");
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-2 0, -1 0)", 'a', 1, "POINT(-1 0)");
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-2 0, 0 0)", 'c', 2, "POINT(-1 0)", "POINT(0 0)");
// ------
// ---------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-2 0, 1 0)", 'c', 2, "POINT(-1 0)", "POINT(1 0)");
// ------
// ------------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-2 0, 2 0)", 'c', 2, "POINT(-1 0)", "POINT(1 0)");
// starting at s1
// /
// //
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-1 -1, 0 0)", 'c', 2, "POINT(-1 -1)", "POINT(0 0)");
// //
// //
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-1 -1, 1 1)", 'e', 2, "POINT(-1 -1)", "POINT(1 1)");
// | /
// |/
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-1 -1, 2 2)", 'f', 1, "POINT(-1 -1)");
// ------
// ---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-1 0, 0 0)", 'c', 2, "POINT(-1 0)", "POINT(0 0)");
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-1 0, 1 0)", 'e', 2, "POINT(-1 0)", "POINT(1 0)");
// ------
// ---------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-1 0, 2 0)", 'c', 2, "POINT(-1 0)", "POINT(1 0)");
// starting inside
// //
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(0 0, 1 1)", 'c', 2, "POINT(0 0)", "POINT(1 1)");
// |/
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(0 0, 2 2)", 's', 1, "POINT(0 0)");
// ------
// ---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(0 0, 1 0)", 'c', 2, "POINT(0 0)", "POINT(1 0)");
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(0 0, 2 0)", 'c', 2, "POINT(0 0)", "POINT(1 0)");
// starting at p2
// |
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(1 1, 2 2)", 'a', 1, "POINT(1 1)");
// ------
// ---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(1 0, 2 0)", 'a', 1, "POINT(1 0)");
// disjoint, crossing
// /
// |
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-3 -3, -2 -2)", 'd', 0);
// |
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(2 2, 3 3)", 'd', 0);
// disjoint, collinear
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-3 0, -2 0)", 'd', 0);
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(2 0, 3 0)", 'd', 0);
// degenerated
// /
// *
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-2 -2, -2 -2)", 'd', 0);
// /
// *
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(-1 -1, -1 -1)", '0', 1, "POINT(-1 -1)");
// /
// *
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(0 0, 0 0)", '0', 1, "POINT(0 0)");
// *
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(1 1, 1 1)", '0', 1, "POINT(1 1)");
// *
// /
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 -1, 1 1)", "SEGMENT(2 2, 2 2)", 'd', 0);
// similar to above, collinear
// * ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-2 0, -2 0)", 'd', 0);
// *------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(-1 0, -1 0)", '0', 1, "POINT(-1 0)");
// ---*---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(0 0, 0 0)", '0', 1, "POINT(0 0)");
// ------*
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(1 0, 1 0)", '0', 1, "POINT(1 0)");
// ------ *
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 0, 1 0)", "SEGMENT(2 0, 2 0)", 'd', 0);
// Northern hemisphere
// --- ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-3 50, -2 50)", 'd', 0);
// ------
// ---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-2 50, -1 50)", 'a', 1, "POINT(-1 50)");
// \/
// /\ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-2 50, 0 50)", 'i', 1, "POINT(-0.5 50.0032229484023)");
// ________
// / _____\ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-2 50, 1 50)", 't', 1, "POINT(1 50)");
// _________
// / _____ \ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-2 50, 2 50)", 'd', 0);
// ______
// /___ \ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-1 50, 0 50)", 'f', 1, "POINT(-1 50)");
// ------
// ------
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-1 50, 1 50)", 'e', 2, "POINT(-1 50)", "POINT(1 50)");
// ________
// /_____ \ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(-1 50, 2 50)", 'f', 1, "POINT(-1 50)");
// ______
// / ___\ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(0 50, 1 50)", 't', 1, "POINT(1 50)");
// \/
// /\ (avoid multi-line comment)
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(0 50, 2 50)", 'i', 1, "POINT(0.5 50.0032229484023)");
// ------
// ---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(1 50, 2 50)", 'a', 1, "POINT(1 50)");
// ------ ---
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(-1 50, 1 50)", "SEGMENT(2 50, 3 50)", 'd', 0);
// ___|
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(0 0, 1 0)", "SEGMENT(1 0, 1 1)", 'a', 1, "POINT(1 0)");
// ___|
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(1 0, 1 1)", "SEGMENT(0 0, 1 0)", 'a', 1, "POINT(1 0)");
// |/
// /|
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(10 -1, 20 1)", "SEGMENT(12.5 -1, 12.5 1)", 'i', 1, "POINT(12.5 -0.50051443471392)");
// |/
// /|
test_spherical_strategy<segment_t, point_t>(
"SEGMENT(10 -1, 20 1)", "SEGMENT(17.5 -1, 17.5 1)", 'i', 1, "POINT(17.5 0.50051443471392)");
}
int test_main(int, char* [])
{
//test_spherical<float>();
test_spherical<double>();
return 0;
}
Reference in New Issue View Git Blame Copy Permalink