改进弧长法 判断点是否在多边形内 c#代码

public class Coordinate
{

/// <summary>
/// 经度X
/// </summary>
public double X;

/// <summary>
/// 纬度Y
/// </summary>
public double Y;

/// <summary>
/// 所在象限，原点，坐标轴1,2,3,4
/// 为0时当正号处理
/// </summary>
public int Quadrant
{
get
{
if (X >= 0)
{
if (Y >= 0) return 1;
if (Y < 0) return 4;
}
if (X < 0)
{
if (Y >= 0) return 2;
if (Y < 0) return 3;

}

return 0;
}
}

}

/// <summary>
/// 弧长法判断点是否在多边形内
/// 第一个是若P[i]的某个坐标为0时，一律当正号处理；第二点是若被测点和多边形的顶点重合时要特殊处理。
/// </summary>
/// <param name="p">多边形顶点集合，有序</param>
/// <param name="corrdinate">要判断的点</param>
/// <returns>是否在多边形内</returns>
/// <remarks>
/// </remarks>
public bool InPolygon(IList<Coordinate> p, Coordinate corrdinate)
{
foreach (var item in p)
{
item.Y = item.Y - corrdinate.Y;
item.X = item.X - corrdinate.X;
}
//移除相邻相同象限的
for (var i = p.Count - 2; i > 0; i--)
{
if (p[i - 1].Quadrant == p[i].Quadrant)
p.RemoveAt(i);
}

//最后一个和第一个点连，将第一个点加到最后
p.Add(p[0]);

var count = 0;
for (var i = 0; i <= p.Count - 2; i++)
{
if (p[i].X == 0 && p[i].Y == 0) return true;//顶点重合

var f = p[i].X * p[i + 1].Y - p[i].Y * p[i + 1].X;
if (f == 0) return true;//在边上,特殊情况，如边(3,0)-(-3,0)

var t = p[i + 1].Quadrant - p[i].Quadrant;
if (Math.Abs(t) == 2) t = f > 0 ? 2 : -2;
if (t == -3) t = 1;//4->1（1-4=-3） =1
if (t == 3) t = -1;//1->4（4-1=3）=-1
count += t;

}
if (Math.Abs(count) == 2) return true;//在边上
if (Math.Abs(count) == 4) return true;//在内部
if (count == 0) return false;//在外部
return false;
}