Brep GetLineContainment

csharp

利用 Brep.GetLineContainment() 判断点和多边形关系

1. 
   
2.         public static bool IsInside(this Point3d pt, Point3dCollection pts)
   
3.         {
   
4.             DBObjectCollection dbCol = new DBObjectCollection();
   
5.             Polyline pl = new Polyline();
   
6.             Plane plan = pl.GetPlane();
   
7.             for (int i = 0; i < pts.Count; i++)
   
8.             {
   
9.                 pl.AddVertexAt(i, pts[i].Convert2d(plan), 0, 0, 0);
   
10.             }
    
11.             pl.Closed = true;
    
12.             dbCol.Add(pl);
    
13.             Region region = Region.CreateFromCurves(dbCol)[0] as Region;
    
14.             Brep brep = new Brep(region);
    
15.             LineSegment3d ln = new LineSegment3d(pt, pt);
    
16.             Hit[] hit = brep.GetLineContainment(ln, 1);
    
17.             if (hit == null)
    
18.             {
    
19.                 return false;
    
20.             }
    
21.             else
    
22.             {
    
23.                 return true;
    
24.             }
    
25.         }
    

复制代码

雪山飞狐(lzh)

为什么不用GetPointContainment

csharp

雪山飞狐(lzh) 发表于 2014-9-2 18:46
为什么不用GetPointContainment

有这个函数？只在Arx手册中看到Point,没注意到net中这个函数
再修改修改

csharp

这样改一改

1. 
   
2.         public static int IsInside(this Point3d pt, IEnumerable<Point3d> pts)
   
3.         {
   
4.             double pi = Math.PI;
   
5.             double fpi = -pi;
   
6.             Point3dCollection npts = new Point3dCollection(pts.ToArray());
   
7.             if (pts.First( ) != pts.Last())
   
8.             {
   
9.                 npts.Add(pts.First( ));
   
10.             }
    
11.             Plane plane = new Plane(new Point3d(0, 0, 0), new Vector3d(0, 0, 1));
    
12.             double tol = 0.0;
    
13.             for (int i = 0; i < npts.Count - 1; i++)
    
14.             {
    
15.                 Point3d p1 = npts[i];
    
16.                 Point3d p2 = npts[i + 1];
    
17.                 LineSegment3d ln = new LineSegment3d(p1, p2);
    
18.                 if (ln.IsOn(pt))
    
19.                 {
    
20.                     return 0;
    
21.                 }
    
22.                 Vector3d vec1 = pt - p1;
    
23.                 Vector3d vec2 = pt - p2;
    
24.                 double an = vec1.AngleOnPlane(plane) - vec2.AngleOnPlane(plane);
    
25.                 if (an > pi)
    
26.                 {
    
27.                     an -= pi;
    
28.                 }
    
29.                 else if (an < fpi)
    
30.                 {
    
31.                     an += pi;
    
32.                 }
    
33.                 tol += an;
    
34.             }
    
35.             if (Math.Abs(tol).Equals(pi))
    
36.             {
    
37.                 return 1;
    
38.             }
    
39.             else
    
40.             {
    
41.                 return -1;
    
42.             }
    
43.         }
    

复制代码

1. 
   
2.         public static int IsInside(this Point3d pt, Curve cv)
   
3.         {
   
4.             if (!cv.Closed)
   
5.             {
   
6.                 return -1;
   
7.             }
   
8.             else if (cv is Circle)
   
9.             {
   
10.                 Circle circle = cv as Circle;
    
11.                 Point3d pc = circle.Center;
    
12.                 double radius = circle.Radius;
    
13.                 double d = pt.DistanceTo(pc);
    
14.                 if (Math .Abs( d-radius) < 1e-6)
    
15.                 {
    
16.                     return 0;
    
17.                 }
    
18.                 else if (d > radius)
    
19.                 {
    
20.                     return -1;
    
21.                 }
    
22.                 else
    
23.                 {
    
24.                     return 1;
    
25.                 }
    
26.             }
    
27.             else if (cv is Ellipse)
    
28.             {
    
29.                 Curve3d c3d = cv.GetGeCurve();
    
30.                 EllipticalArc3d e3d = c3d as EllipticalArc3d;
    
31.                 if (e3d.IsOn(pt))
    
32.                 {
    
33.                     return 0;
    
34.                 }
    
35.                 else if (e3d.IsInside(pt))
    
36.                 {
    
37.                     return 1;
    
38.                 }
    
39.                 else
    
40.                 {
    
41.                     return -1;
    
42.                 }
    
43.             }
    
44.             else
    
45.             {
    
46.                 DBObjectCollection dCol = new DBObjectCollection();
    
47.                 dCol.Add(cv);
    
48.                 Region region = Region.CreateFromCurves(dCol )[0] as Region;
    
49.                 Brep brep = new Brep(region);
    
50.                 PointContainment pc;
    
51.                 brep.GetPointContainment(pt, out pc);
    
52.                 switch (pc)
    
53.                 {
    
54.                     case PointContainment.Inside://内环内时?
    
55.                         return 1;
    
56.                     case PointContainment.OnBoundary://区域内或边上
    
57.                         return 0;
    
58.                     case PointContainment.Outside:
    
59.                         return -1;
    
60.                 }
    
61.             }
    
62.             return -1;
    
63.         }

复制代码

雪山飞狐(lzh)

直接用这个枚举 inside.outside.onboundary作为返回值更直观些

csharp

雪山飞狐(lzh) 发表于 2014-9-2 21:28
直接用这个枚举 inside.outside.onboundary作为返回值更直观些

直接return pc;  ?