c++ - 在boost::geometry中获取凸包的交集

经过多次尝试使用浮点数导致了由于舍入问题而导致的异常，我认为使用整数算术作为解决方法可以解决问题。但是，现在我遇到了完全相同的问题。

我正在尝试计算各种点集的凸包的交集:

#include <iostream>
#include <boost/geometry/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/multi/geometries/multi_point.hpp>
#include <boost/geometry/geometries/polygon.hpp>

int main()
{
    typedef boost::geometry::model::d2::point_xy<int> Point;
    typedef boost::geometry::model::multi_point<Point> MultiPoint;
    typedef boost::geometry::model::ring<Point> Polygon;

    MultiPoint mp0, mp1;
    boost::geometry::read_wkt("MULTIPOINT((54 74),(54 75),(54 75),(62 75),(86 75),(94 75),(118 75),(124 75),(13 50),(13 51),(147 130),(281 51),(281 50))", mp0);
    boost::geometry::read_wkt("MULTIPOINT((52 74),(54 75),(135 90),(175 74),(54 74),(52 74))", mp1);

    Polygon hull0, hull1;
    boost::geometry::convex_hull(mp0, hull0);
    boost::geometry::convex_hull(mp1, hull1);

    std::vector<Polygon> results;
    boost::geometry::intersection(hull0, hull1, results);

    assert(results.size() == 1);

    // This results in the exception.
    assert(!boost::geometry::detail::overlay::has_self_intersections(results[0]));

    return EXIT_SUCCESS;
}

使用boost::geometry::overlay_invalid_input_exception失败。

凸包hull0和hull1看起来像这样:

我做错什么了吗？我真的很想不必自己计算凸包和交集，这似乎是很多不必要的容易出错的工作。

最佳答案

用例看起来还可以。

我们进行了一些数字健壮性升级，但尚未发布(增强版1.55)。如果您想测试它们或提出一些更详细的问题，建议您通过Boost.Geometry邮件列表http://lists.boost.org/mailman/listinfo.cgi/geometry与我们联系。

关于c++ - 在boost::geometry中获取凸包的交集，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/23157391/