1// Boost.Geometry Index
2// Unit Test
3
4// Copyright (c) 2011-2013 Adam Wulkiewicz, Lodz, Poland.
5
6// Use, modification and distribution is subject to the Boost Software License,
7// Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
8// http://www.boost.org/LICENSE_1_0.txt)
9
10#include <geometry_index_test_common.hpp>
11
12#include <boost/geometry/index/detail/algorithms/segment_intersection.hpp>
13
14//#include <boost/geometry/io/wkt/read.hpp>
15
16template <typename Box, typename Point, typename RelativeDistance>
17void test_segment_intersection(Box const& box, Point const& p0, Point const& p1,
18 bool expected_result,
19 RelativeDistance expected_rel_dist)
20{
21 RelativeDistance rel_dist;
22 bool value = bgi::detail::segment_intersection(box, p0, p1, rel_dist);
23 BOOST_CHECK(value == expected_result);
24 if ( value && expected_result )
25 BOOST_CHECK_CLOSE(rel_dist, expected_rel_dist, 0.0001);
26}
27
28template <typename Box, typename Point, typename RelativeDistance>
29void test_geometry(std::string const& wkt_g, std::string const& wkt_p0, std::string const& wkt_p1,
30 bool expected_result,
31 RelativeDistance expected_rel_dist)
32{
33 Box box;
34 bg::read_wkt(wkt_g, box);
35 Point p0, p1;
36 bg::read_wkt(wkt_p0, p0);
37 bg::read_wkt(wkt_p1, p1);
38 test_segment_intersection(box, p0, p1, expected_result, expected_rel_dist);
39}
40
41#include <boost/geometry/geometries/point_xy.hpp>
42#include <boost/geometry/geometries/point.hpp>
43#include <boost/geometry/geometries/box.hpp>
44
45void test_large_integers()
46{
47 typedef bg::model::point<int, 2, bg::cs::cartesian> int_point_type;
48 typedef bg::model::point<double, 2, bg::cs::cartesian> double_point_type;
49
50 bg::model::box<int_point_type> int_box;
51 bg::model::box<double_point_type> double_box;
52 int_point_type int_p0, int_p1;
53 double_point_type double_p0, double_p1;
54
55 std::string const str_box = "POLYGON((1536119 192000, 1872000 528000))";
56 std::string const str_p0 = "POINT(1535000 191000)";
57 std::string const str_p1 = "POINT(1873000 529000)";
58 bg::read_wkt(wkt: str_box, geometry&: int_box);
59 bg::read_wkt(wkt: str_box, geometry&: double_box);
60 bg::read_wkt(wkt: str_p0, geometry&: int_p0);
61 bg::read_wkt(wkt: str_p1, geometry&: int_p1);
62 bg::read_wkt(wkt: str_p0, geometry&: double_p0);
63 bg::read_wkt(wkt: str_p1, geometry&: double_p1);
64
65 float int_value;
66 bool int_result = bgi::detail::segment_intersection(b: int_box, p0: int_p0, p1: int_p1, relative_distance&: int_value);
67 double double_value;
68 bool double_result = bgi::detail::segment_intersection(b: double_box, p0: double_p0, p1: double_p1, relative_distance&: double_value);
69 BOOST_CHECK(int_result == double_result);
70 if ( int_result && double_result )
71 BOOST_CHECK_CLOSE(int_value, double_value, 0.0001);
72}
73
74int test_main(int, char* [])
75{
76 typedef bg::model::point<int, 2, bg::cs::cartesian> P2ic;
77 typedef bg::model::point<float, 2, bg::cs::cartesian> P2fc;
78 typedef bg::model::point<double, 2, bg::cs::cartesian> P2dc;
79
80 typedef bg::model::point<int, 3, bg::cs::cartesian> P3ic;
81 typedef bg::model::point<float, 3, bg::cs::cartesian> P3fc;
82 typedef bg::model::point<double, 3, bg::cs::cartesian> P3dc;
83
84 test_geometry<bg::model::box<P2ic>, P2ic>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0 0)", wkt_p1: "POINT(2 5)", expected_result: true, expected_rel_dist: 1.0f/5);
85 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0 0)", wkt_p1: "POINT(2 5)", expected_result: true, expected_rel_dist: 1.0f/5);
86 test_geometry<bg::model::box<P2dc>, P2dc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0 0)", wkt_p1: "POINT(2 5)", expected_result: true, expected_rel_dist: 1.0/5);
87 test_geometry<bg::model::box<P3ic>, P3ic>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(0 0 0)", wkt_p1: "POINT(2 5 7)", expected_result: true, expected_rel_dist: 2.0f/7);
88 test_geometry<bg::model::box<P3fc>, P3fc>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(0 0 0)", wkt_p1: "POINT(2 5 7)", expected_result: true, expected_rel_dist: 2.0f/7);
89 test_geometry<bg::model::box<P3dc>, P3dc>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(0 0 0)", wkt_p1: "POINT(2 5 7)", expected_result: true, expected_rel_dist: 2.0/7);
90
91 test_geometry<bg::model::box<P2ic>, P2ic>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(3 4)", wkt_p1: "POINT(0 0)", expected_result: true, expected_rel_dist: 1.0f/3);
92 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(3 4)", wkt_p1: "POINT(0 2)", expected_result: true, expected_rel_dist: 1.0f/3);
93 test_geometry<bg::model::box<P2dc>, P2dc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(3 4)", wkt_p1: "POINT(0 2)", expected_result: true, expected_rel_dist: 1.0/3);
94 test_geometry<bg::model::box<P3ic>, P3ic>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(3 5 6)", wkt_p1: "POINT(0 3 3)", expected_result: true, expected_rel_dist: 1.0f/2);
95 test_geometry<bg::model::box<P3fc>, P3fc>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(3 5 6)", wkt_p1: "POINT(0 3 3)", expected_result: true, expected_rel_dist: 1.0f/2);
96 test_geometry<bg::model::box<P3dc>, P3dc>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(3 5 6)", wkt_p1: "POINT(0 3 3)", expected_result: true, expected_rel_dist: 1.0/2);
97
98 test_geometry<bg::model::box<P2ic>, P2ic>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(1 0)", wkt_p1: "POINT(1 5)", expected_result: true, expected_rel_dist: 1.0f/5);
99 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(1 5)", wkt_p1: "POINT(1 0)", expected_result: true, expected_rel_dist: 1.0f/5);
100 test_geometry<bg::model::box<P2dc>, P2dc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(1 0)", wkt_p1: "POINT(1 5)", expected_result: true, expected_rel_dist: 1.0/5);
101 test_geometry<bg::model::box<P3ic>, P3ic>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(1 3 0)", wkt_p1: "POINT(1 3 7)", expected_result: true, expected_rel_dist: 2.0f/7);
102 test_geometry<bg::model::box<P3fc>, P3fc>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(1 3 7)", wkt_p1: "POINT(1 3 0)", expected_result: true, expected_rel_dist: 1.0f/7);
103 test_geometry<bg::model::box<P3dc>, P3dc>(wkt_g: "POLYGON((0 1 2,2 4 6))", wkt_p0: "POINT(1 3 0)", wkt_p1: "POINT(1 3 7)", expected_result: true, expected_rel_dist: 2.0/7);
104
105 test_geometry<bg::model::box<P2ic>, P2ic>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0 0)", wkt_p1: "POINT(0 5)", expected_result: true, expected_rel_dist: 0.2f);
106 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0 5)", wkt_p1: "POINT(0 0)", expected_result: true, expected_rel_dist: 0.2f);
107 test_geometry<bg::model::box<P2dc>, P2dc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0 0)", wkt_p1: "POINT(0 5)", expected_result: true, expected_rel_dist: 0.2);
108
109 test_geometry<bg::model::box<P2ic>, P2ic>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(3 0)", wkt_p1: "POINT(3 5)", expected_result: false, expected_rel_dist: 0.0f);
110 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(3 5)", wkt_p1: "POINT(3 0)", expected_result: false, expected_rel_dist: 0.0f);
111 test_geometry<bg::model::box<P2dc>, P2dc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(3 0)", wkt_p1: "POINT(3 5)", expected_result: false, expected_rel_dist: 0.0);
112
113 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(1 0)", wkt_p1: "POINT(1 1)", expected_result: true, expected_rel_dist: 1.0f);
114 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(1 4)", wkt_p1: "POINT(1 5)", expected_result: true, expected_rel_dist: 0.0f);
115
116 test_geometry<bg::model::box<P2fc>, P2fc>(wkt_g: "POLYGON((0 1,2 4))", wkt_p0: "POINT(0.5 2)", wkt_p1: "POINT(1.5 3)", expected_result: true, expected_rel_dist: 0.0f);
117
118#ifdef HAVE_TTMATH
119 typedef bg::model::point<ttmath_big, 2, bg::cs::cartesian> P2ttmc;
120 typedef bg::model::point<ttmath_big, 3, bg::cs::cartesian> P3ttmc;
121
122 test_geometry<bg::model::box<P2ttmc>, P2ttmc>("POLYGON((0 1,2 4))", "POINT(0 0)", "POINT(2 5)", true, 1.0f/5);
123 test_geometry<bg::model::box<P3ttmc>, P3ttmc>("POLYGON((0 1 2,2 4 6))", "POINT(0 0 0)", "POINT(2 5 7)", true, 2.0f/7);
124#endif
125
126 test_large_integers();
127
128 return 0;
129}
130

source code of boost/libs/geometry/index/test/algorithms/segment_intersection.cpp