取消

计算点是否在多边形内部

使用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)

登录 GitHub 账号进行评论