ObjectArx 使用 CGAL ，50ms 十万级离散点生成三角网

枫叶棋语

2. void createDT() {
   
3.   setFocusDoc();
   
4.   // 使用最快的Epick内核（现代C++语法）
   
5.   using Kernel = CGAL::Epick;
   
6.   using CGAL_Point = Kernel::Point_2;
   
7.   using Delaunay = CGAL::Delaunay_triangulation_2<Kernel>;
   
8. 
   
9.   DbTrans tr;
   
10.   auto btr = tr.currentSpace();
    
11.   int count = 1000;
    
12.   auto res = MapleGetInt();
    
13.   if (res.Ok) {
    
14.     count = res.Value;
    
15.   }
    
16. 
    
17.   // 生成随机点集
    
18.   std::vector<CGAL_Point> points;
    
19.   points.reserve(count + 1);  // 预分配内存
    
20.   std::random_device rd;
    
21.   std::mt19937 gen(rd());
    
22.   std::uniform_real_distribution<> dis(0.0, 10000.0);
    
23. 
    
24.   for (int i = 0; i <= count; ++i) {
    
25.     points.emplace_back(dis(gen), dis(gen));
    
26.   }
    
27. 
    
28.   // 在CAD中绘制点
    
29.   for (const auto& pt : points) {
    
30.     AcDbPoint* dbpoint = new AcDbPoint({ pt.x(), pt.y(), 0 });
    
31.     dbpoint->setColorIndex(1);
    
32.     tr.addEntity(btr, dbpoint);
    
33.   }
    
34. 
    
35.   // CGAL三角剖分（带空间排序优化）
    
36.   TimeSpan span; //计时
    
37.   CGAL::spatial_sort(points.begin(), points.end());   // 加速关键步骤
    
38.   Delaunay dt(points.begin(), points.end());
    
39.   span.hasRuntime(L"CGAL 三角剖分耗时:"); //打印计时
    
40. 
    
41.   // 绘制三角形
    
42.   for (auto face = dt.finite_faces_begin(); face != dt.finite_faces_end(); ++face) {
    
43.     AcDbPolyline* poly = new AcDbPolyline();
    
44.     poly->setColorIndex(3);
    
45.     poly->setClosed(true);
    
46. 
    
47.     // 添加三角形三个顶点
    
48.     for (int i = 0; i < 3; ++i) {
    
49.       auto& pt = face->vertex(i)->point();
    
50.       poly->addVertexAt(i, { pt.x(), pt.y() });
    
51.     }
    
52.     tr.addEntity(btr, poly);
    
53.   }
    
54. }
    

void createDT() {<br />   setFocusDoc();<br />   // 使用最快的Epick内核（现代C++语法）<br />   using Kernel = CGAL::Epick;<br />   using CGAL_Point = Kernel::Point_2;<br />   using Delaunay = CGAL::Delaunay_triangulation_2<Kernel>;<br /> <br />   DbTrans tr;<br />   auto btr = tr.currentSpace();<br />   int count = 1000;<br />   auto res = MapleGetInt();<br />   if (res.Ok) {<br />     count = res.Value;<br />   }<br /> <br />   // 生成随机点集 <br />   std::vector<CGAL_Point> points;<br />   points.reserve(count + 1);  // 预分配内存 <br />   std::random_device rd;<br />   std::mt19937 gen(rd());<br />   std::uniform_real_distribution<> dis(0.0, 10000.0);<br /> <br />   for (int i = 0; i <= count; ++i) {<br />     points.emplace_back(dis(gen), dis(gen));<br />   }<br /> <br />   // 在CAD中绘制点 <br />   for (const auto& pt : points) {<br />     AcDbPoint* dbpoint = new AcDbPoint({ pt.x(), pt.y(), 0 });<br />     dbpoint->setColorIndex(1);<br />     tr.addEntity(btr, dbpoint);<br />   }<br /> <br />   // CGAL三角剖分（带空间排序优化）<br />   TimeSpan span; //计时<br />   CGAL::spatial_sort(points.begin(), points.end());   // 加速关键步骤 <br />   Delaunay dt(points.begin(), points.end());<br />   span.hasRuntime(L"CGAL 三角剖分耗时:"); //打印计时<br /> <br />   // 绘制三角形 <br />   for (auto face = dt.finite_faces_begin(); face != dt.finite_faces_end(); ++face) {<br />     AcDbPolyline* poly = new AcDbPolyline();<br />     poly->setColorIndex(3);<br />     poly->setClosed(true);<br /> <br />     // 添加三角形三个顶点 <br />     for (int i = 0; i < 3; ++i) {<br />       auto& pt = face->vertex(i)->point();<br />       poly->addVertexAt(i, { pt.x(), pt.y() });<br />     }<br />     tr.addEntity(btr, poly);<br />   }<br /> }<br />

枫叶棋语

绘制很慢，代码耗时全部花在绘制上

bskidtf

速度不如倪尧翻译的那个C库，依赖性也小

枫叶棋语

bskidtf 发表于 2025-6-30 01:13
速度不如倪尧翻译的那个C库，依赖性也小

这个速度已经非常快了啊

gzxl

只需平面（不考虑Z值的），可以做到比CGAL还快。
当然楼主又提供了另一种实现的方式，值得借鉴。