petr@… — Fri, 10 Oct 2014 08:27:41 GMT

I do not have permission to add attachments. Code is below in this comment.

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>
#include <assert.h>
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);
    assert(res==true);
    // res is false :-(
}
int main(int argc, char* argv[])
{
  test_point_on_surface();
  return 0;
}

status, milestone changed; resolution set

awulkiew — Fri, 10 Oct 2014 18:49:04 GMT

status new → closed
resolution → fixed
milestone To Be Determined → Boost 1.57.0

Boost C++ Libraries: Ticket #10643: point_on_surface bug, computed point isn't on the surface of polygon

status, milestone changed; resolution set