Hi Terry,


Op 24 sep. 2013 om 10:43 heeft Terry Hole <hujie.eagle@gmail.com> het volgende geschreven:

> Hi,
>
> The intersects does not return the correct result for point in polygon with spherical equatorial.
> These are two cases:
> 1. a polygon with pole (north/south) in it
>     assume polygon PL: POLYGON((0 80,-90 80, -180 80, 90 80, 0 80))
>     a point P: POINT(0 90)
>     boost::geometry::intersects(P, PL) should return TRUE, but return false in boost 1.54.
> 2. point in the opposite side of the polygon
>     assume polygon PL: POLYGON((30 0,30 30,90 30, 90 0, 30 0))
>     a point P: POINT(-120 21)
>     boost::geometry::intersects(P, PL) should return FALSE, but return true in boost 1.54.
>
> The issue may be caused by the spherical_side_formula can not handle this two cases.

Yes, polygons over the pole are not handled correctly. I kind of knew this, it has to be fixed once. Can you create  Boost Ticket with this program? If not, I will do it.


Thanks for your report,
Barend


_______________________________________________
Geometry mailing list
Geometry@lists.boost.org
http://lists.boost.org/mailman/listinfo.cgi/geometry