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

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 #include #include

void test_point_on_surface() { typedef boost::geometry::model::point 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@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users