您好,登錄后才能下訂單哦!
本文實例為大家分享了python3射線法判斷點是否在多邊形內的具體代碼,供大家參考,具體內容如下
#!/usr/bin/python3.4 # -*- coding:utf-8 -*- def isPointinPolygon(point, rangelist): #[[0,0],[1,1],[0,1],[0,0]] [1,0.8] # 判斷是否在外包矩形內,如果不在,直接返回false lnglist = [] latlist = [] for i in range(len(rangelist)-1): lnglist.append(rangelist[i][0]) latlist.append(rangelist[i][1]) print(lnglist, latlist) maxlng = max(lnglist) minlng = min(lnglist) maxlat = max(latlist) minlat = min(latlist) print(maxlng, minlng, maxlat, minlat) if (point[0] > maxlng or point[0] < minlng or point[1] > maxlat or point[1] < minlat): return False count = 0 point1 = rangelist[0] for i in range(1, len(rangelist)): point2 = rangelist[i] # 點與多邊形頂點重合 if (point[0] == point1[0] and point[1] == point1[1]) or (point[0] == point2[0] and point[1] == point2[1]): print("在頂點上") return False # 判斷線段兩端點是否在射線兩側 不在肯定不相交 射線(-∞,lat)(lng,lat) if (point1[1] < point[1] and point2[1] >= point[1]) or (point1[1] >= point[1] and point2[1] < point[1]): # 求線段與射線交點 再和lat比較 point12lng = point2[0] - (point2[1] - point[1]) * (point2[0] - point1[0])/(point2[1] - point1[1]) print(point12lng) # 點在多邊形邊上 if (point12lng == point[0]): print("點在多邊形邊上") return False if (point12lng < point[0]): count +=1 point1 = point2 print(count) if count%2 == 0: return False else: return True if __name__ == '__main__': print(isPointinPolygon([0.8,0.8], [[0,0],[1,1],[0,1],[0,0]]))
以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持億速云。
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。