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@lists.boost.org http://lists.boost.org/mailman/listinfo.cgi/boost-users