Получение точек пересечения двух окружностей в CGAL - PullRequest
Новая
       51

Получение точек пересечения двух окружностей в CGAL

0 голосов
/ 27 ноября 2018

Я пытаюсь найти точки пересечения двух кругов единичного радиуса в CGAL.Вдохновленный примерами кодов CGAL и учебными пособиями, мне удалось создать следующий код.Но, к сожалению, я не могу распечатать пункты.Я заметил, что размер вектора - это число точек пересечения.Любая помощь будет оценена, так как я новичок в CGAL.

    #include <CGAL/Circular_kernel_intersections.h>
    #include <CGAL/Exact_circular_kernel_2.h>

    typedef CGAL::Exact_circular_kernel_2             Circular_k;
    typedef CGAL::Point_2<Circular_k>                 Point_2;
    typedef CGAL::Circle_2<Circular_k>                Circle_2;
    typedef CGAL::Circular_arc_2<Circular_k>          Circular_arc_2;
    typedef CGAL::CK2_Intersection_traits<Circular_k, Circle_2, Circle_2>::type Intersection_result;

    using namespace std;

    int main() {

        Point_2 p(2,2), r(5.5,2);
        Circle_2 c1(p,1), c2(r,1);

        vector<Intersection_result> res;
        intersection(c1,c2,back_inserter(res));

        cout << res.size() << endl;
    }

1 Ответ

0 голосов
/ 27 ноября 2018

Я проверил ваш код, и он, кажется, работает.Обратите внимание, что точки пересечения хранятся как Circular_arc_point_2.Используя ваш код, я добавил только следующие строки: 

using boostRetVal = std::pair<CGAL::Circular_arc_point_2<CGAL::Filtered_bbox_circular_kernel_2<CGAL::Circular_kernel_2<CGAL::Cartesian<CGAL::Gmpq>, CGAL::Algebraic_kernel_for_circles_2_2<CGAL::Gmpq> > > > , unsigned>;

for(const auto& element : res) {
  auto algPoint = std::get<0>( boost::get< boostRetVal >(element) );
  auto point = Point_2(to_double(algPoint.x()), to_double(algPoint.y()));
  std::cout << point << std::endl;
}

В качестве точек я использовал p(0,0), r(2,0), в качестве кружков c1(p,4), c2(r,1), тогда полученный результат был: 

2
7/4 -4360591588697965/4503599627370496
7/4 4360591588697965/4503599627370496
...