Point IsInside Points

csharp

Lisp http://bbs.xdcad.net/forum.php?mod=viewthread&tid=668873
-1 外部，0 边上或者顶点，1 内部

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

复制代码

雪山飞狐(lzh)

GetAngleTo可以直接获取两向量的夹角

public double GetAngleTo(
    Vector3d vector
);
获取两向量间的角度，在区间[0,PI]内

public double GetAngleTo(
    Vector3d vector,
    Vector3d referenceVector
);
以referenceVector为观察方向，相对于vector的角度，在区间[0,2 * PI]内
double angle = vec.GetAngleTo(Vector3d.XAxis, -Vector3d.ZAxis);

雪山飞狐(lzh)

ln.IsOn(pt) || pt == p1？后面的判断会不会多余

csharp

雪山飞狐(lzh) 发表于 2014-9-2 18:53
GetAngleTo可以直接获取两向量的夹角

public double GetAngleTo(

angleto都是正值，这个判断不能用，需要负值
后面判断是多余