void test_spherical()
{
// Test spherical boxes
// See also http://www.gcmap.com/mapui?P=1E45N-19E45N-19E55N-1E55N-1E45N,10E55.1N,10E45.1N
bg::model::box<Point> box;
bg::read_wkt("POLYGON((1 45,19 55))", box);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.1), box), true);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.2), box), true);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.3), box), true);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.4), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.1), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.2), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.3), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.4), box), true);
// Crossing the dateline (Near Tuvalu)
// http://www.gcmap.com/mapui?P=178E10S-178W10S-178W6S-178E6S-178E10S,180W5.999S,180E9.999S
// http://en.wikipedia.org/wiki/Tuvalu
bg::model::box<Point> tuvalu(Point(178, -10), Point(-178, -6));
BOOST_CHECK_EQUAL(bg::within(Point(180, -8), tuvalu), true);
BOOST_CHECK_EQUAL(bg::within(Point(-180, -8), tuvalu), true);
BOOST_CHECK_EQUAL(bg::within(Point(180, -5.999), tuvalu), false);
BOOST_CHECK_EQUAL(bg::within(Point(180, -10.001), tuvalu), true);
}
void test_spherical()
{
typedef typename bg::coordinate_type<Point>::type ct;
bg::model::polygon<Point> wrangel;
// SQL Server check (no geography::STWithin, so check with intersection trick)
/*
with q as (
select geography::STGeomFromText('POLYGON((-178.569 71.5641,-179.034 71.5977,-179.305 71.5514,-179.629 71.5772,-180 71.5358,179.53 71.4383,178.872 71.2175,178.618 71.0355,178.791 70.7964,179.273 70.8886,179.678 70.8955,-180 70.9972,-179.274 70.9078,-178.819 70.98,-177.939 71.0375,-177.62 71.1166,-177.439 71.2269,-177.503 71.2775,-177.833 71.3461,-178.018 71.4497,-178.569 71.5641))',4326) as wrangel
)
select wrangel.STArea()/1000000.0
,geography::STGeomFromText('POINT(-179.3 71.27)',4326).STIntersection(wrangel).STAsText() as workaround_within_1
,geography::STGeomFromText('POINT(-179.9 70.95)',4326).STIntersection(wrangel).STAsText() as workaround_within_2
,geography::STGeomFromText('POINT(179.9 70.95)',4326).STIntersection(wrangel).STAsText() as workaround_within_3
from q
-> 7669.10402181435 POINT (-179.3 71.27) GEOMETRYCOLLECTION EMPTY GEOMETRYCOLLECTION EMPTY
PostGIS knows Within for Geography neither, and the intersection trick gives the same result
*/
bg::read_wkt("POLYGON((-178.568604 71.564148,-178.017548 71.449692,-177.833313 71.3461,-177.502838 71.277466 ,-177.439453 71.226929,-177.620026 71.116638,-177.9389 71.037491,-178.8186 70.979965,-179.274445 70.907761,-180 70.9972,179.678314 70.895538,179.272766 70.888596,178.791016 70.7964,178.617737 71.035538,178.872192 71.217484,179.530273 71.4383 ,-180 71.535843 ,-179.628601 71.577194,-179.305298 71.551361,-179.03421 71.597748,-178.568604 71.564148))", wrangel);
bool within = bg::within(Point(-179.3, 71.27), wrangel);
BOOST_CHECK_EQUAL(within, true);
within = bg::within(Point(-179.9, 70.95), wrangel);
BOOST_CHECK_EQUAL(within, false);
within = bg::within(Point(179.9, 70.95), wrangel);
BOOST_CHECK_EQUAL(within, false);
// Test using great circle mapper
// http://www.gcmap.com/mapui?P=5E52N-9E53N-7E50N-5E52N,7E52.5N,8E51.5N,6E51N
bg::model::polygon<Point> triangle;
bg::read_wkt("POLYGON((5 52,9 53,7 50,5 52))", triangle);
BOOST_CHECK_EQUAL(bg::within(Point(7, 52.5), triangle), true);
BOOST_CHECK_EQUAL(bg::within(Point(8.0, 51.5), triangle), false);
BOOST_CHECK_EQUAL(bg::within(Point(6.0, 51.0), triangle), false);
// Test spherical boxes
// See also http://www.gcmap.com/mapui?P=1E45N-19E45N-19E55N-1E55N-1E45N,10E55.1N,10E45.1N
bg::model::box<Point> box;
bg::read_wkt("POLYGON((1 45,19 55))", box);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.1), box), true);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.2), box), true);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.3), box), true);
BOOST_CHECK_EQUAL(bg::within(Point(10, 55.4), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.1), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.2), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.3), box), false);
BOOST_CHECK_EQUAL(bg::within(Point(10, 45.4), box), true);
// Crossing the dateline (Near Tuvalu)
// http://www.gcmap.com/mapui?P=178E10S-178W10S-178W6S-178E6S-178E10S,180W5.999S,180E9.999S
// http://en.wikipedia.org/wiki/Tuvalu
bg::model::box<Point> tuvalu(Point(178, -10), Point(-178, -6));
BOOST_CHECK_EQUAL(bg::within(Point(180, -8), tuvalu), true);
BOOST_CHECK_EQUAL(bg::within(Point(-180, -8), tuvalu), true);
BOOST_CHECK_EQUAL(bg::within(Point(180, -5.999), tuvalu), false);
BOOST_CHECK_EQUAL(bg::within(Point(180, -10.001), tuvalu), true);
}
