Geometry Boost List: Re: [geometry] intersects returns incorrect result in spherical equatorial

Subject: Re: [geometry] intersects returns incorrect result in spherical equatorial
From: Barend Gehrels (barend_at_[hidden])
Date: 2013-09-24 06:02:23

Hi Terry,

Op 24 sep. 2013 om 10:43 heeft Terry Hole <hujie.eagle_at_[hidden]> 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