
Geometry : 
Subject: Re: [geometry] intersects returns incorrect result in spherical equatorial
From: Barend Gehrels (barend_at_[hidden])
Date: 20130924 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
Geometry list run by mateusz at loskot.net