CGAL包中是否有一个函数按斜率对行(Line_2)进行排序? 或者任何人都可以推荐一种排序算法来考虑退化情况,如垂直线?
答案 0 :(得分:3)
您可以使用CGAL::Direction_2
。 Line_2
有一个构造函数,operator <
之间有一个Direction_2
,允许对它们进行排序。
您可以使用std::map<CGAL::Direction_2, CGAL::Line_2>
,在此地图中插入of Direction_2
对及其对应的Line_2
。您直接获取排序行,std::map
已排序。
答案 1 :(得分:3)
CGAL::compare_slopes
个对象的Line_2
免费功能,以及相应的K::Compare_slope_2
仿函数(可以从k
的给定实例K
中提取k.compare_slope_2_object()
)。
您可以在以下示例中看到它:
#include <algorithm>
#include <iostream>
#include <vector>
#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef typename K::Line_2 Line_2;
struct Slope_comparator {
// K k;
// Slope_comparator (const K &k = K()) : k(k) {}
bool operator() (const Line_2 &l1, const Line_2 &l2) const {
return (CGAL::compare_slopes (l1, l2) < 0);
// return (k.compare_slope_2_object()(l1, l2) < 0);
}
};
int main () {
std::vector< Line_2 > l;
l.push_back (Line_2 ( 1., 1., 0.));
l.push_back (Line_2 ( 1., -1., 0.));
l.push_back (Line_2 ( 0., 1., 0.)); // vertical
l.push_back (Line_2 (1e-100, 1., 0.)); // almost vertical
l.push_back (Line_2 (1e-100, -1., 0.)); // almost vertical
l.push_back (Line_2 ( 0., -1., 1.)); // also vertical
l.push_back (Line_2 (1e-100, 1., 2.)); // also almost vertical
l.push_back (Line_2 (1e-100, -1., 3.)); // also almost vertical
l.push_back (Line_2 ( 1., 0., 0.)); // horizontal
l.push_back (Line_2 ( 1., 1e-100, 0.)); // almost horizontal
l.push_back (Line_2 ( -1., 1e-100, 0.)); // almost horizontal
l.push_back (Line_2 ( -1., 0., 4.)); // also horizontal
l.push_back (Line_2 ( -1., 1e-100, 5.)); // also almost horizontal
l.push_back (Line_2 ( 1., 1e-100, 6.)); // also almost horizontal
std::cout << "insertion order:" << std::endl;
for (int i = 0; i < l.size(); ++i)
std::cout << " " << l[i] << std::endl;
std::sort (l.begin(), l.end(), Slope_comparator());
std::cout << "sorted order:" << std::endl;
for (int i = 0; i < l.size(); ++i)
std::cout << " " << l[i] << std::endl;
}
请注意,compare_slopes
会比较面向的行,主要是方向;如果这不是你想要的,你必须规范你的线(例如确保所有 y - 坐标为正)。
答案 2 :(得分:1)
事实上,确实没有必要定义struct Slope_comparator
并使其调用全局函数。您可以直接将仿函数CGAL::Compare_slope_2
传递给std::sort()
函数。