我试图找出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;
}
答案 0 :(得分:3)
我测试了您的代码,它似乎可以工作。请注意,交点存储为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