Re: [Boost-users] Boost.Geometry, probable bug in point_on_surface

Subject: Re: [Boost-users] Boost.Geometry, probable bug in point_on_surface
From: Adam Wulkiewicz (adam.wulkiewicz_at_[hidden])
Date: 2014-10-09 10:44:04

Hi Petr,

I'll top-post because there is only code below.

You're right, something's wrong with point_on_surface().
Could you please create a ticket: 
http://www.boost.org/development/bugs.html ?

FYI, your WKT POLYGON is CCW but you're creating a CW one using default 
template parameters of model::polygon<>.

Regards,
Adam

Petr Doubrava wrote:
> #include <boost/geometry.hpp>
> #include <boost/geometry/geometries/point_xy.hpp>
> #include <boost/geometry/geometries/polygon.hpp>
>
>    void test_point_on_surface()
>    {
>        typedef boost::geometry::model::point<double, 2, boost::geometry::cs::cartesian> Point;
>        boost::geometry::model::polygon<Point> poly;
>
>        boost::geometry::read_wkt("POLYGON((1074699.93 703064.65, 1074702.10 703054.62, 1074704.53 703061.40, 1074703.90 703064.58, 1074699.93 703064.65))",poly);
>        boost::geometry::correct(poly);
>
>        Point point;
>        boost::geometry::point_on_surface(poly,point);
>
>        //computed point isn't on the surface of polygon
>        //[0]    1074839.1013086310 double
>        //[1]    703154.77900308778 double
>
>        bool res = boost::geometry::within(point,poly);
>        // res is false
>    }
> _______________________________________________
> Boost-users mailing list
> Boost-users_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/boost-users
>