Boost::geometry::covered_by错误的结果

boost::geometry::covered_by wrong results

本文关键字:错误 结果 by geometry covered Boost      更新时间:2023-10-16

我有一个boost::geometry::covered_by方法的问题,以确定点是否在多边形内。在下面的例子中,我有一个点正好在多边形之外。点距离多边形的x坐标(0.1377)。多边形的第三个点的y坐标为4.999999999999999982。这就引出了一个问题。当我将4.999999999999999982更改为5.00时,一切正常。该怎么办?

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>
#include <boost/geometry/io/wkt/wkt.hpp>

int main()
{
typedef boost::geometry::model::d2::point_xy<double> bg_point;
typedef boost::geometry::model::polygon< bg_point, false, true > bg_polygon;
// point is tested (out of polygon)
bg_point    p(0.1377, 5.00);
// polygon
bg_polygon  poly;
boost::geometry::read_wkt("POLYGON((0.1277 4.97,  0.1277 5.00, 0.1278 4.9999999999999982, 0.1278 4.97, 0.1277 4.97))", poly);
bool inside;
// inside = true here (ERROR!)
inside = boost::geometry::covered_by(p, poly);
// change 4.9999999999999982 to 5.00
boost::geometry::read_wkt("POLYGON((0.1277 4.97,  0.1277 5.00, 0.1278 5.00, 0.1278 4.97, 0.1277 4.97))", poly);
// inside = false here (OK)
inside = boost::geometry::covered_by(p, poly);
return 0;
}

这听起来像是一个bug。您应该将其发布到boost跟踪:

https://svn.boost.org/trac/boost/