Files
geometry/test/algorithms/set_operations/intersection/intersection_pl_l.cpp
2015-03-10 16:50:38 +02:00

703 lines
19 KiB
C++

// Boost.Geometry (aka GGL, Generic Geometry Library)
// Copyright (c) 2015, 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 Menelaos Karavelas, on behalf of Oracle
#ifndef BOOST_TEST_MODULE
#define BOOST_TEST_MODULE test_intersection_pointlike_linear
#endif
#include <iostream>
#include <algorithm>
#include <boost/test/included/unit_test.hpp>
#include "../test_set_ops_pointlike.hpp"
#include <boost/geometry/geometries/point.hpp>
#include <boost/geometry/geometries/linestring.hpp>
#include <boost/geometry/geometries/segment.hpp>
#include <boost/geometry/geometries/multi_point.hpp>
#include <boost/geometry/geometries/multi_linestring.hpp>
typedef bg::model::point<double, 2, bg::cs::cartesian> point_type;
typedef bg::model::segment<point_type> segment_type;
typedef bg::model::linestring<point_type> linestring_type;
typedef bg::model::multi_point<point_type> multi_point_type;
typedef bg::model::multi_linestring<linestring_type> multi_linestring_type;
//===========================================================================
//===========================================================================
//===========================================================================
BOOST_AUTO_TEST_CASE( test_intersection_point_segment )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** POINT / SEGMENT INTERSECTION ***" << std::endl;
std::cout << std::endl;
#endif
typedef point_type P;
typedef segment_type S;
typedef multi_point_type MP;
typedef test_set_op_of_pointlike_geometries
<
P, S, MP, bg::overlay_intersection
> tester;
tester::apply
("psi01",
from_wkt<P>("POINT(0 0)"),
from_wkt<S>("SEGMENT(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("psi02",
from_wkt<P>("POINT(0 0)"),
from_wkt<S>("SEGMENT(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("psi03",
from_wkt<P>("POINT(1 1)"),
from_wkt<S>("SEGMENT(0 0,2 2)"),
from_wkt<MP>("MULTIPOINT(1 1)")
);
tester::apply
("psi04",
from_wkt<P>("POINT(3 3)"),
from_wkt<S>("SEGMENT(0 0,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
}
BOOST_AUTO_TEST_CASE( test_intersection_point_linestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** POINT / LINESTRING INTERSECTION ***" << std::endl;
std::cout << std::endl;
#endif
typedef point_type P;
typedef linestring_type L;
typedef multi_point_type MP;
typedef test_set_op_of_pointlike_geometries
<
P, L, MP, bg::overlay_intersection
> tester;
tester::apply
("pli01",
from_wkt<P>("POINT(0 0)"),
from_wkt<L>("LINESTRING(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pli02",
from_wkt<P>("POINT(0 0)"),
from_wkt<L>("LINESTRING(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("pli03",
from_wkt<P>("POINT(1 1)"),
from_wkt<L>("LINESTRING(0 0,2 2)"),
from_wkt<MP>("MULTIPOINT(1 1)")
);
tester::apply
("pli04",
from_wkt<P>("POINT(3 3)"),
from_wkt<L>("LINESTRING(0 0,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
// linestrings with more than two points
tester::apply
("pli05",
from_wkt<P>("POINT(1 1)"),
from_wkt<L>("LINESTRING(0 0,1 1,2 2)"),
from_wkt<MP>("MULTIPOINT(1 1)")
);
tester::apply
("pli06",
from_wkt<P>("POINT(2 2)"),
from_wkt<L>("LINESTRING(0 0,1 1,4 4)"),
from_wkt<MP>("MULTIPOINT(2 2)")
);
tester::apply
("pli07",
from_wkt<P>("POINT(10 10)"),
from_wkt<L>("LINESTRING(0 0,1 1,4 4)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pli08",
from_wkt<P>("POINT(0 1)"),
from_wkt<L>("LINESTRING(0 0,1 1,4 4)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pli09",
from_wkt<P>("POINT(4 4)"),
from_wkt<L>("LINESTRING(0 0,1 1,4 4)"),
from_wkt<MP>("MULTIPOINT(4 4)")
);
tester::apply
("pli10",
from_wkt<P>("POINT(0 0)"),
from_wkt<L>("LINESTRING(0 0,1 1,4 4)"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
}
BOOST_AUTO_TEST_CASE( test_intersection_point_multilinestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** POINT / MULTILINESTRING INTERSECTION ***" << std::endl;
std::cout << std::endl;
#endif
typedef point_type P;
typedef multi_linestring_type ML;
typedef multi_point_type MP;
typedef test_set_op_of_pointlike_geometries
<
P, ML, MP, bg::overlay_intersection
> tester;
tester::apply
("pmli01",
from_wkt<P>("POINT(0 0)"),
from_wkt<ML>("MULTILINESTRING((1 1,2 2))"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pmli02",
from_wkt<P>("POINT(0 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1))"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("pmli03",
from_wkt<P>("POINT(1 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 2))"),
from_wkt<MP>("MULTIPOINT(1 1)")
);
tester::apply
("pmli04",
from_wkt<P>("POINT(3 3)"),
from_wkt<ML>("MULTILINESTRING((0 0,2 2))"),
from_wkt<MP>("MULTIPOINT()")
);
// linestrings with more than two points
tester::apply
("pmli05",
from_wkt<P>("POINT(1 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,2 2))"),
from_wkt<MP>("MULTIPOINT(1 1)")
);
tester::apply
("pmli06",
from_wkt<P>("POINT(2 2)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(2 2)")
);
tester::apply
("pmli07",
from_wkt<P>("POINT(10 10)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pmli08",
from_wkt<P>("POINT(0 1)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pmli09",
from_wkt<P>("POINT(4 4)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(4 4)")
);
tester::apply
("pmli10",
from_wkt<P>("POINT(0 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
// multilinestrings with more than one linestring
tester::apply
("pmli11",
from_wkt<P>("POINT(0 0)"),
from_wkt<ML>("MULTILINESTRING((-10,-10),(0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("pmli12",
from_wkt<P>("POINT(0 0)"),
from_wkt<ML>("MULTILINESTRING((-10 0,0 0,10 0),(0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("pmli13",
from_wkt<P>("POINT(0 0)"),
from_wkt<ML>("MULTILINESTRING((-10 0,10 0),(0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("pmli14",
from_wkt<P>("POINT(-20 0)"),
from_wkt<ML>("MULTILINESTRING((-10 0,10 0),(0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pmli15",
from_wkt<P>("POINT(0 1)"),
from_wkt<ML>("MULTILINESTRING((-10 0,10 0),(0 0,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("pmli16",
from_wkt<P>("POINT(1 0)"),
from_wkt<ML>("MULTILINESTRING((-10 0,10 0),(-1 0,2 0,20 0))"),
from_wkt<MP>("MULTIPOINT(1 0)")
);
tester::apply
("pmli17",
from_wkt<P>("POINT(1 0)"),
from_wkt<ML>("MULTILINESTRING((-10 0,10 0),(0 -1,0 2,0 4))"),
from_wkt<MP>("MULTIPOINT(1 0)")
);
}
BOOST_AUTO_TEST_CASE( test_intersection_multipoint_segment )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTIPOINT / SEGMENT INTERSECTION ***" << std::endl;
std::cout << std::endl;
#endif
typedef multi_point_type MP;
typedef segment_type S;
typedef test_set_op_of_pointlike_geometries
<
MP, S, MP, bg::overlay_intersection
> tester;
tester::apply
("mpsi01",
from_wkt<MP>("MULTIPOINT(0 0)"),
from_wkt<S>("SEGMENT(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi02",
from_wkt<MP>("MULTIPOINT(0 0)"),
from_wkt<S>("SEGMENT(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("mpsi03",
from_wkt<MP>("MULTIPOINT(0 0,0 0)"),
from_wkt<S>("SEGMENT(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi04",
from_wkt<MP>("MULTIPOINT(0 0,0 0)"),
from_wkt<S>("SEGMENT(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpsi05",
from_wkt<MP>("MULTIPOINT(0 0,0 0,1 0)"),
from_wkt<S>("SEGMENT(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsf06",
from_wkt<MP>("MULTIPOINT(0 0,0 0,1 0)"),
from_wkt<S>("SEGMENT(1 0,2 0)"),
from_wkt<MP>("MULTIPOINT(1 0)")
);
tester::apply
("mpsi07",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(0 0,1 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpsi07a",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(0 0,3 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)")
);
tester::apply
("mpsi08",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<S>("SEGMENT(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi09",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(-1 0,1 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpsi10",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(1 0,3 0)"),
from_wkt<MP>("MULTIPOINT(2 0)")
);
tester::apply
("mpsi11",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(-1 0,3 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)")
);
tester::apply
("mpsi12",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(3 0,3 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi12a",
from_wkt<MP>("MULTIPOINT(0 0,2 0,0 0)"),
from_wkt<S>("SEGMENT(3 0,3 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi12b",
from_wkt<MP>("MULTIPOINT(0 0,2 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(3 0,3 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi12c",
from_wkt<MP>("MULTIPOINT(0 0,2 0,0 0,2 0,0 0)"),
from_wkt<S>("SEGMENT(3 0,3 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi13",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(2 0,2 0)"),
from_wkt<MP>("MULTIPOINT(2 0)")
);
tester::apply
("mpsi14",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<S>("SEGMENT(0 0,0 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpsi15",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<S>("SEGMENT(0 0,1 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpsi16",
from_wkt<MP>("MULTIPOINT(0 0,2 0,0 0,2 0,0 0)"),
from_wkt<S>("SEGMENT(-1 0,10 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0,0 0,2 0,2 0)")
);
}
BOOST_AUTO_TEST_CASE( test_intersection_multipoint_linestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTIPOINT / LINESTRING INTERSECTION ***" << std::endl;
std::cout << std::endl;
#endif
typedef multi_point_type MP;
typedef linestring_type L;
typedef test_set_op_of_pointlike_geometries
<
MP, L, MP, bg::overlay_intersection
> tester;
tester::apply
("mpli01",
from_wkt<MP>("MULTIPOINT(0 0)"),
from_wkt<L>("LINESTRING(1 1,2 2,3 3,4 4,5 5)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpli02",
from_wkt<MP>("MULTIPOINT(0 0)"),
from_wkt<L>("LINESTRING(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT(0 0)")
);
tester::apply
("mpli03",
from_wkt<MP>("MULTIPOINT(0 0,0 0)"),
from_wkt<L>("LINESTRING(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpli04",
from_wkt<MP>("MULTIPOINT(0 0,0 0)"),
from_wkt<L>("LINESTRING(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpli05",
from_wkt<MP>("MULTIPOINT(0 0,0 0,1 0)"),
from_wkt<L>("LINESTRING(1 1,2 2)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpli06",
from_wkt<MP>("MULTIPOINT(0 0,0 0,1 0)"),
from_wkt<L>("LINESTRING(1 0,2 0)"),
from_wkt<MP>("MULTIPOINT(1 0)")
);
tester::apply
("mpli07",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(0 0,1 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpli08",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<L>("LINESTRING(0 0,1 1)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpli09",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(-1 0,1 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpli10",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(1 0,3 0)"),
from_wkt<MP>("MULTIPOINT(2 0)")
);
tester::apply
("mpli10a",
from_wkt<MP>("MULTIPOINT(2 0,0 0,2 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(1 0,3 0)"),
from_wkt<MP>("MULTIPOINT(2 0,2 0,2 0)")
);
tester::apply
("mpli11",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(-1 0,3 0)"),
from_wkt<MP>("MULTIPOINT(2 0,0 0,0 0)")
);
tester::apply
("mpli12",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(3 0,3 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpli13",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(2 0,2 0)"),
from_wkt<MP>("MULTIPOINT(2 0)")
);
tester::apply
("mpli14",
from_wkt<MP>("MULTIPOINT(0 0,0 0,2 0)"),
from_wkt<L>("LINESTRING(0 0,0 0)"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
tester::apply
("mpli15",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<L>("LINESTRING(0 0,1 0,2 0,3 0)"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpli16",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<L>("LINESTRING()"),
from_wkt<MP>("MULTIPOINT()")
);
}
BOOST_AUTO_TEST_CASE( test_intersection_multipoint_multilinestring )
{
#ifdef BOOST_GEOMETRY_TEST_DEBUG
std::cout << std::endl << std::endl << std::endl;
std::cout << "*** MULTIPOINT / MULTILINESTRING INTERSECTION ***"
<< std::endl;
std::cout << std::endl;
#endif
typedef multi_point_type MP;
typedef multi_linestring_type ML;
typedef test_set_op_of_pointlike_geometries
<
MP, ML, MP, bg::overlay_intersection
> tester;
tester::apply
("mpmli01",
from_wkt<MP>("MULTIPOINT(0 0,1 0,2 0)"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(1 0)")
);
tester::apply
("mpmli02",
from_wkt<MP>("MULTIPOINT(2 2,3 3,0 0,0 0,2 2,1 1,1 1,1 0,1 0)"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(2 2,2 2,3 3,1 1,1 1,1 0,1 0)")
);
tester::apply
("mpmli03",
from_wkt<MP>("MULTIPOINT(5 5,3 3,0 0,0 0,5 5,1 1,1 1,1 0,1 0)"),
from_wkt<ML>("MULTILINESTRING((1 0,1 1,1 1,4 4))"),
from_wkt<MP>("MULTIPOINT(1 0,1 0,3 3,1 1,1 1)")
);
tester::apply
("mpmli04",
from_wkt<MP>("MULTIPOINT(0 0,1 1,1 0,1 1)"),
from_wkt<ML>("MULTILINESTRING((1 0,0 0,1 1,0 0))"),
from_wkt<MP>("MULTIPOINT(1 0,0 0,1 1,1 1)")
);
tester::apply
("mpmli05",
from_wkt<MP>("MULTIPOINT(0 0,1 0,2 0,5 0)"),
from_wkt<ML>("MULTILINESTRING((0 1,0 0,1 2,1 0),\
(0 1,2 0,0 10,2 0,2 10,5 10,5 0,5 10))"),
from_wkt<MP>("MULTIPOINT(0 0,1 0,2 0,5 0)")
);
tester::apply
("mpmli06",
from_wkt<MP>("MULTIPOINT(0 0,1 1,1 0,1 1)"),
from_wkt<ML>("MULTILINESTRING(())"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpmli07",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<ML>("MULTILINESTRING(())"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpmli08",
from_wkt<MP>("MULTIPOINT()"),
from_wkt<ML>("MULTILINESTRING((0 0,1 0),(9 0,10 0))"),
from_wkt<MP>("MULTIPOINT()")
);
tester::apply
("mpmli09",
from_wkt<MP>("MULTIPOINT(0 0,1 1,0 0,0 0,0 0,0 0,0 0,0 0,\
0 0,0 0,0 0,0 0,0 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,0 1),(0 0,2 0),(0 0,0 3),\
(0 0,4 0),(0 0,0 5),(0 0,6 0),(0 0,7 0),(0 0,8 0))"),
from_wkt<MP>("MULTIPOINT(0 0,0 0,0 0,0 0,0 0,0 0,0 0,\
0 0,0 0,0 0,0 0,0 0)")
);
tester::apply
("mpmli09a",
from_wkt<MP>("MULTIPOINT(0 0,1 1,0 0)"),
from_wkt<ML>("MULTILINESTRING((0 0,0 1),(0 0,2 0),(0 0,0 3),\
(0 0,4 0),(0 0,0 5),(0 0,6 0),(0 0,7 0),(0 0,8 0))"),
from_wkt<MP>("MULTIPOINT(0 0,0 0)")
);
}