Triangle.NET

st788796

一个开源的高效 三角网 库 https://github.com/garykac/triangle.net

一个简单测试，没有使用 USE_Z 生成凸包 和 测试拾取点所在的三角形

2. 
   
3.         [CommandMethod("tt")]
   
4.         public void CmdTinTriangle()
   
5.         {
   
6.             var doc = Application.DocumentManager.MdiActiveDocument;
   
7.             var ed = doc.Editor;
   
8.             var filter = new SelectionFilter(new[] { new TypedValue(0, "POINT") });
   
9.             var pso = new PromptSelectionOptions() { MessageForAdding = "\nselect points: " };
   
10.             var psr = ed.GetSelection(pso, filter);
    
11.             if (psr.Status != PromptStatus.OK) return;
    
12.             using (var tr = doc.Database.TransactionManager.StartTransaction())
    
13.             {
    
14.                 try
    
15.                 {
    
16.                     var polygon = new TriangleNet.Geometry.Polygon();
    
17.                     foreach (ObjectId id in psr.Value.GetObjectIds())
    
18.                     {
    
19.                         var pnt = (DBPoint)tr.GetObject(id, OpenMode.ForRead);
    
20.                         polygon.Add(new TriangleNet.Geometry.Vertex(pnt.Position.X, pnt.Position.Y));
    
21.                     }
    
22. 
    
23.                     var options = new ConstraintOptions { Convex = true };
    
24.                     var mesh = (TriangleNet.Mesh)polygon.Triangulate(options);
    
25. 
    
26. 
    
27.                     var btr = (BlockTableRecord)tr.GetObject(doc.Database.CurrentSpaceId, OpenMode.ForWrite);
    
28.                     var pts = new List<Point2d>();
    
29.                     foreach (var sg in mesh.Segments)
    
30.                     {
    
31.                         if (sg.Label != 1) continue;
    
32. 
    
33.                         var v1 = mesh.Vertices.First(t => t.ID == sg.P0);
    
34.                         pts.Add(new Point2d(v1.X, v1.Y));
    
35. 
    
36.                     }
    
37. 
    
38.                     var pline = new Polyline();
    
39.                     for (int i = 0; i < pts.Count; i++)
    
40.                     {
    
41.                         pline.AddVertexAt(i,pts[i], 0.0, 0.0, 0.0);
    
42.                     }
    
43.                     pline.Closed = true;
    
44. 
    
45.                     btr.AppendEntity(pline);
    
46.                     tr.AddNewlyCreatedDBObject(pline, true);
    
47. 
    
48.                     var ppr = doc.Editor.GetPoint("\ntest point: ");
    
49.                     if (ppr.Status == PromptStatus.OK)
    
50.                     {
    
51.                         var qs = new TriangleNet.Tools.TriangleQuadTree(mesh);
    
52. 
    
53.                         var tri =(Triangle) qs.Query(ppr.Value.X, ppr.Value.Y);
    
54. 
    
55.                         var pnts = mesh.Vertices;
    
56.                         var p1 = pnts.First(t => t.ID == tri.GetVertexID(0));
    
57.                         var p2 = pnts.First(t => t.ID == tri.GetVertexID(1));
    
58.                         var p3 = pnts.First(t => t.ID == tri.GetVertexID(2));
    
59. 
    
60.                         Addpline(tr, btr,
    
61.                             new List<Point2d>()
    
62.                                 {new Point2d(p1.X, p1.Y), new Point2d(p2.X, p2.Y), new Point2d(p3.X, p3.Y)});
    
63. 
    
64.                     }
    
65. 
    
66.                     tr.Commit();
    
67.                 }
    
68.                 catch
    
69.                 {
    
70.                     tr.Abort();
    
71.                 }
    
72.             }
    
73.         }
    
74. 
    
75.         private void Addpline(Transaction tr, BlockTableRecord btr, List<Point2d> pts)
    
76.         {
    
77.             var pline = new Polyline();
    
78.             for (var i = 0; i < pts.Count; i++)
    
79.             {
    
80.                 pline.AddVertexAt(i, pts[i], 0.0, 0.0, 0.0);
    
81.             }
    
82.             pline.Closed = true;
    
83. 
    
84.             btr.AppendEntity(pline);
    
85.             tr.AddNewlyCreatedDBObject(pline, true);
    
86.         }
    
87. 
    

<br />         [CommandMethod("tt")]<br />         public void CmdTinTriangle()<br />         {<br />             var doc = Application.DocumentManager.MdiActiveDocument;<br />             var ed = doc.Editor;<br />             var filter = new SelectionFilter(new[] { new TypedValue(0, "POINT") });<br />             var pso = new PromptSelectionOptions() { MessageForAdding = "\nselect points: " };<br />             var psr = ed.GetSelection(pso, filter);<br />             if (psr.Status != PromptStatus.OK) return;<br />             using (var tr = doc.Database.TransactionManager.StartTransaction())<br />             {<br />                 try<br />                 {<br />                     var polygon = new TriangleNet.Geometry.Polygon();<br />                     foreach (ObjectId id in psr.Value.GetObjectIds())<br />                     {<br />                         var pnt = (DBPoint)tr.GetObject(id, OpenMode.ForRead);<br />                         polygon.Add(new TriangleNet.Geometry.Vertex(pnt.Position.X, pnt.Position.Y));<br />                     }<br /> <br />                     var options = new ConstraintOptions { Convex = true };<br />                     var mesh = (TriangleNet.Mesh)polygon.Triangulate(options);<br /> <br /> <br />                     var btr = (BlockTableRecord)tr.GetObject(doc.Database.CurrentSpaceId, OpenMode.ForWrite);<br />                     var pts = new List<Point2d>();<br />                     foreach (var sg in mesh.Segments)<br />                     {<br />                         if (sg.Label != 1) continue;<br /> <br />                         var v1 = mesh.Vertices.First(t => t.ID == sg.P0);<br />                         pts.Add(new Point2d(v1.X, v1.Y));<br /> <br />                     }<br /> <br />                     var pline = new Polyline();<br />                     for (int i = 0; i < pts.Count; i++)<br />                     {<br />                         pline.AddVertexAt(i,pts[i], 0.0, 0.0, 0.0);<br />                     }<br />                     pline.Closed = true;<br /> <br />                     btr.AppendEntity(pline);<br />                     tr.AddNewlyCreatedDBObject(pline, true);<br /> <br />                     var ppr = doc.Editor.GetPoint("\ntest point: ");<br />                     if (ppr.Status == PromptStatus.OK)<br />                     {<br />                         var qs = new TriangleNet.Tools.TriangleQuadTree(mesh);<br /> <br />                         var tri =(Triangle) qs.Query(ppr.Value.X, ppr.Value.Y);<br /> <br />                         var pnts = mesh.Vertices;<br />                         var p1 = pnts.First(t => t.ID == tri.GetVertexID(0));<br />                         var p2 = pnts.First(t => t.ID == tri.GetVertexID(1));<br />                         var p3 = pnts.First(t => t.ID == tri.GetVertexID(2));<br /> <br />                         Addpline(tr, btr,<br />                             new List<Point2d>()<br />                                 {new Point2d(p1.X, p1.Y), new Point2d(p2.X, p2.Y), new Point2d(p3.X, p3.Y)});<br /> <br />                     }<br /> <br />                     tr.Commit();<br />                 }<br />                 catch<br />                 {<br />                     tr.Abort();<br />                 }<br />             }<br />         }<br /> <br />         private void Addpline(Transaction tr, BlockTableRecord btr, List<Point2d> pts)<br />         {<br />             var pline = new Polyline();<br />             for (var i = 0; i < pts.Count; i++)<br />             {<br />                 pline.AddVertexAt(i, pts[i], 0.0, 0.0, 0.0);<br />             }<br />             pline.Closed = true;<br /> <br />             btr.AppendEntity(pline);<br />             tr.AddNewlyCreatedDBObject(pline, true);<br />         }<br /> <br />

CAD56390

三角函数库