交叉 number/Winding 数多边形测试表示为二进制 (Python)

Question

我正在尝试实现绕组数或交叉数测试，主要围绕布尔运算构建。布尔值要求是由于底层数据集的方法和效率，使得委托变量进行计数不是布尔值的最佳选择。

交叉数似乎最容易实现（我认为），因为它本质上是二进制的（偶数 (0) 与奇数 (1)），其中每侧交叉数测试的结果可以是 xor-ed 使用之前的结果，例如在下面给出的代码中，其中 xyz 是我们评估的坐标。代码改编自最后的http://geomalgorithms.com/a03-_inclusion.html。

#Original points:
pts=[[100,100],[200,200],[300,100],[400,300],[300,400],[200,300],[100,100]]

#Extremas:
min=[pts[0][0],pts[0][1]]
max=[pts[0][0],pts[0][1]]

for i in pts:
  for j in range(2):
    if i[j]<min[j]:
      min[j]=i[j]
    if i[j]>max[j]:
      max[j]=i[j]

#New dimensions:
w=max[0]-min[0]
h=max[1]-min[1]
if len(sys.argv) > 2:
  xyz=[int(sys.argv[1]),int(sys.argv[2])]
else:
  xyz=[200,100]

#Normalize by cutting off lower than minimum, higher than maximum:
for i,p in enumerate(pts):
  pts[i]=[p[0]-min[0],p[1]-min[1]]

x=0
y=1
logic=None
counting=0
for i in range(len(pts)-1):
  test=(  ( (pts[i][y] <= xyz[y]) and (pts[i+1][y] > xyz[y]) ) or \
           ( (pts[i][y] > xyz[y]) and (pts[i+1][y] <= xyz[y]) ) ) and \
             (xyz[x] < pts[i][x] + ( (xyz[y]-pts[i][y])/(pts[i+1][y]-pts[i][y]) ) * (pts[i+1][x] - pts[i][x]))
  if logic is None:
    logic=test
  else:
    logic^=test
  if test:
    counting+=1
print logic
print counting

结果：对于整个图像，二进制流会生成这些图像，其中每个方块都是一个步骤。显然有些不对劲，但我似乎无法弄清楚为什么它在绕过右下角后变得一团糟......有什么想法吗？

Answer 1

啊哈！

&!==and 和 |!==or。通过更改它起作用的运算符。

交叉 number/Winding 数多边形测试表示为二进制 (Python)

Crossing number/Winding number polygon test expressed as binary (Python)

python

vector

geometry

point-in-polygon