CGALでポリゴン交差点を取得

Question

2つの交差するポリゴンの交差線（最初の交差点から最後の交差点まで）を簡単にどのように取得できますか？明確にするために画像を参照してください、緑色の線は私が望むものです。

現在、私は私が交差ポリゴンを取得し、両方のポリゴンの境界上にあるポイントを見つけ、これらは交点である必要があり、次のアルゴリズムを使用します。ここにコード内にあります：

Polygon_2 P,Q; 
Pwh_list_2     intR; 
Pwh_list_2::const_iterator it; 
CGAL::intersection(P, Q, std::back_inserter(intR)); 

//Loop through intersection polygons 
for (it = intR.begin(); it != intR.end(); ++it) { 
    boost::numeric::ublas::vector<double> firstIntersectPoint(3), lastIntersectPoint(3); 
    Polygon_2 Overlap = it->outer_boundary(); 
    typename CGAL::Polygon_2<Kernel>::Vertex_iterator vit; 
    int pointNr = 1; 

    //Loop through points of intersection polygon to find first and last intersection point. 
    for (vit = Overlap.vertices_begin(); vit != Overlap.vertices_end(); ++vit) { 
     CGAL::Bounded_side bsideThis = P.bounded_side(*vit); 
     CGAL::Bounded_side bsideArg = Q.bounded_side(*vit); 
     if (bsideThis == CGAL::ON_BOUNDARY && bsideArg == CGAL::ON_BOUNDARY && pointNr == 1) { 
      firstIntersectPoint <<= 0, CGAL::to_double(vit->x()), CGAL::to_double(vit->y()); 
      pointNr = 2; 
     } 
     else if (bsideThis == CGAL::ON_BOUNDARY && bsideArg == CGAL::ON_BOUNDARY && pointNr == 2) { 
      lastIntersectPoint <<= 0, CGAL::to_double(vit->x()), CGAL::to_double(vit->y()); 
      pointNr = 2; 
     } 
    } 

    //RESULT 
    std::cout << firstIntersectPoint << std::endl; 
    std::cout << lastIntersectPoint << std::endl; 
} 

これは動作しますが、正しい方法ではないと思います。誰かが私にこれを行う正しい方法であるかどうかを示すことができますか、それをより良くする方法を私の指針に教えてください。

Answer 1

2つのポリゴンのセグメントを2D配列に挿入します。次に、次数が4の頂点を見つけます。一般的なケースでは、2以上の頂点があることに注意してください。何もないかもしれない、とあなたは「セグメントのリストの構築を避け、代わりに直接イテレータアダプタを使用して配列にポリゴンのセグメントを挿入することができます。1.

std::list<X_monotone_curve_2> segments; 
for (const auto& pgn : polygons) { 
    for (auto it = pgn.curves_begin(); it != pgn.curves_end(); ++it) 
    segments.push_back(*it); 
} 
Arrangement_2 arr; 
insert(arr, segments.begin(), segments.end()); 
for (auto it = arr.begin_vertices(); it != arr.end_vertices(); ++it) { 
    if (4 == it->degree()) 
    ... 
} 

があるかもしれません。 （これは純粋な汎用プログラミングであり、CGALとは関係ありません）。