交叉数/绕组数多边形测试表示为二进制（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。通过改变它运作的运算符。