提前抱歉,目前在手机上!
所以我基本上有一个包含大约50,000个lat/long元组的列表(列表1)和另一个包含大约1,800个lat/long元组的列表(列表2)。
我想要做的是:对于列表1中的每个列表元素,我想要找到列表2中列表元素中最近的点,这样我基本上就得到了一个大约50,000个值的列表,这些值代表了最小距离。
我在使用geopy.distance计算单个元素的距离时没有遇到任何问题,但是,我被for循环实现所困,非常感谢任何帮助!
多谢。
from math import sin, cos, sqrt, atan2
def distanceCheck(lat1, lat2, lon1, lon2):
R = 6373.0
dlon = lon2 - lon1
dlat = lat2 - lat1
a = (sin(dlat/2))**2 + cos(lat1) * cos(lat2) * (sin(dlon/2))**2
c = 2 * atan2(sqrt(a), sqrt(1-a))
distance = R * c
return distance
distarr = []
for p1 in list1:
minDist = None
point = None
for p2 in list2:
#DISTANCE CHECK HERE -
check = distanceCheck(p1.lat, p2.lat, p1.lon, p2.lon)
if not minDist:
minDist = check
point = p2
else:
if check < minDist:
minDist = check
point = p2
distarr.append({'min': minDist, 'to': point, 'from': p1})
print("{}".format(distarr))
清单1和清单2是带有lat和lon的列表。 希望这有帮助