]>
Commit | Line | Data |
---|---|---|
7c673cae FG |
1 | // Boost.Geometry (aka GGL, Generic Geometry Library) |
2 | // Unit Test | |
3 | ||
4 | // Copyright (c) 2010-2012 Barend Gehrels, Amsterdam, the Netherlands. | |
5 | ||
6 | // This file was modified by Oracle on 2014, 2016. | |
7 | // Modifications copyright (c) 2014-2016 Oracle and/or its affiliates. | |
8 | // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle | |
9 | ||
10 | // Use, modification and distribution is subject to the Boost Software License, | |
11 | // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at | |
12 | // http://www.boost.org/LICENSE_1_0.txt) | |
13 | ||
14 | ||
15 | #include <strategies/test_within.hpp> | |
16 | ||
17 | ||
18 | template <typename Point> | |
19 | void test_cartesian() | |
20 | { | |
21 | typedef bg::model::polygon<Point> polygon; | |
22 | ||
23 | std::string const box = "POLYGON((0 0,0 2,2 2,2 0,0 0))"; | |
24 | std::string const triangle = "POLYGON((0 0,0 4,6 0,0 0))"; | |
25 | std::string const with_hole = "POLYGON((0 0,0 3,3 3,3 0,0 0),(1 1,2 1,2 2,1 2,1 1))"; | |
26 | ||
27 | bg::strategy::within::winding<Point> s; | |
28 | ||
29 | ||
30 | test_geometry<Point, polygon>("b1", "POINT(1 1)", box, s, true); | |
31 | test_geometry<Point, polygon>("b2", "POINT(3 3)", box, s, false); | |
32 | ||
33 | // Test ALL corners (officialy false but some strategies might answer true) | |
34 | test_geometry<Point, polygon>("b3a", "POINT(0 0)", box, s, false); | |
35 | test_geometry<Point, polygon>("b3b", "POINT(0 2)", box, s, false); | |
36 | test_geometry<Point, polygon>("b3c", "POINT(2 2)", box, s, false); | |
37 | test_geometry<Point, polygon>("b3d", "POINT(2 0)", box, s, false); | |
38 | ||
39 | // Test ALL sides (officialy false but some strategies might answer true) | |
40 | test_geometry<Point, polygon>("b4a", "POINT(0 1)", box, s, false); | |
41 | test_geometry<Point, polygon>("b4b", "POINT(1 2)", box, s, false); | |
42 | test_geometry<Point, polygon>("b4c", "POINT(2 1)", box, s, false); | |
43 | test_geometry<Point, polygon>("b4d", "POINT(1 0)", box, s, false); | |
44 | ||
45 | ||
46 | test_geometry<Point, polygon>("t1", "POINT(1 1)", triangle, s, true); | |
47 | test_geometry<Point, polygon>("t2", "POINT(3 3)", triangle, s, false); | |
48 | ||
49 | test_geometry<Point, polygon>("t3a", "POINT(0 0)", triangle, s, false); | |
50 | test_geometry<Point, polygon>("t3b", "POINT(0 4)", triangle, s, false); | |
51 | test_geometry<Point, polygon>("t3c", "POINT(5 0)", triangle, s, false); | |
52 | ||
53 | test_geometry<Point, polygon>("t4a", "POINT(0 2)", triangle, s, false); | |
54 | test_geometry<Point, polygon>("t4b", "POINT(3 2)", triangle, s, false); | |
55 | test_geometry<Point, polygon>("t4c", "POINT(2 0)", triangle, s, false); | |
56 | ||
57 | ||
58 | test_geometry<Point, polygon>("h1", "POINT(0.5 0.5)", with_hole, s, true); | |
59 | test_geometry<Point, polygon>("h2a", "POINT(1.5 1.5)", with_hole, s, false); | |
60 | test_geometry<Point, polygon>("h2b", "POINT(5 5)", with_hole, s, false); | |
61 | ||
62 | test_geometry<Point, polygon>("h3a", "POINT(1 1)", with_hole, s, false); | |
63 | test_geometry<Point, polygon>("h3b", "POINT(2 2)", with_hole, s, false); | |
64 | test_geometry<Point, polygon>("h3c", "POINT(0 0)", with_hole, s, false); | |
65 | ||
66 | test_geometry<Point, polygon>("h4a", "POINT(1 1.5)", with_hole, s, false); | |
67 | test_geometry<Point, polygon>("h4b", "POINT(1.5 2)", with_hole, s, false); | |
68 | ||
69 | // Lying ON (one of the sides of) interior ring | |
70 | test_geometry<Point, polygon>("#77-1", "POINT(6 3.5)", | |
71 | "POLYGON((5 3,5 4,4 4,4 5,3 5,3 6,5 6,5 5,7 5,7 6,8 6,8 5,9 5,9 2,8 2,8 1,7 1,7 2,5 2,5 3),(6 3,8 3,8 4,6 4,6 3))", | |
72 | s, false); | |
73 | } | |
74 | ||
75 | template <typename T> | |
76 | void test_spherical() | |
77 | { | |
78 | typedef bg::model::point<T, 2, bg::cs::spherical_equatorial<bg::degree> > point; | |
79 | typedef bg::model::polygon<point> polygon; | |
80 | ||
81 | bg::strategy::within::winding<point> s; | |
82 | ||
83 | ||
84 | // Ticket #9354 | |
85 | test_geometry<point, polygon>( | |
86 | "#9354", | |
87 | "POINT(-78.1239 25.9556)", | |
88 | "POLYGON((-97.08466667 25.95683333, -97.13683333 25.954, -97.1 26, -97.08466667 25.95683333))", | |
89 | s, | |
90 | false); | |
91 | ||
92 | test_geometry<point, polygon>( | |
93 | "sph1N", | |
94 | "POINT(0 10.001)", | |
95 | "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))", | |
96 | s, | |
97 | bg::strategy::side::spherical_side_formula<>::apply( | |
98 | point(-10, 10), | |
99 | point(10, 10), | |
100 | point(0, (T)10.001)) == -1 // right side | |
101 | /*true*/); | |
102 | test_geometry<point, polygon>( | |
103 | "sph1S", | |
104 | "POINT(0 -10.001)", | |
105 | "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))", | |
106 | s, | |
107 | bg::strategy::side::spherical_side_formula<>::apply( | |
108 | point(10, -10), | |
109 | point(-10, -10), | |
110 | point(0, (T)-10.001)) == -1 // right side | |
111 | /*true*/); | |
112 | ||
113 | test_geometry<point, polygon>( | |
114 | "sph2S", | |
115 | "POINT(0 10.001)", | |
116 | "POLYGON((-10 20, 10 20, 10 10, -10 10, -10 20))", | |
117 | s, | |
118 | bg::strategy::side::spherical_side_formula<>::apply( | |
119 | point(10, 10), | |
120 | point(-10, 10), | |
121 | point(0, (T)10.001)) == -1 // right side | |
122 | /*false*/); | |
123 | ||
124 | test_geometry<point, polygon>( | |
125 | "sph3N", | |
126 | "POINT(0 10)", | |
127 | "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))", | |
128 | s, | |
129 | bg::strategy::side::spherical_side_formula<>::apply( | |
130 | point(-10, 10), | |
131 | point(10, 10), | |
132 | point(0, (T)10.001)) == -1 // right side | |
133 | /*true*/); | |
134 | test_geometry<point, polygon>( | |
135 | "sph3S", | |
136 | "POINT(0 -10)", | |
137 | "POLYGON((-10 10, 10 10, 10 -10, -10 -10, -10 10))", | |
138 | s, | |
139 | bg::strategy::side::spherical_side_formula<>::apply( | |
140 | point(10, -10), | |
141 | point(-10, -10), | |
142 | point(0, (T)-10.001)) == -1 // right side | |
143 | /*true*/); | |
144 | ||
145 | test_geometry<point, polygon>( | |
146 | "sphEq1", | |
147 | "POINT(179 10)", | |
148 | "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))", | |
149 | s, | |
150 | true, | |
151 | false); | |
152 | test_geometry<point, polygon>( | |
153 | "sphEq2", | |
154 | "POINT(179 10)", | |
155 | "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))", | |
156 | s, | |
157 | false, | |
158 | false); | |
159 | test_geometry<point, polygon>( | |
160 | "sphEq3", | |
161 | "POINT(-179 10)", | |
162 | "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))", | |
163 | s, | |
164 | true, | |
165 | false); | |
166 | test_geometry<point, polygon>( | |
167 | "sphEq4", | |
168 | "POINT(-179 10)", | |
169 | "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))", | |
170 | s, | |
171 | false, | |
172 | false); | |
173 | ||
174 | test_geometry<point, polygon>( | |
175 | "sphEq5", | |
176 | "POINT(169 10)", | |
177 | "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))", | |
178 | s, | |
179 | false, | |
180 | false); | |
181 | test_geometry<point, polygon>( | |
182 | "sphEq6", | |
183 | "POINT(-169 10)", | |
184 | "POLYGON((170 20, -170 20, -170 10, 170 10, 170 20))", | |
185 | s, | |
186 | false, | |
187 | false); | |
188 | test_geometry<point, polygon>( | |
189 | "sphEq7", | |
190 | "POINT(169 10)", | |
191 | "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))", | |
192 | s, | |
193 | false, | |
194 | false); | |
195 | test_geometry<point, polygon>( | |
196 | "sphEq8", | |
197 | "POINT(-169 10)", | |
198 | "POLYGON((170 10, -170 10, -170 0, 170 0, 170 10))", | |
199 | s, | |
200 | false, | |
201 | false); | |
202 | } | |
203 | ||
204 | int test_main(int, char* []) | |
205 | { | |
206 | test_cartesian<bg::model::point<float, 2, bg::cs::cartesian> >(); | |
207 | test_cartesian<bg::model::point<double, 2, bg::cs::cartesian> >(); | |
208 | ||
209 | test_spherical<float>(); | |
210 | test_spherical<double>(); | |
211 | ||
7c673cae FG |
212 | return 0; |
213 | } |