如何在boost几何中适应自定义多边形类型

Question

我试图使用boost几何算法与我自己的自定义多边形类型。但是我得到了编译器错误(在Visual 2019 Windows 10中)。

我已经将我想要做的事情简化为下面的代码。

在my_custom_polygon.hpp中

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>

using point = boost::geometry::model::d2::point_xy<double>;
using ring = boost::geometry::model::ring<point>;

struct MyPoly
{
    ring exteriorRing;
    std::vector<ring> interiorRings;
};

using polygon = MyPoly;
using multipolygon = std::vector<MyPoly>;
//using polygon = boost::geometry::model::polygon<point>;
//using multipolygon = boost::geometry::model::multi_polygon<polygon>;

namespace boost::geometry::traits
{
    template<> struct tag<MyPoly> { using type = polygon_tag; };
    template<> struct tag<std::vector<MyPoly>> { using type = multi_polygon_tag; };
    template<> struct ring_const_type<MyPoly> { using type = const ring; };
    template<> struct ring_mutable_type<MyPoly> { using type = ring; };
    template<> struct interior_const_type<MyPoly> { using type = const std::vector<ring>; };
    template<> struct interior_mutable_type<MyPoly> { using type = std::vector<ring>; };

    template<> struct exterior_ring<MyPoly>
    {
        static ring& get(MyPoly& poly) { return poly.exteriorRing; }
        static const ring& get(const MyPoly& poly) { return poly.exteriorRing; }
    };

    template<> struct interior_rings<MyPoly>
    {
        static std::vector<ring>& get(MyPoly& poly) { return poly.interiorRings; }
        static const std::vector<ring>& get(const MyPoly& poly) { return poly.interiorRings; }
    };
}

在my_custom_polygon.cpp中

int main(int argc, const char** argv)
{
    const double buffer_distance = 1.0;
    const int points_per_circle = 36;
    boost::geometry::strategy::buffer::distance_symmetric<double> distance_strategy(buffer_distance);
    boost::geometry::strategy::buffer::join_round join_strategy(points_per_circle);
    boost::geometry::strategy::buffer::end_round end_strategy(points_per_circle);
    boost::geometry::strategy::buffer::point_circle circle_strategy(points_per_circle);
    boost::geometry::strategy::buffer::side_straight side_strategy;

    multipolygon result;

    point p{0.0, 0.0};

    boost::geometry::buffer(p, result,
            distance_strategy, side_strategy,
            join_strategy, end_strategy, circle_strategy);

    return 0
}

这无法使用第70行的boost/geometry/algorithms/detail/overlay/convert_ring.hpp中的错误C2664进行编译。

它说它不能将参数2从'boost::geometry::model::ring‘转换为'Geometry2 &’。

但是，如果我对多边形和多多边形类型使用.hpp文件中的注释行，它就会编译和运行得很好。

我显然没有正确地调整多边形。

有人有什么想法吗？

谢谢

Answer 1

第一想到读标题..。哦，孩子:又来了:)

花了一段时间(清理一下)，我发现原因是ring_{mutable,const}_type::type需要作为这个算法的参考。

这使我认为我以前做的修改不是最优的，导致不必要的环复制，例如应用boost几何学调整几何对象模型的进一步问题和带着它的多边形类型？。

因此，没有进一步的担心：

住在Coliru

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <iostream>

namespace bg = boost::geometry;

using Point = bg::model::d2::point_xy<double>;
using Ring  = bg::model::ring<Point>;
using Rings = std::vector<Ring>;

struct MyPoly {
    Ring  exteriorRing;
    Rings interiorRings;
};

using MyMultiPoly = std::vector<MyPoly>;

namespace boost::geometry::traits {
    template <> struct tag<MyPoly>                   { using type = polygon_tag;  };
    template <> struct ring_mutable_type<MyPoly>     { using type = Ring&;        };
    template <> struct ring_const_type<MyPoly>       { using type = const Ring&;  };
    template <> struct interior_mutable_type<MyPoly> { using type = Rings;        };
    template <> struct interior_const_type<MyPoly>   { using type = const Rings;  };

    template<> struct exterior_ring<MyPoly> {
        static auto& get(MyPoly& poly)       { return poly.exteriorRing; } 
        static auto& get(const MyPoly& poly) { return poly.exteriorRing; } 
    };

    template<> struct interior_rings<MyPoly> {
        static auto& get(MyPoly& poly)       { return poly.interiorRings; } 
        static auto& get(const MyPoly& poly) { return poly.interiorRings; } 
    };
} // namespace boost::geometry::traits

namespace boost::geometry::traits {
    template <> struct tag<MyMultiPoly> { using type = multi_polygon_tag; };
} // namespace boost::geometry::traits

int main() {
    MyMultiPoly result;
    Point p{0.0, 0.0};

    namespace bs = bg::strategy::buffer;
    const double buffer_distance   = 1.0;
    const int    points_per_circle = 36;

    bs::distance_symmetric<double> distance(buffer_distance);
    bs::join_round                 join(points_per_circle);
    bs::end_round                  end(points_per_circle);
    bs::point_circle               circle(points_per_circle);
    bs::side_straight              side;

    bg::buffer(p, result, distance, side, join, end, circle);

    std::cout << "result: " << bg::wkt(result) << "\n";
}

打印

result: MULTIPOLYGON(((1 0,0.984808 -0.173648,0.939693 -0.34202,0.866025 -0.5,0.766044 -0.642788,0.642788 -0.766044,0.5 -0.866025,0.34202 -0.939693,0.173648 -0.984808,6.12323e-17 -1,-0.173648 -0.984808,-0.34202 -0.939693,-0.5 -0.866025,-0.642788 -0.766044,-0.766044 -0.642788,-0.866025 -0.5,-0.939693 -0.34202,-0.984808 -0.173648,-1 -1.45473e-15,-0.984808 0.173648,-0.939693 0.34202,-0.866025 0.5,-0.766044 0.642788,-0.642788 0.766044,-0.5 0.866025,-0.34202 0.939693,-0.173648 0.984808,-2.84823e-15 1,0.173648 0.984808,0.34202 0.939693,0.5 0.866025,0.642788 0.766044,0.766044 0.642788,0.866025 0.5,0.939693 0.34202,0.984808 0.173648,1 0)))