我正在尝试编写一个 Ray Casting python 代码,它将返回一个True
or列表False
,它将显示点是否在多边形中。我有一个包含多个多边形的纬度/经度的列表,以及两个附加列表,一个是纬度,另一个是点的经度。所以基本上我要做的是使用列表,看看多个点是否在多个多边形中。任何建议都会非常有帮助,因为我是新的 Python,并且边走边学!
下面是创建光线投射以查看一个点是否是一个多边形的函数,因此需要查看多个点是否在多个多边形中,为每个点返回 True 或 False,因此返回[True, False, True, False]
多边形中一点的光线投射:
def point_in_poly(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:
xints = (y-p1y)*(p2x-p1x)/float(p2y-p1y)+p1x
if p1x == p2x or x <= xints:
inside = not inside
p1x,p1y = p2x,p2y
return inside
## Test
polygon = [(0,11),(11,10),(11,0),(0,1)]
point_x = 12
point_y = 12
## Call the function with the points and the polygon
print (point_in_poly(point_x,point_y,polygon))