我正在尝试检测一个给定的点(x,y)是否在一个n*2数组的多边形内。但是似乎一些位于多边形边界上的点会返回不包括在内的结果。
def point_inside_polygon(x,y,poly):
n = len(poly)
inside =False
p1x,p1y = poly[0]
for i in range(n+1):
p2x,p2y = poly[i % n]
if y > min(p1y,p2y):
if y <= max(p1y,p2y):
if x <= max(p1x,p2x):
if p1y != p2y:
xinters = (y-p1y)*(p2x-p1x)/float((p2y-p1y))+p1x
if p1x == p2x or x <= xinters:
inside = not inside
p1x,p1y = p2x,p2y
return inside
p1x == p2x
是不好的:它可能相等也可能不相等,存在精度损失问题。 - Jean-François Fabre