ブースト/ポリゴンを使用してポリゴン内の自己交差を見つける方法

Question

ポリゴンで自己交差を見つける必要があります。 私はブーストがその能力を持っていることを知っています。 しかし、交差点に関する情報を得るためにturn_infoを使用する方法を理解できません。セグメントがどのように交差しているかのように。 誰でも助けることができますか？ ありがとう

Answer 1

Boost Geometryで定義された概念は自己交差を許可しないため、正しくはできません。 （私は1.59だと思うので、新しい）

ただし、間接的にあなたがして、自己交差点に関するいくつかの情報を取得するために、検証機能を使用することができます。

std::string reason; 
poly p; 
bg::read_wkt("POLYGON((0 0, 0 4, 2 4, 2 2, 6 2, 6 6, 2 6, 2 4, 0 4, 0 8, 8 8, 8 0, 0 0))", expected); 
bool ok = bg::is_valid(p, reason); 
std::cout << "Expected: " << bg::dsv(p) << (ok?" valid":" invalid: '" + reason + "'") << "\n"; 

プリント：期待

：（（（0 （0,4）、（2,4）、（2,2）、（6,2）、（6,6）、（2,6）、（2,4）、（0,4） ）、（0,8）、（8,8）、（8,0）、（0、0）））無効： 'ジオメトリに無効な自己交差があります。自己交差点は（0、4）にあった。メソッド：t;操作：x/u;セグメントID {source、multi、ring、segment}：{0、-1、-1、0}/{0、-1、-1、7} '

Answer 2

これは自己交差点。情報を取得するための

namespace bg = boost::geometry; 
    using namespace std; 
    typedef bg::model::d2::point_xy<double> point_2d; 
    typedef bg::model::polygon<boost::geometry::model::d2::point_xy<double> > Polygon; 

    Polygon poly { { { 10, 10 }, { 20, 10 }, { 20, 5 }, { 25, 5 }, { 25, 7 }, { 30, 7 }, { 30, 3 }, { 25, 3 }, { 25, 5 }, { 20, 5 }, { 20, 0 }, { 10, 0 }, { 10, 10 } }}; 


    typedef bg::point_type<Polygon>::type point_type; 
    typedef boost::geometry::detail::overlay::turn_info<point_type, boost::geometry::segment_ratio<double> > TurnInfoType; 

    bg::detail::no_rescale_policy robust_policy; 
    bg::detail::self_get_turn_points::no_interrupt_policy interrupt_policy; 
    std::vector<TurnInfoType> turns; 

boost::geometry::self_turns<boost::geometry::detail::overlay::assign_null_policy>(poly.outer(), robust_policy, turns, interrupt_policy); 

だけのようなものを使用します。私は突然、私は[：ブースト・ジオメトリタグ]を想定し実現

turns[i].operations[0].seg_id.segment_index