Re: [Boost-users] How to use geometry::within with a polygon that encompasses a pole

8 Jan 2013


      On Mon, Jan 7, 2013 at 5:20 PM, Cameron Sabol  wrote:
...
Hi,
I'm have an issue with polygons from the geometry library using the spherical_equatorial coordinates when the polygon is meant to wrap a pole.
e.g.
namespace bg = boost::geometry;
typedef bg::model::point
    <
        double, 2, bg::cs::spherical_equatorialbg::degree
    > LonLatPoint;
bg::model::polygon < LonLatPoint> northPoleBoundary;
double lon = -180;
while (lon <= 180)
{
    bg::append(northPoleBoundary, LonLatPoint(lon, 88));
    lon += 5;
}
LonLatPoint northPole(90, 0);
bool isInNorthPole = bg::within(northPole, northPoleBoundary);
This code is giving me isInNorthPole == false. Am I constructing my polygon incorrectly? Is this a bug? Or is this behavior simply not implemented? Thanks for your help!
I think you meant:

LonLatPoint northPole(0, 90);

--Mike

Re: [Boost-users] How to use geometry::within with a polygon that encompasses a pole

Michael Fawcett