boost/geometry
2
0
Fork 0
mirror of https://github.com/boostorg/geometry.git synced 2026-01-20 16:42:12 +00:00
Code Issues Projects Releases Wiki Activity
Files
develop
geometry/test/algorithms/set_operations/difference/difference_linear_linear.cpp
Vissarion Fisikopoulos 75c7dfded7 chore: Replace typedefs by using
2025-03-18 13:20:09 +02:00

1464 lines
54 KiB
C++
Raw Permalink Blame History

// Boost.Geometry (aka GGL, Generic Geometry Library)
// Copyright (c) 2014-2025, Oracle and/or its affiliates.
// Licensed under the Boost Software License version 1.0.
// http://www.boost.org/users/license.html
// Contributed and/or modified by Vissarion Fysikopoulos, on behalf of Oracle
// Contributed and/or modified by Menelaos Karavelas, on behalf of Oracle
#include <iostream>
#ifndef BOOST_TEST_MODULE
#define BOOST_TEST_MODULE test_difference_linear_linear
#endif
#ifdef BOOST_GEOMETRY_TEST_DEBUG
#define BOOST_GEOMETRY_DEBUG_TURNS
#define BOOST_GEOMETRY_DEBUG_SEGMENT_IDENTIFIER
#endif
#include <boost/test/included/unit_test.hpp>
#include "test_difference_linear_linear.hpp"
#include <boost/geometry/geometries/linestring.hpp>
#include <boost/geometry/geometries/multi_linestring.hpp>
#include <boost/geometry/algorithms/difference.hpp>
using point_type = bg::model::point<double,2,bg::cs::cartesian>;
using segment_type = bg::model::segment<point_type>;
using linestring_type = bg::model::linestring<point_type>;
using multi_linestring_type = bg::model::multi_linestring<linestring_type>;
using L = linestring_type;
using ML = multi_linestring_type;
//===========================================================================
//===========================================================================
//===========================================================================
BOOST_AUTO_TEST_CASE( test_difference_linestring_linestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** LINESTRING / LINESTRING DIFFERENCE ***" << std::endl;
std::cout << std::endl;
#endif
using tester = test_difference_of_geometries<L, L, ML>;
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 1,2 1,3 2)"),
from_wkt<L>("LINESTRING(0 2,1 1,2 1,3 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1),(2 1,3 2))"),
"lldf00");
tester::apply
(from_wkt<L>("LINESTRING(0 2,1 1,2 1,3 0)"),
from_wkt<L>("LINESTRING(0 0,1 1,2 1,3 2)"),
from_wkt<ML>("MULTILINESTRING((0 2,1 1),(2 1,3 0))"),
"lldf00-1");
tester::apply
(from_wkt<L>("LINESTRING(0 0,5 0)"),
from_wkt<L>("LINESTRING(3 0,4 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,3 0),(4 0,5 0))"),
"lldf01");
tester::apply
(from_wkt<L>("LINESTRING(3 0,4 0)"),
from_wkt<L>("LINESTRING(0 0,5 0)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf01-1");
tester::apply
(from_wkt<L>("LINESTRING(0 0,4 0)"),
from_wkt<L>("LINESTRING(3 0,6 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,3 0))"),
"lldf01-2");
tester::apply
(from_wkt<L>("LINESTRING(3 0,6 0)"),
from_wkt<L>("LINESTRING(0 0,4 0)"),
from_wkt<ML>("MULTILINESTRING((4 0,6 0))"),
"lldf01-3");
tester::apply
(from_wkt<L>("LINESTRING(0 0,6 0)"),
from_wkt<L>("LINESTRING(0 0,4 0)"),
from_wkt<ML>("MULTILINESTRING((4 0,6 0))"),
"lldf01-4");
tester::apply
(from_wkt<L>("LINESTRING(0 0,4 0)"),
from_wkt<L>("LINESTRING(0 0,6 0)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf01-5");
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<L>("LINESTRING(0 0,1 1,2 0,3 1,4 0,5 0,6 1,7 -1,8 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,4 0),(5 0,20 0))"),
"lldf01-6");
tester::apply
(from_wkt<L>("LINESTRING(-20 0,20 0)"),
from_wkt<L>("LINESTRING(0 0,1 1,2 0,3 1,4 0,5 0,6 1,7 -1,8 0)"),
from_wkt<ML>("MULTILINESTRING((-20 0,4 0),(5 0,20 0))"),
"lldf01-7");
tester::apply
(from_wkt<L>("LINESTRING(0 0,4 0)"),
from_wkt<L>("LINESTRING(2 0,4 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0))"),
"lldf01-8");
tester::apply
(from_wkt<L>("LINESTRING(2 0,4 0)"),
from_wkt<L>("LINESTRING(0 0,4 0)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf01-9");
tester::apply
(from_wkt<L>("LINESTRING(0 0,2 0)"),
from_wkt<L>("LINESTRING(4 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0))"),
"lldf01-10");
tester::apply
(from_wkt<L>("LINESTRING(0 0,2 0)"),
from_wkt<L>("LINESTRING(2 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0))"),
"lldf01-11");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,4 0)"),
from_wkt<L>("LINESTRING(3 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0,3 0))"),
"lldf01-11a");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,4 0)"),
from_wkt<L>("LINESTRING(3 0,4 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0,3 0))"),
"lldf01-11b");
tester::apply
(from_wkt<L>("LINESTRING(0 0,6 0)"),
from_wkt<L>("LINESTRING(2 0,4 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(5 0,6 0))"),
"lldf01-12");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,5 5,10 5,15 0)"),
from_wkt<L>("LINESTRING(-1 6,0 5,15 5)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0,5 5),(10 5,15 0))"),
"lldf02");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0)"),
from_wkt<L>("LINESTRING(-1 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0))"),
"lldf03");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,5 5,10 5,15 0,20 0)"),
from_wkt<L>("LINESTRING(-1 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((1 0,5 5,10 5,15 0))"),
"lldf04");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,25 1)"),
from_wkt<L>("LINESTRING(-1 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0),\
(20 0,25 1))"),
"lldf05");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,30 0)"),
from_wkt<L>("LINESTRING(-1 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0))"),
"lldf05-1");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,31 0)"),
from_wkt<L>("LINESTRING(-1 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0),\
(30 0,31 0))"),
"lldf06");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,31 0)"),
from_wkt<L>("LINESTRING(-1 0,25 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0),\
(30 0,31 0))"),
"lldf07");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,31 0)"),
from_wkt<L>("LINESTRING(-1 0,19 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0),\
(30 0,31 0))"),
"lldf08");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,30 0,31 1)"),
from_wkt<L>("LINESTRING(-1 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0),\
(30 0,31 1))"),
"lldf09");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,5 5,10 5,15 0,20 0,30 0,31 1)"),
from_wkt<L>("LINESTRING(-1 -1,0 0,1 0,2 1,3 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(1 0,5 5,10 5,15 0),\
(30 0,31 1))"),
"lldf10");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,4 0,5 5,10 5,15 0,20 0,\
30 0,31 1)"),
from_wkt<L>("LINESTRING(-1 -1,0 0,1 0,2 0,2.5 1,3 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(2 0,3 0),\
(4 0,5 5,10 5,15 0),(30 0,31 1))"),
"lldf11");
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,1 0,4 0,5 5,10 5,15 0,31 0)"),
from_wkt<L>("LINESTRING(-1 -1,0 0,1 0,2 0,2.5 1,3 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0),(2 0,3 0),\
(4 0,5 5,10 5,15 0),(30 0,31 0))"),
"lldf11-1");
tester::apply
(from_wkt<L>("LINESTRING(0 0,2 0,3 1)"),
from_wkt<L>("LINESTRING(0 0,2 0,3 1)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf12");
tester::apply
(from_wkt<L>("LINESTRING(0 0,2 0,3 1)"),
from_wkt<L>("LINESTRING(3 1,2 0,0 0)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf12-1");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,2 1,3 5,4 0)"),
from_wkt<L>("LINESTRING(1 0,2 1,3 5,4 0,5 10)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"lldf13");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,2 0,2.5 0,3 1)"),
from_wkt<L>("LINESTRING(0 0,2 0,2.5 0,3 1)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf14");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,2 1,3 5,4 0)"),
from_wkt<L>("LINESTRING(1 0,2 1,3 5)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(3 5,4 0))"),
"lldf15");
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,2 1,3 2)"),
from_wkt<L>("LINESTRING(0.5 0,1 0,3 2,4 5)"),
from_wkt<ML>("MULTILINESTRING((0 0,0.5 0))"),
"lldf16");
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0,20 1,30 1)"),
from_wkt<L>("LINESTRING(1 1,2 0,3 1,20 1,25 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(25 1,30 1))"),
"lldf17");
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0,20 1,21 0,30 0)"),
from_wkt<L>("LINESTRING(1 1,2 0,3 1,20 1,25 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1,21 0,30 0))"),
"lldf18");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0,5 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19a");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19b");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0,5 0,6 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19c");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0,3 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19d");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0,3 0,3 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19e");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,4 0,5 0,5 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19f");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<L>("LINESTRING(1 0,5 0,20 1,4 1,5 0,5 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,30 0))"),
"lldf19g");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0,30 30,10 30,10 -10,15 0,40 0)"),
from_wkt<L>("LINESTRING(5 5,10 0,10 30,20 0,25 0,25 25,50 0,35 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,20 0),(25 0,30 0,30 30,10 30),\
(10 0,10 -10,15 0,20 0),(25 0,35 0))"),
"lldf20");
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0,30 30,10 30,10 -10,15 0,40 0)"),
from_wkt<L>("LINESTRING(5 5,10 0,10 30,20 0,25 0,25 25,50 0,15 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,15 0),(30 0,30 30,10 30),\
(10 0,10 -10,15 0))"),
"lldf20a");
tester::apply
(from_wkt<L>("LINESTRING(0 0,18 0,19 0,30 0)"),
from_wkt<L>("LINESTRING(2 2,5 -1,15 2,18 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,18 0),(20 0,30 0))"),
"lldf21"
);
tester::apply
(from_wkt<L>("LINESTRING(2 2,5 -1,15 2,18 0,20 0)"),
from_wkt<L>("LINESTRING(0 0,18 0,19 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((2 2,5 -1,15 2,18 0))"),
"lldf21a"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0)"),
from_wkt<L>("LINESTRING(1 0,4 0,2 1,5 1,4 0,8 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(8 0,10 0))"),
"lldf22"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0)"),
from_wkt<L>("LINESTRING(4 0,5 0,5 1,1 1,1 0,4 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(5 0,10 0))"),
"lldf23"
);
// the following two tests have been discussed with by Adam
tester::apply
(from_wkt<L>("LINESTRING(1 0,1 1,2 1)"),
from_wkt<L>("LINESTRING(2 1,1 1,1 0)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf24"
);
tester::apply
(from_wkt<L>("LINESTRING(1 0,1 1,2 1)"),
from_wkt<L>("LINESTRING(1 2,1 1,1 0)"),
from_wkt<ML>("MULTILINESTRING((1 1,2 1))"),
"lldf25"
);
tester::apply
(from_wkt<L>("LINESTRING(6 0,0 0,5 0)"),
from_wkt<L>("LINESTRING(2 0,-10 0)"),
from_wkt<ML>("MULTILINESTRING((6 0,2 0),(2 0,5 0))"),
"lldf27a"
);
tester::apply
(from_wkt<L>("LINESTRING(6 0,0 0,5 0)"),
from_wkt<L>("LINESTRING(2 0,-1 0,-10 0)"),
from_wkt<ML>("MULTILINESTRING((6 0,2 0),(2 0,5 0))"),
"lldf27b"
);
tester::apply
(from_wkt<L>("LINESTRING(6 0,0 0,5 0)"),
from_wkt<L>("LINESTRING(2 0,0 0,-10 0)"),
from_wkt<ML>("MULTILINESTRING((6 0,2 0),(2 0,5 0))"),
"lldf27c"
);
tester::apply
(from_wkt<L>("LINESTRING(2 0,0 0,-10 0)"),
from_wkt<L>("LINESTRING(6 0,0 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,-10 0))"),
"lldf27d"
);
tester::apply
(from_wkt<L>("LINESTRING(-3 6,-3 0,-3 5)"),
from_wkt<L>("LINESTRING(-3 2,-3 0,-3 -10)"),
from_wkt<ML>("MULTILINESTRING((-3 6,-3 2),(-3 2,-3 5))"),
"lldf28a"
);
tester::apply
(from_wkt<L>("LINESTRING(-3 2,-3 0,-3 -10)"),
from_wkt<L>("LINESTRING(-3 6,-3 0,-3 5)"),
from_wkt<ML>("MULTILINESTRING((-3 0,-3 -10))"),
"lldf28b"
);
tester::apply
(from_wkt<L>("LINESTRING(-3 6,-3 0,-3 5)"),
from_wkt<L>("LINESTRING(-3 2,-3 0,-3 -10)"),
from_wkt<ML>("MULTILINESTRING((-3 6,-3 2),(-3 2,-3 5))"),
"lldf28c"
);
tester::apply
(from_wkt<L>("LINESTRING(-7 -8,3 0,4 -1)"),
from_wkt<L>("LINESTRING(-5 -4,3 0,4 -1,7 -4)"),
from_wkt<ML>("MULTILINESTRING((-7 -8,3 0))"),
"lldf29a"
);
tester::apply
(from_wkt<L>("LINESTRING(-7 -8,3 0,4 -1,-7 10)"),
from_wkt<L>("LINESTRING(-5 -4,3 0,4 -1,2 -1)"),
from_wkt<ML>("MULTILINESTRING((-7 -8,3 0),(3 0,-7 10))"),
"lldf29b"
);
tester::apply
(from_wkt<L>("LINESTRING(-7 -8,3 0,4 -1,-7 10)"),
from_wkt<L>("LINESTRING(-5 -4,3 0,4 -1,7 -4,2 -1)"),
from_wkt<ML>("MULTILINESTRING((-7 -8,3 0),(3 0,-7 10))"),
"lldf29c"
);
tester::apply
(from_wkt<L>("LINESTRING(-2 -2,-4 0,1 -8,-2 6,8 5,-7 -8,\
3 0,4 -1,-7 10,-4 10)"),
from_wkt<L>("LINESTRING(-5 -4,3 0,4 -1,7 -4,2 -1,-4 -1,-2 6)"),
from_wkt<ML>("MULTILINESTRING((-2 -2,-4 0,1 -8,-2 6,8 5,-7 -8,\
3 0),(3 0,-7 10,-4 10))"),
"lldf29d"
);
#ifdef BOOST_GEOMETRY_TEST_FAILURES
tester::apply
(from_wkt<L>("LINESTRING(8 5,5 1,-2 3,1 10)"),
from_wkt<L>("LINESTRING(1.9375 1.875,\
1.7441860465116283 1.9302325581395348,\
-0.7692307692307692 2.6483516483516487,\
-2 3,-1.0071942446043165 5.316546762589928)"),
from_wkt<ML>("MULTILINESTRING((8 5,5 1,-2 3,1 10))"),
"lldf30a"
);
tester::apply
(from_wkt<L>("LINESTRING(1.9375 1.875,\
1.7441860465116283 1.9302325581395348,\
-0.7692307692307692 2.6483516483516487,\
-2 3,-1.0071942446043165 5.316546762589928)"),
from_wkt<L>("LINESTRING(8 5,5 1,-2 3,1 10)"),
from_wkt<ML>("MULTILINESTRING((1.9375 1.875,\
1.7441860465116283 1.9302325581395348,\
-0.7692307692307692 2.6483516483516487,\
-2 3,-1.0071942446043165 5.316546762589928))"),
"lldf30b"
);
tester::apply
(from_wkt<L>("LINESTRING(5 -8,-7 -6,-3 6,-3 1,-5 4,-1 0,8 5,\
5 1,-2 3,1 10,8 5,6 2,7 4)"),
from_wkt<L>("LINESTRING(1.9375 1.875,\
1.7441860465116283 1.9302325581395348,\
-0.7692307692307692 2.6483516483516487,\
-2 3,-1.0071942446043165 5.316546762589928)"),
from_wkt<ML>("MULTILINESTRING((5 -8,-7 -6,-3 6,-3 1,-5 4,-1 0,8 5,\
5 1,-2 3,1 10,8 5,6 2,7 4))"),
"lldf30c"
);
#endif
tester::apply
(from_wkt<L>("LINESTRING(8 1, 4 .4)"),
from_wkt<L>("LINESTRING(0 -.2, 8 1)"),
from_wkt<ML>("MULTILINESTRING()"),
"lldf31s"
);
{
ut_settings settings{1e-10};
tester::apply
(from_wkt<L>("LINESTRING(8 1, 4 .4,2 8)"),
from_wkt<L>("LINESTRING(0 -.2, 8 1)"),
from_wkt<ML>("MULTILINESTRING((4 .4,2 8))"),
"lldf31x",
settings
);
tester::apply
(from_wkt<L>("LINESTRING(0 5, 8 1, 4 .4, 2 8)"),
from_wkt<L>("LINESTRING(0 -.2, 8 1, -.5 7)"),
from_wkt<ML>("MULTILINESTRING((0 5,8 1),(4 .4,2 8))"),
"lldf31y",
settings
);
tester::apply
(from_wkt<L>("LINESTRING(0 5, 8 1, 4 .4, 2 8)"),
from_wkt<L>("LINESTRING(0 -.2, 8 1, -.5 7, 6 +.2)"),
from_wkt<ML>("MULTILINESTRING((0 5,8 1),(4 .4,2 8))"),
"lldf31y+",
settings
);
tester::apply
(from_wkt<L>("LINESTRING(10.0002 2,9 -1032.34324, .3 8, 0 5, 8 1, 4 .4, 2 8)"),
from_wkt<L>("LINESTRING(0 -.2, 8 1, -.5 7, 6 +.2)"),
from_wkt<ML>("MULTILINESTRING((10.0002 2,9 -1032.34324,.3 8,0 5,8 1),(4 .4,2 8))"),
"lldf31z",
settings
);
}
tester::apply
(from_wkt<L>("LINESTRING(0 0, 8 1, -.5 7)"),
from_wkt<L>("LINESTRING(0 5, 8 1, 4 .5, 2 8)"),
from_wkt<ML>("MULTILINESTRING((0 0,4 .5),(8 1,-.5 7))"),
"lldf32"
);
}
BOOST_AUTO_TEST_CASE( test_difference_linestring_multilinestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** LINESTRING / MULTILINESTRING DIFFERENCE ***"
<< std::endl;
std::cout << std::endl;
#endif
using tester = test_difference_of_geometries<L, ML, ML>;
// disjoint linestrings
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0,20 1)"),
from_wkt<ML>("MULTILINESTRING((1 1,2 2,4 3),(1 1,2 2,5 3))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1))"),
"lmldf01"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0,20 1)"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),(1 1,3 0,4 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(4 0,10 0,20 1))"),
"lmldf02"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0,20 1)"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),(1 1,3 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(5 0,10 0,20 1))"),
"lmldf03"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0,20 1)"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(4 0,10 0,20 1))"),
"lmldf04"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,101 0)"),
from_wkt<ML>("MULTILINESTRING((-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"lmldf07"
);
tester::apply
(from_wkt<L>("LINESTRING(-1 1,0 0,101 0)"),
from_wkt<ML>("MULTILINESTRING((-1 -1,0 0,50 0),\
(19 -1,20 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0))"),
"lmldf07a"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,101 0)"),
from_wkt<ML>("MULTILINESTRING((-1 -1,0 0,50 0),\
(19 -1,20 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING()"),
"lmldf07b"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,101 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 1,2 0),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"lmldf08"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,2 0.5,3 0,101 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 1,2 0.5),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0,2 0.5,3 0))"),
"lmldf09"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,1 0,1.5 0,2 0.5,3 0,101 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 1,1 0,2 0.5),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(1.5 0,2 0.5,3 0))"),
"lmldf10"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(1 1,2 0,18 0,19 1),(2 1,3 0,17 0,18 1),\
(3 1,4 0,16 0,17 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf12"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((1 0,19 0,20 1),\
(2 0,18 0,19 1),(3 0,17 0,18 1),\
(4 0,16 0,17 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf13"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1,19 1,18 0,2 0,\
1 1,2 1,3 0,17 0,18 1,17 1,16 0,4 0,3 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf14"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 2,6 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf15"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(6 0,4 2,2 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf15a"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 2,5 0,6 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf16"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,20 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(6 0,5 0,4 2,2 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"lmldf16a"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 0,5 2,20 2,25 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"lmldf17"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 0,5 2,20 2,25 0,26 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"lmldf17a"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,18 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"lmldf18"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,18 0,19 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,18 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"lmldf18a"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,18 0,19 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,18 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(20 0,30 0))"),
"lmldf18b"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,18 0,19 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,25 0,26 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,25 0),(26 0,30 0))"),
"lmldf18c"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,18 0,19 0,30 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,25 0,21 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,21 0),(25 0,30 0))"),
"lmldf18d"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0)"),
from_wkt<ML>("MULTILINESTRING((0 5,1 0,9 0,10 5),(0 1,2 0,3 1),\
(0 -2,3 0,4 4),(0 -5,4 0,5 0,6 3))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10))"),
"lmldf19"
);
tester::apply
(from_wkt<L>("LINESTRING(0 0,10 0)"),
from_wkt<ML>("MULTILINESTRING((-1 0,0 0),(10 0,12 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
"lmldf20"
);
}
BOOST_AUTO_TEST_CASE( test_difference_multilinestring_linestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTILINESTRING / LINESTRING DIFFERENCE ***"
<< std::endl;
std::cout << std::endl;
#endif
using tester = test_difference_of_geometries<ML, L, ML>;
// disjoint linestrings
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
from_wkt<L>("LINESTRING(1 1,2 2,4 3)"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
"mlldf01"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
from_wkt<L>("LINESTRING(1 1,2 0,4 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(4 0,10 0,20 1),\
(1 0,2 0),(4 0,7 0))"),
"mlldf02"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,101 0))"),
from_wkt<L>("LINESTRING(-1 -1,1 0,101 0,200 -1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"mlldf03"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<L>("LINESTRING(0 1,1 0,19 0,20 1,19 1,18 0,2 0,\
1 1,2 1,3 0,17 0,18 1,17 1,16 0,4 0,3 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlldf04"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((-1 0,-1 10),(0 0,20 0),(25 0,30 0))"),
from_wkt<L>("LINESTRING(0 1,1 0,19 0,20 1,19 1,18 0,2 0,\
1 1,2 1,3 0,17 0,18 1,17 1,16 0,4 0,3 1)"),
from_wkt<ML>("MULTILINESTRING((-1 0,-1 10),(0 0,1 0),(19 0,20 0),(25 0,30 0))"),
"mlldf05"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((-3 2,-3 0,-3 -10))"),
from_wkt<L>("LINESTRING(-3 6,-3 0,-3 5)"),
from_wkt<ML>("MULTILINESTRING((-3 0,-3 -10))"),
"mlldf06a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((6 6,-3 2,-3 0,-3 -10,9 -2))"),
from_wkt<L>("LINESTRING(-3 6,-3 0,-3 5,2 -3,-6 10,5 0,2 8,\
-6 1,10 -6)"),
from_wkt<ML>("MULTILINESTRING((6 6,-3 2),(-3 0,-3 -10,9 -2))"),
"mlldf06b"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 -3,5 4,6 6,-3 2,-3 0,-3 -10,\
9 -2,9 5,5 -5,-4 -8,9 0))"),
from_wkt<L>("LINESTRING(-3 6,-3 0,-3 5,2 -3,-6 10,5 0,2 8,\
-6 1,10 -6)"),
from_wkt<ML>("MULTILINESTRING((0 -3,5 4,6 6,-3 2),\
(-3 0,-3 -10,9 -2,9 5,5 -5,-4 -8,9 0))"),
"mlldf06c"
);
}
BOOST_AUTO_TEST_CASE( test_difference_multilinestring_multilinestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTILINESTRING / MULTILINESTRING DIFFERENCE ***"
<< std::endl;
std::cout << std::endl;
#endif
using tester = test_difference_of_geometries<ML, ML, ML>;
// disjoint linestrings
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 2,4 3),(1 1,2 2,5 3))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
"mlmldf01"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),(1 1,3 0,4 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(4 0,10 0,20 1),\
(1 0,2 0),(4 0,7 0))"),
"mlmldf02"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),(1 1,3 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(5 0,10 0,20 1),\
(1 0,2 0),(5 0,7 0))"),
"mlmldf03"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(4 0,10 0,20 1),\
(1 0,2 0),(4 0,7 0))"),
"mlmldf04"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 1),(1 0,7 0),\
(10 10,20 10,30 20))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),\
(10 20,15 10,25 10,30 15))"),
from_wkt<ML>("MULTILINESTRING((0 0,2 0),(4 0,10 0,20 1),\
(1 0,2 0),(4 0,7 0),(10 10,15 10),(20 10,30 20))"),
"mlmldf05"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 10),(1 0,7 0),\
(10 10,20 10,30 20))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),\
(-1 -1,0 0,9 0,11 10,12 10,13 3,14 4,15 5),\
(10 20,15 10,25 10,30 15))"),
from_wkt<ML>("MULTILINESTRING((9 0,10 0,13 3),(15 5,20 10),\
(10 10,11 10),(12 10,15 10),(20 10,30 20))"),
"mlmldf06"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((1 1,2 0,4 0),\
(-1 -1,0 0,9 0,11 10,12 10,13 3,14 4,15 5),\
(10 20,15 10,25 10,30 15))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,20 10),(1 0,7 0),\
(10 10,20 10,30 20))"),
from_wkt<ML>("MULTILINESTRING((1 1,2 0),(-1 -1,0 0),(9 0,11 10),\
(12 10,13 3),(10 20,15 10),(20 10,25 10,30 15))"),
"mlmldf06a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"mlmldf07"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((-1 1,0 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((-1 -1,0 0,50 0),\
(19 -1,20 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((-1 1,0 0))"),
"mlmldf07a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((-1 -1,0 0,50 0),\
(19 -1,20 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING()"),
"mlmldf07b"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 1,2 0),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"mlmldf08"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,1 0,2 0.5,3 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 1,2 0.5),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0,2 0.5,3 0))"),
"mlmldf09"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,1 0,1.5 0,2 0.5,3 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 1,1 0,2 0.5),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(1.5 0,2 0.5,3 0))"),
"mlmldf10"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,1 1,100 1,101 0),\
(0 0,101 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1,2 1,3 0,4 0,5 1,6 1,\
7 0,8 0,9 1,10 1,11 0,12 0,13 1,14 1,15 0),\
(-1 -1,1 0,101 0,200 -1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1),(2 1,5 1),(6 1,9 1),\
(10 1,13 1),(14 1,100 1,101 0),(0 0,1 0))"),
"mlmldf11"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(1 1,2 0,18 0,19 1),(2 1,3 0,17 0,18 1),\
(3 1,4 0,16 0,17 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf12"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,19 0,20 1),\
(2 0,18 0,19 1),(3 0,17 0,18 1),\
(4 0,16 0,17 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf13"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1,19 1,18 0,2 0,\
1 1,2 1,3 0,17 0,18 1,17 1,16 0,4 0,3 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf14"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 2,6 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf15"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(6 0,4 2,2 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf15a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 2,5 0,6 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf16"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,20 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(6 0,5 0,4 2,2 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,20 0))"),
"mlmldf16a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,30 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 0,5 2,20 2,25 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"mlmldf17"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,30 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,4 0,5 2,20 2,25 0,26 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"mlmldf17a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,30 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,18 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"mlmldf18"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,18 0,19 0,30 0))"),
from_wkt<ML>("MULTILINESTRING((0 1,1 0,19 0,20 1),\
(2 2,5 -1,15 2,18 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(19 0,30 0))"),
"mlmldf18a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((-1 0,0 0),(10 0,12 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
"mlmldf19"
);
{
ut_settings settings{1e-10};
tester::apply
(from_wkt<ML>("MULTILINESTRING((1 5, -4.3 -.1), (0 6, 8.6 6, 189.7654 5, 1 3, 6 3, 3 5, 6 2.232432, 0 4), (-6 5, 1 2.232432), (3 -1032.34324, 9 0, 189.7654 1, -1.4 3, 3 189.7654, +.3 10.0002, 1 5, 6 3, 5 1, 9 1, 10.0002 -1032.34324, -0.7654 0, 5 3, 3 4), (2.232432 2.232432, 8.6 +.4, 0.0 2.232432, 4 0, -8.8 10.0002), (1 0, 6 6, 7 2, -0 8.4), (-0.7654 3, +.6 8, 4 -1032.34324, 1 6, 0 4), (0 7, 2 1, 8 -7, 7 -.7, -1032.34324 9), (5 0, 10.0002 4, 8 7, 3 3, -8.1 5))"),
from_wkt<ML>("MULTILINESTRING((5 10.0002, 2 7, -0.7654 0, 5 3), (0 -0.7654, 4 10.0002, 4 +.1, -.8 3, -.1 8, 10.0002 2, +.9 -1032.34324))"),
from_wkt<ML>("MULTILINESTRING((1 5,-4.3 -0.1),(0 6,8.6 6,189.7654 5,1 3,6 3,3 5,6 2.232432,0 4),(-6 5,1 2.232432),(5 3,3 4),(3 -1032.34324,9 0,189.7654 1,-1.4 3,3 189.7654,0.3 10.0002,1 5,6 3,5 1,9 1,10.0002 -1032.34324,-0.7654 0),(2.232432 2.232432,8.6 0.4,0 2.232432,4 0,-8.8 10.0002),(1 0,6 6,7 2,-0 8.4),(-0.7654 3,0.6 8,4 -1032.34324,1 6,0 4),(0 7,2 1,8 -7,7 -0.7,-1032.34324 9),(5 0,10.0002 4,8 7,3 3,-8.1 5))"),
"mlmldf24",
settings
);
}
{
ut_settings settings{1e-10};
tester::apply
(from_wkt<ML>("MULTILINESTRING((-.4 2, 2.232432 3, 6 9, 8 189.7654, -1032.34324 5.4, 2.232432 9), (-1032.34324 3, 8 -1.6), (0 -.2, 8 1, -.5 7, 6 +.2))"),
from_wkt<ML>("MULTILINESTRING((-8 1, 4.8 6, 2 +.5), (10.0002 2,9 -1032.34324, .3 8, 0 5, 8 1, 4 .4, 2 8), (6 7, +.1 7, 0 -.5))"),
from_wkt<ML>("MULTILINESTRING((-0.4 2,2.232432 3,6 9,8 189.7654,-1032.34324 5.4,2.232432 9),(-1032.34324 3,8 -1.6),(0 -0.2,4 0.4),(8 1,-0.5 7,6 0.2))"),
"mlmldf25",
settings
);
}
{
std::string wkt1 = "MULTILINESTRING((0 -18,17 -9,8.422018348623853 -6.477064220183486),(4.0893470790378 -5.202749140893471,0 -4,6.049382716049383 -7.4567901234567895),(9.21951219512195 -9.26829268292683,14 -12),(-7 3,-5.903614457831325 6.289156626506024),(-5.277777777777778 8.166666666666668,-5 9,-5.445859872611465 8.26751592356688),(-6.342710997442455 6.794117647058823,-19 -14,-6.2338779019776425 6.668959587274291),(-5.3004291845493565 8.180257510729614,-3.579710144927536 10.966183574879228),(-3.390625 11.272321428571429,2 20),(20 -19,14 -3,3.258064516129032 1.774193548387096),(-5.275590551181104 5.5669291338582685,-13 9,-8.678414096916299 9.480176211453745),(-7.65625 9.59375,-4 10),(-2 -4,-0.15873015873015883 -0.31746031746031766),(1.875 3.75,9 18,6 8,13 17),(12 -11,6 4))";
std::string wkt2 = "MULTILINESTRING((1 16,-16 18),(17 10,15 4),(0 -14,-7 3),(-16 3,-6.471885336273427 6.942668136714442),(-4.557823129251701 7.73469387755102,-2.220994475138122 8.701657458563536),(-1.7326203208556137 8.90374331550802,13 15,15 -16,11 -11),(-6 11,-4.285714285714287 7.571428571428573),(-0.5882352941176467 0.1764705882352935,4 -9,-4 -8,-4 4.1),(-4 7.4,-4 11.666666666666666),(-4 12.142857142857142,-4 19),(-17 -1,-3.5989717223650377 3.638817480719794),(0.8699690402476783 5.185758513931889,9 8),(18 10,2.814814814814815 2.4074074074074057),(-0.9090909090909101 0.545454545454545,-4 -1,-3 -16,1.272141706924316 -1.9629629629629637),(2.666666666666667 2.6190476190476204,4 7),(-1 18,11 -6))";
std::string wkt_out = "MULTILINESTRING((0 -18,17 -9,8.4220183486238529 -6.47706422018348604),(4.08934707903780037 -5.20274914089347096,0 -4,6.04938271604938294 -7.45679012345678949),(9.21951219512195053 -9.26829268292683039,14 -12),(-7 3,-5.90361445783132499 6.28915662650602414),(-5.27777777777777768 8.16666666666666785,-5 9,-5.44585987261146531 8.26751592356687937),(-6.34271099744245515 6.79411764705882337,-19 -14,-6.23387790197764247 6.66895958727429061),(-5.30042918454935652 8.18025751072961427,-3.57971014492753614 10.9661835748792278),(-3.390625 11.2723214285714288,2 20),(20 -19,14 -3,3.25806451612903203 1.77419354838709609),(-5.27559055118110365 5.56692913385826849,-13 9,-8.67841409691629906 9.48017621145374534),(-7.65625 9.59375,-4 10),(-2 -4,-0.158730158730158832 -0.317460317460317665),(1.875 3.75,9 18,6 8,13 17),(12 -11,6 4))";
ML mls1, mls2, mls_out;
bg::read_wkt(wkt1, mls1);
bg::read_wkt(wkt2, mls2);
bg::read_wkt(wkt_out, mls_out);
tester::apply(mls1, mls2, mls_out, "mlmldf26", 1e-10);
}
{
std::string wkt1 = "MULTILINESTRING((12 -11,6 4))";
std::string wkt2 = "MULTILINESTRING((18 10,2.814814814814815 2.4074074074074057),(-1 18,11 -6))";
std::string wkt_out = "MULTILINESTRING((12 -11,6 4))";
ML mls1, mls2, mls_out;
bg::read_wkt(wkt1, mls1);
bg::read_wkt(wkt2, mls2);
bg::read_wkt(wkt_out, mls_out);
tester::apply(mls1, mls2, mls_out, "mlmldf27");
}
}
#ifndef BOOST_GEOMETRY_TEST_NO_DEGENERATE
BOOST_AUTO_TEST_CASE( test_difference_ml_ml_degenerate )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTILINESTRING / MULTILINESTRING DIFFERENCE"
<< " (DEGENERATE) ***"
<< std::endl;
std::cout << std::endl;
#endif
using tester = test_difference_of_geometries<ML, ML, ML>;
// the following test cases concern linestrings with duplicate
// points and possibly linestrings with zero length.
tester::apply
(from_wkt<ML>("MULTILINESTRING((5 5,5 5),(0 0,18 0,18 0,\
19 0,19 0,19 0,30 0),(2 0,2 0),(4 10,4 10))"),
from_wkt<ML>("MULTILINESTRING((-10 0,-9 0),(0 10,5 0,20 0,20 0,30 10),\
(1 1,2 2),(1 10,1 10,1 0,1 0,1 -10),\
(2 0,2 0),(3 0,3 0,3 0),(0 0,0 10),\
(4 0,4 10),(5 5,5 5))"),
from_wkt<ML>("MULTILINESTRING((0 0,5 0),(20 0,30 0))"),
"mlmldf20a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((5 5,5 5),(0 0,0 0,18 0,18 0,\
19 0,19 0,19 0,30 0,30 0),(2 0,2 0),(4 10,4 10))"),
from_wkt<ML>("MULTILINESTRING((-10 0,-9 0),(0 10,5 0,20 0,20 0,30 10),\
(1 1,1 1,2 2,2 2),(1 10,1 10,1 0,1 0,1 -10),\
(2 0,2 0),(3 0,3 0,3 0),(0 0,0 0,0 10,0 10),\
(4 0,4 10,4 10),(5 5,5 5))"),
from_wkt<ML>("MULTILINESTRING((0 0,5 0),(20 0,30 0))"),
"mlmldf20aa"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((10 10,10 10),(0 0,0 0,18 0,18 0,\
19 0,19 0,19 0,30 0,30 0),(2 0,2 0),(4 10,4 10))"),
from_wkt<ML>("MULTILINESTRING((-10 0,-9 0),(0 10,5 0,20 0,20 0,30 10),\
(1 1,1 1,2 2,2 2),(1 10,1 10,1 0,1 0,1 -10),\
(2 0,2 0),(3 0,3 0,3 0),(0 0,0 0,0 10,0 10),\
(4 0,4 10,4 10),(5 5,5 5))"),
from_wkt<ML>("MULTILINESTRING((10 10,10 10),(0 0,5 0),(20 0,30 0))"),
"mlmldf20aaa"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((10 10),(0 0,0 0,18 0,18 0,\
19 0,19 0,19 0,30 0,30 0),(2 0,2 0),(4 10,4 10))"),
from_wkt<ML>("MULTILINESTRING((-10 0,-9 0),(0 10,5 0,20 0,20 0,30 10),\
(1 1,1 1,2 2,2 2),(1 10,1 10,1 0,1 0,1 -10),\
(2 0,2 0),(3 0,3 0,3 0),(0 0,0 0,0 10,0 10),\
(4 0,4 10,4 10),(5 5,5 5))"),
from_wkt<ML>("MULTILINESTRING((10 10,10 10),(0 0,5 0),(20 0,30 0))"),
"mlmldf20aaaa"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,0 0),(1 1,1 1))"),
from_wkt<ML>("MULTILINESTRING((1 1,1 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,0 0))"),
"mlmldf21"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,0 0),(2 2,2 2),(1 1,1 1))"),
from_wkt<ML>("MULTILINESTRING((1 1,1 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,0 0),(2 2,2 2))"),
"mlmldf22"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,0 0),(1 1,1 1),(2 2,2 2))"),
from_wkt<ML>("MULTILINESTRING((1 1,1 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,0 0),(2 2,2 2))"),
"mlmldf23"
);
}
#endif // BOOST_GEOMETRY_TEST_NO_DEGENERATE
BOOST_AUTO_TEST_CASE( test_difference_ml_ml_spikes )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTILINESTRING / MULTILINESTRING DIFFERENCE"
<< " (WITH SPIKES) ***"
<< std::endl;
std::cout << std::endl;
#endif
using tester = test_difference_of_geometries<ML, ML, ML>;
// the following test cases concern linestrings with spikes
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,9 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0))"),
"mlmldf-spikes-01"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((9 0,1 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0))"),
"mlmldf-spikes-02"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,9 0,2 0,8 0,3 0,7 0,4 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0))"),
"mlmldf-spikes-03"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,3 0,2 0,4 0,3 0,5 0,4 0,6 0,\
5 0,7 0,6 0,8 0,7 0,9 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0))"),
"mlmldf-spikes-04"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,6 0,5 0),(7 0,8 0,7 0),\
(9 1,9 0,9 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(6 0,7 0),(8 0,10 0))"),
"mlmldf-spikes-05"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,6 0,5 0),(7 0,8 0,7 0),\
(9 0,9 2,9 1))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(6 0,7 0),(8 0,10 0))"),
"mlmldf-spikes-05a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,6 0,5 0),(9 0,6 0,8 0),\
(11 0,8 0,12 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0))"),
"mlmldf-spikes-06"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((-1 0,0 0,-2 0),(11 0,10 0,12 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
"mlmldf-spikes-07"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((-1 -1,0 0,-2 -2),(11 1,10 0,12 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
"mlmldf-spikes-07a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,6 0,5 0),(11 0,10 0,12 0),\
(7 5,7 0,8 0,6.5 0,8.5 0,8.5 5))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(6 0,6.5 0),(8.5 0,10 0))"),
"mlmldf-spikes-08"
);
// now the first geometry has a spike
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,7 0,4 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,8 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(8 0,10 0))"),
"mlmldf-spikes-09"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,7 0,4 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,7 0,4 0,9 0))"),
"mlmldf-spikes-09a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,7 0,4 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,5 0),(9 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((5 0,7 0,5 0),(5 0,9 0))"),
"mlmldf-spikes-09b"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,7 0,4 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,5 0),(6 0,10 0))"),
from_wkt<ML>("MULTILINESTRING((5 0,6 0),(6 0,5 0),(5 0,6 0))"),
"mlmldf-spikes-09c"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,8 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(8 0,10 0,8 0))"),
"mlmldf-spikes-10"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,8 0,4 0),(2 0,9 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0,9 0))"),
"mlmldf-spikes-11"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((11 1,10 0,12 2))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
"mlmldf-spikes-12"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((11 -1,10 0,12 -2))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
"mlmldf-spikes-12a"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
from_wkt<ML>("MULTILINESTRING((11 0,10 0,12 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,10 0,5 0))"),
"mlmldf-spikes-13"
);
// the following three tests have been discussed with Adam
tester::apply
(from_wkt<ML>("MULTILINESTRING((1 0,1 1,2 1))"),
from_wkt<ML>("MULTILINESTRING((1 2,1 1,1 2))"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1,2 1))"),
"mlmldf-spikes-14"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,1 0,0 0))"),
from_wkt<ML>("MULTILINESTRING((2 0,1 0,2 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0,0 0))"),
"mlmldf-spikes-15"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((1 0,1 1,2 1))"),
from_wkt<ML>("MULTILINESTRING((2 0,1 1,2 0))"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1,2 1))"),
"mlmldf-spikes-16"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((1 0,1 1,2 1))"),
from_wkt<ML>("MULTILINESTRING((2 1,1 1,2 1))"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1))"),
"mlmldf-spikes-17"
);
// test cases sent by Adam on the mailing list (equal spikes)
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,1 1,0 0))"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,0 0))"),
from_wkt<ML>("MULTILINESTRING()"),
"mlmldf-spikes-18"
);
tester::apply
(from_wkt<ML>("MULTILINESTRING((0 0,1 1,0 0))"),
from_wkt<ML>("MULTILINESTRING((1 1,0 0,1 1))"),
from_wkt<ML>("MULTILINESTRING()"),
"mlmldf-spikes-19"
);
}
BOOST_AUTO_TEST_CASE( test_difference_ls_mls_geo_rad )
{
using pt = bg::model::point<double, 2, bg::cs::geographic<bg::radian>>;
using ls = bg::model::linestring<pt>;
using mls = bg::model::multi_linestring<ls>;
bg::srs::spheroid<double> sph_wgs84(6378137.0, 6356752.3142451793);
boost::geometry::strategy::intersection::geographic_segments<> wgs84(sph_wgs84);
ls g1 = from_wkt<ls>("LINESTRING(0 0,-3.14159265358979 0)");
mls g2 = from_wkt<mls>("MULTILINESTRING((-2.1467549799530232 -0.12217304763960295,"
"-2.5481807079117185 -0.90757121103705041,"
"-2.6529004630313784 0.85521133347722067,"
" 0.92502450355699373 0.62831853071795796,"
"-2.5307274153917754 0,"
" 2.8099800957108676 1.0646508437165401,"
"-1.6057029118347816 -1.5009831567151219,"
" 0.2268928027592626 1.0646508437165401,"
"-2.199114857512853 -0.017453292519943278,"
" 0 0.31415926535897898,"
" 0 0.57595865315812822,"
" 1.0471975511965967 -0.73303828583761765,"
" 2.1118483949131366 -0.54105206811824158))");
mls out;
bg::difference(g1, g2, out, wgs84);
check_result(g1, g2, out,
from_wkt<mls>("MULTILINESTRING((0 0,0 0.31415926535897897853),"
"(0 0.57595865315812821983,-3.1415926535897900074 0))"),
"geo_lmldf-1");
}
Reference in New Issue View Git Blame Copy Permalink