1 // Boost.Geometry (aka GGL, Generic Geometry Library)
3 // Copyright (c) 2007-2012 Barend Gehrels, Amsterdam, the Netherlands.
4 // Copyright (c) 2018 Adam Wulkiewicz, Lodz, Poland.
6 // This file was modified by Oracle on 2013, 2014, 2015.
7 // Modifications copyright (c) 2013-2015 Oracle and/or its affiliates.
9 // Contributed and/or modified by Adam Wulkiewicz, on behalf of Oracle
11 // Use, modification and distribution is subject to the Boost Software License,
12 // Version 1.0. (See accompanying file LICENSE_1_0.txt or copy at
13 // http://www.boost.org/LICENSE_1_0.txt)
15 #include "test_relate.hpp"
16 #include "nan_cases.hpp"
20 void test_linestring_polygon()
22 typedef bg::model::linestring
<P
> ls
;
23 typedef bg::model::polygon
<P
> poly
;
24 typedef bg::model::ring
<P
> ring
;
27 test_geometry
<ls
, poly
>("LINESTRING(11 0,11 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "FF1FF0212");
28 test_geometry
<ls
, ring
>("LINESTRING(11 0,11 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "FF1FF0212");
31 test_geometry
<ls
, poly
>("LINESTRING(0 0,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1FFF0F212");
32 test_geometry
<ls
, poly
>("LINESTRING(5 0,5 5,10 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1FFF0F212");
33 test_geometry
<ls
, poly
>("LINESTRING(5 1,5 5,9 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1FF0FF212");
36 test_geometry
<ls
, poly
>("LINESTRING(11 1,11 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "FF1FF0212");
38 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "FF1F00212");
40 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 5,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F00212");
41 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 10,0 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F00212");
42 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 0,0 0)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F00212");
43 test_geometry
<ls
, poly
>("LINESTRING(0 -1,1 0,2 0)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F00212");
45 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 5,5 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1010F0212");
47 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 5,11 5,11 1)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F01FFF212");
49 test_geometry
<ls
, poly
>("LINESTRING(11 1,10 5,10 10,11 5,11 1)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11FFF212");
52 test_geometry
<ls
, poly
>("LINESTRING(0 0,5 0,5 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "11F00F212");
54 test_geometry
<ls
, poly
>("LINESTRING(5 0,5 5,10 5,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "11FF0F212");
57 test_geometry
<ls
, poly
>("LINESTRING(1 0,2 0,3 1,4 0,5 0)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "11FF0F212");
59 test_geometry
<ls
, poly
>("LINESTRING(1 0,2 0,3 -1,4 0,5 0)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F0F212");
62 test_geometry
<ls
, poly
>("LINESTRING(5 5,10 5,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "11F00F212");
64 test_geometry
<ls
, poly
>("LINESTRING(10 10,10 5,5 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "11F00F212");
66 test_geometry
<ls
, poly
>("LINESTRING(15 5,10 5,10 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F00212");
68 test_geometry
<ls
, poly
>("LINESTRING(10 10,10 5,15 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "F11F00212");
73 test_geometry
<ls
, poly
>("LINESTRING(5 5,10 5,15 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1010F0212");
74 test_geometry
<ls
, poly
>("LINESTRING(5 5,5 5,5 5,10 5,10 5,10 5,15 10,15 10,15 10)",
75 "POLYGON((0 0,0 0,0 0,0 10,0 10,0 10,10 10,10 10,10 10,10 0,10 0,10 0,0 0,0 0,0 0))",
77 test_geometry
<ls
, poly
>("LINESTRING(5 5,5 5,5 5,10 0,10 0,10 0,15 10,15 10,15 10)",
78 "POLYGON((0 0,0 0,0 0,0 10,0 10,0 10,10 10,10 10,10 10,10 0,10 0,10 0,0 0,0 0,0 0))",
81 test_geometry
<ls
, poly
>("LINESTRING(15 10,15 10,15 10,10 5,10 5,10 5,5 5,5 5,5 5)",
82 "POLYGON((0 0,0 0,0 0,0 10,0 10,0 10,10 10,10 10,10 10,10 0,10 0,10 0,0 0,0 0,0 0))",
84 test_geometry
<ls
, poly
>("LINESTRING(15 10,15 10,15 10,10 0,10 0,10 0,5 5,5 5,5 5)",
85 "POLYGON((0 0,0 0,0 0,0 10,0 10,0 10,10 10,10 10,10 10,10 0,10 0,10 0,0 0,0 0,0 0))",
89 //test_geometry<ls, poly>("LINESTRING(5 5,5 5,5 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1010F0212");
90 test_geometry
<ls
, poly
>("LINESTRING(5 5,5 5,5 5,15 5,15 5,15 5)", "POLYGON((0 0,0 10,10 10,10 0,0 0))", "1010F0212");
93 // non-simple polygon with hole
94 test_geometry
<ls
, poly
>("LINESTRING(9 1,10 5,9 9)",
95 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
97 test_geometry
<ls
, poly
>("LINESTRING(10 1,10 5,10 9)",
98 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
100 test_geometry
<ls
, poly
>("LINESTRING(2 8,10 5,2 2)",
101 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
103 test_geometry
<ls
, poly
>("LINESTRING(10 1,10 5,2 2)",
104 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
106 test_geometry
<ls
, poly
>("LINESTRING(10 1,10 5,2 8)",
107 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
110 // non-simple polygon with hole, linear ring
111 test_geometry
<ls
, poly
>("LINESTRING(9 1,10 5,9 9,1 9,1 1,9 1)",
112 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
114 test_geometry
<ls
, poly
>("LINESTRING(10 5,10 9,11 5,10 1,10 5)",
115 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
117 test_geometry
<ls
, poly
>("LINESTRING(11 5,10 1,10 5,10 9,11 5)",
118 "POLYGON((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5))",
121 // non-simple polygon with self-touching holes
122 test_geometry
<ls
, poly
>("LINESTRING(7 1,8 5,7 9)",
123 "POLYGON((0 0,0 10,10 10,10 0,0 0),(8 1,9 1,9 9,8 9,8 1),(2 2,8 5,2 8,2 2))",
125 test_geometry
<ls
, poly
>("LINESTRING(8 2,8 5,8 8)",
126 "POLYGON((0 0,0 10,10 10,10 0,0 0),(8 1,9 1,9 9,8 9,8 1),(2 2,8 5,2 8,2 2))",
128 test_geometry
<ls
, poly
>("LINESTRING(2 8,8 5,2 2)",
129 "POLYGON((0 0,0 10,10 10,10 0,0 0),(8 1,9 1,9 9,8 9,8 1),(2 2,8 5,2 8,2 2))",
132 // non-simple polygon self-touching
133 test_geometry
<ls
, poly
>("LINESTRING(9 1,10 5,9 9)",
134 "POLYGON((0 0,0 10,10 10,10 5,2 8,2 2,10 5,10 0,0 0))",
136 test_geometry
<ls
, poly
>("LINESTRING(10 1,10 5,10 9)",
137 "POLYGON((0 0,0 10,10 10,10 5,2 8,2 2,10 5,10 0,0 0))",
139 test_geometry
<ls
, poly
>("LINESTRING(2 8,10 5,2 2)",
140 "POLYGON((0 0,0 10,10 10,10 5,2 8,2 2,10 5,10 0,0 0))",
143 // non-simple polygon self-touching, linear ring
144 test_geometry
<ls
, poly
>("LINESTRING(9 1,10 5,9 9,1 9,1 1,9 1)",
145 "POLYGON((0 0,0 10,10 10,10 5,2 8,2 2,10 5,10 0,0 0))",
147 test_geometry
<ls
, poly
>("LINESTRING(10 5,10 9,11 5,10 1,10 5)",
148 "POLYGON((0 0,0 10,10 10,10 5,2 8,2 2,10 5,10 0,0 0))",
150 test_geometry
<ls
, poly
>("LINESTRING(11 5,10 1,10 5,10 9,11 5)",
151 "POLYGON((0 0,0 10,10 10,10 5,2 8,2 2,10 5,10 0,0 0))",
154 // polygons with some ring equal to the linestring
155 test_geometry
<ls
, poly
>("LINESTRING(0 0,10 0,10 10,0 10,0 0)",
156 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
158 test_geometry
<ls
, poly
>("LINESTRING(0 0,10 0,10 10,0 10,0 0)",
159 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,5 5,2 8,2 2))",
161 test_geometry
<ls
, poly
>("LINESTRING(2 2,5 5,2 8,2 2)",
162 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,5 5,2 8,2 2))",
165 // self-IP going on the boundary then into the exterior and to the boundary again
166 test_geometry
<ls
, poly
>("LINESTRING(2 10,5 10,5 15,6 15,5 10,8 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))",
168 // self-IP going on the boundary then into the interior and to the boundary again
169 test_geometry
<ls
, poly
>("LINESTRING(2 10,5 10,5 5,6 5,5 10,8 10)", "POLYGON((0 0,0 10,10 10,10 0,0 0))",
172 // self-IP with a hole -> B to I to B to E
173 test_geometry
<ls
, poly
>("LINESTRING(0 0,3 3)", "POLYGON((0 0,0 10,10 10,10 0,0 0),(0 0,9 1,9 9,1 9,0 0))",
178 typedef bg::model::polygon
<P
, false> ccwpoly
;
181 test_geometry
<ls
, ccwpoly
>("LINESTRING(11 1,10 5,5 5)", "POLYGON((0 0,10 0,10 10,0 10,0 0))", "1010F0212");
183 test_geometry
<ls
, ccwpoly
>("LINESTRING(11 1,10 1,10 5,5 5)", "POLYGON((0 0,10 0,10 10,0 10,0 0))", "1110F0212");
184 test_geometry
<ls
, ccwpoly
>("LINESTRING(11 1,10 5,10 1,5 5)", "POLYGON((0 0,10 0,10 10,0 10,0 0))", "1110F0212");
186 test_geometry
<ls
, ccwpoly
>("LINESTRING(5 1,10 5,11 1)", "POLYGON((0 0,10 0,10 10,0 10,0 0))", "1010F0212");
188 test_geometry
<ls
, ccwpoly
>("LINESTRING(5 5,10 1,10 5,11 5)", "POLYGON((0 0,10 0,10 10,0 10,0 0))", "1110F0212");
189 test_geometry
<ls
, ccwpoly
>("LINESTRING(5 5,10 5,10 1,11 5)", "POLYGON((0 0,10 0,10 10,0 10,0 0))", "1110F0212");
196 test_geometry
<ls
, poly
>("LINESTRING(0 0,2 2,3 3,1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
197 test_geometry
<ls
, poly
>("LINESTRING(0 0,3 3,1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
198 test_geometry
<ls
, poly
>("LINESTRING(0 0,2 2,1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
199 test_geometry
<ls
, poly
>("LINESTRING(1 1,3 3,2 2)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
200 test_geometry
<ls
, poly
>("LINESTRING(1 1,2 2,1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FFFF212");
202 test_geometry
<ls
, poly
>("LINESTRING(3 3,1 1,0 0,2 2)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
203 test_geometry
<ls
, poly
>("LINESTRING(3 3,0 0,2 2)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
204 test_geometry
<ls
, poly
>("LINESTRING(2 2,0 0,1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
205 test_geometry
<ls
, poly
>("LINESTRING(3 3,1 1,2 2)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FF0F212");
206 test_geometry
<ls
, poly
>("LINESTRING(2 2,1 1,2 2)", "POLYGON((0 0,3 3,3 0,0 0))", "F1FFFF212");
208 test_geometry
<ls
, poly
>("LINESTRING(0 0,2 2,4 4,1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F11F0F212");
210 test_geometry
<ls
, poly
>("LINESTRING(0 1,1 1,0 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F01FFF212");
211 test_geometry
<ls
, poly
>("LINESTRING(0 1,3 3,0 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F01FFF212");
212 test_geometry
<ls
, poly
>("LINESTRING(0 1,0 0,0 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F01FFF212");
214 test_geometry
<ls
, poly
>("LINESTRING(0 1,1 1,-1 1)", "POLYGON((0 0,3 3,3 0,0 0))", "F01FF0212");
219 // note that those are invalid Polygons and in general the result is undefined
221 // MySQL bug 15.12.2014 (assertion for invalid polygon)
222 test_geometry
<ls
, poly
>("LINESTRING(6 3,9 0)",
223 "POLYGON((0 0,5 8,6 1,6 3,8 1,5 4,9 6,2 5,7 4,1 7,0 0))",
225 test_geometry
<ls
, poly
>("LINESTRING(9 0,2 7)",
226 "POLYGON((4 1,3 2,5 9,8 4,4 5,3 6,8 1,6 2,2 4,6 0,4 1))",
229 test_geometry
<ls
, poly
>("LINESTRING(6 3,9 0)",
230 "POLYGON((6 1,6 3,8 1,5 4,6 1))",
232 test_geometry
<ls
, poly
>("LINESTRING(6 3,8 1,9 0)",
233 "POLYGON((6 1,6 3,8 1,5 4,6 1))",
236 // entry-point spikes are not handled
237 test_geometry
<ls
, poly
>("LINESTRING(9 0,6 3)",
238 "POLYGON((6 1,6 3,8 1,5 4,6 1))",
240 test_geometry
<ls
, poly
>("LINESTRING(9 0,8 1,6 3)",
241 "POLYGON((6 1,6 3,8 1,5 4,6 1))",
244 // Ls going out-in on 2 collinear spikes touching each other
245 // this works partially, other invalid effect is noticeable here,
246 // the bottom of a spike shouldn't be there so the algorithm
247 // thinks the linestring is going into the interior
248 test_geometry
<ls
, poly
>("LINESTRING(2 2,6 6)",
249 "POLYGON((0 0,0 2,2 2,4 4,3 3,5 3,5 5,4 4,7 7,7 0,0 0))",
251 // Ls going out-in on 2 collinear spikes NOT touching each other
252 // This is not supported, plus the same as above
253 test_geometry
<ls
, poly
>("LINESTRING(2 2,6 6)",
254 "POLYGON((0 0,0 2,2 2,3 3,2 2,5 3,5 5,4 4,7 7,7 0,0 0))",
259 test_geometry
<ls
, poly
>("LINESTRING(3 7, 9 1)",
260 "POLYGON((0 0,0 10,10 10,10 0,5 5,0 0))",
263 // inspired by the bug for mls/poly
264 test_geometry
<ls
, poly
>("LINESTRING(0 0, 7 7)",
265 "POLYGON((5 5,5 15,15 15,15 5,5 5))",
269 test_geometry
<ls
, poly
>("LINESTRING(4 5, 0 6, 5 6)",
270 "POLYGON((0 0,0 10,10 10,10 0,0 0),(4 4,6 4,6 6,4 6,4 4))",
272 test_geometry
<ls
, poly
>("LINESTRING(2 3, 4 5, 0 6, 5 6)",
273 "POLYGON((0 0,0 10,10 10,10 0,0 0),(4 4,6 4,6 6,4 6,4 4))",
275 test_geometry
<ls
, poly
>("LINESTRING(2 9, 1 1, 10 1, 10 10, 1 10, 0 6, 5 6)",
276 "POLYGON((0 0,0 10,10 10,10 0,0 0),(4 4,4 6,6 6,6 4,4 4))",
280 template <typename P
>
281 void test_linestring_multi_polygon()
283 typedef bg::model::linestring
<P
> ls
;
284 typedef bg::model::polygon
<P
> poly
;
285 typedef bg::model::multi_polygon
<poly
> mpoly
;
287 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,10 9)",
288 "MULTIPOLYGON(((0 20,0 30,10 30,10 20,0 20)),((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)))",
290 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,10 9)",
291 "MULTIPOLYGON(((0 20,0 30,10 30,10 20,0 20)),((0 0,0 10,10 10,10 0,0 0)))",
294 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,2 2)",
295 "MULTIPOLYGON(((0 20,0 30,10 30,10 20,0 20)),((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)))",
297 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,2 2)",
298 "MULTIPOLYGON(((0 20,0 30,10 30,10 20,0 20)),((0 0,0 10,10 10,10 0,0 0)))",
301 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,2 2)",
302 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)),((10 5,3 3,3 7,10 5)))",
304 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,2 8)",
305 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)),((10 5,3 3,3 7,10 5)))",
307 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,3 3)",
308 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)),((10 5,3 3,3 7,10 5)))",
310 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,3 7)",
311 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)),((10 5,3 3,3 7,10 5)))",
313 test_geometry
<ls
, mpoly
>("LINESTRING(10 1,10 5,5 5)",
314 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(10 5,2 8,2 2,10 5)),((10 5,3 3,3 7,10 5)))",
317 test_geometry
<ls
, mpoly
>("LINESTRING(0 0,10 0,10 10,0 10,0 0)",
318 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((20 0,20 10,30 20,20 0)))",
321 // degenerated points
322 test_geometry
<ls
, mpoly
>("LINESTRING(5 5,10 10,10 10,10 10,15 15)",
323 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((10 10,10 20,20 20,20 10,10 10)))",
326 // self-IP polygon with a hole and second polygon with a hole -> B to I to B to B to I to B to E
327 test_geometry
<ls
, mpoly
>("LINESTRING(0 0,3 3)",
328 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(0 0,9 1,9 9,1 9,0 0)),((0 0,2 8,8 8,8 2,0 0),(0 0,7 3,7 7,3 7,0 0)))",
330 // self-IP polygon with a hole and second polygon -> B to I to B to B to I
331 test_geometry
<ls
, mpoly
>("LINESTRING(0 0,3 3)",
332 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(0 0,9 1,9 9,1 9,0 0)),((0 0,2 8,8 8,8 2,0 0)))",
334 test_geometry
<ls
, mpoly
>("LINESTRING(0 0,3 3)",
335 "MULTIPOLYGON(((0 0,2 8,8 8,8 2,0 0)),((0 0,0 10,10 10,10 0,0 0),(0 0,9 1,9 9,1 9,0 0)))",
338 // MySQL report 18.12.2014 (https://svn.boost.org/trac/boost/ticket/10887)
339 test_geometry
<ls
, mpoly
>("LINESTRING(5 -2,5 2)",
340 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
342 test_geometry
<ls
, mpoly
>("LINESTRING(5 -2,5 5)",
343 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
345 test_geometry
<ls
, mpoly
>("LINESTRING(5 -2,5 0)",
346 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
348 // MySQL report 18.12.2014 - extended
349 test_geometry
<ls
, mpoly
>("LINESTRING(5 -2,5 0)",
350 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)),((5 0,7 1,7 -1,5 0)))",
352 test_geometry
<ls
, mpoly
>("LINESTRING(0 0,5 0)",
353 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)),((5 0,7 1,7 -1,5 0)))",
357 test_geometry
<ls
, mpoly
>("LINESTRING(5 5,0 0,10 0)",
358 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
360 test_geometry
<ls
, mpoly
>("LINESTRING(5 5,0 0,0 10)",
361 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
364 test_geometry
<ls
, mpoly
>("LINESTRING(5 5,0 0,2 1)",
365 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
367 test_geometry
<ls
, mpoly
>("LINESTRING(5 5,0 0,5 -5)",
368 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
370 test_geometry
<ls
, mpoly
>("LINESTRING(5 5,0 0,5 -5,5 1)",
371 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
373 test_geometry
<ls
, mpoly
>("LINESTRING(-5 5,0 0,5 -5)",
374 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
378 template <typename P
>
379 void test_multi_linestring_polygon()
381 typedef bg::model::linestring
<P
> ls
;
382 typedef bg::model::multi_linestring
<ls
> mls
;
383 typedef bg::model::polygon
<P
> poly
;
384 typedef typename
bg::coordinate_type
<P
>::type coord_t
;
386 test_geometry
<mls
, poly
>("MULTILINESTRING((11 11, 20 20),(5 7, 4 1))",
387 "POLYGON((0 0,0 10,10 10,10 0,0 0),(2 2,4 2,4 4,2 4,2 2))",
389 test_geometry
<mls
, poly
>("MULTILINESTRING((10 0, 18 12),(2 2,2 1))",
390 "POLYGON((5 0,0 -5,-5 0,0 5,5 0))",
393 if ( BOOST_GEOMETRY_CONDITION((
394 boost::is_same
<coord_t
, double>::value
)) )
396 // assertion failure in 1.57
397 test_geometry
<mls
, poly
>("MULTILINESTRING((-0.59322033898305082 -8.0508474576271194,-2.882352941176471 -7.7647058823529411,-2.8823529411764706 -7.7647058823529411,-3.7361111111111112 -6.5694444444444446,-3.4404145077720205 -5.766839378238342,-4.1864406779661012 -3.6779661016949152,-7.5252525252525251 -5.5858585858585865,-7.5862068965517242 -5.1896551724137936,-4.47887323943662 -2.859154929577465,-4.5789473684210531 -2.5789473684210527,-3 -1,-2.9310344827586206 -0.86206896551724144,-3.1764705882352944 -0.70588235294117663,-4.7401960784313726 -2.1274509803921577,-5.3255813953488369 -0.48837209302325502,-4.7872340425531918 0.31914893617021284,-5.8571428571428577 1.0000000000000007,-5.3255813953488369 -0.48837209302325502,-5.9473684210526319 -1.4210526315789465,-8 2,-7.7333333333333334 2.1939393939393939,-8.8294573643410867 2.891472868217055,-8.8556701030927822 3.061855670103093,-7.5999999999999996 3.6000000000000001,-7.7999999999999998 3.7999999999999998,-7.75 3.7916666666666665,-7.5471698113207548 3.6226415094339623,-7.3200000000000003 3.7200000000000002,-3.473684210526315 3.0789473684210527,-3.2549019607843133 3.2156862745098036,-2.9999999999999982 3.1428571428571423,-3.1733333333333325 3.2666666666666666,-2.9180327868852456 3.4262295081967209,-2.8723404255319145 3.1063829787234041,-2.1111111111111112 2.8888888888888888,-2.1428571428571428 2.8571428571428572,-1.8433734939759043 2.8072289156626509,-1.8396226415094346 2.8113207547169816,-1.6486486486486487 2.756756756756757,-1.76510067114094 2.8926174496644301,-0.53846153846153855 4.2307692307692308,1.8148148148148147 5.4074074074074074,1.588235294117647 2.2352941176470589,1.819672131147541 2.1967213114754101,2 4,2 2.1666666666666665,2.3538461538461544 2.1076923076923078,2 1.6875000000000004,2 -2,1.2173913043478262 -3.8260869565217392,1.7375886524822697 1.3758865248226959,1.5073170731707317 1.1024390243902444,1.1428571428571428 -4,-0.59322033898305082 -8.0508474576271194),(1.666666666666667 1.5999999999999988,1.5675675675675675 1.8378378378378377,1.4374999999999991 1.8750000000000002,1.0487804878048776 2.3414634146341466,0.46666666666666712 2.6060606060606055,0.086956521739131043 2.2608695652173911,1.4374999999999991 1.8750000000000002,1.666666666666667 1.5999999999999988))",
398 "POLYGON((-2.333333333333333 -8.6666666666666661,-4.3253012048192767 -8.168674698795181,-4.1194968553459113 -7.6100628930817606,-2.8823529411764706 -7.7647058823529411,-2.882352941176471 -7.7647058823529411,-2.263157894736842 -8.6315789473684212,-2.333333333333333 -8.6666666666666661))",
400 test_geometry
<mls
, poly
>("MULTILINESTRING((-2.333333333333333 -8.6666666666666661,-4.3253012048192767 -8.168674698795181,-4.1194968553459113 -7.6100628930817606,-2.8823529411764706 -7.7647058823529411,-2.882352941176471 -7.7647058823529411,-2.263157894736842 -8.6315789473684212,-2.333333333333333 -8.6666666666666661))",
401 "POLYGON((-0.59322033898305082 -8.0508474576271194,-2.882352941176471 -7.7647058823529411,-2.8823529411764706 -7.7647058823529411,-3.7361111111111112 -6.5694444444444446,-3.4404145077720205 -5.766839378238342,-4.1864406779661012 -3.6779661016949152,-7.5252525252525251 -5.5858585858585865,-7.5862068965517242 -5.1896551724137936,-4.47887323943662 -2.859154929577465,-4.5789473684210531 -2.5789473684210527,-3 -1,-2.9310344827586206 -0.86206896551724144,-3.1764705882352944 -0.70588235294117663,-4.7401960784313726 -2.1274509803921577,-5.3255813953488369 -0.48837209302325502,-4.7872340425531918 0.31914893617021284,-5.8571428571428577 1.0000000000000007,-5.3255813953488369 -0.48837209302325502,-5.9473684210526319 -1.4210526315789465,-8 2,-7.7333333333333334 2.1939393939393939,-8.8294573643410867 2.891472868217055,-8.8556701030927822 3.061855670103093,-7.5999999999999996 3.6000000000000001,-7.7999999999999998 3.7999999999999998,-7.75 3.7916666666666665,-7.5471698113207548 3.6226415094339623,-7.3200000000000003 3.7200000000000002,-3.473684210526315 3.0789473684210527,-3.2549019607843133 3.2156862745098036,-2.9999999999999982 3.1428571428571423,-3.1733333333333325 3.2666666666666666,-2.9180327868852456 3.4262295081967209,-2.8723404255319145 3.1063829787234041,-2.1111111111111112 2.8888888888888888,-2.1428571428571428 2.8571428571428572,-1.8433734939759043 2.8072289156626509,-1.8396226415094346 2.8113207547169816,-1.6486486486486487 2.756756756756757,-1.76510067114094 2.8926174496644301,-0.53846153846153855 4.2307692307692308,1.8148148148148147 5.4074074074074074,1.588235294117647 2.2352941176470589,1.819672131147541 2.1967213114754101,2 4,2 2.1666666666666665,2.3538461538461544 2.1076923076923078,2 1.6875000000000004,2 -2,1.2173913043478262 -3.8260869565217392,1.7375886524822697 1.3758865248226959,1.5073170731707317 1.1024390243902444,1.1428571428571428 -4,-0.59322033898305082 -8.0508474576271194),(1.666666666666667 1.5999999999999988,1.5675675675675675 1.8378378378378377,1.4374999999999991 1.8750000000000002,1.0487804878048776 2.3414634146341466,0.46666666666666712 2.6060606060606055,0.086956521739131043 2.2608695652173911,1.4374999999999991 1.8750000000000002,1.666666666666667 1.5999999999999988))",
406 test_geometry
<mls
, poly
>("MULTILINESTRING((6 6,15 15),(0 0, 7 7))",
407 "POLYGON((5 5,5 15,15 15,15 5,5 5))",
409 test_geometry
<mls
, poly
>("MULTILINESTRING((15 15,6 6),(0 0, 7 7))",
410 "POLYGON((5 5,5 15,15 15,15 5,5 5))",
413 test_geometry
<mls
, poly
>("MULTILINESTRING((15 15,6 6),(4 14,6 16))",
414 "POLYGON((5 5,5 15,15 15,15 5,5 5))",
418 test_geometry
<mls
, poly
>("MULTILINESTRING((4 10, 3 10, 10 6),(5 0, 7 5, 9 10))",
419 "POLYGON((0 0,0 10,10 10,10 0,5 5,0 0))",
423 test_geometry
<mls
, poly
>("MULTILINESTRING((3 10, 1 5, 1 10, 3 4, 7 8, 6 10, 10 2))",
424 "POLYGON((0 0,0 10,10 10,10 0,0 0))",
428 // assertion failure in relate->boundary_checker->std::equal_range with msvc
429 if (BOOST_GEOMETRY_CONDITION(is_nan_case_supported
<mls
>::value
))
435 std::string wkt2
= "POLYGON((1.1e+308 1.2e+308,-1 -9,1 1e+12,1.1e+308 7.8e+307,1.1e+308 1.2e+308),(3 2,1 1,8e+307 1e+308,3 2),(258 2049,1 -3,1 1,-6 9,258 2049))";
437 bg::read_wkt(wkt2
, g2
);
439 check_geometry(g1
, g2
, wkt1
, wkt2
, "*********");
443 template <typename P
>
444 void test_multi_linestring_multi_polygon()
446 typedef bg::model::linestring
<P
> ls
;
447 typedef bg::model::polygon
<P
> poly
;
448 typedef bg::model::multi_linestring
<ls
> mls
;
449 typedef bg::model::multi_polygon
<poly
> mpoly
;
451 // polygons with some ring equal to the linestrings
452 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,10 0,10 10,0 10,0 0),(20 20,50 50,20 80,20 20))",
453 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
456 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,10 0,10 10,0 10,0 0),(2 2,5 5,2 8,2 2))",
457 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0),(2 2,5 5,2 8,2 2)))",
461 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,10 0,10 10),(10 10,0 10,0 0))",
462 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
464 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,10 0,10 10),(10 10,0 10,0 0),(20 20,50 50,20 80,20 20))",
465 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
469 test_geometry
<mls
, mpoly
>("MULTILINESTRING((20 20,30 30),(30 30,40 40))",
470 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)))",
473 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 5,0 5),(5 5,5 0),(10 10,10 5,5 5,5 10,10 10))",
474 "MULTIPOLYGON(((0 0,0 5,5 5,5 0,0 0)),((5 5,5 10,10 10,10 5,5 5)),((5 5,10 1,10 0,5 5)))",
476 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 5,0 5),(5 5,5 0),(0 5,0 0,5 0),(10 10,10 5,5 5,5 10,10 10))",
477 "MULTIPOLYGON(((0 0,0 5,5 5,5 0,0 0)),((5 5,5 10,10 10,10 5,5 5)),((5 5,10 1,10 0,5 5)))",
479 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 5,0 0),(5 5,5 0),(10 10,10 5,5 5,5 10,10 10))",
480 "MULTIPOLYGON(((0 0,0 5,5 5,5 0,0 0)),((5 5,5 10,10 10,10 5,5 5)),((5 5,10 1,10 0,5 5)))",
483 // MySQL report 18.12.2014 - extended
484 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 -2,4 -2,5 0),(5 -2,6 -2,5 0))",
485 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
487 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,0 1,5 0),(0 0,0 -1,5 0))",
488 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
490 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 -2,4 -2,5 0),(6 -2,5 0))",
491 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
493 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,0 1,5 0),(0 -1,5 0))",
494 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
496 test_geometry
<mls
, mpoly
>("MULTILINESTRING((0 0,5 0),(5 -2,5 0))",
497 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
499 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 -2,5 0),(0 0,5 0))",
500 "MULTIPOLYGON(((5 0,0 5,10 5,5 0)),((5 0,10 -5,0 -5,5 0)))",
503 // 22.01.2015 - extended
504 test_geometry
<mls
, mpoly
>("MULTILINESTRING((10 10,0 10),(5 5,0 0,10 0))",
505 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
507 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 5,0 0,5 -5),(0 0,9 1))",
508 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
510 test_geometry
<mls
, mpoly
>("MULTILINESTRING((5 -5,0 0,5 5),(0 0,5 -1))",
511 "MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((0 0,0 -10,-10 -10,-10 0,0 0)))",
515 template <typename P
>
518 test_linestring_polygon
<P
>();
519 test_linestring_multi_polygon
<P
>();
520 test_multi_linestring_polygon
<P
>();
521 test_multi_linestring_multi_polygon
<P
>();
524 int test_main( int , char* [] )
526 test_all
<bg::model::d2::point_xy
<int> >();
527 test_all
<bg::model::d2::point_xy
<double> >();