| 1 |
|
|---|
| 2 | #include <iostream>
|
|---|
| 3 | #include <boost/geometry.hpp>
|
|---|
| 4 | #include <boost/geometry/io/wkt/wkt.hpp>
|
|---|
| 5 |
|
|---|
| 6 |
|
|---|
| 7 |
|
|---|
| 8 | namespace bg = boost::geometry;
|
|---|
| 9 |
|
|---|
| 10 | void intersectsTest(size_t maxtimes) {
|
|---|
| 11 | typedef bg::model::point<double, 2, bg::cs::spherical_equatorial<bg::degree> > SEP;
|
|---|
| 12 | typedef boost::geometry::model::polygon<SEP> polygon_type;
|
|---|
| 13 | polygon_type testpl;
|
|---|
| 14 | SEP testp;
|
|---|
| 15 | bool r = false;
|
|---|
| 16 | /// wkt for spherical equatorial, format is "longtitude latitude"
|
|---|
| 17 | boost::geometry::read_wkt("POLYGON((30 0,30 30,90 30, 90 0, 30 0))", testpl);
|
|---|
| 18 | boost::geometry::assign_values(testp, -120,21);
|
|---|
| 19 | std::cout << "Result 1: " << boost::geometry::intersects(testpl, testp) << std::endl;
|
|---|
| 20 | boost::geometry::read_wkt("POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))", testpl);
|
|---|
| 21 | boost::geometry::assign_values(testp, 0,90);
|
|---|
| 22 | std::cout << "Result 2: " << boost::geometry::intersects(testpl, testp) << std::endl;
|
|---|
| 23 | }
|
|---|
