使用PNPoly算法计算指定点是否在多边形内部
代码
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
void Main()
{
var points = new List<Point>{
new Point{ X = 819.8061, Y = 615.9394 },
new Point{ X = 652.1153 ,Y = 469.6320 },
new Point{ X = 894.0853 ,Y = 554.0401 },
new Point{ X = 1042.6436, Y = 396.4782 },
new Point{ X = 981.8698 ,Y = 615.9394 },
new Point{ X = 819.8061 ,Y = 615.9394 },
};
var point = new Point { X=882.8309 ,Y=579.9252 };
var polygon = Polygon(points.Select(i => i.X).ToArray(), points.Select(i => i.Y).ToArray(), point.X, point.Y);
Console.WriteLine(polygon);
point = new Point { X=878.1040 ,Y=530.4058 };
polygon = Polygon(points.Select(i => i.X).ToArray(), points.Select(i => i.Y).ToArray(), point.X, point.Y);
Console.WriteLine(polygon);
}
private static bool Polygon(double[] xs, double[] ys, double x, double y)
{
bool contained = false; // 点是否包含在多边形内
double xMin = xs.Min();
double xMax = xs.Max();
double yMin = ys.Min();
double yMax = ys.Max();
if (x > xMax || x < xMin || y > yMax || y < yMin)
{
contained = false;
}
// 核心算法部分
int N = xs.Length;
for (int i = 0, j = N - 1; i < N; j = i++)
{
if (((ys[j] > y) != (ys[i] > y))
&& (x < (xs[j] - xs[i]) * (y - ys[i]) / (ys[j] - ys[i]) + xs[i]))
{
contained = !contained;
}
}
return contained;
}
struct Point
{
public double X { get; set; }
public double Y { get; set; }
}
参考资料
本文会经常更新,请阅读原文: https://dashenxian.github.io/post/%E8%AE%A1%E7%AE%97%E7%82%B9%E6%98%AF%E5%90%A6%E5%9C%A8%E5%A4%9A%E8%BE%B9%E5%BD%A2%E5%86%85%E9%83%A8 ,以避免陈旧错误知识的误导,同时有更好的阅读体验。
本作品采用 知识共享署名-非商业性使用-相同方式共享 4.0 国际许可协议 进行许可。欢迎转载、使用、重新发布,但务必保留文章署名 小神仙 (包含链接: https://dashenxian.github.io ),不得用于商业目的,基于本文修改后的作品务必以相同的许可发布。如有任何疑问,请 与我联系 (125880321@qq.com) 。