关于delaunay的源码

一大波蓝猫

在DelaunayTriangulator.cpp里 有如下代码

1. // begin triangulation   
   
2.     GLuint pidx = 0;
   
3.     osg::Vec3Array::const_iterator i;   
   
4.    
   
5.     OSG_INFO << "DelaunayTriangulator: triangulating vertex grid (" << (points->size()-3) <<" points)\n";   
   
6. 
   
7.     for (i=points->begin(); i!=points->end(); ++i, ++pidx)
   
8.     {
   
9. 
   
10.         // don't process supertriangle vertices
    
11.         if (pidx > last_valid_index) break;
    
12. 
    
13.         Edge_set edges;
    
14. 
    
15.         // iterate through triangles
    
16.         Triangle_list::iterator j, next_j;
    
17.         for (j=triangles.begin(); j!=triangles.end(); j = next_j)
    
18.         {
    
19. 
    
20.             next_j = j;
    
21.             ++next_j;
    
22. 
    
23.             // get the circumcircle (x,y centre & radius)
    
24.             osg::Vec3 cc = j->get_circumcircle();
    
25. 
    
26.             // OPTIMIZATION: since points are pre-sorted by the X component,
    
27.             // check whether we can discard this triangle for future operations
    
28.             float xdist = i->x() - cc.x();
    
29.             // this is where the circumcircles radius rather than R^2 is faster.
    
30.             // original code used r^2 and needed to test xdist*xdist>cc.z && i->x()>cc.x().
    
31.             if ((xdist ) > cc.z() )
    
32.             {
    
33.                 discarded_tris.push_back(*j); // these are not needed for further tests as no more
    
34.                 // points will ever lie inside this triangle.
    
35.                 triangles.erase(j);
    
36.             }
    
37.             else
    
38.             {
    
39. 
    
40.                 // if the point lies in the triangle's circumcircle then add
    
41.                 // its edges to the edge list and remove the triangle
    
42.                 if (point_in_circle(*i, cc))
    
43.                 {
    
44.                     for (int ei=0; ei<3; ++ei)
    
45.                     {
    
46. 
    
47.                         std::pair<Edge_set::iterator, bool> result = edges.insert(j->get_edge(ei));
    
48.                         if (!result.second)
    
49.                         {
    
50.                             // cast away constness of a set element, which is
    
51.                             // safe in this case since the set_duplicate is
    
52.                             // not used as part of the Less operator.
    
53.                             Edge& edge = const_cast<Edge&>(*(result.first));
    
54.                             // not clear why this change is needed? But prevents removal of twice referenced edges??
    
55.                       //      edge.set_duplicate(true);
    
56.                             edge.set_duplicate(!edge.get_duplicate());
    
57.                         }
    
58.                     }                    
    
59.                     triangles.erase(j);
    
60.                 }
    
61.                                 //add my own code for remove the triangles where at least one edge is above a threshold
    
62.             }
    
63.                        
    
64.         }

复制代码

我想在//add my own code for remove the triangles where at least one edge is above a threshold 这里添加一些处理 就是去掉那些边长大于一个阈值的三角形 但是不知道他的数据结构是什么？ 应该去哪里看？ 或者有知道的告诉下  怎么获得三角形的边长？感激不尽！