判断点是否在多边形内（C#版本）

菜鸟Liu

   /// <summary>
        /// 判断点是否在闭合多线段内
        /// </summary>
        /// <param name="pt">待判断点</param>
        /// <param name="poly">目标多边形</param>
        /// <returns>在多边形内则返回1</returns>
        
        private double min(double x,double y)
        {
            if (x > y)
                return y;
            else
                return x;
        }
        private double max(double x, double y)
        {
            if (x > y)
                return x;
            else
                return y;
        }
        public int PtinPolygon(Point3d pt,Point3d[] ptPolygon)
        {
            int nCount = ptPolygon.Length;
            // 记录是否在多边形的边上
            bool isBeside = false;
            // 多边形外接矩形
            double maxx,maxy,minx,miny;
            if (nCount > 0)
            {
                maxx = ptPolygon[0].X;
                minx = ptPolygon[0].X;
                maxy = ptPolygon[0].Y;
                miny = ptPolygon[0].Y;
                for (int j = 1; j < nCount;j++)
                {
                    if (ptPolygon[j].X >= maxx)
                        maxx = ptPolygon[j].X;
                    else if (ptPolygon[j].X <= minx)
                        minx = ptPolygon[j].X;
                    if (ptPolygon[j].Y >= maxy)
                        maxy = ptPolygon[j].Y;
                    else if (ptPolygon[j].Y <= miny)
                        miny = ptPolygon[j].Y;
                }
                if ((pt.X > maxx) || (pt.X < minx) || (pt.Y > maxy) || (pt.Y < miny))
                    return -1;
            }
            // 射线法判断
            int nCross = 0;
            for (int i = 0; i < nCount;i++)
            {
                Point3d p1 = ptPolygon;
                Point3d p2 = ptPolygon[(i + 1) % nCount];
                if (p1.Y==p2.Y)   
                {
                    if (pt.Y == p1.Y && pt.X >= min(p1.X,p2.X) && pt.X <= max(p1.X,p2.X))
                    {
                        isBeside = true;
                        continue;
                    }
                }
                // 交点在p1p2延长线上
                if (pt.Y < min(p1.Y, p2.Y) || pt.Y > max(p1.Y, p2.Y))
                    continue;
                // 求交点的X坐标
                double x = (double)(pt.Y - p1.Y) * (double)(p2.X - p1.X) / (double)(p2.Y - p1.Y) + p1.X;
                if (x > pt.X)
                    nCross++;  // 只统计单边交点
                else if (x == pt.X)
                    isBeside = true;
            }
            if (isBeside)
                return 0;  //多边形上
            else if (nCross % 2 == 1)
                return 1;  // 多边形内
            return -1;     // 多边形外
        }


该贴已经同步到 菜鸟Liu的微博

zz0147

，你直接用系统api吧。

gyl

zz0147 发表于 2012-3-25 20:48
，你直接用系统api吧。

能具体说一下吗？

jiangfei200809

用lz的方法 在有些时候判别出现了问题
你看下这样代码应该好点

2. public int ptInPolygon1(Point3d pt, Polyline pPolyline)
   
3.         {
   
4.             int count = pPolyline.NumberOfVertices;
   
5.             //构建多边形外接矩形
   
6.             double maxx = double.MinValue, maxy = double.MinValue, minx = double.MaxValue, miny = double.MaxValue;
   
7.             for (int i = 0; i < count; i++)
   
8.             {
   
9.                 if (pPolyline.GetPoint3dAt(i).X > maxx)
   
10.                 {
    
11.                     maxx = pPolyline.GetPoint3dAt(i).X;
    
12.                 }
    
13.                 if (pPolyline.GetPoint3dAt(i).Y > maxy)
    
14.                 {
    
15.                     maxy = pPolyline.GetPoint3dAt(i).Y;
    
16.                 }
    
17.                 if (pPolyline.GetPoint3dAt(i).X < minx)
    
18.                 {
    
19.                     minx = pPolyline.GetPoint3dAt(i).X;
    
20.                 }
    
21.                 if (pPolyline.GetPoint3dAt(i).Y < miny)
    
22.                 {
    
23.                     miny = pPolyline.GetPoint3dAt(i).Y;
    
24.                 }
    
25.             }
    
26.             if (pt.X > maxx || pt.Y > maxy || pt.X < minx || pt.Y < miny)
    
27.             {
    
28.                 return -1;
    
29.             }
    
30. 
    
31.             Line line1 = new Line(new Point3d(maxx, pt.Y, 0), new Point3d(minx, pt.Y, 0));
    
32. 
    
33.             int crossCount = 0;
    
34. 
    
35.             using (Transaction tran = pDatabase.TransactionManager.StartTransaction())
    
36.             {
    
37.                 BlockTableRecord pBlockTableRecord = tran.GetObject(pDatabase.CurrentSpaceId, OpenMode.ForWrite) as BlockTableRecord;
    
38. 
    
39.                 Point3dCollection crossPoint = new Point3dCollection();
    
40. 
    
41.                 line1.IntersectWith(pPolyline, Intersect.OnBothOperands, crossPoint, 0, 0);
    
42.                 if (crossPoint.Count >= 1)
    
43.                 {
    
44.                     for (int n = 0; n < crossPoint.Count; n++)
    
45.                     {
    
46.                         Point3d crossPt = crossPoint[n];
    
47.                         if (crossPt.X > pt.X)
    
48.                         {
    
49.                             crossCount++;
    
50. 
    
51.                             Circle pCircle = new Circle(crossPt, Vector3d.ZAxis, 2);
    
52. 
    
53.                             pBlockTableRecord.AppendEntity(pCircle);
    
54.                             tran.AddNewlyCreatedDBObject(pCircle, true);
    
55.                         }
    
56.                     }
    
57. 
    
58.                 }
    
59.                 Circle circle = new Circle(pt, Vector3d.ZAxis, 2);
    
60.                 pBlockTableRecord.AppendEntity(circle);
    
61.                 tran.AddNewlyCreatedDBObject(circle, true);
    
62. 
    
63.                 pBlockTableRecord.AppendEntity(line1);
    
64.                 tran.AddNewlyCreatedDBObject(line1, true);
    
65.                 tran.Commit();
    
66.             }
    
67.             return crossCount % 2;
    
68. 
    
69.         }

public int ptInPolygon1(Point3d pt, Polyline pPolyline)<br />         {<br />             int count = pPolyline.NumberOfVertices;<br />             //构建多边形外接矩形<br />             double maxx = double.MinValue, maxy = double.MinValue, minx = double.MaxValue, miny = double.MaxValue;<br />             for (int i = 0; i < count; i++)<br />             {<br />                 if (pPolyline.GetPoint3dAt(i).X > maxx)<br />                 {<br />                     maxx = pPolyline.GetPoint3dAt(i).X;<br />                 }<br />                 if (pPolyline.GetPoint3dAt(i).Y > maxy)<br />                 {<br />                     maxy = pPolyline.GetPoint3dAt(i).Y;<br />                 }<br />                 if (pPolyline.GetPoint3dAt(i).X < minx)<br />                 {<br />                     minx = pPolyline.GetPoint3dAt(i).X;<br />                 }<br />                 if (pPolyline.GetPoint3dAt(i).Y < miny)<br />                 {<br />                     miny = pPolyline.GetPoint3dAt(i).Y;<br />                 }<br />             }<br />             if (pt.X > maxx || pt.Y > maxy || pt.X < minx || pt.Y < miny)<br />             {<br />                 return -1;<br />             }<br /> <br />             Line line1 = new Line(new Point3d(maxx, pt.Y, 0), new Point3d(minx, pt.Y, 0));<br /> <br />             int crossCount = 0;<br /> <br />             using (Transaction tran = pDatabase.TransactionManager.StartTransaction())<br />             {<br />                 BlockTableRecord pBlockTableRecord = tran.GetObject(pDatabase.CurrentSpaceId, OpenMode.ForWrite) as BlockTableRecord;<br /> <br />                 Point3dCollection crossPoint = new Point3dCollection();<br /> <br />                 line1.IntersectWith(pPolyline, Intersect.OnBothOperands, crossPoint, 0, 0);<br />                 if (crossPoint.Count >= 1)<br />                 {<br />                     for (int n = 0; n < crossPoint.Count; n++)<br />                     {<br />                         Point3d crossPt = crossPoint[n];<br />                         if (crossPt.X > pt.X)<br />                         {<br />                             crossCount++;<br /> <br />                             Circle pCircle = new Circle(crossPt, Vector3d.ZAxis, 2);<br /> <br />                             pBlockTableRecord.AppendEntity(pCircle);<br />                             tran.AddNewlyCreatedDBObject(pCircle, true);<br />                         }<br />                     }<br /> <br />                 }<br />                 Circle circle = new Circle(pt, Vector3d.ZAxis, 2);<br />                 pBlockTableRecord.AppendEntity(circle);<br />                 tran.AddNewlyCreatedDBObject(circle, true);<br /> <br />                 pBlockTableRecord.AppendEntity(line1);<br />                 tran.AddNewlyCreatedDBObject(line1, true);<br />                 tran.Commit();<br />             }<br />             return crossCount % 2;<br /> <br />         }

sdaulj

东西不错，顶了，感谢分享