/// <summary> /// 点P是否在多边形区域内 /// </summary> /// <param name="PG">多边形PL</param> /// <param name="P">点P</param> /// <returns>如果点P在区域内返回True,否则返回False.</returns> /// <remarks>此处不考虑平行边的情况</remarks> public static Boolean InPolygon(PolygonI PG, PointI P) { //count ← 0; //以P为端点,作从右向左的射线L; //for 多边形的每条边s // do if P在边s上 // then return true; // if s不是水平的 // then if s的一个端点在L上 // if 该端点是s两端点中纵坐标较大的端点 // then count ← count+1 // else if s和L相交 // then count ← count+1; //if count mod 2 = 1 // then return true; //else return false; Int32 count = 0; LineI L = new LineI(P, new PointI(Int32.MaxValue, P.Y)); foreach (LineI S in PG) { if (LineAlgorithm.OnLine(S, P) == true) { return(true); } if (LineAlgorithm.Gradient(S) != 0)//不是水平线段 { if (LineAlgorithm.OnLine(L, S.Starting) || LineAlgorithm.OnLine(L, S.End)) { if (LineAlgorithm.OnLine(L, S.Starting) && S.Starting.Y > S.End.Y) { ++count; } if (LineAlgorithm.OnLine(L, S.End) && S.End.Y > S.Starting.Y) { ++count; } } else if (LineAlgorithm.HasIntersection(S, L) > 0) { ++count; } } } //如果X的水平射线和多边形的交点数是奇数个则点在多边形内,否则不在区域内。 if (count % 2 == 0) { return(false); } return(true); }
/// <summary> /// 线段L是否在多边形区域内 /// </summary> /// <param name="PG">多边形PG</param> /// <param name="L">线段L</param> /// <returns>如果线段L在区域内返回True,否则返回False.</returns> public static Boolean InPolygon(PolygonI PG, LineI L) { //if 线端PQ的端点不都在多边形内 // then return false; //点集pointSet初始化为空; //for 多边形的每条边s // do if 线段的某个端点在s上 // then 将该端点加入pointSet; // else if s的某个端点在线段PQ上 // then 将该端点加入pointSet; // else if s和线段PQ相交 // 这时候已经可以肯定是内交了 // then return false; //将pointSet中的点按照X-Y坐标排序; //for pointSet中每两个相邻点 pointSet[i] , pointSet[ i+1] // do if pointSet[i] , pointSet[ i+1] 的中点不在多边形中 // then return false; //return true; List <PointI> PointList = new List <PointI>(); if (InPolygon(PG, L.Starting) == false) { return(false); } foreach (LineI S in PG) { if (LineAlgorithm.OnLine(S, L.Starting) == true) { PointList.Add(L.Starting); } else if (LineAlgorithm.OnLine(S, L.End) == true) { PointList.Add(L.End); } else if (LineAlgorithm.OnLine(L, S.Starting) == true) { PointList.Add(S.Starting); } else if (LineAlgorithm.OnLine(L, S.End) == true) { PointList.Add(S.End); } else if (LineAlgorithm.HasIntersection(L, S) > 0) { return(false); } } PointI[] OrderedPointList = PointList.ToArray(); for (Int32 i = 0; i < (OrderedPointList.Length - 1); ++i) { for (Int32 j = 0; j < (OrderedPointList.Length - i - 1); j++) { MinMax <PointI> MM = new MinMax <PointI>(OrderedPointList[j], OrderedPointList[j + 1]); OrderedPointList[j] = MM.Min; OrderedPointList[j + 1] = MM.Max; } } for (Int32 i = 0; i < (OrderedPointList.Length - 1); ++i) { if (false == InPolygon(PG, PointAlgorithm.MidPoint(OrderedPointList[i], OrderedPointList[i + 1]))) { return(false); } } return(true); }